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.

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.

We present results from a multiwavelength study of the blazar PKS 1954-388 at radio, UV, X-ray, and gamma-ray energies. A RadioAstron observation at 1.66 GHz in June 2012 resulted in the detection of interferometric fringes on baselines of 6.2 Earth-diameters. This suggests a source frame brightness temperature of greater than 2x10^12 K, well in excess of both equipartition and inverse Compton limits and implying the existence of Doppler boosting in the core. An 8.4 GHz TANAMI VLBI image, made less than a month after the RadioAstron observations, is consistent with a previously reported superluminal motion for a jet component. Flux density monitoring with the Australia Telescope Compact Array confirms previous evidence for long-term variability that increases with observing frequency. A search for more rapid variability revealed no evidence for significant day-scale flux density variation. The ATCA light-curve reveals a strong radio flare beginning in late 2013 which peaks higher, and earlier, at higher frequencies. Comparison with the Fermi gamma-ray light-curve indicates this followed ~9 months after the start of a prolonged gamma-ray high-state -- a radio lag comparable to that seen in other blazars. The multiwavelength data are combined to derive a Spectral Energy Distribution, which is fitted by a one-zone synchrotron-self-Compton (SSC) model with the addition of external Compton (EC) emission.

Quasi-simultaneous observations of the Flat Spectrum Radio Quasar PKS 2326-502 were carried out in the gamma-ray, X-ray, UV, optical, near-infrared, and radio bands. Thanks to these observations we are able to characterize the spectral energy distribution of the source during two flaring and one quiescent gamma-ray states. These data were used to constrain one-zone leptonic models of the spectral energy distributions of each flare and investigate the physical conditions giving rise to them. While modeling one flare only required changes to the electron spectrum, the other flare needed changes in both the electron spectrum and the size of the emitting region with respect to the quiescent state. These results are consistent with an emerging pattern of two broad classes of flaring states seen in blazars.

Simultaneous broadband spectral and temporal studies of blazars are an important tool for investigating active galactic nuclei (AGN) jet physics. We study the spectral evolution between quiescent and flaring periods of 22 radio-loud AGN through multi-epoch, quasi-simultaneous broadband spectra. For many of these sources these are the first broadband studies. We use a Bayesian block analysis of \Fermi/LAT light curves in order to determine time ranges of constant flux for constructing quasi-simultaneous SEDs. The shapes of the resulting 81 SEDs are described by two logarithmic parabolas and a blackbody spectrum where needed. For low states the peak frequencies and luminosities agree well with the blazar sequence, higher luminosity implying lower peak frequencies. This is not true for sources in a high state. The $\gamma$-ray photon index in Fermi/LAT correlates with the synchrotron peak frequency in low and intermediate states. No correlation is present in high states. The black hole mass cannot be determined from the SEDs. Surprisingly, the thermal excess often found in FSRQs at optical/UV wavelengths can be described by blackbody emission and not an accretion disk spectrum. The "harder-when-brighter" trend, typically seen in X-ray spectra of flaring blazars, is visible in the blazar sequence. Our results for low and intermediate states, as well as the Compton dominance, are in agreement with previous results. Black hole mass estimates using the parameters from Bonchi (2013) are in agreement with some of the more direct measurements. For two sources, estimates disagree by more than four orders of magnitude, possibly due to boosting effects. The shapes of the thermal excess seen predominantly in flat spectrum radio quasars are inconsistent with a direct accretion disk origin.

We report on Fermi Large Area Telescope (LAT) and multi-wavelength results on the recently-discovered very-high-energy (VHE, $E>$ 100 GeV) blazar S4 0954+65 ($z=0.368$) during an exceptionally bright optical flare in 2015 February. During the time period (2015 February, 13/14, or MJD 57067) when the MAGIC telescope detected VHE $\gamma$-ray emission from the source, the Fermi-LAT data indicated a significant spectral hardening at GeV energies, with a power-law photon index of $1.8 \pm 0.1$---compared with the 3FGL value (averaged over four years of observation) of $2.34 \pm 0.04$. In contrast, Swift/XRT data showed a softening of the X-ray spectrum, with a photon index of $1.72 \pm 0.08$ (compared with $1.38 \pm 0.03$ averaged during the flare from MJD 57066 to 57077), possibly indicating a modest contribution of synchrotron photons by the highest-energy electrons superposed on the inverse Compton component. Fitting of the quasi-simultaneous ($<1$ day) broadband spectrum with a one-zone synchrotron plus inverse-Compton model revealed that GeV/TeV emission could be produced by inverse-Compton scattering of external photons from the dust torus. We emphasize that a flaring blazar showing high flux of $\gtrsim 1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ ($E>$ 100 MeV) and a hard spectral index of $\Gamma_{\rm GeV} < 2.0$ detected by Fermi-LAT on daily time scales is a promising target for TeV follow-up by ground-based Cherenkov telescopes to discover high-redshift blazars, investigate their temporal variability and spectral features in the VHE band, and also constrain the intensity of the extragalactic background light.

