Showing 1–50 of 186 results for author: Hovatta, T

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.

Comments: 23 pages, 17 figures, 5 tables, submitted to ApJ

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.

Comments: Accepted for publication in A&A, 24 pages, 10 figures, 2 tables, 3 electronic tables

Distinguishing radiation mechanisms and particle populations in blazar jets through long-term multi-band monitoring with RINGO3 and Fermi

Authors: Callum McCall, Helen Jermak, Iain A. Steele, Iván Agudo, Ulisses Barres de Almeida, Talvikki Hovatta, Gavin P. Lamb, Elina Lindfors, Carole Mundell

Abstract: We present the results of seven years of multicolour photometric monitoring of a sample of 31 $γ$-ray bright blazars using the RINGO3 polarimeter on the Liverpool Telescope from 2013--2020. We explore the relationships between simultaneous observations of flux in three optical wavebands along with Fermi $γ$-ray data in order to explore the radiation mechanisms and particle populations in blazar je… ▽ More We present the results of seven years of multicolour photometric monitoring of a sample of 31 $γ$-ray bright blazars using the RINGO3 polarimeter on the Liverpool Telescope from 2013--2020. We explore the relationships between simultaneous observations of flux in three optical wavebands along with Fermi $γ$-ray data in order to explore the radiation mechanisms and particle populations in blazar jets. We find significant correlations between optical and $γ$-ray flux with no detectable time lag, suggesting leptonic emission processes in the jets of these sources. Furthermore, we find the spectral behaviour against optical and $γ$-ray flux for many sources is best fit logarithmically. This is suggestive of a transition between bluer-/redder-when-brighter into stable-when-brighter behaviour during high activity states; a behaviour that might be missed in poorly sampled data, resulting in apparent linear relationships. △ Less

Submitted 11 July, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

Comments: 32 pages, 35 figures

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 23 July, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 44 pages JCAP style (excluding author list and references), 19 figures; minor changes to match published version

Journal ref: JCAP 07 (2024) 047

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.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

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.

Comments: Accepted for publication in Astronomy & Astrophysics. 19 pages, 9 figures. Corresponding authors: Lea Heckmann, Axel Arbet Engels, David Paneque

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.

Comments: 31 pages, 37 figures, submitted to MNRAS

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.

Comments: Accepted for publication in MNRAS. 16 pages, 18 figures, 12 tables

A Kpc-scale Radio Polarization Study of PG BL Lacs with the uGMRT

Authors: Janhavi Baghel, P. Kharb, T. Hovatta, S. Gulati, E. Lindfors, Silpa S

Abstract: We present here uGMRT band 4 (~650MHz) polarization images of 8 BL~Lac objects belonging to the Palomar-Green (PG) `blazar' sample. A large fraction of the sources (~63%) reveal core-halo radio structures with most of the polarization detected in the inner core-jet regions. PG1101+385 and PG2254+075 exhibit a `spine-sheath structure' in polarization. The core-halo and `spine-sheath' structures are… ▽ More We present here uGMRT band 4 (~650MHz) polarization images of 8 BL~Lac objects belonging to the Palomar-Green (PG) `blazar' sample. A large fraction of the sources (~63%) reveal core-halo radio structures with most of the polarization detected in the inner core-jet regions. PG1101+385 and PG2254+075 exhibit a `spine-sheath structure' in polarization. The core-halo and `spine-sheath' structures are consistent with the Unified Scheme suggestion that BL~Lacs are the pole-on beamed counterparts of Fanaroff-Riley (FR) type I radio galaxies. PG1418+546 and PG0851+203 (OJ287) show the presence of terminal hotspots similar to FR type II radio galaxies. They were also found to be low-spectrally peaked BL Lacs, supportive of the `blazar envelope' scenario for BL~Lacs and quasars. Fractional polarization ranges from 1-13% in the cores and 2-26% in the inner jets/lobes of the sample BL Lacs. Compared to the varied radio morphology of quasars from the PG `blazar' sample, the BL~Lacs appear to be less diverse. A comparison of the inferred core magnetic (B-) field structures on arcsec- (kpc-) scales w.r.t. the Very Long Baseline Interferometry (VLBI) jet direction does not reveal any preferred orientation, suggesting that if large-scale ordered B-fields exist, they do so on scales smaller than probed by the current observations. However, the presence of polarized emission on arcsec-scales suggests that any mixing of thermal plasma with the synchrotron emitting plasma is insufficient to fully depolarize the emission via the internal depolarization process. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: 17 pages; 12 figures; MNRAS accepted

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

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.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 682, A114 (2024)

Variability studies of active galactic nuclei from the long-term monitoring program with the Cherenkov Telescope Array

Authors: Guillaume Grolleron, Josefa Becerra González, Jonathan Biteau, Matteo Cerruti, Roger Grau, Lucas Gréaux, Talvikki Hovatta, Jean-Philippe Lenain, Elina Lindfors, Walter Max-Moerbeck, Davide Miceli, Abelardo Moralejo, Kari Nilsson, Elisa Pueschel, Arkadipta Sarkar, Sofia Suutarinen

