Motivated by the identification of the blazar TXS 0506+056 as the first promising high-energy neutrino counterpart candidate, we search for additional neutrino blazars candidates among the Fermi-LAT detected blazars. We investigate the multi-wavelength behavior from radio to GeV gamma rays of blazars found to be in spatial coincidence with single high-energy neutrinos and lower-energy neutrino flare candidates. In addition, we compare the average gamma-ray emission of the potential neutrino-emitting sources to the entire sample of gamma-ray blazars. We find that neutrino-emitting blazar candidates are statistically compatible with both hypothesis of a linear correlation and of no correlation between neutrino and gamma-ray energy flux.

The $\gamma$-ray production mechanism and its localization in blazars are still a matter of debate. The main goal of this paper is to constrain the location of the high-energy emission in the blazar TXS 2013+370 and to study the physical and geometrical properties of the inner jet region on sub-pc scales. VLBI observations at 86 GHz and space-VLBI at 22 GHz allowed us to image the jet base with an angular resolution of $\sim$0.4 pc. By employing CLEAN imaging and Gaussian model-fitting, we performed a thorough kinematic analysis, which provided estimates of the jet speed, orientation, and component ejection times. Additionally, we studied the jet expansion profile and used the information on the jet geometry to estimate the location of the jet apex. VLBI data were combined with single-dish measurements to search for correlated activity between the radio and $\gamma$-ray emission. The high-resolution VLBI imaging revealed the existence of a spatially bent jet, described by moving and stationary features. New jet features are observed to emerge from the core, accompanied by flaring activity in radio bands and $\gamma$ rays. The analysis of the transverse jet width profile constrains the location of the mm core to lie $\leq$ 2 pc downstream of the jet apex, and also reveals the existence of a transition from parabolic to conical jet expansion at a distance of $\sim$54 pc from the core, corresponding to $\sim$1.5$\times$10$^{\rm 6}$ Schwarzschild radii. The cross-correlation analysis reveals a strong correlation between the radio and $\gamma$-ray data, with the 1 mm emission lagging $\sim$49 days behind the $\gamma$ rays. Based on this, we infer that the high energy emission is produced at a distance of $\sim$1 pc from the VLBI core, suggesting that the seed photon fields for the external Compton mechanism originate either in the dusty torus or in the broad-line region.

Radio galaxies host relativistic jets oriented away from our line of sight, making them challenging targets for Very High Energy (VHE, E$>$100 GeV) $\gamma$-ray detectors. Indeed, out of $\sim100$ extragalactic sources detected at $E>100$ GeV, only six are radio galaxies. The new Cherenkov Telescope Array (CTA) will provide order-of-magnitude improvements in sensitivity with respect to the present generation of ground-based $\gamma$-ray observatories, opening new frontiers for high-energy studies of radio galaxies. In a previous paper, we studied the CTA detection prospects of misaligned jets from the Third Fermi-LAT catalog (3FGL) sample. In this work, we take advantage of the expanded sample from the Fourth Fermi-LAT catalog (4FGL), including roughly double the number of sources. We simulate CTA observations of 41 $\gamma$-ray radio galaxies, extrapolating their GeV spectrum into the TeV range assuming different spectral shapes. We predict that the CTA will detect eleven new TeV radio galaxies with an observational campaign of 50 hours per source, under the realistic assumption of a spectral cutoff at 0.5 TeV. This would increase the sample of VHE radio galaxies by a factor of three, and result in the first detection of FR II radio galaxies at these energies. By simulating CTA observations with 5 hours exposure, we predict that two Fermi-LAT radio galaxies should already be within reach of current TeV observatories. Finally we have investigated the prospects for a CTA detection of extended VHE emission from the lobe-dominated FR I Fornax A, and predict that this will be possible for integration times $\gtrsim\!50$ hours. We conclude that the CTA will significantly impact our understanding of misaligned jets at TeV energies, allowing us to perform population studies and compare the two main radio galaxy subclasses for the first time in this energy band.

