As the largest radio telescope in the world, the Square Kilometre Array (SKA) will lead the next generation of radio astronomy. The feats of engineering required to construct the telescope array will be matched only by the techniques developed to exploit the rich scientific value of the data. To drive forward the development of efficient and accurate analysis methods, we are designing a series of data challenges that will provide the scientific community with high-quality data sets for testing and evaluating new techniques. In this paper, we present a description and results from the first such Science Data Challenge 1 (SDC1). Based on SKA MID continuum simulated observations and covering three frequencies (560, 1400, and 9200 MHz) at three depths (8, 100, and 1000 h), SDC1 asked participants to apply source detection, characterization, and classification methods to simulated data. The challenge opened in 2018 November, with nine teams submitting results by the deadline of 2019 April. In this work, we analyse the results for eight of those teams, showcasing the variety of approaches that can be successfully used to find, characterize, and classify sources in a deep, crowded field. The results also demonstrate the importance of building domain knowledge and expertise on this kind of analysis to obtain the best performance. As high-resolution observations begin revealing the true complexity of the sky, one of the outstanding challenges emerging from this analysis is the ability to deal with highly resolved and complex sources as effectively as the unresolved source population.

We monitored the position of the close-by (about 370 Mpc) gamma-ray burst GRB 190829A, which originated from a massive star collapse, through very long baseline interferometry (VLBI) observations with the European VLBI Network and the Very Long Baseline Array, carrying out a total of nine observations between 9 and 117 days after the gamma-ray burst at 5 and 15 GHz, with a typical resolution of a few milliarcseconds. From a state-of-the art analysis of these data, we obtained valuable limits on the source size and expansion rate. The limits are in agreement with the size evolution entailed by a detailed modeling of the multiwavelength light curves with a forward-plus-reverse shock model, which agrees with the observations across almost 18 orders of magnitude in frequency (including the HESS data at TeV photon energies) and more than 4 orders of magnitude in time. Thanks to the multiwavelength, high-cadence coverage of the afterglow, inherent degeneracies in the afterglow model are broken to a large extent, allowing us to capture some unique physical insights; we find a low prompt emission efficiency of ≲10⁻³, a low fraction of relativistic electrons in the forward shock downstream χ _e < 13% (90% credible level), and a rapid decay of the magnetic field in the reverse shock downstream after the shock crossing. While our model assumes an on-axis jet, our VLBI astrometry is not sufficiently tight as to exclude any off-axis viewing angle, but we can exclude the line of sight to have been more than ∼2° away from the border of the gamma-ray-producing region based on compactness arguments.

Blazars are a sub-class of quasars with Doppler boosted jets oriented close to the line of sight, and thus efficient probes of supermassive black hole growth and their environment, especially at high redshifts. Here we report on Very Long Baseline Interferometry observations of a blazar J0906 + 6930 at z = 5.47, which enabled the detection of polarised emission and measurement of jet proper motion at parsec scales. The observations suggest a less powerful jet compared with the general blazar population, including lower proper motion and bulk Lorentz factor. This coupled with a previously inferred high accretion rate indicate a transition from an accretion radiative power to a jet mechanical power based transfer of energy and momentum to the surrounding gas. While alternative scenarios could not be fully ruled out, our results indicate a possibly nascent jet embedded in and interacting with a dense medium resulting in a jet bending.

The fast-rising blue optical transient AT2018cow indicated unusual early-phase characteristics unlike relatively better studied explosive transients. Its afterglow may be produced by either a relativistically beamed (jetted) or intrinsically luminous (nonjetted) ejecta and carries observational signatures of the progenitor and environment. High-resolution monitoring can distinguish between these scenarios and clarify the progenitor nature. We present very long baseline interferometry (VLBI) observations of AT2018cow at 5 GHz involving 21 radio telescopes from the European VLBI Network with five sessions spanning ≈1 yr. With an astrometric precision up to 25 μas per epoch, the rapidly fading compact mas-scale source is found to be nonjetted with a proper motion of ≤0.15 mas yr⁻¹ (0.14 c). This and a dense (number density ≈10⁴–10⁵ cm⁻³) magnetized environment (magnetic field strength ≥0.84 G) are characteristic of a newly formed magnetar driven central engine, originating in the successful explosion of a low-mass star.

We present the results from the 43 GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas × 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43 GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milliarcsecond and submilliarcsecond scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from the literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ∼7 GHz, and decreases between ∼7 and 240 GHz but increases again above 240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and γ-ray flux densities. Our correlation is tighter than those in the literature derived from lower-frequency VLBI data, suggesting that the γ-ray emission is produced more cospatially with the 43 GHz VLBA core emission. This correlation can also be extrapolated to the unbeamed AGN population, implying that a universal γ-ray production mechanism might be at work for all types of AGNs.

Dual active galactic nuclei and supermassive black hole binaries at kiloparsec and parsec‐scale separations, respectively, are expected during stages of galaxy merger and evolution. Their observational identification can address a range of areas of current astrophysics frontiers including the final parsec problem and their contribution toward the emission of low‐frequency gravitational waves. This has however been difficult to achieve with current spectroscopy and time domain strategies. Very long baseline interferometry as a method of directly imaging radio structures with milli‐arcsecond and submilli‐arcsecond resolutions is introduced as a possible means of detecting dual active galactic nuclei and supermassive black hole binaries. We motivate its usage with expected observational signatures and cite some studies from literature to illustrate its current status and present an updated list of candidates imaged with high‐resolution radio observations. We then recall some shortcomings of the method with possible solutions and discuss future directions, relevant to large surveys with the upcoming Square Kilometer Array and future space very long baseline interferometry missions.