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.

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.

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.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

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

Very-high-energy $\gamma$-ray observations of the active galaxy IC 310 with the MAGIC telescopes have revealed fast variability with doubling time scales of less than 4.8min. This implies that the emission region in IC 310 is smaller than 20% of the gravitational radius of the central supermassive black hole with a mass of $3\times 10^8 M_\odot$, which poses serious questions on the emission mechanism and classification of this enigmatic object. We report on the first quasi-simultaneous multi-frequency VLBI observations of IC 310 conducted with the EVN. We find a blazar-like one-sided core-jet structure on parsec scales, constraining the inclination angle to be less than $\sim 20^\circ$ but very small angles are excluded to limit the de-projected length of the large-scale radio jet.

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.