The discovery of extraterrestrial very-high-energy neutrinos by the IceCube collaboration has launched a quest for the identification of their astrophysical sources. Gamma-ray blazars have been predicted to yield a cumulative neutrino signal exceeding the atmospheric background above energies of 100 TeV, assuming that both the neutrinos and the gamma-ray photons are produced by accelerated protons in relativistic jets. Since the background spectrum falls steeply with increasing energy, the individual events with the clearest signature of being of an extraterrestrial origin are those at PeV energies. Inside the large positional-uncertainty fields of the first two PeV neutrinos detected by IceCube, the integrated emission of the blazar population has a sufficiently high electromagnetic flux to explain the detected IceCube events, but fluences of individual objects are too low to make an unambiguous source association. Here, we report that a major outburst of the blazar PKS B1424-418 occurred in temporal and positional coincidence with the third PeV-energy neutrino event (IC35) detected by IceCube. Based on an analysis of the full sample of gamma-ray blazars in the IC35 field and assuming a photo-hadronic emission model, we show that the long-term average gamma-ray emission of blazars as a class is in agreement with both the measured all-sky flux of PeV neutrinos and the spectral slope of the IceCube signal. The outburst of PKS B1424-418 has provided an energy output high enough to explain the observed PeV event, indicative of a direct physical association.

Using high-resolution radio imaging with VLBI techniques, the TANAMI program has been observing the parsec-scale radio jets of southern (declination south of -30{\deg}) gamma-ray bright AGN simultaneously with Fermi/LAT monitoring of their gamma-ray emission. We present the radio and gamma-ray properties of the TANAMI sources based on one year of contemporaneous TANAMI and Fermi/LAT data. A large fraction (72%) of the TANAMI sample can be associated with bright gamma-ray sources for this time range. Association rates differ for different optical classes with all BL Lacs, 76% of quasars and just 17% of galaxies detected by the LAT. Upper limits were established on the gamma-ray flux from TANAMI sources not detected by LAT. This analysis led to the identification of three new Fermi sources whose detection was later confirmed. The gamma-ray and radio luminosities are related by $L_\gamma \propto L_r^{0.89+-0.04}$. The brightness temperatures of the radio cores increase with the average gamma-ray luminosity, and the presence of brightness temperatures above the inverse Compton limit implies strong Doppler boosting in those sources. The undetected sources have lower gamma/radio luminosity ratios and lower contemporaneous brightness temperatures. Unless the Fermi/LAT-undetected blazars are strongly gamma-ray-fainter than the Fermi/LAT-detected ones, their gamma-ray luminosity should not be significantly lower than the upper limits calculated here.

As part of the TANAMI multiwavelength progam, we discuss new X-ray observations of the $\gamma$-ray and radio-loud Narrow Line Seyfert galaxy ($\gamma$-NLS1) PKS 2004-447. The active galaxy is a member of a small sample of radio-loud NLS1s detected in $\gamma$-rays by the Fermi Large Area Telescope. It is the radio-loudest and only southern-hemisphere source in this sample. We present results from our X-ray monitoring program comprised of Swift snapshot observations from 2012 through 2014 and two new X-ray observations with XMM-Newton in 2012. We analyze the X-ray spectrum and variability of this peculiar source using supplementary archival data from 2004 and 2011. The (0.5-10) keV spectrum is well described by a flat power law, which can be interpreted as non-thermal emission from a relativistic jet. The source exhibits moderate flux variability on timescales of both months and years. Correlated brightness variations in the (0.5-2) keV and (2-10) keV bands are explained by a single variable spectral component, such as the jet. A possible soft excess seen in the data from 2004 cannot be confirmed by the new \xmm{} observations in 2012. Any contribution to the total flux in 2004 is less than $20\%$ of the power-law component. The (0.5-10) keV luminosities of PKS 2004-447 are in the range of (0.5--2.7)$\times10^{44}\,\mathrm{erg\,s}^{-1}$. A comparison of the X-ray properties among the known $\gamma$-NLS1 galaxies shows that X-ray spectrum is typically dominated by a flat power law without intrinsic absorption. These objects are moderately variable in their brightness, while spectral variability is observed in at least two sources. The major difference across the X-ray spectra of $\gamma$-NLS1s is the luminosity, which spans a range of almost two orders of magnitude, from $10^{44}\,\mathrm{erg\,s}^{-1}$ to $10^{46}\,\mathrm{erg\,s}^{-1}$ in the (0.5-10) keV band.

