Search

We present a study of the radio and optical properties of the high-frequency peaker (HFP) blazar PKS 1614+051 at $z=3.21$ based on the data covering the time period of 1997-2024. The radio data are represented by the almost instantaneous 1-22 GHz measurements from the SAO RAS RATAN-600 radio telescope, the 5 and 8 GHz data from the IAA RAS RT-32 telescopes, and the 37 GHz data from the RT-22 telescope of CrAO RAS. The optical measurements in the $R$ band were collected with the SAO RAS 1-m Zeiss-1000 and 0.5-m AS-500/2 telescopes and the ZTF archive data. We have found low overall variability indices (0.1-0.2) and a spectral peak around 4.6 GHz, which is stable during the long-term period of monitoring. An analysis of the radio light curves reveals significant time delays (0.6 to 6.4 years) between the radio frequencies along with variability timescales ranging from 0.2 to 1.8 years in the source's rest frame, which is similar to the blazars at lower redshifts. Spectral modeling suggests the presence of both synchrotron-self absorption (SSA) and free-free absorption (FFA) processes. Based on the SSA model, we provide estimates of the magnetic field strength, which peaks at $\sim\!30$ mG. A spectroscopic study with the BTA SCORPIO-I spectrograph has found signs of the regular motion of a neutral hydrogen envelope around the blazar center, which confirms the presence of enough gaseous matter to form an external FFA screen. The results highlight the importance of multiwavelength and long-term monitoring to understand the physical mechanisms driving the variability in high-redshift blazars.

In the era of precision measurements of neutrino oscillation parameters, it is necessary for experiments to disentangle discrepancies that may indicate physics beyond the Standard Model in the neutrino sector. KM3NeT/ORCA is a water Cherenkov neutrino detector under construction and anchored at the bottom of the Mediterranean Sea. The detector is designed to study the oscillations of atmospheric neutrinos and determine the neutrino mass ordering. This paper focuses on the initial configuration of ORCA, referred to as ORCA6, which comprises six out of the foreseen 115 detection units of photosensors. A high-purity neutrino sample was extracted during 2020 and 2021, corresponding to an exposure of 433 kton-years. This sample is analysed following a binned log-likelihood approach to search for invisible neutrino decay, in a three-flavour neutrino oscillation scenario, where the third neutrino mass state $\nu_3$ decays into an invisible state, e.g. a sterile neutrino. The resulting best fit of the invisible neutrino decay parameter is $\alpha_3 = 0.92^{+1.08}_{-0.57}\times 10^{-4}~\mathrm{eV^2}$, corresponding to a scenario with $\theta_{23}$ in the second octant and normal neutrino mass ordering. The results are consistent with the Standard Model, within a $2.1\,\sigma$ interval., Comment: 20 pages, 8 figures

KM3NeT/ORCA is an underwater neutrino telescope under construction in the Mediterranean Sea. Its primary scientific goal is to measure the atmospheric neutrino oscillation parameters and to determine the neutrino mass ordering. ORCA can constrain the oscillation parameters $\Delta m^{2}_{31}$ and $\theta_{23}$ by reconstructing the arrival direction and energy of multi-GeV neutrinos crossing the Earth. Searches for deviations from the Standard Model of particle physics in the forward scattering of neutrinos inside Earth matter, produced by Non-Standard Interactions, can be conducted by investigating distortions of the standard oscillation pattern of neutrinos of all flavours. This work reports on the results of the search for non-standard neutrino interactions using the first six detection units of ORCA and 433 kton-years of exposure. No significant deviation from standard interactions was found in a sample of 5828 events reconstructed in the 1 GeV$-$1 TeV energy range. The flavour structure of the non-standard coupling was constrained at 90\% confidence level to be $|\varepsilon_{\mu\tau} | \leq 5.4 \times 10^{-3}$, $|\varepsilon_{e\tau} | \leq 7.4 \times 10^{-2}$, $|\varepsilon_{e\mu} | \leq 5.6 \times 10^{-2}$ and $-0.015 \leq \varepsilon_{\tau\tau} - \varepsilon_{\mu\mu} \leq 0.017$. The results are comparable to the current most stringent limits placed on the parameters by other experiments.

Indirect dark matter detection methods are used to observe the products of dark matter annihilations or decays originating from astrophysical objects where large amounts of dark matter are thought to accumulate. With neutrino telescopes, an excess of neutrinos is searched for in nearby dark matter reservoirs, such as the Sun and the Galactic Centre, which could potentially produce a sizeable flux of Standard Model particles. The KM3NeT infrastructure, currently under construction, comprises the ARCA and ORCA undersea \v{C}erenkov neutrino detectors located at two different sites in the Mediterranean Sea, offshore of Italy and France, respectively. The two detector configurations are optimised for the detection of neutrinos of different energies, enabling the search for dark matter particles with masses ranging from a few GeV/c$^2$ to hundreds of TeV/c$^2$. In this work, searches for dark matter annihilations in the Galactic Centre and the Sun with data samples taken with the first configurations of both detectors are presented. No significant excess over the expected background was found in either of the two analyses. Limits on the velocity-averaged self-annihilation cross section of dark matter particles are computed for five different primary annihilation channels in the Galactic Centre. For the Sun, limits on the spin-dependent and spin-independent scattering cross sections of dark matter with nucleons are given for three annihilation channels.

Recent observations of the Galactic component of the high-energy neutrino flux, together with the detection of the diffuse Galactic gamma-ray emission up to sub-PeV energies, open new possibilities to study the acceleration and propagation of cosmic rays in the Milky Way. At the same time, both large non-astrophysical backgrounds at TeV energies and scarcity of neutrino events in the sub-PeV band currently limit these analyses. Here we use the sample of cascade events with estimated neutrino energies above 200 TeV, detected by the partially deployed Baikal Gigaton Volume Detector (GVD) in six years of operation, to test the continuation of the Galactic neutrino spectrum to sub-PeV energies. We find that the distribution of the arrival directions of Baikal-GVD cascades above 200 TeV in the sky suggests an excess of neutrinos from low Galactic latitudes. We find the excess above 200 TeV also in the most recent IceCube public data sets, both of cascades and tracks. The significant (3.6 sigma in the combined analysis) flux of Galactic neutrinos above 200 TeV challenges often-used templates for neutrino search based on cosmic-ray simulations., Comment: 7 pages, 5 figures, revtex 4.2

The KM3NeT Collaboration has tackled a common challenge faced by the astroparticle physics community, namely adapting the experiment-specific simulation software to work with the CORSIKA air shower simulation output. The proposed solution is an extension of the open-source code gSeaGen, allowing for the transport of muons generated by CORSIKA to a detector of any size at an arbitrary depth. The gSeaGen code was not only extended in terms of functionalities but also underwent a thorough redesign of the muon propagation routine, resulting in a more accurate and efficient simulation. This paper presents the capabilities of the new gSeaGen code as well as prospects for further developments., Comment: 27 pages, 13 figures, submitted to Computer Physics Communications