In this paper, we characterize the first $\gamma$-ray flaring episode of the FSRQ PKS 0346-27 (z=0.991), as revealed by Fermi-LAT monitoring data, and the concurrent multi-wavelength variability observed from radio through X-rays. The quasi-simultaneous multi-wavelength coverage allowed us to construct time-resolved spectral energy distributions (SEDs). PKS 0346-27 entered an elevated $\gamma$-ray activity state starting from the beginning of 2018. The high-state continued throughout the year, displaying the highest fluxes in May 2018. We find evidence of short-time scale variability down to $\sim$1.5 hours, which constrains the $\gamma$-ray emission region to be compact. The extended flaring period was characterized by a persistently harder spectrum with respect to the quiescent state, indicating changes in the broadband spectral properties of the source. This was confirmed by the multi-wavelength observations, which show a shift in the position of the two SED peaks by $\sim$2 orders of magnitude in energy and peak flux value. As a result, during the high state the non-thermal jet emission completely outshines the thermal contribution from the dust torus and accretion disk. The broadband SED of PKS 0346-27 transitions from a typical Low-Synchrotron-Peaked (LSP) to the Intermediate-Synchrotron-Peaked (ISP) class, a behavior previously observed in other flaring $\gamma$-ray sources. Our one-zone leptonic emission model of the high-state SEDs constrains the $\gamma$-ray emission region to have a lower magnetic field, larger radius, and higher maximum electron Lorentz factors with respect to the quiescent SED. Finally, we note that the bright and hard $\gamma$-ray spectrum observed during the peak of flaring activity in May 2018 implies that PKS 0346-27 could be a promising target for future ground-based Cherenkov observatories such as the CTA.

In the framework of the TANAMI multi-wavelength and VLBI monitoring, we study the evolution of the parsec-scale radio emission in radio galaxies in the southern hemisphere and their relationship to the $\gamma$-ray properties. In this first paper, we focus on Fermi-LAT-detected sources. We perform a kinematic analysis for five $\gamma$-ray detected radio galaxies using multi-epoch 8.4 GHz VLBI images, deriving limits on intrinsic jet parameters. We analyzed Fermi-LAT data in order to study possible connections between the $\gamma$-ray properties and the pc-scale jets of Fermi-LAT-detected radio galaxies, both in terms of variability and average properties. We discuss the individual source results and draw preliminary conclusions on sample properties including published VLBI results from the MOJAVE survey, with a total of fifteen sources. We find that the first $\gamma$-ray detection of Pictor A might be associated with the passage of a new VLBI component through the radio core. For the peculiar AGN PKS 0521-36, we detect subluminal parsec-scale jet motions, and we confirm the presence of fast $\gamma$-ray variability in the source down to timescales of 6 hours. We robustly confirm the presence of significant superluminal motion, up to $\beta_{app}\sim$3, in the jet of the TeV radio galaxy PKS 0625-35. Finally, we place a lower limit on the age of the Compact Symmetric Object (CSO) PKS 1718-649. We draw some preliminary conclusions on the relationship between pc-scale jets and $\gamma$-ray emission in radio galaxies. We find that the VLBI core flux density correlates with the $\gamma$-ray flux, as seen in blazars. On the other hand, the $\gamma$-ray luminosity does not show any dependence on the core brightness temperature and core dominance, two indicators of Doppler boosting, suggesting that $\gamma$-ray emission in radio galaxies is not driven by orientation-dependent effects.

In March 2018, Mukherjee (2018) reported the detection by VERITAS of very-high-energy emission (VHE; > 100 GeV) from 3C 264. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this gamma-ray emitting jet. The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and cover a time period of ~2 years. Through a segmented wavelet decomposition method we estimate the apparent displacement of individual plasma features; we then perform a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the SED. Proper motion is mostly detected along the edges of the flow, which appears strongly limb-brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ~$11.5$ c. This constrains the jet viewing angle to assume relatively small values ($\theta\lesssim10^{\circ}$). In the acceleration region, extending up to a de-projected distance of ~$4.8\times10^4$ Schwarzschild radii (~$11$ pc), the jet is collimating ($r\propto z^{0.40\pm 0.04}$), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated ($U_B/U_e>1$), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling.

High-redshift ($z>2$) blazars are the most powerful members of the blazar family. Yet, only a handful of them have both X-ray and $\gamma$-ray detection, thereby making it difficult to characterize the energetics of the most luminous jets. Here, we report, for the first time, the Fermi-Large Area Telescope detection of the significant $\gamma$-ray emission from the high-redshift blazar DA 193 ($z=2.363$). Its time-averaged $\gamma$-ray spectrum is soft ($\gamma$-ray photon index = $2.9\pm0.1$) and together with a relatively flat hard X-ray spectrum (14$-$195 keV photon index = $1.5\pm0.4$), DA 193 presents a case to study a typical high-redshift blazar with inverse Compton peak being located at MeV energies. An intense GeV flare was observed from this object in the first week of 2018 January, a phenomenon rarely observed from high-redshift sources. What makes this event a rare one is the observation of an extremely hard $\gamma$-ray spectrum (photon index = $1.7\pm0.2$), which is somewhat unexpected since high-redshift blazars typically exhibit a steep falling spectrum at GeV energies. The results of our multi-frequency campaign, including both space- (Fermi, NuSTAR, and Swift) and ground-based (Steward and Nordic Optical Telescope) observatories, are presented and this peculiar $\gamma$-ray flare is studied within the framework of a single-zone leptonic emission scenario.