Abstract: Blazars are active galactic nuclei (AGN) with a relativistic jet oriented toward the observer. This jet is composed of accelerated particles which can display emission over the entire electromagnetic spectrum. Spectral variability has been observed on short- and long-time scales in AGN, with a power spectral density (PSD) that can show a break at frequencies below the well-known red-noise process.… ▽ More Blazars are active galactic nuclei (AGN) with a relativistic jet oriented toward the observer. This jet is composed of accelerated particles which can display emission over the entire electromagnetic spectrum. Spectral variability has been observed on short- and long-time scales in AGN, with a power spectral density (PSD) that can show a break at frequencies below the well-known red-noise process. This break frequency in the PSD has been observed in X-rays to scale with the accretion regime and the mass of the central black hole. It is expected that a break could also be seen in the very-high-energy gamma rays, but constraining the shape of the PSD in these wavelengths has not been possible with the current instruments. The Cherenkov Telescope Array (CTA) will be more sensitive by a factor of five to ten depending on energy than the current generation of imaging atmospheric Cherenkov telescopes, therefore it will be possible with CTA to reconstruct the PSD with a high accuracy, bringing new information about AGN variability. In this work, we focus on the AGN long-term monitoring program planned with CTA. The program is proposed to begin with early-start observing campaigns with CTA precursors. This would allow us to probe longer time scales on the AGN PSD. △ Less

Submitted 21 September, 2023; originally announced September 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:submit/5125530)

Report number: CTA-ICRC/2023/19

Bright blazar flares with CTA

Authors: M. Cerruti, J. Finke, G. Grolleron, J. P. Lenain, T. Hovatta, M. Joshi, E. Lindfors, P. Morris, M. Petropoulou, P. Romano, S. Vercellone, M. Zacharias

Abstract: The TeV extragalactic sky is dominated by blazars, radio-loud active galactic nuclei with a relativistic jet pointing towards the Earth. Blazars show variability that can be quite exceptional both in terms of flux (orders of magnitude of brightening) and time (down to the minute timescale). This bright flaring activity contains key information on the physics of particle acceleration and photon pro… ▽ More The TeV extragalactic sky is dominated by blazars, radio-loud active galactic nuclei with a relativistic jet pointing towards the Earth. Blazars show variability that can be quite exceptional both in terms of flux (orders of magnitude of brightening) and time (down to the minute timescale). This bright flaring activity contains key information on the physics of particle acceleration and photon production in the emitting region, as well as the structure and physical properties of the jet itself. The TeV band is accessed from the ground by Cherenkov telescopes that image the pair cascade triggered by the interaction of the gamma ray with the Earth's atmosphere. The Cherenkov Telescope Array (CTA) represents the upcoming generation of imaging atmospheric Cherenkov telescopes, with a significantly higher sensitivity and larger energy coverage with respect to current instruments. It will thus provide us with unprecedented statistics on blazar light-curves and spectra. In this contribution we present the results from realistic simulations of CTA observations of bright blazar flares, taking as input state-of-the-art numerical simulations of blazar emission models and including all relevant observational constraints. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/27

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.

Journal ref: ApJ 964 (2024) 3

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.

Comments: 17 pages, 12 figures

Quantitative comparisons of VHE gamma-ray blazar flares with relativistic reconnection models

Authors: J. Jormanainen, T. Hovatta, I. M. Christie, E. Lindfors, M. Petropoulou, I. Liodakis

Abstract: The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute time scales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection… ▽ More The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute time scales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection layer, is a candidate solution to this problem. In this work, we employ state-of-the-art particle-in-cell simulations in a statistical comparison with observations of a flaring episode of a well-known blazar, Mrk 421, at very high energy (VHE, E > 100 GeV). We tested the predictions of our model by generating simulated VHE light curves that we compared quantitatively with methods that we have developed for a precise evaluation of theoretical and observed data. With our analysis, we can constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma, viewing angle and the reconnection layer orientation in the blazar jet. Our analysis favours parameter spaces with magnetic field strength 0.1 G, rather large viewing angles (6-8 degrees), and misaligned layer angles, offering a strong candidate explanation for the Doppler crisis often observed in the jets of high synchrotron peaking blazars. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 678, A140 (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.

Journal ref: A&A 677, A144 (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.

Comments: 18 pages, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

Probing AGN variability with the Cherenkov Telescope Array

Authors: F. Cangemi, T. Hovatta, E. Lindfors, M. Cerruti, J. Becerra-Gonzalez, J. Biteau, C. Boisson, M. Böttcher, E. de Gouveia Dal Pino, D. Dorner, G. Grolleron, J. -P. Lenain, M. Manganaro, W. Max-Moerbeck, P. Morris, K. Nilsson, L. Passos Reis, P. Romano, O. Sergijenko, F. Tavecchio, S. Vercellone, S. Wagner, M. Zacharias

Abstract: Relativistic jets launched by Active Galactic Nuclei are among the most powerful particle accelerators in the Universe. The emission over the entire electromagnetic spectrum of these relativistic jets can be extremely variable with scales of variability from less than few minutes up to several years. These variability patterns, which can be very complex, contain information about the acceleration… ▽ More Relativistic jets launched by Active Galactic Nuclei are among the most powerful particle accelerators in the Universe. The emission over the entire electromagnetic spectrum of these relativistic jets can be extremely variable with scales of variability from less than few minutes up to several years. These variability patterns, which can be very complex, contain information about the acceleration processes of the particles and the area(s) of emission. Thanks to its sensitivity, five-to twenty-times better than the current generation of Imaging Atmospheric Cherenkov Telescopes depending on energy, the Cherenkov Telescope Array will be able to follow the emission from these objects with a very accurate time sampling and over a wide spectral coverage from 20 GeV to > 20 TeV and thus reveal the nature of the acceleration processes at work in these objects. We will show the first results of our lightcurve simulations and long-term behavior of AGN as will be observed by CTA, based on state-of-art particle acceleration models. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: Proceedings for Gamma 2022, 8 pages, 6 figures

Interactions between the jet and disk wind in a nearby radio intermediate quasar III Zw 2

Authors: Ailing Wang, Tao An, Shaoguang Guo, Prashanth Mohan, Wara Chamani, Willem A. Baan, Talvikki Hovatta, Heino Falcke, Tim J. Galvin, Natasha Hurley-Walker, Sumit Jaiswal, Anne Lahteenmaki, Baoqiang Lao, Weijia Lv, Merja Tornikoski, Yingkang Zhang

Abstract: Disk winds and jets are ubiquitous in active galactic nuclei (AGN), and how these two components interact remains an open question. We study the radio properties of a radio-intermediate quasar III Zw 2. We detect two jet knots J1 and J2 on parsec scales, which move at a mildly apparent superluminal speed of $1.35\,c$. Two $γ$-ray flares were detected in III Zw 2 in 2009--2010, corresponding to the… ▽ More Disk winds and jets are ubiquitous in active galactic nuclei (AGN), and how these two components interact remains an open question. We study the radio properties of a radio-intermediate quasar III Zw 2. We detect two jet knots J1 and J2 on parsec scales, which move at a mildly apparent superluminal speed of $1.35\,c$. Two $γ$-ray flares were detected in III Zw 2 in 2009--2010, corresponding to the primary radio flare in late 2009 and the secondary radio flare in early 2010. The primary 2009 flare was found to be associated with the ejection of J2. The secondary 2010 flare occurred at a distance of $\sim$0.3 parsec from the central engine, probably resulting from the collision of the jet with the accretion disk wind. The variability characteristics of III Zw 2 (periodic radio flares, unstable periodicity, multiple quasi-periodic signals and possible harmonic relations between them) can be explained by the global instabilities of the accretion disk. These instabilities originating from the outer part of the warped disk propagate inwards and can lead to modulation of the accretion rate and consequent jet ejection. At the same time, the wobbling of the outer disk may also lead to oscillations of the boundary between the disk wind and the jet tunnel, resulting in changes in the jet-wind collision site. III Zw 2 is one of the few cases observed with jet-wind interactions, and the study in this paper is of general interest for gaining insight into the dynamic processes in the nuclear regions of AGN. △ Less

Submitted 28 December, 2022; originally announced December 2022.

Comments: accepted by ApJ

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.

Comments: 20 pages, 7 figures. Accepted in A&A. Corresponding authors: Jorge Otero-Santos; Daniel Morcuende; Vandad Fallah Ramazani; Daniela Dorner; David Paneque (mailto: [email protected])

Journal ref: A&A 670, A49 (2023)

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.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

Optical polarization from colliding stellar stream shocks in a tidal disruption event

Authors: I. Liodakis, K. I. I. Koljonen, D. Blinov, E. Lindfors, K. D. Alexander, T. Hovatta, M. Berton, A. Hajela, J. Jormanainen, K. Kouroumpatzakis, N. Mandarakas, K. Nilsson

Abstract: A tidal disruption event (TDE) occurs when a supermassive black hole rips apart a passing star. Part of the stellar material falls toward the black hole, forming an accretion disk that in some cases launches a relativistic jet. We performed optical polarimetry observations of a TDE, AT 2020mot. We find a peak linear polarization degree of $25\pm4$%, consistent with highly polarized synchrotron rad… ▽ More A tidal disruption event (TDE) occurs when a supermassive black hole rips apart a passing star. Part of the stellar material falls toward the black hole, forming an accretion disk that in some cases launches a relativistic jet. We performed optical polarimetry observations of a TDE, AT 2020mot. We find a peak linear polarization degree of $25\pm4$%, consistent with highly polarized synchrotron radiation, as is typically observed from relativistic jets. However, our radio observations, taken up to 8 months after the optical peak, do not detect the corresponding radio emission expected from a relativistic jet. We suggest that the linearly polarized optical emission instead arises from shocks that occur during accretion disk formation, as the stream of stellar material collides with itself. △ Less

Submitted 27 November, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 41 pages, 10 figures, 4 tables, published in Science

Journal ref: Volume 380, Issue 6645, pp. 656-658 (2023)

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.

Comments: 8 pages, 1 figure, 1 table, accepted for publication in A&A

Journal ref: A&A 666, A36 (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.

Comments: 18 pages, 6 figures; V2 to match published version, typo fixed

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.

Comments: Accepted to ApJ

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.

Comments: 10 pages, 3 figures, European VLBI Network Mini-Symposium and Users' Meeting (EVN2021), to be published in Proceedings of Science

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.

Comments: ApJ accepted. Comments are welcome

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.

Comments: 21 pages, 12 figures, 5 tables, accepted for publication in MNRAS

The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

Authors: S. O'Neill, S. Kiehlmann, A. C. S. Readhead, M. F. Aller, R. D. Blandford, I. Liodakis, M. L. Lister, P. Mróz, C. P. O'Dea, T. J. Pearson, V. Ravi, M. Vallisneri, K. A. Cleary, M. J. Graham, K. J. B. Grainge, M. W. Hodges, T. Hovatta, A. Lähteenmäki, J. W. Lamb, T. J. W. Lazio, W. Max-Moerbeck, V. Pavlidou, T. A. Prince, R. A. Reeves, M. Tornikoski , et al. (2 additional authors not shown)

Abstract: Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio lig… ▽ More Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field. △ Less

Submitted 24 January, 2022; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 24 pages, 13 figure, 3 Tables, accepted for publication in APJL

Confronting observations of VHE gamma-ray blazar flares with reconnection models

Authors: J. Jormanainen, T. Hovatta, E. Lindfors, I. Christie, M. Petropoulou, I. Liodakis

Abstract: Several models have been suggested to explain the fast gamma-ray variability observed in blazars, but its origin is still debated. One scenario is magnetic reconnection, a process that can efficiently convert magnetic energy to energy of relativistic particles accelerated in the reconnection layer. In our study, we compare results from state-of-the-art particle-in-cell simulations with observation… ▽ More Several models have been suggested to explain the fast gamma-ray variability observed in blazars, but its origin is still debated. One scenario is magnetic reconnection, a process that can efficiently convert magnetic energy to energy of relativistic particles accelerated in the reconnection layer. In our study, we compare results from state-of-the-art particle-in-cell simulations with observations of blazars at Very High Energy (VHE, E > 100 GeV) gamma-rays. Our goal is to test our model predictions on fast gamma-ray variability with data and to constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma and the reconnection layer orientation in the blazar jet. For this first comparison, we used the remarkably well-sampled VHE gamma-ray light curve of Mrk 421 observed with the MAGIC and VERITAS telescopes in 2013.The simulated VHE light curves were generated using the observable parameters of Mrk 421, such as the jet power, bulk Lorentz factor, and the jet viewing angle, and sampled as real data. Our results pave the way for future model-to-data comparison with next-generation Cherenkov telescopes, which will help further constrain the different variability models. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: To be published in the proceedings of 37th International Cosmic Ray Conference

MOJAVE XIX: Brightness Temperatures and Intrinsic Properties of Blazar Jets

Authors: D. C. Homan, M. H. Cohen, T. Hovatta, K. I. Kellermann, Y. Y. Kovalev, M. L. Lister, A. V. Popkov, A. B. Pushkarev, E. Ros, T. Savolainen

Abstract: We present multi-epoch, parsec-scale core brightness temperature observations of 447 AGN jets from the MOJAVE and 2cm Survey programs at 15 GHz from 1994 to 2019. The brightness temperature of each jet over time is characterized by its median value and variability. We find that the range of median brightness temperatures for AGN jets in our sample is much larger than the variations within individu… ▽ More We present multi-epoch, parsec-scale core brightness temperature observations of 447 AGN jets from the MOJAVE and 2cm Survey programs at 15 GHz from 1994 to 2019. The brightness temperature of each jet over time is characterized by its median value and variability. We find that the range of median brightness temperatures for AGN jets in our sample is much larger than the variations within individual jets, consistent with Doppler boosting being the primary difference between the brightness temperatures of jets in their median state. We combine the observed median brightness temperatures with apparent jet speed measurements to find the typical intrinsic Gaussian brightness temperature of (4.1 +- 0.6)*10^10 K, suggesting that jet cores are at or below equipartition between particle and magnetic field energy in their median state. We use this value to derive estimates for the Doppler factor for every source in our sample. For the 309 jets with both apparent speed and brightness temperature data, we estimate their Lorentz factors and viewing angles to the line of sight. Within the BL Lac optical class, we find that high-synchrotron-peaked (HSP) BL Lacs have smaller Doppler factors, lower Lorentz factors, and larger angles to the line of sight than intermediate and low-synchrotron-peaked (LSP) BL Lacs. We confirm that AGN jets with larger Doppler factors measured in their parsec-scale radio cores are more likely to be detected in gamma rays, and we find a strong correlation between gamma-ray luminosity and Doppler factor for the detected sources. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 29 pages, 16 figures, 4 tables; accepted by ApJ; full electronic tables and all Figure Set plots are available as ancillary files

Journal ref: ApJ 923 (2021) 67

Identifying changing jets through their radio variability

Authors: I. Liodakis, T. Hovatta, M. F. Aller, H. D. Aller, M. A. Gurwell, A. Lähteenmäki, M. Tornikoski

Abstract: Supermassive black holes launch highly relativistic jets with velocities reaching Lorentz factors as high as $Γ>50$. How the jets accelerate to such high velocities and where along the jet do they reach terminal velocity are open questions that are tightly linked to their structure, launching and dissipation mechanisms. Changes in the beaming factor along the jets could potentially reveal jet acce… ▽ More Supermassive black holes launch highly relativistic jets with velocities reaching Lorentz factors as high as $Γ>50$. How the jets accelerate to such high velocities and where along the jet do they reach terminal velocity are open questions that are tightly linked to their structure, launching and dissipation mechanisms. Changes in the beaming factor along the jets could potentially reveal jet acceleration, deceleration, or bending. We aim to (1) quantify the relativistic effects in multiple radio frequencies and (2) study possible jet velocity--viewing angle variations at parsec scales. We used the state-of-the-art code Magnetron to model light curves from the University of Michigan Radio Observatory and the Metsähovi Radio Observatory's monitoring programs in five frequencies covering about 25 years of observations in the 4.8-37~GHz range for 61 sources. We supplement our data set with high-frequency radio observations in the 100-340~GHz range from ALMA, CARMA, and SMA. For each frequency we estimate the Doppler factor which we use to quantify possible changes in the relativistic effects along the jets. The majority of our sources do not show any statistically significant difference in their Doppler factor across frequencies. This is consistent with constant velocity in a conical jet, as expected at parsec scales. However, our analysis reveals 17 sources where relativistic beaming changes as a function of frequency. In the majority of cases the Doppler factor increases towards lower frequencies. Only 1253-053 shows the opposite behavior. By exploring their jet properties we find that the jet of 0420-014 is likely bent across the 4.8-340~GHz range. For 0212+735 the jet is likely parabolic, and still accelerating in the 4.8-37~GHz range. We discuss possible interpretations for the trends found in the remaining sources. △ Less

Submitted 26 August, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: 15 pages, 5 figure, 3 Tables, accepted for publication in A&A

Journal ref: A&A 654, A169 (2021)

SMILE: Search for MIlli-LEnses

Authors: C. Casadio, D. Blinov, A. C. S. Readhead, I. W. A. Browne, P. N. Wilkinson, T. Hovatta, N. Mandarakas, V. Pavlidou, K. Tassis, H. K. Vedantham, J. A. Zensus, V. Diamantopoulos, K. E. Dolapsaki, K. Gkimisi, G. Kalaitzidakis, M. Mastorakis, K. Nikolaou, E. Ntormousi, V. Pelgrims, K. Psarras

Abstract: Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact ob… ▽ More Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact objects in the mass range $\sim$ 10$^{6}$ -- 10$^{9}$$M_{\odot}$, we search for strong gravitational lenses on milli (mas)-arcseconds scales (< 150 mas). For our search, we used the Astrogeo VLBI FITS image database -- the largest publicly available database, containing multi-frequency VLBI data of 13828 individual sources. We used the citizen science approach to visually inspect all sources in all available frequencies in search for images with multiple compact components on mas-scales. At the final stage, sources were excluded based on the surface brightness preservation criterion. We obtained a sample of 40 sources that passed all steps and therefore are judged to be milli-arcsecond lens candidates. These sources are currently followed-up with on-going European VLBI Network (EVN) observations at 5 and 22 GHz. Based on spectral index measurements, we suggest that two of our candidates have a higher probability to be associated with gravitational lenses. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 6 pages, 1 figure; accepted for publication

A decade of joint MOJAVE-Fermi AGN monitoring: localisation of the gamma-ray emission region

Authors: I. G. Kramarenko, A. B. Pushkarev, Y. Y. Kovalev, M. L. Lister, T. Hovatta, T. Savolainen

Abstract: Within the MOJAVE VLBA program (Monitoring of Jets in AGN with VLBA Experiments), we have accumulated observational data at 15 GHz for hundreds of jets in $γ$-ray bright active galactic nuclei since the beginning of the Fermi scientific observations in August 2008. We investigated a time delay between the flux density of AGN parsec-scale radio emission at 15 GHz and 0.1$-$300 GeV Fermi LAT photon… ▽ More Within the MOJAVE VLBA program (Monitoring of Jets in AGN with VLBA Experiments), we have accumulated observational data at 15 GHz for hundreds of jets in $γ$-ray bright active galactic nuclei since the beginning of the Fermi scientific observations in August 2008. We investigated a time delay between the flux density of AGN parsec-scale radio emission at 15 GHz and 0.1$-$300 GeV Fermi LAT photon flux, taken from constructed light curves using weekly and adaptive binning. The correlation analysis shows that radio is lagging $γ$-ray radiation by up to 8 months in the observer's frame, while in the source frame, the typical delay is about 2-3 months. If the jet radio emission, excluding the opaque core, is considered, significant correlation is found at greater time lags. We supplement these results with VLBI kinematics and core shift data to conclude that the dominant high-energy production zone is typically located at a distance of several parsecs from the central nucleus. We also found that quasars have on average more significant correlation peak, more distant $γ$-ray emission region from the central engine and shorter variability time scale compared to those of BL Lacertae objects. △ Less

Submitted 16 November, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 12 pages, 12 figures, accepted for publication in MNRAS

Journal ref: MNRAS 510 (2022) 469

Probing extreme environments with the Cherenkov Telescope Array

Authors: C. Boisson, A. M. Brown, A. Burtovoi, M. Cerruti, M. Chernyakova, T. Hassan, J. -P. Lenain, M. Manganaro, P. Romano, H. Sol, F. Tavecchio, S. Vercellone, L. Zampieri, R. Zanin, A. Zech, I. Agudo, R. Alves Batista, E. O. Anguner, L. A. Antonelli, M. Backes, C. Balazs, J. Becerra González, C. Bigongiari, E. Bissaldi, J. Bolmont , et al. (105 additional authors not shown)

Abstract: The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and ex… ▽ More The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and extreme environments are also the most likely sources of multi-messenger emission. The most energetic part of the electromagnetic spectrum corresponds to the very-high-energy (VHE, E>100 GeV) gamma-ray regime, which can be extensively studied with ground based Imaging Atmospheric Cherenkov Telescopes (IACTs). The results obtained by the current generation of IACTs, such as H.E.S.S., MAGIC, and VERITAS, demonstrate the crucial importance of the VHE band in understanding the non-thermal emission of extreme environments in our Universe. In some objects, the energy output in gamma rays can even outshine the rest of the broadband spectrum. The Cherenkov Telescope Array (CTA) is the next generation of IACTs, which, with cutting edge technology and a strategic configuration of ~100 telescopes distributed in two observing sites, in the northern and southern hemispheres, will reach better sensitivity, angular and energy resolution, and broader energy coverage than currently operational IACTs. With CTA we can probe the most extreme environments and considerably boost our knowledge of the non-thermal Universe. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted as input to ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium

Multi-messenger and transient astrophysics with the Cherenkov Telescope Array

Authors: Ž. Bošnjak, A. M. Brown, A. Carosi, M. Chernyakova, P. Cristofari, F. Longo, A. López-Oramas, M. Santander, K. Satalecka, F. Schüssler, O. Sergijenko, A. Stamerra, I. Agudo, R. Alves Batista, E. Amato, E. O. Anguner, L. A. Antonelli, M. Backes, Csaba Balazs, L. Baroncelli, J. Becker Tjus, C. Bigongiari, E. Bissaldi, C. Boisson, J. Bolmont , et al. (120 additional authors not shown)

Abstract: The discovery of gravitational waves, high-energy neutrinos or the very-high-energy counterpart of gamma-ray bursts has revolutionized the high-energy and transient astrophysics community. The development of new instruments and analysis techniques will allow the discovery and/or follow-up of new transient sources. We describe the prospects for the Cherenkov Telescope Array (CTA), the next-generati… ▽ More The discovery of gravitational waves, high-energy neutrinos or the very-high-energy counterpart of gamma-ray bursts has revolutionized the high-energy and transient astrophysics community. The development of new instruments and analysis techniques will allow the discovery and/or follow-up of new transient sources. We describe the prospects for the Cherenkov Telescope Array (CTA), the next-generation ground-based gamma-ray observatory, for multi-messenger and transient astrophysics in the decade ahead. CTA will explore the most extreme environments via very-high-energy observations of compact objects, stellar collapse events, mergers and cosmic-ray accelerators. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted to ASTRONET roadmap on behalf of the CTA consortium

The Time-Dependent Distribution of Optical Polarization Angle Changes in Blazars

Authors: S. Kiehlmann, D. Blinov, I. Liodakis, V. Pavlidou, A. C. S. Readhead, E. Angelakis, C. Casadio, T. Hovatta, N. Kylafis, A. Mahabal, N. Mandarakas, I. Myserlis, G. V. Panopoulou, T. J. Pearson, A. Ramaprakash, P. Reig, R. Skalidis, A. Slowikowska, K. Tassis, J. A. Zensus

Abstract: At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-inte… ▽ More At optical wavelengths, blazar Electric Vector Position Angle (EVPA) rotations linked with gamma-ray activity have been the subject of intense interest and systematic investigation for over a decade. One difficulty in the interpretation of EVPA rotations is the inherent 180° ambiguity in the measurements. It is therefore essential, when studying EVPA rotations, to ensure that the typical time-interval between successive observations -- i.e. the cadence -- is short enough to ensure that the correct modulo 180° value is selected. This optimal cadence depends on the maximum intrinsic EVPA rotation speed in blazars, which is currently not known. In this paper we address the following questions for the RoboPol sample: What range of rotation speeds for rotations greater than 90° can we expect? What observation cadence is required to detect such rotations? Have rapid rotations been missed in EVPA rotation studies thus far? What fraction of data is affected by the ambiguity? And how likely are detected rotations affected by the ambiguity? We answer these questions with three seasons of optical polarimetric observations of a statistical sample of blazars sampled weekly with the RoboPol instrument and an additional season with daily observations. We model the distribution of EVPA changes on time scales from 1-30 days and estimate the fraction of changes exceeding 90°. We show that at least daily observations are necessary to measure >96% of optical EVPA variability in the RoboPol sample of blazars correctly and that intra-day observations are needed to measure the fastest rotations that have been seen thus far. △ Less

Submitted 14 July, 2021; v1 submitted 6 April, 2021; originally announced April 2021.

Comments: accepted to MNRAS, 20 pages, 15 figures

ALMA/NICER observations of GRS 1915+105 indicate a return to a hard state

Authors: Karri I. I. Koljonen, Talvikki Hovatta

Abstract: GRS 1915+105 is a transient black hole X-ray binary consistently emitting 10-100% of the Eddington luminosity in the X-ray band during the last three decades until mid-2018 when the source luminosity suddenly decreased by an order of magnitude. This phase was followed by a change to a state with even lower average X-ray fluxes never seen before during the outburst but presenting renewed flaring ac… ▽ More GRS 1915+105 is a transient black hole X-ray binary consistently emitting 10-100% of the Eddington luminosity in the X-ray band during the last three decades until mid-2018 when the source luminosity suddenly decreased by an order of magnitude. This phase was followed by a change to a state with even lower average X-ray fluxes never seen before during the outburst but presenting renewed flaring activity at different wavelengths. Nevertheless, the mean fluxes were still in decline. GRS 1915+105 has the longest orbital period known among low-mass X-ray binaries, the largest accretion disk size, and therefore the largest mass supply for accretion. The high inclination of the disk allows the study of geometrical effects of the accretion flow such as changes in the height-to-radius ratio or the effect of accretion disk winds on the intrinsic emission that are expected during the outburst decay. In addition, the transient jet is expected to change to a compact, self-absorbed, steady jet. We conducted two full polarization Atacama Large Millimeter Array observations to study the jet properties during the outburst decay by analyzing the spectral, polarization, and intra-epoch variability for both observation epochs. In addition, we analyzed nearly daily Neutron Star Interior Composition Explorer pointing observations consisting of modeling X-ray power spectral densities, spectral energy distributions, and lightcurves with a physically motivated model to follow the changing accretion disk properties throughout the outburst decay and relating them to the jet emission. We show that the X-ray and mm spectral, timing, and polarization properties are consistent with those of a typical decaying X-ray binary outburst and that GRS 1915+105 has descended into the low-luminosity hard X-ray state. △ Less

Submitted 3 February, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 18 pages, 11 figures, 5 tables. Accepted for publication in A&A. Abstract abridged to meet ArXiv size limit

VHE gamma-ray detection of FSRQ QSO B1420+326 and modeling of its enhanced broadband state in 2020

Authors: 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, 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, Ž. Bošnjak, G. Busetto , et al. (209 additional authors not shown)

Abstract: Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over di… ▽ More Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over different phases of the flare. Methods. The source was observed with a number of instruments in radio, near infrared, optical (including polarimetry and spectroscopy), ultra-violet, X-ray and gamma-ray bands. We use dedicated optical spectroscopy results to estimate the accretion disk and the dust torus luminosity. We perform spectral energy distribution modeling in the framework of combined Synchrotron-Self-Compton and External Compton scenario in which the electron energy distribution is partially determined from acceleration and cooling processes. Results. During the enhanced state the flux of both SED components drastically increased and the peaks were shifted to higher energies. Follow up observations with the MAGIC telescopes led to the detection of very-high-energy gamma-ray emission from this source, making it one of only a handful of FSRQs known in this energy range. Modeling allows us to constrain the evolution of the magnetic field and electron energy distribution in the emission region. The gamma-ray flare was accompanied by a rotation of the optical polarization vector during a low polarization state. Also, a new, superluminal radio knot contemporaneously appeared in the radio image of the jet. The optical spectroscopy shows a prominent FeII bump with flux evolving together with the continuum emission and a MgII line with varying equivalent width. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 20 pages, 15 figures, 7 tables, accepted for publication in A&A

Journal ref: A&A 647, A163 (2021)

The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135

Authors: A. C. S. Readhead, V. Ravi, I. Liodakis, M. L. Lister, V. Singh, M. F. Aller, R. D. Blandford, I. W. A. Browne, V. Gorjian, K. J. B. Grainge, M. A. Gurwell, M. W. Hodges, T. Hovatta, S. Kiehlmann, A. Lähteenmäki, T. McAloone, W. Max-Moerbeck, V. Pavlidou, T. J. Pearson, A. L. Peirson, E. S. Perlman, R. A. Reeves, B. T. Soifer, G. B. Taylor, M. Tornikoski , et al. (4 additional authors not shown)

Abstract: PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and i… ▽ More PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: Accepted for publication in The Astrophysical Journal. 28 pages, 15 figures

RoboPol: AGN polarimetric monitoring data

Authors: D. Blinov, S. Kiehlmann, V. Pavlidou, G. V. Panopoulou, R. Skalidis, E. Angelakis, C. Casadio, E. N. Einoder, T. Hovatta, K. Kokolakis, A. Kougentakis, A. Kus, N. Kylafis, E. Kyritsis, A. Lalakos, I. Liodakis, S. Maharana, E. Makrydopoulou, N. Mandarakas, G. M. Maragkakis, I. Myserlis, I. Papadakis, G. Paterakis, T. J. Pearson, A. N. Ramaprakash , et al. (7 additional authors not shown)

Abstract: We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the dataset presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was… ▽ More We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGN), covering observations between 2013, when the instrument was commissioned, and 2017. In total, the dataset presented in this paper includes 5068 observations of 222 AGN with Dec > -25 deg. We describe the current version of the RoboPol pipeline that was used to process and calibrate the entire dataset, and we make the data publicly available for use by the astronomical community. Average quantities summarising optopolarimetric behaviour (average degree of polarization, polarization variability index) are also provided for each source we have observed and for the time interval we have followed it. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: Accepted to MNRAS

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

Association of IceCube neutrinos with radio sources observed at Owens Valley and Metsähovi Radio Observatories

Authors: T. Hovatta, E. Lindfors, S. Kiehlmann, W. Max-Moerbeck, M. Hodges, I. Liodakis, A. Lahteenmaki, T. J. Pearson, A. C. S. Readhead, R. A. Reeves, S. Suutarinen, J. Tammi, M. Tornikoski

Abstract: Identifying the most likely sources for high-energy neutrino emission has been one of the main topics in high-energy astrophysics ever since the first observation of high-energy neutrinos by the IceCube Neutrino Observatory. Active galactic nuclei with relativistic jets, also known as blazars, have been considered to be one of the main candidates because of their ability to accelerate particles to… ▽ More Identifying the most likely sources for high-energy neutrino emission has been one of the main topics in high-energy astrophysics ever since the first observation of high-energy neutrinos by the IceCube Neutrino Observatory. Active galactic nuclei with relativistic jets, also known as blazars, have been considered to be one of the main candidates because of their ability to accelerate particles to high energies. We study the connection between radio emission and IceCube neutrino events using data from the Owens Valley Radio Observatory (OVRO) and Metsähovi Radio Observatory blazar monitoring programs. We identify sources in our radio monitoring sample that are positionally consistent with IceCube high-energy neutrino events. We estimate their mean flux density and variability amplitudes around the neutrino arrival time, and compare these with values from random samples to establish the significance of our results. We find radio source associations within our samples with 15 high-energy neutrino events detected by IceCube. Nearly half of the associated sources are not detected in the $γ$-ray energies, but their radio variability properties and Doppler boosting factors are similar to the $γ$-ray detected objects in our sample, meaning that they could still be potential neutrino emitters. We find that the number of strongly flaring objects in our statistically complete OVRO samples is unlikely to be a random coincidence (at $2σ$ level). Based on our results, we conclude that although it is clear that not all neutrino events are associated with strong radio flaring blazars, observations of large-amplitude radio flares in a blazar at the same time as a neutrino event are unlikely to be a random coincidence. △ Less

Submitted 15 April, 2021; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A (18 pages, 6 figures). Accepted version (v2) includes changes to methodology and associations, but the main conclusions are unaffected. V3 minor fix to numbers in Sect. 3.3.1. and some language edits

Journal ref: A&A 650, A83 (2021)

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella, M. Cerruti , et al. (196 additional authors not shown)

Abstract: Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $ν_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In Augu… ▽ More Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $ν_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $γ$-ray telescope FACT during a high $γ$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $γ$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $γ$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $Γ$ index of the intrinsic spectrum in the VHE $γ$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz △ Less

Submitted 1 August, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: Accepted in MNRAS. Corresponding authors: M. Manganaro, A. Arbet Engels, D. Dorner

Journal ref: MNRAS 496, 3912 (2020)

Testing two-component models on very-high-energy gamma-ray emitting BL Lac objects

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (179 additional authors not shown)

Abstract: Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline… ▽ More Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Journal ref: A&A 640, A132 (2020)

Relativistic Jets of Blazars

Authors: Talvikki Hovatta, Elina Lindfors

Abstract: Relativistic jets of active galactic nuclei have been known to exist for 100 years. Blazars with their jet pointing close to our line of sight are some of the most variable and extreme objects in the universe, showing emission from radio to very-high-energy gamma rays. In this review, we cover relativistic jets of blazars from an observational perspective with the main goal of discussing how obser… ▽ More Relativistic jets of active galactic nuclei have been known to exist for 100 years. Blazars with their jet pointing close to our line of sight are some of the most variable and extreme objects in the universe, showing emission from radio to very-high-energy gamma rays. In this review, we cover relativistic jets of blazars from an observational perspective with the main goal of discussing how observations can be used to constrain theoretical models. We cover a range of topics from multiwavelength observations to imaging of jets with a special emphasis on current open questions in the field. △ Less

Submitted 15 August, 2020; v1 submitted 13 March, 2020; originally announced March 2020.

Comments: To appear in New Astronomy Reviews special volume '100 Years of Astrophysical Jets' (eds. Rob Fender and Ralph Wijers). Added a few more references and fixed some typos

A decade of multi-wavelength observations of the TeV blazar 1ES 1215+303: Extreme shift of the synchrotron peak frequency and long-term optical-gamma-ray flux increase

Authors: Janeth Valverde, Deirdre Horan, Denis Bernard, Stephen Fegan, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, W. Cui, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan , et al. (64 additional authors not shown)

Abstract: Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radi… ▽ More Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed. △ Less

Submitted 12 February, 2020; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: ApJ 891 (2020) 170

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

Authors: MAGIC collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babic, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra Gonzalez, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Z. Bosnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (215 additional authors not shown)

Abstract: We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on… ▽ More We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 46 pages, 17 figures, 9 tables, accepted for publication in ApJ Supplements