Search

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

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

We analyze the variability of the parsec-scale jet directions in active galactic nuclei (AGNs). Our analysis involves 317 AGNs at frequencies ranging from 2 to 43 GHz, and is made possible by developing an automatic inner jet direction measurement procedure. We find strong significant variations in a one quarter of these AGNs; the effect is likely ubiquitous, and not detected in the rest due to a limited sensitivity and observations epoch coverage. Average apparent jet rotation speeds range from 0.21 deg/yr at 2 GHz to 1.04 deg/yr at 43 GHz. This strong frequency dependence indicates that the variability cannot be explained by jet components propagating ballistically without acceleration: more complex jet shapes or patterns are required. Still, we demonstrate that the apparent direction changes are predominantly caused by the jet nozzle rotations, and not by individual components propagating transversely to the jet. In this work, we focus on variability scales much longer than the times of observations, that is > 50 years. Using our measurements, we bound potential periods to less than 1000 years in the source rest frame for 90% AGNs in the sample. This allows us to constrain mechanisms causing these variations if they are periodic, such as instabilities, disk-driven precession, or binary black hole effects., Comment: 13 pages, 11 figures, 2 tables, submitted to MNRAS

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

We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery will expose universal features of a black hole's spacetime that are distinct from the complex astrophysics of the emitting plasma, allowing the first direct measurements of a supermassive black hole's spin. In addition to studying the properties of the nearby supermassive black holes M87* and Sgr A*, BHEX will measure the properties of dozens of additional supermassive black holes, providing crucial insights into the processes that drive their creation and growth. BHEX will also connect these supermassive black holes to their relativistic jets, elucidating the power source for the brightest and most efficient engines in the universe. BHEX will address fundamental open questions in the physics and astrophysics of black holes that cannot be answered without submillimeter space VLBI. The mission is enabled by recent technological breakthroughs, including the development of ultra-high-speed downlink using laser communications, and it leverages billions of dollars of existing ground infrastructure. We present the motivation for BHEX, its science goals and associated requirements, and the pathway to launch within the next decade., Comment: Proceedings for SPIE Astronomical Telescopes and Instrumentation

We present a baseline science operations plan for the Black Hole Explorer (BHEX), a space mission concept aiming to confirm the existence of the predicted sharp ``photon ring" resulting from strongly lensed photon trajectories around black holes, as predicted by general relativity, and to measure its size and shape to determine the black hole's spin. BHEX will co-observe with a ground-based very long baseline interferometric (VLBI) array at high-frequency radio wavelengths, providing unprecedented high resolution with the extension to space that will enable photon ring detection and studies of active galactic nuclei. Science operations require a simultaneous coordination between BHEX and a ground array of large and small radio apertures to provide opportunities for surveys and imaging of radio sources, while coordination with a growing network of optical downlink terminals provides the data rates necessary to build sensitivity on long baselines to space. Here we outline the concept of operations for the hybrid observatory, the available observing modes, the observation planning process, and data delivery to achieve the mission goals and meet mission requirements., Comment: 12 pages, 10 figures, submitted to SPIE Astronomical Telescopes + Instrumentation (AS24)

We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251\deg $\le l \le$ 358\deg and 2\deg $\le l \le$ 61\deg at $|b| \le 1.5\deg $). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $\mu$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance., Comment: Accepted for publication in MNRAS. The data release is live and links can be found in the Data Availability Statement in the paper

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

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

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

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

Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies., Comment: The paper on 29 pages contains 12 figures, accepted for publication in Acta Astronautica, based on the presentation at the 72nd IAC, Dubai, UAE, October 2021

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 study milliarcsecond-scale properties of the persistent radio counterpart to FRB 121102 and investigate the spectro-polarimetric properties of a bright burst. For the former, we use European VLBI Network (EVN) observations in 2017 at 1.7 and 4.8 GHz. For the latter, we re-analyse the 1.7-GHz data from the 100-m Effelseberg telescope taken in 2016. These observations predate other polarimetric studies of FRB 121102, and yield the highest burst Faraday rotation measure (RM) to date, RM = 1.27*10^5 rad m^-2, consistent with the decreasing RM trend. The fractional polarization of the burst emission is 15% at 1.7 GHz. This can be reconciled with the high fractional polarization at higher frequencies if the Faraday width of the burst environment is 150 rad m^-2 - a bare 0.1% of the total Faraday rotation. The width may originate from minor non-uniformities in the Faraday screen, or from effects in the emitting region itself. The upper limit on the persistent source size is 1 pc, barely consistent with a young supernova (SN) scenario. The flux variability limit of <10% is not in favor of the young SN scenario, and challenges other interpretations as well. The fractional polarization of the faint persistent source is constrained at <25% at 4.8 GHz ruling out a common origin with the highly polarized individual bursts., Comment: 9 pages, 8 figures; accepted to MNRAS

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

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

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

Observational studies of collimation in jets in active galactic nuclei (AGN) are a key to understanding their formation and acceleration processes. We have performed an automated search for jet shape transitions in a sample of 367 AGN using VLBA data at 15 GHz and 1.4 GHz. This search has found ten out of 29 nearby jets at redshifts z<0.07 with a transition from a parabolic to conical shape, while the full analyzed sample is dominated by distant AGN with a typical z about 1. The ten AGN are UGC00773, NGC1052, 3C111, 3C120, TXS0815-094, Mrk180, PKS1514+00, NGC6251, 3C371, and BL Lac. We conclude that the geometry transition may be a common effect in AGN jets. It can be observed only when sufficient linear resolution is obtained. Supplementing these results with previously reported shape breaks in the nearby AGN 1H0323+342 and M87, we estimate that the break occurs at 10^5-10^6 gravitational radii from the nucleus. We suggest that the jet shape transition happens when the bulk plasma kinetic energy flux becomes equal to the Poynting energy flux, while the ambient medium pressure is assumed to be governed by Bondi accretion. In general, the break point may not coincide with the Bondi radius. The observational data supports our model predictions on the jet acceleration and properties of the break point., Comment: 18 pages; electronic table 3 is available from the preprint source; a typo in the abstract is corrected; accepted by MNRAS

Physical parameters of AGN jets observed with Very Long Baseline Interferometry (VLBI) are usually inferred from the core shift measurements or flux and size measured at a peak frequency of the synchrotron spectrum. Both are preceded by modelling of the observed VLBI jet structure with a simple Gaussian templates. We propose to infer the jets parameters using the inhomogeneous jet model directly - bypassing the modelling of the source structure with a Gaussian templates or image deconvolution. We applied Bayesian analysis to multi-frequency VLBA observations of radio galaxy NGC 315 and found that its parsec-scale jet is well described by the inhomogeneous conical model. Our results favour electron-positron jet. We also detected a component in a counter jet. Its position implies the presence of an external absorber with a steep density gradient at close ($r=0.1$ pc) distance from the central engine., Comment: 15 pages, 14 figures. Revised version published in MNRAS

The apparent position of jet base (core) in radio-loud active galactic nuclei changes with frequency because of synchrotron self-absorption. Studying this `core shift` effect enables us to reconstruct properties of the jet regions close to the central engine. We report here results from core shift measurements in AGNs observed with global VLBI at 2 and 8 GHz at epochs from 1994 to 2016. Our sample contains 40 objects observed at least 10 times during that period. The core shift is determined using a new automatic procedure introduced to minimize possible biases. The resulting multiple epoch measurements of the core position are employed for examining temporal variability of the core shift. We argue that the core shift variability is a common phenomenon, as established for 33 of 40 AGNs we study. Our analysis shows that the typical offsets between the core positions at 2 and 8 GHz are about 0.5 mas and they vary in time. Typical variability of the individual core positions is about 0.3 mas. The measurements show a strong dependence between the core position and its flux density, suggesting that changes in both are likely related to the nuclear flares injecting denser plasma into the flow. We determine that density of emitting relativistic particles significantly increases during these flares, while relative magnetic field changes less and in the opposite direction., Comment: 22 pages, 15 figures, 2 tables; accepted to MNRAS

