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.

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.

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.

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$.

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.

Centaurus A is the closest active galactic nucleus. High resolution imaging using Very Long Baseline Interferometry (VLBI) enables us to study the spectral and kinematic behavior of the radio jet-counterjet system on sub-parsec scales, providing essential information for jet emission and formation models. Our aim is to study the structure and spectral shape of the emission from the central-parsec region of Cen A. As a target of the Southern Hemisphere VLBI monitoring program TANAMI (Tracking Active Galactic Nuclei with Milliarcsecond Interferometry), VLBI observations of Cen A are made regularly at 8.4 and 22.3 GHz with the Australian Long Baseline Array (LBA) and associated telescopes in Antarctica, Chile, and South Africa. The first dual-frequency images of this source are presented along with the resulting spectral index map. An angular resolution of 0.4 mas x 0.7 mas is achieved at 8.4 GHz, corresponding to a linear scale of less than 0.013 pc. Hence, we obtain the highest resolution VLBI image of Cen A, comparable to previous space-VLBI observations. By combining with the 22.3 GHz image, which has been taken without contributing transoceanic baselines at somewhat lower resolution, we present the corresponding dual-frequency spectral index distribution along the sub-parsec scale jet revealing the putative emission regions for recently detected gamma-rays from the core region by Fermi/LAT. We resolve the innermost structure of the milliarcsecond scale jet and counterjet system of Cen A into discrete components. The simultaneous observations at two frequencies provide the highest resolved spectral index map of an AGN jet allowing us to identify multiple possible sites as the origin of the high energy emission.

We present gamma-ray observations with the LAT on board the Fermi Gamma-Ray Telescope of the nearby radio galaxy Centaurus~A. The previous EGRET detection is confirmed, and the localization is improved using data from the first 10 months of Fermi science operation. In previous work, we presented the detection of the lobes by the LAT; in this work, we concentrate on the gamma-ray core of Cen~A. Flux levels as seen by the LAT are not significantly different from that found by EGRET, nor is the extremely soft LAT spectrum ($\G=2.67\pm0.10_{stat}\pm0.08_{sys}$ where the photon flux is $\Phi\propto E^{-\G}$). The LAT core spectrum, extrapolated to higher energies, is marginally consistent with the non-simultaneous HESS spectrum of the source. The LAT observations are complemented by simultaneous observations from Suzaku, the Swift Burst Alert Telescope and X-ray Telescope, and radio observations with the Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry (TANAMI) program, along with a variety of non-simultaneous archival data from a variety of instruments and wavelengths to produce a spectral energy distribution (SED). We fit this broadband data set with a single-zone synchrotron/synchrotron self-Compton model, which describes the radio through GeV emission well, but fails to account for the non-simultaneous higher energy TeV emission observed by HESS from 2004-2008. The fit requires a low Doppler factor, in contrast to BL Lacs which generally require larger values to fit their broadband SEDs. This indicates the $\g$-ray emission originates from a slower region than that from BL Lacs, consistent with previous modeling results from Cen~A. This slower region could be a slower moving layer around a fast spine, or a slower region farther out from the black hole in a decelerating flow.

We introduce the TANAMI program (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) which is monitoring an initial sample of 43 extragalactic jets located south of -30 degrees declination at 8.4 GHz and 22 GHz since 2007. All aspects of the program are discussed. First epoch results at 8.4 GHz are presented along with physical parameters derived therefrom. We present first epoch images for 43 sources, some observed for the first time at milliarcsecond resolution. Parameters of these images as well as physical parameters derived from them are also presented and discussed. These and subsequent images from the TANAMI survey are available at http://pulsar.sternwarte.uni-erlangen.de/tanami/ We obtain reliable, high dynamic range images of the southern hemisphere AGN. All the quasars and BL Lac objects in the sample have a single-sided radio morphology. Galaxies are either double-sided, single-sided or irregular. About 28% of the TANAMI sample has been detected by LAT during its first three months of operations. Initial analysis suggests that when galaxies are excluded, sources detected by LAT have larger opening angles than those not detected by LAT. Brightness temperatures of LAT detections and non-detections seem to have similar distributions. The redshift distributions of the TANAMI sample and sub-samples are similar to those seen for the bright gamma-ray AGN seen by LAT and EGRET but none of the sources with a redshift above 1.8 have been detected by LAT.

We present a summary of the observation strategy of TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry), a monitoring program to study the parsec-scale structure and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Australian Long Baseline Array (LBA) and the trans-oceanic antennas Hartebeesthoek, TIGO, and O'Higgins. TANAMI is focusing on extragalactic sources south of -30 degrees declination with observations at 8.4 GHz and 22 GHz every ~2 months at milliarcsecond resolution. The initial TANAMI sample of 43 sources has been defined before the launch of the Fermi Gamma Ray Space Telescope to include the most promising candidates for bright gamma-ray emission to be detected with its Large Area Telescope (LAT). Since November 2008, we have been adding new sources to the sample, which now includes all known radio- and gamma-ray bright AGN of the Southern Hemisphere. The combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN and the nature of relativistic jets.

The TANAMI (Tracking AGN with Austral Milliarcsecond Interferometry) and associated programs provide comprehensive radio monitoring of extragalactic gamma-ray sources south of declination -30 degrees. Joint quasi-simultaneous observations between the Fermi Gamma-ray Space Telescope and ground based observatories allow us to discriminate between competing theoretical blazar emission models. High resolution VLBI observations are the only way to spatially resolve the sub-parsec level emission regions where the high-energy radiation originates. The gap from radio to gamma-ray energies is spanned with near simultaneous data from the Swift satellite and ground based optical observatories. We present early results from the TANAMI program in the context of this panchromatic suite of observations.

TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) is a monitoring program to study the parsec-scale structures and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Long Baseline Array and associated telescopes. Extragalactic jets south of -30 degrees declination are observed at 8.4 GHz and 22 GHz every two months at milliarcsecond resolution. The initial TANAMI sample is a hybrid radio and gamma-ray selected sample since the combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN.