High-mass protostars (M$_{\star} >$ 8 M$_{\odot}$) are thought to gain the majority of their mass via short, intense bursts of growth. This episodic accretion is thought to be facilitated by gravitationally unstable and subsequently inhomogeneous accretion disks. Limitations of observational capabilities, paired with a lack of observed accretion burst events has withheld affirmative confirmation of the association between disk accretion, instability and the accretion burst phenomenon in high-mass protostars. Following its 2019 accretion burst, a heat-wave driven by a burst of radiation propagated outward from the high-mass protostar G358.93-0.03-MM1. Six VLBI (very long baseline interferometry) observations of the raditively pumped 6.7 GHz methanol maser were conducted during this period, tracing ever increasing disk radii as the heat-wave propagated outward. Concatenating the VLBI maps provided a sparsely sampled, milliarcsecond view of the spatio-kinematics of the accretion disk covering a physical range of $\sim$ 50 - 900 AU. We term this observational approach `heat-wave mapping'. We report the discovery of a Keplerian accretion disk with a spatially resolved four-arm spiral pattern around G358.93-0.03-MM1. This result positively implicates disk accretion and spiral arm instabilities into the episodic accretion high-mass star formation paradigm.

This paper reports observations of a 22 GHz water maser `superburst' in the G25.65+1.05 massive star forming region, conducted in response to an alert from the Maser Monitoring Organisation (M2O). Very long baseline interferometry (VLBI) observations using the European VLBI Network (EVN) recorded a maser flux density of $1.2 \times 10^{4}$ Jy. The superburst was investigated in the spectral, structural and temporal domains and its cause was determined to be an increase in maser path length generated by the superposition of multiple maser emitting regions aligning in the line of sight to the observer. This conclusion was based on the location of the bursting maser in the context of the star forming region, its complex structure, and its rapid onset and decay.

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High-spatial-resolution measurements of the source size and displacement can discriminate between these scenarios. We present very-long-baseline interferometry observations, performed 207.4 days after the merger by using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milli-arc seconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet.

The first binary neutron star merger, detected through both electromagnetic radiation and gravitational waves on the 17th of August 2017, raised the question whether a narrow relativistic jet or a more isotropic outflow was launched as a consequence of the merger. High resolution measurements of the source size and position can provide the answer. Very Long Baseline Interferometry observations, performed 207.4 days after the binary merger through a global network of 32 radio telescopes spread over five continents, constrain the apparent source size to be smaller than 2 milliarcseconds at the 90% confidence level. This excludes the possibility that a nearly isotropic, mildly relativistic outflow is responsible for the emission, as in this case its apparent size, after more than six months of expansion, should have been significantly larger and resolved by the VLBI observation. Our size measurement proves that in at least 10% of neutron star mergers a structured relativistic jet should be produced.

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.

Context. The binary system AR Sco contains an M star and the only known radio-pulsing white dwarf. The system shows emission from radio to X-rays, likely dominated by synchrotron radiation. The mechanism that produces most of this emission remains unclear. Two competing scenarios have been proposed: collimated outflows, and direct interaction between the magnetospheres of the white dwarf and the M star.Aims. The two proposed scenarios can be tested via very long baseline interferometric radio observations.Methods. We conducted a radio observation with the Australian Long Baseline Array (LBA) on 20 Oct. 2016 at 8.5 GHz to study the compactness of the radio emission. Simultaneous data with the Australian Telescope Compact Array (ATCA) were also recorded for a direct comparison of the obtained flux densities.Results. AR Sco shows radio emission compact on milliarcsecond angular scales (less than or similar to 0.02 AU, or 4 R-circle dot). The emission is orbitally modulated, with an average flux density of approximate to 6.5 mJy. A comparison with the simultaneous ATCA data shows that no flux is resolved out on mas scales, implying that the radio emission is produced in this compact region. Additionally, the obtained radio light curves on hour timescales are consistent with the optical light curve.Conclusions. The radio emission in AR Sco is likely produced in the magnetosphere of the M star or the white dwarf, and we see no evidence for a radio outflow or collimated jets significantly contributing to the radio emission.

The TANAMI program has been observing parsec-scale radio jets of southern (declination south of − 30°) γ -ray bright AGN, simultaneously with Fermi /LAT monitoring of their γ -ray emission, via high-resolution radio imaging with Very Long Baseline Interferometry techniques. We present the radio and γ -rayproperties 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 γ -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 γ -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 γ -ray and radio luminosities are related by L γ ∝ L r 0.89±0.04 . The brightness temperatures of the radio cores increase with the average γ -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 γ /radio luminosity ratios and lower contemporaneous brightness temperatures. Unless the Fermi /LAT-undetected blazars are much γ -ray-fainter than the Fermi /LAT-detected sources, their γ -ray luminosity should not be significantly lower than the upper limits calculated here.

Context. Simultaneous broadband spectral and temporal studies of blazars are an important tool for investigating active galactic nuclei (AGN) jet physics.Aims. We study the spectral evolution between quiescent and flaring periods of 22 radio-loud AGN through multiepoch, quasisimultaneous broadband spectra. For many of these sources these are the first broadband studies.Methods. We use a Bayesian block analysis of Fermi/LAT light curves to determine time ranges of constant flux for constructing quasi-simultaneous spectral energy distributions (SEDs). The shapes of the resulting 81 SEDs are described by two logarithmic parabolas and a blackbody spectrum where needed.Results. The peak frequencies and luminosities agree well with the blazar sequence for low states with higher luminosity implying lower peak frequencies. This is not true for sources in high states. 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.Conclusions. The so-called 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 recently published parameters are in agreement with some of the more direct measurements. For two sources, estimates disagree by more than four orders of magnitude, possibly owing to boosting effects. The shapes of the thermal excess seen predominantly in flat spectrum radio quasars are inconsistent with a direct accretion disk origin.