Blazars are beamed active galactic nuclei known for their strong multi-wavelength variability on timescales from years down to minutes. We aim to investigate the suitability of the twisting jet model presented in previous works to explain the multi-wavelength behaviour of BL Lacertae, the prototype of one of the blazar classes. According to this model, the jet is inhomogeneous, curved, and twisting, and the long-term variability is due to changes in the Doppler factor due to variations in the orientation of the jet-emitting regions. We analysed optical data of the source obtained during monitoring campaigns organised by the Whole Earth Blazar Telescope (WEBT) in 2019-2022, together with radio data from the WEBT and other teams, and gamma-ray data from the Fermi satellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical and gamma-ray brightness maxima. The application of the twisting jet model to the source light curves allows us to infer the wiggling motion of the optical, radio, and gamma-ray jet-emitting regions. The optical-radio correlation shows that the changes in the radio viewing angle follow those in the optical viewing angle by about 120 days, and it suggests that the jet is composed of plasma filaments, which is in agreement with some radio high-resolution observations of other sources. The gamma-ray emitting region is found to be co-spatial with the optical one, and the analysis of the gamma-optical correlation is consistent with both the geometric interpretation and a synchrotron self-Compton (SSC) origin of the high-energy photons. We propose a geometric scenario where the jet is made up of a pair of emitting plasma filaments in a sort of double-helix curved rotating structure, whose wiggling motion produces changes in the Doppler beaming and can thus explain the observed multi-wavelength long-term variability., Comment: In press for A&A

Neutrinos described as an open quantum system may interact with the environment which introduces stochastic perturbations to their quantum phase. This mechanism leads to a loss of coherence along the propagation of the neutrino $-$ a phenomenon commonly referred to as decoherence $-$ and ultimately, to a modification of the oscillation probabilities. Fluctuations in space-time, as envisaged by various theories of quantum gravity, are a potential candidate for a decoherence-inducing environment. Consequently, the search for decoherence provides a rare opportunity to investigate quantum gravitational effects which are usually beyond the reach of current experiments. In this work, quantum decoherence effects are searched for in neutrino data collected by the KM3NeT/ORCA detector from January 2020 to November 2021. The analysis focuses on atmospheric neutrinos within the energy range of a few GeV to $100\,\mathrm{GeV}$. Adopting the open quantum system framework, decoherence is described in a phenomenological manner with the strength of the effect given by the parameters $\Gamma_{21}$ and $\Gamma_{31}$. Following previous studies, a dependence of the type $\Gamma_{ij} \propto (E/E_0)^n$ on the neutrino energy is assumed and the cases $n = -2,-1$ are explored. No significant deviation with respect to the standard oscillation hypothesis is observed. Therefore, $90\,\%$ CL upper limits are estimated as $\Gamma_{21} < 4.6\cdot 10^{-21}\,$GeV and $\Gamma_{31} < 8.4\cdot 10^{-21}\,$GeV for $n = -2$, and $\Gamma_{21} < 1.9\cdot 10^{-22}\,$GeV and $\Gamma_{31} < 2.7\cdot 10^{-22}\,$GeV for $n = -1$, respectively., Comment: 17 pages, 5 figures

KM3NeT/ORCA is a water Cherenkov neutrino detector under construction and anchored at the bottom of the Mediterranean Sea. The detector is designed to study oscillations of atmospheric neutrinos and determine the neutrino mass ordering. This paper focuses on an initial configuration of ORCA, referred to as ORCA6, which comprises six out of the foreseen 115 detection units of photo-sensors. A high-purity neutrino sample was extracted, corresponding to an exposure of 433 kton-years. The sample of 5828 neutrino candidates is analysed following a binned log-likelihood method in the reconstructed energy and cosine of the zenith angle. The atmospheric oscillation parameters are measured to be $\sin^2\theta_{23}= 0.51^{+0.04}_{-0.05}$, and $ \Delta m^2_{31} = 2.18^{+0.25}_{-0.35}\times 10^{-3}~\mathrm{eV^2} \cup \{-2.25,-1.76\}\times 10^{-3}~\mathrm{eV^2}$ at 68\% CL. The inverted neutrino mass ordering hypothesis is disfavoured with a p-value of 0.25., Comment: 29 pages, 12 figures

High-redshift blazars provide valuable input to studies of the evolution of active galactic nuclei (AGN) jets and provide constraints on cosmological models. Detections at high energies ($0.1<\mathrm{E}<100$ GeV) of these distant sources are rare, but when they exhibit bright gamma-ray flares, we are able to study them. However, contemporaneous multi-wavelength observations of high-redshift objects ($z>4$) during their different periods of activity have not been carried out so far. An excellent opportunity for such a study arose when the blazar TXS 1508+572 ($z=4.31$) exhibited a $\gamma$-ray flare in 2022 February in the $0.1-300$ GeV range with a flux 25 times brighter than the one reported in the in the fourth catalog of the \textit{Fermi} Large Area Telescope. Our goal is to monitor the morphological changes, spectral index and opacity variations that could be associated with the preceding $\gamma$-ray flare in TXS 1508+572 to find the origin of the high-energy emission in this source. We also plan to compare the source characteristics in the radio band to the blazars in the local Universe ($z<0.1$). In addition, we aim to collect quasi-simultaneous data to our multi-wavelength observations of the object, making TXS 1508+572 the first blazar in the early Universe ($z>4$) with contemporaneous multi-frequency data available in its high state. In order to study the parsec-scale structure of the source, we performed three epochs of very-long-baseline interferometry (VLBI) follow-up observations with the Very Long Baseline Array (VLBA) supplemented with the Effelsberg 100-m Telescope at 15, 22, and 43 GHz, which corresponds to 80, 117, and 228 GHz in the rest frame of TXS 1508+572. In addition, one 86 GHz (456 GHz) measurement was performed by the VLBA and the Green Bank Telescope during the first epoch.

Context : Understanding the mechanisms that launch and shape powerful relativistic jets from supermassive black holes (SMBHs) in high-redshift active galactic nuclei (AGN) is crucial for probing the co-evolution of SMBHs and galaxies over cosmic time. Aims :We study the high-redshift ($z=3.396$) blazar OH~471 to explore the jet launching mechanism in the early Universe. Methods : Using multi-frequency radio monitoring observations and high-resolution Very Long Baseline Interferometry imaging over three decades, we study the milliarcsecond structure and long-term variability of OH~471. Results : Spectral modelling of the radio flux densities reveals a synchrotron self-absorbed spectrum indicating strong magnetic fields within the compact core. By applying the flux freezing approximation, we estimate the magnetic flux carried by the jet and find that it reaches or exceeds theoretical predictions for jets powered by black hole spin energy via the Blandford-Znajek mechanism. This implies that OH~471 was in a magnetically arrested disk (MAD) state where the magnetic flux accumulated near the horizon regulates the accretion flow, allowing efficient extraction of black hole rotational energy. Conclusions : Our study demonstrates the dominance of MAD accretion in powering the prominent radio flares and relativistic jets observed in the radio-loud AGN OH~471 and statistical studies of large samples of high-redshift AGN will shed light on the role of MAD accretion in launching and accelerating the earliest relativistic jets., Comment: 16 pages, 7 figures, 3 tables

We present an optical-to-radio study of the BL Lac object S4 0954+658 observations during 1998-2023. The measurements were obtained with the SAO RAS Zeiss-1000 1-m and AS-500/2 0.5-m telescopes in 2003-2023, with the RATAN-600 radio telescope at 1.25 (0.96, 1.1), 2.3, 4.7 (3.7, 3.9), 8.2 (7.7), 11.2, 22.3 (21.7) GHz in 1998-2023, with the IAA RAS RT-32 Zelenchukskaya and Badary telescopes at 5.05 and 8.63 GHz in 2020--2023, and with the RT-22 single-dish telescope of CrAO RAS at 36.8 GHz in 2009-2023. In this period the blazar had been showing extremely high broadband activity with the variability amplitude of flux densities up to 70-100% both in the optical and radio domains. In the period of 2014-2023 the blazar had been showing the historically highest activity in the radio wavelengths, and we detected multiple radio flares of varying amplitude and duration. The large flares last on average from 0.3 to 1 year at 22-36.8 GHz and slightly longer at 5-11.2 GHz. The optical flares are shorter and last 7-50 days. In the most active epoch of 2018-2023 the characteristic time scale $\tau$ of variation at 5-22 GHz is about 100 days and about 1000 days for the state with lower activity in 2009-2014. We found a general correlation between the optical, radio, and $\gamma$-ray flux variations, which suggests that we observe the same photon population from different emission regions. We estimated linear size of this region as 0.5-2 pc for different epochs. A broadband two components radio spectrum of S4 0954+658 jet was modelled by using both electrons and protons as emitting particles. It is shown that the synchrotron radio waves in this AGN may be generated by relativistic protons., Comment: Published in Astrophysical Bulletin Vol.78, N4 (2023)

We conducted multi-epoch, multi-frequency parsec-scale studies on the gigahertz-peaked spectrum quasar PKS 0858-279 with the Very Long Baseline Array (VLBA). Our observations on 2005-11-26 elucidated a weak core, characterized by an inverted spectrum, and a distinctly bent jet that exhibited a notable bright feature in its Stokes I emission. Through comprehensive analysis of polarization and spectral data, we inferred the formation of a shock wave within this feature, stemming from interactions with a dense cloud in the ambient medium. In this paper, VLBI-Gaia astrometry further reinforces the core identification. With a deep analysis of six additional VLBA epochs spanning from 2007 to 2018, we observed that while the quasar's parsec-scale structure remained largely consistent, there were discernible flux density changes. These variations strongly imply the recurrent ejection of plasma into the jet. Complementing our VLBA data, RATAN-600 observations of the integrated spectra suggested an interaction between standing and travelling shock waves in 2005. Moreover, our multi-epoch polarization analysis revealed a drastic drop in rotation measure values from 6000 rad/m^2 to 1000 rad/m^2 within a single year, attributable to diminishing magnetic fields and particle density in an external cloud. This change is likely instigated by a shock in the cloud, triggered by the cloud's interaction with the jet, subsequently prompting its expansion. Notably, we also observed a significant change in the magnetic field direction of the jet, from being perpendicular post its observed bend to being perpendicular prior to the bend - an alteration possibly induced by the dynamics of shock waves., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

We present the results of our study of cross-correlations between long-term multi-band observations of the radio variability of the blazar 3C 279. More than a decade (2008-2022) of radio data were collected at seven different frequencies ranging from 2 GHz to 230 GHz. The multi-band radio light curves show variations in flux, with the prominent flare features appearing first at higher-frequency and later in lower-frequency bands. This behavior is quantified by cross-correlation analysis, which finds that the emission at lower-frequency bands lags that at higher-frequency bands. Lag versus frequency plots are well fit by straight lines with negative slope, typically ~-30 day/GHz. We discuss these flux variations in conjunction with the evolution of bright moving knots seen in multi-epoch VLBA maps to suggest possible physical changes in the jet that can explain the observational results. Some of the variations are consistent with the predictions of shock models, while others are better explained by a changing Doppler beaming factor as the knot trajectory bends slightly, given a small viewing angle to the jet., Comment: Submitted revised version to MNRAS journal, 11 pages, 6 figures, 4 tables

Baikal-GVD is a next generation, kilometer-scale neutrino telescope currently under construction in Lake Baikal. GVD consists of multi-megaton subarrays (clusters) and is designed for the detection of astrophysical neutrino fluxes at energies from a few TeV up to 100 PeV. The large detector volume and modular design of Baikal-GVD allows for the measurements of the astrophysical diffuse neutrino flux to be performed already at early phases of the array construction. We present here recent results of the measurements on the diffuse cosmic neutrino flux obtained with the Baikal-GVD neutrino telescope using cascade-like events., Comment: Presented at the 38th International Cosmic Ray Conference (ICRC 2023)

Active galaxies, especially blazars, are among the most promising neutrino source candidates. To date, ANTARES searches for these objects considered GeV-TeV $\gamma$-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\sigma$ 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\sigma$ indicates a $p=1.4\%$ ($2.5\sigma$ 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 $\gamma$-ray flares of the J0242+1101 blazar at a $p=0.5\%$ ($2.9\sigma$ 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.

We report on the observation of a rare neutrino event detected by Baikal-GVD in April 2021. The event GVD210418CA is the highest-energy cascade observed by Baikal-GVD so far from the direction below the horizon. The estimated cascade energy is $224\pm75$~TeV. The evaluated signalness parameter of GVD210418CA is 97.1\% using an assumption of the E$^{-2.46}$ spectrum of astrophysical neutrinos. The arrival direction of GVD210418CA is near the position of the well-known radio blazar TXS~0506+056, with the angular distance being within a 90\% directional uncertainty region of the Baikal-GVD measurement. The event was followed by a radio flare observed by the RATAN-600 radio telescope, further strengthening the case for the neutrino-blazar association., Comment: 9 pages, 6 figures, Contribution to the 38th International Cosmic Rays Conference (ICRC2023). arXiv admin note: text overlap with arXiv:2210.01650

We present full-track high-resolution radio observations of the jet of the galaxy M87 at 8 and 15 GHz. These observations were taken over three consecutive days in May 2009 using the Very Long Baseline Array (VLBA), one antenna of the Very Large Array (VLA), and the Effelsberg 100 m telescope. Our produced images have dynamic ranges exceeding 20,000:1 and resolve linear scales down to approximately 100 Schwarzschild radii, revealing a limb-brightened jet and a faint, steep spectrum counter-jet. We performed jet-to-counter-jet analysis, which helped estimate the physical parameters of the flow. The rich internal structure of the jet is dominated by three helical threads, likely produced by the Kelvin-Helmholtz (KH) instability developing in a supersonic flow with a Mach number of approximately 20 and an enthalpy ratio of around 0.3. We produce a CLEAN imaging bias-corrected 8-15GHz spectral index image, which shows spectrum flattening in regions of helical thread intersections. This further supports the KH origin of the observed internal structure of the jet. We detect polarised emission in the jet at distances of approximately 20 milliarcseconds from the core and find Faraday rotation which follows a transverse gradient across the jet. We apply Faraday rotation correction to the polarisation position angle and find that the position angle changes as a function of distance from the jet axis, which suggests the presence of a helical magnetic field., Comment: 15 pages, 10 figures, 1 table, MNRAS: Accepted 2023 October 03. Received 2023 September 27; in original form 2023 March 25

Baikal-GVD has recently published its first measurement of the diffuse astrophysical neutrino flux, performed using high-energy cascade-like events. We further explore the Baikal-GVD cascade dataset collected in 2018-2022, with the aim to identify possible associations between the Baikal-GVD neutrinos and known astrophysical sources. We leverage the relatively high angular resolution of the Baikal-GVD neutrino telescope (2-3 deg.), made possible by the use of liquid water as the detection medium, enabling the study of astrophysical point sources even with cascade events. We estimate the telescope's sensitivity in the cascade channel for high-energy astrophysical sources and refine our analysis prescriptions using Monte-Carlo simulations. We primarily focus on cascades with energies exceeding 100 TeV, which we employ to search for correlation with radio-bright blazars. Although the currently limited neutrino sample size provides no statistically significant effects, our analysis suggests a number of possible associations with both extragalactic and Galactic sources. Specifically, we present an analysis of an observed triplet of neutrino candidate events in the Galactic plane, focusing on its potential connection with certain Galactic sources, and discuss the coincidence of cascades with several bright and flaring blazars., Comment: 10 pages, 3 figures, accepted for publication in MNRAS

The advancement of neutrino observatories has sparked a surge in multi-messenger astronomy. Multiple neutrino associations among blazars are reported while neutrino production site is located within their central (sub)parsecs. Yet many questions remain on the nature of those processes. The next generation Event Horizon Telescope (ngEHT) is uniquely positioned for these studies, as its high frequency and resolution can probe both the accretion disk region and the parsec-scale jet. This opens up new opportunities for connecting the two regions and unraveling the proton acceleration and neutrino production in blazars. We outline observational strategies for ngEHT and highlight what it can contribute to the multi-messenger study of blazars., Comment: 9 pages; accepted to Galaxies Special Issue "From Vision to Instrument: Creating a Next-Generation Event Horizon Telescope for a New Era of Black Hole Science"

Bright blazars were found to be prominent neutrino sources, and a number of IceCube events were associated with them. Evaluating high-energy photon emission of such blazars is crucial for better understanding of the processes and regions where neutrinos are produced. Here, we focus on hard X-ray emission observed by the SRG/ART-XC telescope, by the Swift/BAT imager, and by the INTEGRAL/IBIS telescope. Their energy range ~10 keV is well-suited for probing photons that potentially participate in neutrino production by interacting with ultrarelativistic protons. We find that neutrino-associated blazars tend to demonstrate remarkably strong X-ray emission compared to other VLBI blazars in the sky. Both neutrinos and hard X-rays are found to come from blazars at cosmological distances z ~ 1, and are boosted by relativistic beaming that makes it possible to detect them on Earth. Our results suggest that neutrinos are produced within compact blazar jets, with target X-ray photons emitted from accelerated jet regions., Comment: 19 pages, 8 figures, 2 tables; published in JCAP

In the past years, evidence has started piling up that some high-energy cosmic neutrinos can be associated with blazars in flaring states. On February 26, 2022, a new blazar-neutrino coincidence has been reported: the track-like neutrino event IC220225A detected by IceCube is spatially coincident with the flat-spectrum radio quasar PKS 0215+015. Like previous associations, this source was found to be in a high optical and ${\gamma}$-ray state. Moreover, the source showed a bright radio outburst, which substantially increases the probability of a true physical association. We have performed six observations with the VLBA shortly after the neutrino event with a monthly cadence and are monitoring the source with the Effelsberg 100m-Telescope, and with the Australia Compact Telescope Array. Here, we present first results on the contemporary parsec-scale jet structure of PKS 0215+015 in total intensity and polarization to constrain possible physical processes leading to neutrino emission in blazars., Comment: 4 pages, 2 figures, IAU Symposium 375: The Multimessenger Chakra of Blazar Jets

Very long baseline interferometry (VLBI) observations show that some active galactic nuclei (AGN) jets exhibit bending even at parsec scales. The nature of bending is comprehensively analysed only for a small number of individual AGN, and the overall trends in shape of the substantially curved jets are unclear. In this work, we analyse outflows in AGN on the basis of publicly available multi-frequency VLBI images. Nearly 73000 images of about 11000 AGN are studied. Our research reveals that about 5% of them show a significantly curved jet structure. We characterize the jets geometry by fitting total intensity ridge lines constructed at all available frequencies and epochs with a set of simple models and suggest possible scenarios explaining the observed bending., Comment: 5 pages, 2 figures; presented at the 15th European VLBI Network Mini-Symposium and Users' Meeting; submitted to PoS

Using stacking of images obtained at different epochs, we studied the variability properties of linear polarization of active galactic nucleus (AGN) jets on parsec-scales. Our sample is drawn from the MOJAVE programme, and consists of 436 AGNs manifesting core-jet morphology and having at least five VLBA observing epochs at 15 GHz from January 1996 through August 2019, with some additional archival VLBA data reduced by us. We employed a stacking procedure and constructed maps of (i) standard deviation of fractional polarization and electric vector position angle (EVPA) over epochs as the measure of variability and (ii) median polarization degree to quantify typical values in time. The distributions of these values along and across the jet were analysed for the whole sample for the first time. We found that core EVPA variability is typically higher than that of the jet, presumably due to component blending and outflow bends in the core. The BL Lacertae object cores have lower EVPA variability, compared to that of quasars, possibly due to lower Faraday rotation measure, suggesting a stronger ordered magnetic field component. The EVPA becomes more stable down the jet. Most of the sources showing this trend have a time coverage of more than 12 years and at least 15 epochs. The possible cause could be the increase of stability in the magnetic field direction, reflecting an increase in the fraction of the magnetic field that is ordered. There are no significant optical-class-dependent or spectral-class-dependent relations in the EVPA variability properties in AGN jets., Comment: 15 pages, 11 figures, 5 tables; accepted by MNRAS; full electronic tables are available as ancillary files

Evidence for bright radio blazars being high-energy neutrino sources was found in recent years. However, specifics of how and where these particles get produced still need to be determined. In this paper, we add 14 new IceCube events from 2020-2022 to update our analysis of the neutrino-blazars connection. We test and refine earlier findings by utilising the total of 71 track-like high-energy IceCube events from 2009-2022. We correlate them with the complete sample of 3412 extragalactic radio sources selected by their compact radio emission. We demonstrate that neutrinos are statistically associated with radio-bright blazars with a post-trial p-value of 3*10^-4. In addition to this statistical study, we confirm previous individual neutrino-blazar associations, find and discuss several new ones. Notably, PKS 1741-038 was selected earlier and had a second neutrino detected from its direction in 2022; PKS 0735+168 has experienced a major flare across the whole electromagnetic spectrum coincidently with a neutrino arrival from that direction in 2021., Comment: 10 pages, 5 figures, 2 machine-readable tables; V2: minor changes and updates, accepted by MNRAS

The existence of high-energy astrophysical neutrinos has been unambiguously demonstrated, but their sources remain elusive. IceCube reported an association of a 290-TeV neutrino with a gamma-ray flare of TXS 0506+056, an active galactic nucleus with a compact radio jet pointing to us. Later, radio-bright blazars were shown to be associated with IceCube neutrino events with high statistical significance. These associations remained unconfirmed with the data of independent experiments. Here we report on the detection of a rare neutrino event with the estimated energy of 224+-75 TeV from the direction of TXS 0506+056 by the new Baikal Gigaton Volume Detector (Baikal-GVD) in April 2021. This event is the highest-energy cascade detected so far by the Baikal-GVD neutrino telescope from a direction below horizon. The result supports previous suggestions that radio blazars in general, and TXS 0506+056 in particular, are the sources of high-energy neutrinos, and opens up the cascade channel for the neutrino astronomy., Comment: 9 pages, 6 figures, 3 tables; accepted for publication in MNRAS

We analysed the parsec-scale linear polarization properties of 436 active galactic nuclei (AGN) based on 15 GHz polarimetric Very Long Baseline Array (VLBA) observations. We present polarization and total intensity images averaged over at least five epochs since 1996 January 19 through 2019 August 4. Stacking improves the image sensitivity down to $\sim$30 $\mu$Jy/beam and effectively fills out the jet cross-section both in total intensity and linear polarization. It delineates the long-term persistent magnetic field configuration and its regularity by restoring spatial distributions of the electric vector position angle (EVPA) and fractional polarization, respectively. On average, about ten years of stacking period is needed to reveal the stable and most-complete polarization distribution of a source. We find that the degree of polarization significantly increases down and across the jet towards its edges, typically manifesting U or W-shaped transverse profiles, suggesting a presence of a large-scale helical magnetic field associated with the outflow. In some AGN jets, mainly BL Lacs, we detect quasi-constant fractional polarization profiles across the jet, accompanied by EVPAs that closely follow the outflow. BL Lacs show higher fractional polarization values in their cores and jets than those in quasars up to hectoparsec de-projected scales, while on larger scales, they become comparable. High-synchrotron-peaked BL Lac jets are found to be less polarized than intermediate and low-synchrotron-peaked BL Lacs. The spatial distribution of the EVPAs in BL Lacs tend to align with the local jet direction, while quasars show an excess of orthogonal polarization orientation., Comment: 18 pages, 16 figures, 5 tables; accepted by MNRAS; full electronic tables, animated movie and all Figure Set plots are available as ancillary files

Astrophysical sources of neutrinos detected by large-scale neutrino telescopes remain uncertain. While there exist statistically significant observational indications that a part of the neutrino flux is produced by blazars, numerous theoretical studies suggest also the presence of potential Galactic point sources. Some of them have been observed in gamma rays above 100 TeV. Moreover, cosmic-ray interactions in the Galactic disk guarantee a diffuse neutrino flux. However, these Galactic neutrinos have not been unambiguously detected so far. Here we examine whether such a Galactic component is present among the observed neutrinos of the highest energies. We analyze public track-like IceCube events with estimated neutrino energies above 200 TeV. We examine the distribution of arrival directions of these neutrinos in the Galactic latitude b with the help of a simple unbinned, non-parametric test statistics, the median |b| over the sample. This distribution deviates from that implied by the null hypothesis of the neutrino flux isotropy, and is shifted towards lower b with the p-value of 4*10^{-5}, corresponding to the statistical significance of 4.1 sigma. There exists a significant component of the high-energy neutrino flux of Galactic origin, matching well the multi-messenger expectations from Tibet-ASgamma observations of diffuse Galactic gamma rays at hundreds of TeV. Together with the previously established extragalactic associations, the Galactic component we report here implies that the neutrino sky is rich and is composed of contributions from various classes of sources., Comment: 7 pages, 4 figures. V2: minor changes, accepted by ApJ Letters. Machine-readable list of events is available with this submission

The direction of parsec-scale jets in active galactic nuclei (AGNs) is essential information for many astrophysical and astrometric studies, including linear polarization and magnetic field structure, frequency-dependent synchrotron opacity, proper motion, and reference frame alignment. We developed a rigorous, simple, and completely automated method to measure the directions from calibrated interferometric visibility data at frequencies ranging from 1.4 GHz to 86 GHz. We publish the results for 9220 AGNs with the typical accuracy below 10 degrees. An internal check of the method comparing the directions between different observing frequencies as well as with previous publications verifies the robustness of the measured values., Comment: 11 pages, 10 figures, 2 tables; machine-readable tables available as ancillary files; accepted to ApJS

We used archival very long baseline interferometry (VLBI) data of active galactic nuclei (AGN) observed from 1.4 GHz to 86 GHz to measure the angular size of VLBI radio cores in 8959 AGNs. We analysed their sky distributions, frequency dependencies and created the most densely sampled and complete to date distribution map of large-scale scattering properties of the interstellar medium in our Galaxy. Significant angular broadening of the measured AGN core sizes is detected for the sources seen through the Galactic plane, and this effect is especially strong at low frequencies (e.g., at 2 GHz). The scattering screens containing electron density fluctuations of hot plasma are mainly concentrated in the Galactic plane and manifest clumpy distribution. The region of the strongest scattering is the Galactic centre, where the Galactic bar and the compact radio source Sagittarius A* are located. We have also found the enhancement of scattering strength in regions of the Cygnus constellation, supernova remnants Taurus A, Vela, W78 and Cassiopeia A, and the Orion Nebula. Using multi-frequency observational data of AGN core sizes, we separated the contribution of the intrinsic and scattered sizes to the measured angular diameter for 1411 sources. For the sources observed through the Galactic plane, the contribution of the scattered size component is systematically larger than for those seen outside the Galactic plane. The derived power-law scattering indices are found to be in good agreement with theoretical predictions for the diffractive-dominated scattering of radio emission in a hot plasma with Gaussian distribution of density inhomogeneities., Comment: 15 pages, 14 figures, 7 tables, accepted to MNRAS

High-energy astrophysical neutrinos have been observed by several telescopes in the last decade, but their sources still remained unknown. We address the problem of locating astrophysical neutrinos' sources in a statistical manner. We show that blazars positionally associated with IceCube neutrino detections have stronger parsec-scale radio cores than the rest of the sample. The probability of a chance coincidence is only 4*10^-5 corresponding to a significance of 4.1 sigma. We explicitly list four strong radio blazars as highly probable sources of neutrinos above 200 TeV: 3C 279, NRAO 530, PKS 1741-038, and PKS 2145+067. There are at least 70 more radio-bright blazars that emit neutrinos of lower energies, starting from TeVs. Using continuous RATAN-600 monitoring of VLBI-selected blazars, we find that radio flares at frequencies above 10 GHz coincide with neutrino arrival dates. The most pronounced example of such behavior is PKS 1502+106 that experienced a major flare in 2019. We demonstrate that the majority of IceCube astrophysical neutrino flux derived from muon-track analyses may be explained by blazars, that is, AGNs with bright Doppler-boosted jets. High-energy neutrinos can be produced in photohadronic interactions within parsec-scale relativistic jets. Radio-bright blazars associated with neutrino detections have very diverse gamma-ray properties, which suggests that gamma-rays and neutrinos may be produced in different regions of blazars and not directly related. A narrow jet viewing angle is, however, required to detect either of them. We conclude with discussion of recent independent tests and extensions of our findings., Comment: 8 pages, 5 figures; published in ICRC'2021 proceedings

We present RadioAstron space-based very long baseline interferometry (VLBI) observations of the nearby radio galaxy 3C84 (NGC1275) at the centre of the Perseus cluster. The observations were carried out on September 21-22, 2013 and involved a global array of 24 ground radio telescopes observing at 5 GHz and 22 GHz, together with the Space Radio Telescope (SRT). Furthermore, the Very Long Baseline Array (VLBA) and the phased Very Large Array (VLA) observed the source quasi-simultaneously at 15 GHz and 43 GHz. Fringes between the ground array and the SRT were detected on baseline lengths up to 8.1 times the Earth's diameter, providing unprecedented resolution for 3C 84 at these wavelengths. We note that the corresponding fringe spacing is 125 microarcsec at 5 GHz and 27 microarcsec at 22 GHz. Our space-VLBI images reveal a previously unseen sub-structure inside the compact 1 pc long jet that was ejected about ten years earlier. In the 5 GHz image, we detected, for the first time, low-intensity emission from a cocoon-like structure around the restarted jet. Our results suggest that the increased power of the young jet is inflating a bubble of hot plasma as it carves its way through the ambient medium of the central region of the galaxy. Here, we estimate the minimum energy stored in the mini-cocoon, along with its pressure, volume, expansion speed, and the ratio of heavy particles to relativistic electrons, as well as the density of the ambient medium. About half of the energy delivered by the jet is dumped into the mini-cocoon and the quasi-spherical shape of the bubble suggests that this energy may be transferred to a significantly larger volume of the interstellar medium than what would be accomplished by the well-collimated jet on its own. The pressure of the hot mini-cocoon also provides a natural explanation for the almost cylindrical jet profile seen in the 22 GHz RadioAstron image., Comment: To be published in Astronomy & Astrophysics

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., Comment: 29 pages, 16 figures, 4 tables; accepted by ApJ; full electronic tables and all Figure Set plots are available as ancillary files

We have analyzed the parsec-scale jet kinematics of 447 bright radio-loud AGN, based on 15 GHz VLBA data obtained between 1994 August 31 and 2019 August 4. We present new total intensity and linear polarization maps obtained between 2017 January 1 to 2019 August 4 for 143 of these AGN. We tracked 1923 bright features for five or more epochs in 419 jets. A majority (60%) of the well-sampled jet features show either accelerated or non-radial motion. In 47 jets there is at least one non-accelerating feature with an unusually slow apparent speed. Most of the jets show variations of 10 to 50 deg in their inner jet position angle (PA) over time, although the overall distribution has a continuous tail out to 200 deg. AGN with SEDs peaked at lower frequencies tend to have more variable PAs, with BL Lacs being less variable than quasars. The Fermi LAT gamma-ray associated AGN also tend to have more variable PAs than the non-LAT AGN in our sample. We attribute these trends to smaller viewing angles for the lower spectral peaked and LAT-associated jets. We identified 13 AGN where multiple features emerge over decade-long periods at systematically increasing or decreasing PAs. Since the ejected features do not fill the entire jet cross-section, this behavior is indicative of a precessing flow instability near the jet base. Although some jets show indications of oscillatory PA evolution, we claim no bona fide cases of periodicity since the fitted periods are comparable to the total VLBA time coverage., Comment: 25 pages, 16 figures, 7 tables, accepted by ApJ; full electronic tables are provided as arXiv ancillary files; figure sets are made available from the full paper source https://www.physics.purdue.edu/astro/MOJAVE/paperXVIII.tar. This article is a follow-up to the previous MOJAVE survey paper on blazar jet kinematics (arXiv:1902.09591), in which the same data gathering and calibration procedures were used

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 $\gamma$-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 $\gamma$-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 $\gamma$-ray emission region from the central engine and shorter variability time scale compared to those of BL Lacertae objects., Comment: 12 pages, 12 figures, accepted for publication in MNRAS

We performed multi-frequency studies on the gigahertz-peaked spectrum high-redshift quasar 0858-279. Initially, the source presented itself in the early VLBI images as a very peculiar resolved blob. We observed the quasar with the VLBA at 1.4-24 GHz in a dual-polarization mode. The high spatial resolution and the spectral index maps enabled us to resolve the core-jet structure and locate a weak and compact core by its inverted spectrum. The dominant jet component 20 parsecs away from the core was optically thin above 10 GHz and opaque below it. We also estimated an uncommonly strong magnetic field in the bright jet feature, which turned out to be around 1 G. The Faraday rotation measure maps revealed high RM values over 6000 rad/m^2. Additionally, these maps allowed us to follow the magnetic field direction in the bright jet feature being perpendicular to the propagation direction of the jet. All the results strongly indicated the formation of a shock wave in the dominant component arising from an interaction with the surrounding matter. Using the proposed hypothesis and the core shift approach, we discovered that the magnetic field in the core region is of the order of 0.1 G., Comment: 15 pages, 12 figures, accepted for publication in MNRAS

In recent studies, several AGN have exhibited gradients of the Faraday Rotation Measure (RM) transverse to their parsec-scale jet direction. Faraday rotation likely occurs as a result of a magnetized sheath wrapped around the jet. In the case of 3C 273, using Very Long Baseline Array multi-epoch observations at 5, 8 and 15 GHz in 2009--2010, we observe that the jet RM has changed significantly towards negative values compared with that previously observed. These changes could be explained by a swing of the parsec-scale jet direction which causes synchrotron emission to pass through different portions of the Faraday screen. We develop a model for the jet-sheath system in 3C 273 where the sheath is wider than the single-epoch narrow relativistic jet. We present our oversized sheath model together with a derived wide jet full intrinsic opening angle $\alpha_\mathrm{int}=2.1^\circ$ and magnetic field strength $B_{||}=3$ $\mu$G and thermal particle density $N_\mathrm{e}=125~\mathrm{cm}^{-3}$ at the wide jet--sheath boundary 230 pc downstream (deprojected) from its beginning. Most of the Faraday rotation occurs within the innermost layers of the sheath. The model brings together the jet direction swing and long-term RM evolution and may be applicable to other AGN jets that exhibit changes of their apparent jet direction., Comment: 10 pages, 9 figures. Accepted to the Astrophysical Journal

Context. Supermassive black holes in the centres of radio-loud active galactic nuclei (AGN) can produce collimated relativistic outflows (jets). Magnetic fields are thought to play a key role in the formation and collimation of these jets, but the details are much debated. Aims. We study the innermost jet morphology and magnetic field strength in the AGN 3C 345 with an unprecedented resolution using images obtained within the framework of the key science programme on AGN polarisation of the Space VLBI mission RadioAstron. Methods. We observed the flat spectrum radio quasar 3C 345 at 1.6 GHz on 2016 March 30 with RadioAstron and 18 ground-based radio telescopes in full polarisation mode. Results. Our images, in both total intensity and linear polarisation, reveal a complex jet structure at 300 $\mu$as angular resolution, corresponding to a projected linear scale of about 2 pc or a few thousand gravitational radii. We identify the synchrotron self-absorbed core at the jet base and find the brightest feature in the jet 1.5 mas downstream of the core. Several polarised components appear in the Space VLBI images that cannot be seen from ground array-only images. Except for the core, the electric vector position angles follow the local jet direction, suggesting a magnetic field perpendicular to the jet. This indicates the presence of plane perpendicular shocks in these regions. Additionally, we infer a minimum brightness temperature at the largest $(u,v)$-distances of $1.1\times 10^{12}$ K in the source frame, which is above the inverse Compton limit and an order of magnitude larger than the equipartition value. This indicates locally efficient injection or re-acceleration of particles in the jet to counter the inverse Compton cooling or the geometry of the jet creates significant changes in the Doppler factor, which has to be $>11$ to explain the high brightness temperatures., Comment: 10 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

We present Space-VLBI RadioAstron observations at 1.6 GHz and 4.8 GHz of the flat spectrum radio quasar 3C 273, with detections on baselines up to 4.5 and 3.3 Earth Diameters, respectively. Achieving the best angular resolution at 1.6 GHz to date, we have imaged limb-brightening in the jet, not previously detected in this source. In contrast, at 4.8 GHz, we detected emission from a central stream of plasma, with a spatial distribution complementary to the limb-brightened emission, indicating an origin in the spine of the jet. While a stratification across the jet width in the flow density, internal energy, magnetic field, or bulk flow velocity are usually invoked to explain the limb-brightening, the different jet structure detected at the two frequencies probably requires a stratification in the emitting electron energy distribution. Future dedicated numerical simulations will allow the determination of which combination of physical parameters are needed to reproduce the spine/sheath structure observed by Space-VLBI with RadioAstron in 3C 273, Comment: Accepted by A&A

The Blandford and K\"{o}nigl model of AGN jets predicts that the position of the apparent opaque jet base - the core - changes with frequency. This effect is observed with radio interferometry and is widely used to infer parameters and structure of the innermost jet regions. The position of the radio core is typically estimated by fitting a Gaussian template to the interferometric visibilities. This results in a model approximation error, i.e. a bias that can be detected and evaluated through simulations of observations with a realistic jet model. To assess the bias, we construct an artificial sample of sources based on the AGN jet model evaluated on a grid of the parameters derived from a real VLBI flux-density-limited sample and create simulated VLBI data sets at 2.3, 8.1 and 15.4 GHz. We found that the core position shifts from the true jet apex are generally overestimated. The bias is typically comparable to the core shift random error and can reach a factor of two for jets with large apparent opening angles. This observational bias depends mostly on the ratio between the true core shift and the image resolution. This implies that the magnetic field, the core radial distance and the jet speed inferred from the core shift measurements are overestimated. We present a method to account for the bias., Comment: 11 pages, 9 figures. Published in MNRAS

Recently we have shown that high-energy neutrinos above 200 TeV detected by IceCube are produced within several parsecs in the central regions of radio-bright blazars, that is active galactic nuclei with jets pointing towards us. To independently test this result and extend the analysis to a wider energy range, we use public data for all neutrino energies from seven years of IceCube observations. The IceCube point-source likelihood map is analyzed against the positions of blazars from a statistically complete sample selected by their compact radio flux density. The latter analysis delivers a 3.0 sigma significance with the combined post-trial significance of both studies being 4.1 sigma. The correlation is driven by a large number of blazars. Together with fainter but physically similar sources not included in the sample, they may explain the entire IceCube astrophysical neutrino flux as derived from muon-track analyses. The neutrinos can be produced in interactions of relativistic protons with X-ray self-Compton photons in parsec-scale blazar jets., Comment: 14 pages, 5 figures; accepted to ApJ; v2: extended and clarified discussion, added a diagram; v3: improved numerical estimates

We use the observed jet boundary transition from parabolic to conical shape, which was earlier discovered as possibly a common effect in active galactic nuclei, to estimate a black hole, a jet and an ambient medium parameters. We explained earlier the geometry transition as a consequence of a change in the jet inner properties: a transition from a magnetically dominated to an equipartition regime. This interpretation allows us to estimate a black hole spin, a black hole mass and an ambient pressure amplitude, using the observed jet shape break position and the jet width at the transition point for 11 active galactic nuclei. The black hole spin values obtained using our method are consistent with the lower estimates for the sources with redshift $z<2$ from the spin evolution modelling. We find that the method of a black hole mass determination based on the relation between the broad-line region size and its luminosity may underestimate masses of the sources with large jet viewing angles. We propose a new method for the black hole mass determination, with the obtained masses being in interval $10^8-10^{10}\;M_{\odot}$. The range of the values of the ambient pressure amplitude points to the uniform medium conditions for the sources in our sample, with a tentative indication of higher pressure around FRII sources., Comment: 13 pages, 6 figures, 3 tables; accepted by MNRAS

We have carried out a Chandra X-ray and multi-frequency radio VLBA study of the AGN TXS 0128+554, which is associated with the Fermi gamma-ray source 4FGL J0131.2+5547. The AGN is unresolved in a target 19.3 ks Chandra image, and its spectrum is well fit by a simple absorbed power law model, with no distinguishable spectral features. Its relatively soft X-ray spectrum compared to other CSOs may be indicative of a thermal emission component, for which we were able to obtain an upper temperature limit of kT = 0.08 keV. The compact radio morphology and measured advance speed of 0.32c +- 0.07c indicate a kinematic age of only 82 y +- 17 y, placing TXS 0128+554 among the youngest members of the compact symmetric object (CSO) class. The lack of compact, inverted spectrum hotspots and an emission gap between the bright inner jet and outer radio lobe structure indicate that the jets have undergone episodic activity, and were re-launched a decade ago. The predicted gamma-ray emission from the lobes, based on an inverse Compton-emitting cocoon model, is three orders of magnitude below the observed Fermi LAT flux. A comparison to other Fermi-detected and non-Fermi detected CSOs with redshift z<0.1 indicates that the gamma-ray emission likely originates in the inner jet/core region, and that nearby, recently launched AGN jets are primary candidates for detection by the Fermi LAT instrument., Comment: 22 pages, 9 figures, 7 tables; accepted by ApJ; full electronic table 2 is available from the preprint source

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., Comment: Accepted for publication in the Astrophysical Journal

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet., Comment: 21 pages, 18 figures, accepted for publication in MNRAS

Significant positional offsets of the value from 1 mas to more than 10 mas were found previously between radio (VLBI) and optical (Gaia) positions of active galactic nuclei (AGN). They happen preferentially parallel to the parsec-scale jet direction. AGN with VLBI-to-Gaia offsets pointed downstream the jet are found to have favourably higher optical polarization, as expected if extended optical jets dominate in the emission and shift the Gaia centroid away from the physical nucleus of the source. Upstream offsets with the suggested domination of accretion disks manifest themselves through the observed low optical polarization. Direction of linear optical polarization is confirmed to preferentially align with parsec-scale jets in AGN with dominant jets consistent with a toroidal magnetic field structure. Our findings support the disk-jet interpretation of the observed positional offsets. These results call on an intensification of AGN optical polarization monitoring programs in order to collect precious observational data. Taken together with the continued VLBI and Gaia observations, they will allow researchers to reconstruct detailed models of the disk-jet system in AGN on parsec scales., Comment: 5 pages, 4 figures, 2 tables; accepted by MNRAS Letters; table 2 is electronic only - available from the preprint source; minor updates are introduced

Observational information on high-energy astrophysical neutrinos is being continuously collected by the IceCube observatory. However, the sources of neutrinos are still unknown. In this study, we use radio very-long-baseline interferometry (VLBI) data for a complete VLBI-flux-density limited sample of active galactic nuclei (AGN). We address the problem of the origin of astrophysical neutrinos with energies above 200 TeV in a statistical manner. It is found that AGN positionally associated with IceCube events have typically stronger parsec-scale cores than the rest of the sample. The post-trial probability of a chance coincidence is 0.2%. We select the four strongest AGN as highly probable associations: 3C 279, NRAO 530, PKS 1741-038, and PKS 2145+067. Moreover, we find an increase of radio emission at frequencies above 10 GHz around neutrino arrival times for several other VLBI-selected AGN on the basis of RATAN-600 monitoring. The most pronounced example of such behavior is PKS 1502+106. We conclude that AGN with bright Doppler-boosted jets constitute an important population of neutrino sources. High-energy neutrinos are produced in their central parsec-scale regions, probably in proton-photon interactions at or around the accretion disk. Radio-bright AGN that are likely associated with neutrinos have very diverse gamma-ray properties suggesting that gamma-rays and neutrinos may be produced in different regions of AGN and not directly related. A small viewing angle of the jet-disk axis is, however, required to detect either of them., Comment: 19 pages, 5 figures, 3 tables; accepted to ApJ; v2: extended discussion, added a figure and electronic table; v3: fixed typo in a source name

Detailed studies of relativistic jets in active galactic nuclei (AGN) require high-fidelity imaging at the highest possible resolution. This can be achieved using very long baseline interferometry (VLBI) at radio frequencies, combining worldwide (global) VLBI arrays of radio telescopes with a space-borne antenna on board a satellite. We present multiwavelength images made of the radio emission in the powerful quasar S5 0836+710, obtained using a global VLBI array and the antenna Spektr-R of the RadioAstron mission of the Russian Space Agency, with the goal of studying the internal structure and physics of the relativistic jet in this object. The RadioAstron observations at wavelengths of 18cm, 6cm, and 1.3cm are part of the Key Science Program for imaging radio emission in strong AGN. The internal structure of the jet is studied by analyzing transverse intensity profiles and modeling the structural patterns developing in the flow. The RadioAstron images reveal a wealth of structural detail in the jet of S5 0836+710 on angular scales ranging from 0.02mas to 200mas. Brightness temperatures in excess of $10^{13}$\,K are measured in the jet, requiring Doppler factors of $\ge 100$ for reconciling them with the inverse Compton limit. Several oscillatory patterns are identified in the ridge line of the jet and can be explained in terms of the Kelvin-Helmholtz (KH) instability. The oscillatory patterns are interpreted as the surface and body wavelengths of the helical mode of the KH instability. The interpretation provides estimates of the jet Mach number and of the ratio of the jet to the ambient density, which are found to be $M_\mathrm{j}\approx 12$ and $\eta\approx 0.33$. The ratio of the jet to the ambient density should be conservatively considered an upper limit because its estimate relies on approximations., Comment: 16 pages, 9 figures (9 extra figures in Appendix), in press; Astronomy and Astrophysics