Gamma-ray detected radio-loud narrow-line Seyfert 1 (g-NLS1) galaxies constitute a small but interesting sample of the gamma-ray loud AGN. The radio-loudest g-NLS1 known, PKS 2004-447, is located in the southern hemisphere and is monitored in the radio regime by the multiwavelength monitoring program TANAMI. We aim for the first detailed study of the radio morphology and long-term radio spectral evolution of PKS 2004-447, which are essential to understand the diversity of the radio properties of g-NLS1s. The TANAMI VLBI monitoring program uses the Australian Long Baseline Array (LBA) and telescopes in Antarctica, Chile, New Zealand, and South Africa to monitor the jets of radio-loud active galaxies in the southern hemisphere. Lower resolution radio flux density measurements at multiple radio frequencies over four years of observations were obtained with the Australia Telescope Compact Array (ATCA). The TANAMI VLBI image at 8.4 GHz shows an extended one-sided jet with a dominant compact VLBI core. Its brightness temperature is consistent with equipartition, but it is an order of magnitude below other g-NLS1s with the sample value varying over two orders of magnitude. We find a compact morphology with a projected large-scale size <11 kpc and a persistent steep radio spectrum with moderate flux-density variability. PKS 2004-447 appears to be a unique member of the g-NLS1 sample. It exhibits blazar-like features, such as a flat featureless X-ray spectrum and a core dominated, one-sided parsec-scale jet with indications for relativistic beaming. However, the data also reveal properties atypical for blazars, such as a radio spectrum and large-scale size consistent with Compact-Steep-Spectrum (CSS) objects, which are usually associated with young radio sources. These characteristics are unique among all g-NLS1s and extremely rare among gamma-ray loud AGN.

The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope.The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons --- and hence their neutrino progenitors --- from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for different possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin.Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653$-$329 and 1714$-$336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90\% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than $-2.4$.

Since the discovery of a neutrino flux in excess of the atmospheric background by the IceCube Collaboration, searches for the astrophysical sources have been ongoing. Due to the steeply falling background towards higher energies, the PeV events detected in three years of IceCube data are the most likely ones to be of extraterrestrial origin. Even excluding the PeV events detected so far, the neutrino flux is well above the atmospheric background, so it is likely that a number of sub-PeV events originate from the same astrophysical sources that produce the PeV events. We study the high-energy properties of AGN that are positionally coincident with the neutrino events from three years of IceCube data and show the results for event number 4. IC 4 is a event with a low angular error (7.1$^\circ$) and a large deposited energy of 165 TeV. We use multiwavelength data, including Fermi/LAT and X-ray data, to construct broadband spectra and present parametrizations of the broadband spectral energy distributions with logarithmic parabolas. Assuming the X-ray to {\gamma}-ray emission in blazars originates in the photoproduction of pions by accelerated protons, their predicted neutrino luminosity can be estimated. The measurements of the diffuse extragalactic background by Fermi/LAT gives us an estimate of the flux contributions from faint unresolved blazars. Their contribution increases the number of expected events by a factor of $\sim$2. We conclude that the detection of the IceCube neutrinos IC4, IC14, and IC20 can be explained by the integral emission of blazars, even though no individual source yields a sufficient energy output.

Multiwavelength observations have revealed the highly unusual properties of the gamma-ray source PMN J1603-4904, which are difficult to reconcile with any other well established gamma-ray source class. The object is either a very atypical blazar or compact jet source seen at a larger angle to the line of sight. In order to determine the physical origin of the high-energy emission processes in PMN J1603-4904, we study the X-ray spectrum in detail. We performed quasi-simultaneous X-ray observations with XMM-Newton and Suzaku in 2013 September, resulting in the first high signal-to-noise X-ray spectrum of this source. The 2-10 keV X-ray spectrum can be well described by an absorbed power law with an emission line at 5.44$\pm$0.05 keV (observed frame). Interpreting this feature as a K{\alpha} line from neutral iron, we determine the redshift of PMN J1603-4904 to be z=0.18$\pm$0.01, corresponding to a luminosity distance of 872$\pm$54 Mpc. The detection of a redshifted X-ray emission line further challenges the original BL Lac classification of PMN J1603-4904. This result suggests that the source is observed at a larger angle to the line of sight than expected for blazars, and thus the source would add to the elusive class of gamma-ray loud misaligned-jet objects, possibly a {\gamma}-ray bright young radio galaxy.

Centaurus A is the closest radio-loud active galaxy. Very Long Baseline Interferometry (VLBI) enables us to study the jet-counterjet system on milliarcsecond (mas) scales, providing essential information for jet emission and propagation models. We study the evolution of the central parsec jet structure of Cen A over 3.5 years. The proper motion analysis of individual jet components allows us to constrain jet formation and propagation and to test the proposed correlation of increased high energy flux with jet ejection events. Cen A is an exceptional laboratory for such detailed study as its proximity translates to unrivaled linear resolution, where 1 mas corresponds to 0.018 pc. The first 7 epochs of high-resolution TANAMI VLBI observations at 8 GHz of Cen A are presented, resolving the jet on (sub-)mas scales. They show a differential motion of the sub-pc scale jet with significantly higher component speeds further downstream where the jet becomes optically thin. We determined apparent component speeds within a range of 0.1c to 0.3c, as well as identified long-term stable features. In combination with the jet-to-counterjet ratio we can constrain the angle to the line of sight to ~12{\deg} to 45{\deg}. The high resolution kinematics are best explained by a spine-sheath structure supported by the downstream acceleration occurring where the jet becomes optically thin. On top of the underlying, continuous flow, TANAMI observations clearly resolve individual jet features. The flow appears to be interrupted by an obstacle causing a local decrease in surface brightness and a circumfluent jet behavior. We propose a jet-star interaction scenario to explain this appearance. The comparison of jet ejection times with high X-ray flux phases yields a partial overlap of the onset of the X-ray emission and increasing jet activity, but the limited data do not support a robust correlation.

The IceCube Collaboration has announced the discovery of a neutrino flux in excess of the atmospheric background. Due to the steeply falling atmospheric background spectrum, events at PeV energies are most likely of extraterrestrial origin. We present the multiwavelength properties of the six radio brightest blazars positionally coincident with these events using contemporaneous data of the TANAMI blazar sample, including high-resolution images and spectral energy distributions. Assuming the X-ray to {\gamma}-ray emission originates in the photoproduction of pions by accelerated protons, the integrated predicted neutrino luminosity of these sources is large enough to explain the two detected PeV events.

We investigate the nature and classification of PMNJ1603-4904, a bright radio source close to the Galactic plane, which is associated with one of the brightest hard-spectrum gamma-ray sources detected by Fermi/LAT. It has previously been classified as a low-peaked BL Lac object based on its broadband emission and the absence of optical emission lines. Optical measurements, however, suffer strongly from extinction and the absence of pronounced short-time gamma-ray variability over years of monitoring is unusual for a blazar. We are combining new and archival multiwavelength data in order to reconsider the classification and nature of this unusual gamma-ray source. For the first time, we study the radio morphology at 8.4GHz and 22.3GHz, and its spectral properties on milliarcsecond (mas) scales, based on VLBI observations from the TANAMI program. We combine the resulting images with multiwavelength data in the radio, IR, optical/UV, X-ray, and gamma-ray regimes. PMNJ1603-4904 shows a symmetric brightness distribution at 8.4GHz on mas-scales, with the brightest, and most compact component in the center of the emission region. The morphology is reminiscent of a Compact Symmetric Object (CSO). Such objects have been predicted to produce gamma-ray emission but have not been detected as a class by Fermi/LAT so far. Sparse (u, v)-coverage at 22.3GHz prevents an unambiguous modeling of the source morphology. IR measurements reveal an excess in the spectral energy distribution (SED), which can be modeled with a blackbody with a temperature of about 1600K, and which is usually not present in blazar SEDs. The VLBI data and the shape of the SED challenge the current blazar classification. PMNJ1603-4904 seems to be either a highly peculiar BL Lac object or a misaligned jet source. In the latter case, the intriguing VLBI structure opens room for a possible classification as a gamma-ray bright CSO.

The Southern Hemisphere VLBI monitoring program TANAMI provides dual-frequency (8 GHz and 22 GHz), milliarcsecond monitoring of extragalactic jets south of -30 degrees declination. The TANAMI sample consists of a combined radio and gamma-ray selected subsample of currently 80 AGN jets, with new gamma-ray bright sources being added upon detections by Fermi/LAT. Supporting programs provide simultaneous multiwavelength coverage of all sources, in order to construct broadband spectral energy distributions (SEDs) of flaring and quiescence source states, as well as a rapid follow-up of high-energy flares. This combined setup allows us to continuously study the spectral and structural evolution of highly energetic extragalactic jets and test correlations in different wavebands, providing crucial information on underlying physical mechanisms. Here, we present jet kinematics of Centaurus A and show preliminary VLBI results on PKS 0625-354 and the time-dependent spectral index image of PKS 0537-441.