Context. 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. Aims. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope. Methods. 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. Results. 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 IC 14 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.

Context. 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.Aims. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope.Methods. 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.Results. 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 IC 14 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.

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.

Context. Centaurus A (Cen A) is the closest radio-loud active galactic nucleus. Very Long Baseline Interferometry (VLBI) enables us to study the spectral and kinematic behavior of the radio jet-counterjet system on milliarcsecond scales, providing essential information for jet emission and propagation models.Aims. In the framework of the TANAMI monitoring, we investigate the kinematics and complex structure of Cen A on subparsec scales. We have been studying the evolution of the central parsec jet structure of CenA for 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 a detailed study because its proximity translates to unrivaled linear resolution, where one milliarcsecond corresponds to 0.018 pc.Methods. As a target of the southern-hemisphere VLBI monitoring program TANAMI, observations of CenA are done approximately every six months at 8.4 GHz with the Australian Long Baseline Array (LBA) and associated telescopes in Antarctica, Chile, New Zealand, and South Africa, complemented by quasi-simultaneous 22.3 GHz observations.Results. The first seven epochs of high-resolution TANAMI VLBI observations at 8.4 GHz of CenA are presented, resolving the jet on (sub-) milliarcsecond scales. They show a differential motion of the subparsec scale jet with significantly higher component speeds farther downstream where the jet becomes optically thin. We determined apparent component speeds within a range of 0 : 1 c to 0 : 3 c and identified long-term stable features. In combination with the jet-to-counterjet ratio, we can constrain the angle to the line of sight to theta similar to 12 degrees-45 degrees.Conclusions. 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 circumfluent jet behavior. We propose a jet-star interaction scenario to explain this appearance. The comparison of jet ejection times to 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 Fermi Gamma-Ray Space Telescope has been detecting a wealth of sources where the multiwavelength counterpart is either inconclusive or missing altogether. We present a combination of factors that can be used to identify multiwavelength counterparts to these Fermi unassociated sources. This approach was used to select and investigate seven bright, high-latitude unassociated sources with radio, UV, X-ray, and gamma-ray observations. As a result, four of these sources are candidates to be active galactic nuclei, and one to be a pulsar, while two do not fit easily into these known categories of sources. The latter pair of extraordinary sources might reveal a new category subclass or a new type of gamma-ray emitter. These results altogether demonstrate the power of a multiwavelength approach to illuminate the nature of unassociated Fermi sources.

As a critical part of the Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry (TANAMI) program, in November 2007 the Australia Telescope Compact Array (ATCA) started monitoring the radio spectra of a sample of southern hemisphere active galactic nuclei (AGN) that were selected as likely candidates for detection (as well as a control sample) by the Large Area Telescope (LAT) aboard the Fermi Gamma Ray Space Observatory. The initial sample was chosen based on properties determined from AGN detections by the Energetic Gamma Ray Experiment Telescope (EGRET). Most of the initial sample has been detected by Fermi/LAT and with the addition of new detections the sample has grown to include 226 AGN, 133 of which have data for more than one epoch. For the majority of these AGN, our monitoring program provides the only dynamic radio spectra available. The ATCA receiver suite makes it possible to observe several sources at frequencies between 4.5 and 41 GHz in a few hours, resulting in an excellent measure of spectral index at each epoch. By examining how the spectral index changes over time, we aim to investigate the mechanics of radio and gamma-ray emission from AGN jets.

Rapid variability in the radio flux density of the BL Lac object PKS B1144-379 has been observed at four frequencies, ranging from 1.5 to 15 GHz, with the Very Large Array and the University of Tasmania's Ceduna antenna. Intrinsic and line-of-sight effects were examined as possible causes of this variability, with interstellar scintillation best explaining the frequency dependence of the variability timescales and modulation indices. This scintillation is consistent with a compact source 20-40 mu as or 0.15-0.3 pc in size. The inferred brightness temperature for PKS B1144-379 (assuming that the observed variations are due to scintillation) is 6.2 x 10(12) K at 4.9 GHz, with approximately 10% of the total flux in the scintillating component. We show that scintillation surveys aimed at identifying variability timescales of days to weeks are an effective way to identify the active galactic nuclei with the highest brightness temperatures.

Context. There is growing evidence of relativistic jets in radio-loud narrow-line Seyfert 1 (RL-NLS1) galaxies.Aims. We constrain the observational properties of the radio emission in the first RL-NLS1 galaxy ever detected in gamma-rays, PMN J0948+0022, i.e., its flux density and structure in both total intensity and polarization, its compactness, and variability.Methods. We performed three real-time e-VLBI observations of PMN J0948+0022 at 22 GHz, using a global array including telescopes in Europe, East Asia, and Australia. These are the first e-VLBI science observations ever carried out with a global array, reaching a maximum baseline length of 12 458 km. The observations were part of a large multiwavelength campaign in 2009.Results. The source is detected at all three epochs. The structure is dominated by a bright component, more compact than 55 mu as, with a fainter component at a position angle theta similar to 35 degrees. Relativistic beaming is required by the observed brightness temperature of 3.4 x 10(11) K. Polarization is detected at a level of about 1%.Conclusions. The parameters derived by the VLBI observations, in addition to the broad-band properties, confirm that PMN J0948+0022 is similar to flat spectrum radio quasars. Global e-VLBI is a reliable and promising technique for future studies.