PKS 1718$-$649 is one of the closest and most comprehensively studied candidates of a young active galactic nucleus (AGN) that is still embedded in its optical host galaxy. The compact radio structure, with a maximal extent of a few parsecs, makes it a member of the group of compact symmetric objects (CSO). Its environment imposes a turnover of the radio synchrotron spectrum towards lower frequencies, also classifying PKS 1718$-$649 as gigahertz-peaked radio spectrum (GPS) source. Its close proximity has allowed the first detection of extended X-ray emission in a GPS/CSO source with Chandra that is for the most part unrelated to nuclear feedback. However, not much is known about the nature of this emission. By co-adding all archival Chandra data and complementing these datasets with the large effective area of XMM-Newton, we are able to study the detailed physics of the environment of PKS 1718$-$649. Not only can we confirm that the bulk of the $\lesssim$kiloparsec-scale environment emits in the soft X-rays, but we also identify the emitting gas to form a hot, collisionally ionized medium. While the feedback of the central AGN still seems to be constrained to the inner few parsecs, we argue that supernovae are capable of producing the observed large-scale X-ray emission at a rate inferred from its estimated star formation rate.

TANAMI is a multiwavelength program monitoring active galactic nuclei (AGN) south of -30deg declination including high-resolution Very Long Baseline Interferometry (VLBI) imaging, radio, optical/UV, X-ray and gamma-ray studies. We have previously published first-epoch 8.4GHz VLBI images of the parsec-scale structure of the initial sample. In this paper, we present images of 39 additional sources. The full sample comprises most of the radio- and gamma-ray brightest AGN in the southern quarter of the sky, overlapping with the region from which high-energy (>100TeV) neutrino events have been found. We characterize the parsec-scale radio properties of the jets and compare with the quasi-simultaneous Fermi/LAT gamma-ray data. Furthermore, we study the jet properties of sources which are in positional coincidence with high-energy neutrino events as compared to the full sample. We test the positional agreement of high-energy neutrino events with various AGN samples. Our observations yield the first images of many jets below -30deg declination at milliarcsecond resolution. We find that gamma-ray loud TANAMI sources tend to be more compact on parsec-scales and have higher core brightness temperatures than gamma-ray faint jets, indicating higher Doppler factors. No significant structural difference is found between sources in positional coincidence with high-energy neutrino events and other TANAMI jets. The 22 gamma-ray brightest AGN in the TANAMI sky show only a weak positional agreement with high-energy neutrinos demonstrating that the >100TeV IceCube signal is not simply dominated by a small number of the $\gamma$-ray brightest blazars. Instead, a larger number of sources have to contribute to the signal with each individual source having only a small Poisson probability for producing an event in multi-year integrations of current neutrino detectors.

Misaligned AGN (MAGNs), i.e., radio-loud AGNs with the jet not pointing directly towards us, represent a new class of GeV emitters revealed by the Fermi space telescope. Although they comprise only a small fraction of the high-energy sources, MAGNs are extremely interesting objects offering a different perspective to study high-energy processes with respect to blazars. The aim of this work is to evaluate the impact of the new-generation Cherenkov Telescope Array (CTA) on the MAGN class and propose possible observational strategies to optimize their detection.

Radio galaxies are relatively faint at $\gamma$-ray energies, where they make up only 1-2% of all AGN detected by Fermi-LAT. However, they offer a unique perspective to study the intrinsic properties of AGN jets. For this reason, the combination of $\gamma$-ray and multi-wavelength data with high-resolution VLBI monitoring is a powerful tool to tackle the basic unanswered questions about AGN jets. Here we present preliminary results from a sample study of radio galaxies in the Southern hemisphere observed by the TANAMI VLBI monitoring program. We obtain high-resolution maps at 8.4 and 22.3 GHz, and study the jet kinematics using multi-epoch data. We present a preliminary kinematic analysis for the peculiar $\gamma$-ray AGN PKS 0521$-$36.