Results of a pilot project with the participation of the "Kvazar-KVO" radio interferometry array in observations carried out with the European VLBI Network are presented. The aim of the project was to conduct and analyze multi-frequency (1.7, 2.3, 5.0, 8.4 GHz) observations of the parsec-scale jets of 24 active galactic nuclei. Three observing sessions were successfully carried out in October 2008. Maps of the radio intensity distributions have been constructed in all four frequencies using phase referencing. A method for measuring the frequency-dependent shift of the position of the VLBI core by applying relative astrometry to observations of close triplets of radio sources has been developed. The fundamental possibility of detecting core shifts in ultra-compact sources for which traditional methods based on the achromatic positions of optically thin regions of the jet are not suitable is demonstrated. The conditions for successful measurement of this shift are discussed; these are determined by the closeness of the calibrator used, the effective resolution of the system, the quality of the filling of the $uv$ plane, the relative orientations of the jets in the triplets, and the brightnesses of the sources., Comment: Accepted for publication in Astronomy Reports; 24 pages, 9 figures, 3 tables

We analyze differences in positions of active galactic nuclei between Gaia data release 2 and VLBI and compare the significant VLBI-to-Gaia offsets in more than 1000 objects with their jet directions. Remarkably at least 3/4 of the significant offsets are confirmed to occur downstream or upstream the jet representing a genuine astrophysical effect. Introducing redshift and Gaia color into analysis can help distinguish between the contribution of the host galaxy, jet, and accretion disk emission. We find that strong optical jet emission at least 20-50pc long is required to explain the Gaia positions located downstream from VLBI ones. Offsets in the upstream direction of up to 2 mas are at least partly due to the dominant impact of the accretion disk on the Gaia coordinates and by the effects of parsec-scale radio jet. The host galaxy was found not to play an important role in the detected offsets. BL Lacertae object and Seyfert 2 galaxies are observationally confirmed to have a relatively weak disk and consequently downstream offsets. The disk emission drives upstream offsets in a significant fraction of quasars and Seyfert 1 galaxies when it dominates over the jet in the optical band. The observed behaviour of the different AGN classes is consistent with the unified scheme assuming varying contribution of the obscuring dusty torus and jet beaming., Comment: 7 pages, 6 figures, 1 table; accepted to ApJ; manuscript expanded from letter to article format

We have analyzed the differences in positions of 9081 matched sources between the Gaia DR2 and VLBI catalogues. The median position uncertainty of matched sources in the VLBI catalogue is a factor of two larger than the median position uncertainty in the Gaia DR2. There are 9% matched sources with statistically significant offsets between both catalogues. We found that reported positional errors should be re-scaled by a factor of 1.3 for VLBI and 1.06 for Gaia, and in addition, Gaia errors should be multiplied by the square root of chi square per degree of freedom in order to best fit the normalized position differences to the Rayleigh distribution. We have established that the major contributor to statistically significant position offsets is the presence of optical jets. Among the sources for which the jet direction was determined, the position offsets are parallel to the jet directions for 62% of the outliers. Among the matched sources with significant proper motion, the fraction of objects with proper motion directions parallel to jets is a factor of 3 greater than on average. Such sources have systematically higher chi square per degree of freedom. We explain these proper motions as a manifestation of the source position jitter caused by flares that we have predicted earlier. Therefore, the assumption that quasars are fixed points and therefore, differential proper motions determined with respect to quasar photocenters can be regarded as absolute proper motions, should be treated with a great caution., Comment: 9 pages, 8 figures; accepted by MNRAS

We present RadioAstron space-based very long baseline interferometry (VLBI) observations of the nearby radio galaxy 3C 84 (NGC 1275) at the centre of the Perseus cluster. The observations were carried out during a perigee passage of the Spektr-R spacecraft 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 μas at 5 GHz and 27 μas 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.