Search

In this work we examine the classification accuracy and robustness of a state-of-the-art semi-supervised learning (SSL) algorithm applied to the morphological classification of radio galaxies. We test if SSL with fewer labels can achieve test accuracies comparable to the supervised state-of-the-art and whether this holds when incorporating previously unseen data. We find that for the radio galaxy classification problem considered, SSL provides additional regularisation and outperforms the baseline test accuracy. However, in contrast to model performance metrics reported on computer science benchmarking data-sets, we find that improvement is limited to a narrow range of label volumes, with performance falling off rapidly at low label volumes. Additionally, we show that SSL does not improve model calibration, regardless of whether classification is improved. Moreover, we find that when different underlying catalogues drawn from the same radio survey are used to provide the labelled and unlabelled data-sets required for SSL, a significant drop in classification performance is observered, highlighting the difficulty of applying SSL techniques under dataset shift. We show that a class-imbalanced unlabelled data pool negatively affects performance through prior probability shift, which we suggest may explain this performance drop, and that using the Frechet Distance between labelled and unlabelled data-sets as a measure of data-set shift can provide a prediction of model performance, but that for typical radio galaxy data-sets with labelled sample volumes of O(1000), the sample variance associated with this technique is high and the technique is in general not sufficiently robust to replace a train-test cycle., Comment: Accepted to MNRAS. 14 pages

We present a novel natural language processing (NLP) approach to deriving plain English descriptors for science cases otherwise restricted by obfuscating technical terminology. We address the limitations of common radio galaxy morphology classifications by applying this approach. We experimentally derive a set of semantic tags for the Radio Galaxy Zoo EMU (Evolutionary Map of the Universe) project and the wider astronomical community. We collect 8,486 plain English annotations of radio galaxy morphology, from which we derive a taxonomy of tags. The tags are plain English. The result is an extensible framework which is more flexible, more easily communicated, and more sensitive to rare feature combinations which are indescribable using the current framework of radio astronomy classifications., Comment: 17 pages, 11 Figures, Accepted at MNRAS

We have carried out relativistic three-dimensional simulations of high-power radio sources propagating into asymmetric cluster environments. We offset the environment by 0 or 1 core radii (equal to 144 kpc), and incline the jets by 0, 15, or 45{\deg} away from the environment centre. The different environment encountered by each radio lobe provides a unique opportunity to study the effect of environment on otherwise identical jets. We find that the jets become unstable towards the end of the simulations, even with a Lorentz factor of 5; they nevertheless develop typical FR II radio morphology. The jets propagating into denser environments have consistently shorter lobe lengths and brighter hotspots, while the axial ratio of the two lobes is similar. We reproduce the recently reported observational anti-correlation between lobe length asymmetry and environment asymmetry, corroborating the notion that observed large-scale radio lobe asymmetry can be driven by differences in the underlying environment., Comment: 13 pages, 11 figures. Accepted for publication in MNRAS

We have investigated the evolution of the BL Lac object PKS B1144−379 using the University of Tasmania Ceduna 30-m radio telescope at a frequency of 6.7 GHz and very long baseline interferometry (VLBI) data at 8.6 GHz. Variability time-scales associated with two flares detected in 2005 November and 2008 August were derived from long-term variations in total flux density monitored by Ceduna between 2003 and 2011. A kinematic study of the parsec-scale jet of PKS B1144−379 was performed using VLBI data obtained between 1997 and 2018. Quasi-periodic flarings with a period of ∼3–4 yr were observed. Over the 20-yr interval, the average jet position angle was found to be ~150°.

The highly variable BL Lac object PKS B1144 − 379 was monitored at 6.7 GHz using the Ceduna Radio Telescope with high cadence from 2003 to 2011. Intraday variations due to interstellar scintillation (ISS) were observed throughout the period. To complement our earlier analysis of the ISS and variability of this object, we have investigated the physical origin of changes in the modality of flux density distributions, calculated for ∼14 d observing blocks. Our analysis shows that the flux density distribution is primarily bimodal, but it changes to unimodal during the core brightening and jet expansion phases of the source. The presence of unimodal flux density distributions during these two phases is most likely due to the compactness of the scintillating component and the intrinsic evolution of the source. The existence of unimodality in the flux density distributions associated with specific phases of the source evolution also suggests that changes in the modality are unlikely due to multiple scattering screens. We propose that the physical origin of changes in the modality of the flux density distribution for PKS B1144 − 379 is most likely due to the combination of multiple bright jet features with interstellar scintillation along the line of sight between observer and source. This new approach complements our previous investigations of the temporal evolution of PKS B1144 − 379 that used interstellar scintillation and very long baseline interferometry, and the combination of these techniques provides a crucial starting point for understanding the system.

We present an analytical model for the evolution of extended active galactic nuclei (AGNs) throughout their full lifecycle, including the initial jet expansion, lobe formation, and eventual remnant phases. A particular focus of our contribution is on the early jet expansion phase, which is traditionally not well captured in analytical models. We implement this model within the Radio AGN in Semi-Analytic Environments (RAiSE) framework, and find that the predicted radio source dynamics are in good agreement with hydrodynamic simulations of both low-powered Fanaroff-Riley Type-I and high-powered Type-II radio lobes. We construct synthetic synchrotron surface brightness images by complementing the original RAiSE model with the magnetic field and shock-acceleration histories of a set of Lagrangian tracer particles taken from an existing hydrodynamic simulation. We show that a single set of particles is sufficient for an accurate description of the dynamics and observable features of Fanaroff-Riley Type-II radio lobes with very different jet parameters and ambient density profile normalisations. Our new model predicts that the lobes of young (< 10 Myr) sources will be both longer and brighter than expected at the same age from existing analytical models which lack a jet-dominated expansion phase; this finding has important implications for interpretation of radio galaxy observations. The RAiSE code, written in Python, is publicly available on GitHub and PyPI., 22 pages, 11 figures, 1 table; accepted in MNRAS

Quantifying the energetics and lifetimes of remnant radio-loud active galactic nuclei (AGNs) is much more challenging than for active sources due to the added complexity of accurately determining the time since the central black hole switched off. Independent spectral modelling of remnant lobes enables the derivation of the remnant ratio, $R_\mathrm{rem}$, (i.e. `off-time/source age'), thus reducing this complexity back to that of active sources however, the requirement of high-frequency ($\gtrsim5\,$GHz) coverage makes the application of this technique over large-area radio surveys difficult. In this work we propose a new method, which relies on the observed brightness of backflow of Fanaroff-Riley type~II lobes, combined with the \emph{Radio AGN in Semi-Analytic Environments} (RAiSE) code, to measure the duration of the remnant phase. Sensitive radio observations of the remnant radio galaxy J2253-34 are obtained to provide a robust comparison of this technique with the canonical spectral analysis and modelling methods. We find that the remnant lifetimes modelled by each method are consistent; spectral modelling yields $R_\mathrm{rem} = 0.23\pm0.02$, compared to $R_\mathrm{rem} = 0.26\pm0.02$ from our new method. We examine the viability of applying our proposed technique to low-frequency radio surveys using mock radio source populations, and examine whether the technique is sensitive to any intrinsic properties of radio AGNs. Our results show that the technique can be used to robustly classify active and remnant populations, with the most confident predictions for the remnant ratio, and thus off-time, in the longest-lived radio sources ($>50~$Myr) and those at higher redshifts ($z > 0.1$)., 20 pages, 9 figures, 8 tables. Accepted for publication to MNRAS

Next-generation radio surveys are expected to detect tens of millions of active galactic nuclei (AGNs) with a median redshift of $z$ ≥ 1. Beyond targeted surveys, the vast majority of these objects will not have spectroscopic redshifts, while photometric redshifts for high-redshift AGNs are of limited quality, and even then require optical and infrared photometry. We propose a new approach to measure the redshifts of lobed radio galaxies based exclusively on radio-frequency imaging and broad-band radio photometry. Specifically, our algorithm uses the lobe flux density, angular size and width, and spectral shape to derive probability density functions for the most likely source redshift based on the Radio AGN in Semi-analytic Environments dynamical model. The full physically based model explains 70 per cent of the variation in the spectroscopic redshifts of a high-redshift (2 < $z$ < 4) sample of radio AGNs, compared to at most 27 per cent for any one of the observed attributes in isolation. We find that upper bounds on the angular size, as expected for unresolved sources, are sufficient to yield accurate redshift measurements at $z$ ≥ 2. The error in the model upon calibration using at least nine sources with known spectroscopic redshifts is

The University of Tasmania Ceduna radio telescope has been used to investigate rapid variability in the radio flux density of the BL Lac object PKS B1144-379 at 6.7 GHz. High-cadence monitoring of this extreme scintillator was carried out over a period of approximately nine years, between 2003 and 2011. We have used structure functions created from the intensity time series to determine the characteristic timescale of the variability. The characteristic timescale is consistently observed to increase during certain periods of each year, demonstrating the annual cycle expected for scintillation through an interstellar scattering screen. The best-fitting annual cycle model for each year suggests that the scintillation pattern has an anisotropic structure and that the upper limit of its scattering screen is at a distance of ~0.84 kpc. Higher anisotropy in some of the annual cycle fits suggests that changes in the intrinsic source structure might be influencing the variability timescale. We found a prominent annual cycle is only present in the variability timescale for certain years where other evidence suggests that the core is compact. From our measurements, we calculated that the core angular size varied between 5.65-15.90 $\mu$as (0.05-0.13 pc). The core component was found to be at its most compact during two flares in the total flux density, which were observed in 2005 and 2008. We conclude that the long-term variability in the radio flux density of PKS B1144-379 is due to intrinsic changes in the source and that these affect our ability to measure an annual cycle in its variability time scale., Comment: 24 pages, 22 figures

Feedback from radio jets associated with Active Galactic Nuclei (AGN) plays a profound role in the evolution of galaxies. Kinetic power of these radio jets appears to show temporal variation, but the mechanism(s) responsible for this process are not yet clear. Recently, the LOw Frequency ARray (LOFAR) has uncovered large populations of active, remnant and restarted radio jet populations. By focusing on LOFAR data in the Lockman Hole, in this work we use the Radio AGN in Semi-Analytic Environments (RAiSE) dynamical model to present the first self-consistent modelling analysis of active, remnant and restarted radio source populations. Consistent with other recent work, our models predict that remnant radio lobes fade quickly. Any high (>10 percent) observed fraction of remnant and restarted sources therefore requires a dominant population of short-lived jets. We speculate that this could plausibly be provided by feedback-regulated accretion., 14 pages, 7 figures, accepted for publication in MNRAS. Companion paper to Jurlin et al., also in today's arXiv posting

We model the X-ray surface brightness distribution of emission associated with Fanaroff & Riley type-II radio galaxies. Our approach builds on the RAiSE dynamical model which describes broadband radio-frequency synchrotron evolution of jet-inflated lobes in a wide range of environments. The X-ray version of the model presented here includes: (1) inverse-Compton upscattering of cosmic microwave background radiation; (2) the dynamics of the shocked gas shell and associated bremsstrahlung radiation; and (3) emission from the surrounding ambient medium. We construct X-ray surface brightness maps for a mock catalogue of extended FR-IIs based on the technical characteristics of the eRosita telescope. The integrated X-ray luminosity function at low redshifts ($z\leqslant1$) is found to strongly correlate with the density of the ambient medium in all but the most energetic sources, whilst at high-redshift ($z>1$) the majority of objects are dominated by inverse-Compton lobe emission due to the stronger cosmic microwave background radiation. By inspecting our mock spatial brightness distributions, we conclude that any extended X-ray detection can be attributed to AGN activity at redshifts $z\geqslant1$. We compare the expected detection rates of active and remnant high-redshift radio AGNs for eRosita and LOFAR, and future more sensitive surveys. We find that a factor of ten more remnants can be detected using X-ray wavelengths over radio frequencies at $z>2.2$, increasing to a factor of 100 for redshifts $z>3.1$., 16 pages, 10 figures, 1 table; accepted in MNRAS

We present a method for applying spatially resolved adiabatic and radiative loss processes to synthetic radio emission from hydrodynamic simulations of radio sources from active galactic nuclei (AGN). Lagrangian tracer particles, each representing an ensemble of electrons, are injected into simulations and the position, grid pressure, and time since the last strong shock are recorded. These quantities are used to track the losses of the electron packet through the radio source in a manner similar to the Radio AGN in Semi-analytic Environments (RAiSE) formalism, which uses global source properties to calculate the emissivity of each particle ex-situ. Freedom in the choice of observing parameters, including redshift, is provided through the post-processing nature of this approach. We apply this framework to simulations of jets in different environments, including asymmetric ones. We find a strong dependence of radio source properties on frequency and redshift, in good agreement with observations and previous modelling work. There is a strong evolution of radio spectra with redshift due to the more prominent inverse-Compton losses at high redshift. Radio sources in denser environments have flatter spectral indices, suggesting that spectral index asymmetry may be a useful environment tracer. We simulate intermediate Mach number jets that disrupt before reaching the tip of the lobe, and find that these retain an edge-brightened Fanaroff-Riley Type II morphology, with the most prominent emission remaining near the tip of the lobes for all environments and redshifts we study., 17 pages, 15 figures. Accepted for publication in MNRAS

We present the Cosmological Double Radio Active Galactic Nuclei (CosmoDRAGoN) project: a large suite of simulated AGN jets in cosmological environments. These environments sample the intra-cluster media of galaxy clusters that form in cosmological smooth particle hydrodynamics (SPH) simulations, which we then use as inputs for grid-based hydrodynamic simulations of radio jets. Initially conical jets are injected with a range of jet powers, speeds (both relativistic and non-relativistic), and opening angles; we follow their collimation and propagation on scales of tens to hundreds of kiloparsecs, and calculate spatially-resolved synthetic radio spectra in post-processing. In this paper, we present a technical overview of the project, and key early science results from six representative simulations which produce radio sources with both core- (Fanaroff-Riley Type I) and edge-brightened (Fanaroff-Riley Type II) radio morphologies. Our simulations highlight the importance of accurate representation of both jets and environments for radio morphology, radio spectra, and feedback the jets provide to their surroundings., Comment: 25 pages, 20 figures. Accepted for publication in PASA

Immensely bright quasars and radio-loud active galactic nuclei (AGNs) provide an enticing opportunity to construct standard candles detectable up to the very early universe. An analytic theory is proposed to measure the distance to powerful \citeauthor{FR+1974} type-II radio sources based on their integrated flux density across a broad range of radio frequencies, and the angular size and axis ratio of their synchrotron-emitting lobes. This technique can be used at low-redshift to construct absolute standard candles in conjunction with X-ray observations of the host cluster, or at high-redshift to measure the relative distances of objects and constrain the curvature of our universe. Distances calculated with this method are consistent for dissimilar objects at the same redshift; the two lobes of Cygnus A have flux densities, linear sizes and spectral break frequencies varying by between 15-35\% yet their fitted distances are the same to within 7\%. These distance estimates together yield a transverse comoving distance to Cygnus A of $261_{-55}^{+70}\rm\, Mpc$ corresponding to a Hubble constant of $H_0 = 64_{-13}^{+17}\rm\, km\, s^{-1}\, Mpc^{-1}$. Large samples of suitable FR-II sources could provide a measure of the Hubble constant independent of existing techniques such as the cosmic microwave background, baryon acoustic oscillations, and type 1a supernovae., Comment: 12 pages, 3 figures, 3 tables; accepted in MNRAS

We present the discovery of another Odd Radio Circle (ORC) with the Australian Square Kilometre Array Pathfinder (ASKAP) at 944 MHz. The observed radio ring, ORC J0102-2450, has a diameter of ~70 arcsec or 300 kpc, if associated with the central elliptical galaxy DES J010224.33-245039.5 (z ~ 0.27). Considering the overall radio morphology (circular ring and core) and lack of ring emission at non-radio wavelengths, we investigate if ORC J0102-2450 could be the relic lobe of a giant radio galaxy seen end-on or the result of a giant blast wave. We also explore possible interaction scenarios, for example, with the companion galaxy, DES J010226.15-245104.9, located in or projected onto the south-eastern part of the ring. We encourage the search for further ORCs in radio surveys to study their properties and origin., MNRAS Letters, accepted (5 pages, 3 figures)

Active galactic nuclei (AGN) are typically identified through radio, mid-infrared, or X-ray emission or through the presence of broad and/or narrow emission lines. AGN can also leave an imprint on a galaxy's spectral energy distribution (SED) through the re-processing of photons by the dusty torus. Using the SED fitting code ProSpect with an incorporated AGN component, we fit the far ultraviolet to far-infrared SEDs of $\sim$494,00 galaxies in the D10-COSMOS field and $\sim$230,000 galaxies from the GAMA survey. By combining an AGN component with a flexible star formation and metallicity implementation, we obtain estimates for the AGN luminosities, stellar masses, star formation histories, and metallicity histories for each of our galaxies. We find that ProSpect can identify AGN components in 91 per cent of galaxies pre-selected as containing AGN through narrow-emission line ratios and the presence of broad lines. Our ProSpect-derived AGN luminosities show close agreement with luminosities derived for X-ray selected AGN using both the X-ray flux and previous SED fitting results. We show that incorporating the flexibility of an AGN component when fitting the SEDs of galaxies with no AGN has no significant impact on the derived galaxy properties. However, in order to obtain accurate estimates of the stellar properties of AGN host galaxies, it is crucial to include an AGN component in the SED fitting process. We use our derived AGN luminosities to map the evolution of the AGN luminosity function for $0, Comment: Published in MNRAS

The remnant phase of a radio galaxy begins when the jets launched from an active galactic nucleus are switched off. To study the fraction of radio galaxies in a remnant phase, we take advantage of a$8.31$deg$^2$subregion of the GAMA 23 field which comprises of surveys covering the frequency range 0.1–9 GHz. We present a sample of 104 radio galaxies compiled from observations conducted by the Murchison Widefield Array (216 MHz), the Australia Square Kilometer Array Pathfinder (887 MHz), and the Australia Telescope Compact Array (5.5 GHz). We adopt an ‘absent radio core’ criterion to identify 10 radio galaxies showing no evidence for an active nucleus. We classify these as new candidate remnant radio galaxies. Seven of these objects still display compact emitting regions within the lobes at 5.5 GHz; at this frequency the emission is short-lived, implying a recent jet switch off. On the other hand, only three show evidence of aged lobe plasma by the presence of an ultra-steep-spectrum ($\alpha) and a diffuse, low surface brightness radio morphology. The predominant fraction of young remnants is consistent with a rapid fading during the remnant phase. Within our sample of radio galaxies, our observations constrain the remnant fraction to$4\%\lesssim f_{\mathrm{rem}} \lesssim 10\%$; the lower limit comes from the limiting case in which all remnant candidates with hotspots are simply active radio galaxies with faint, undetected radio cores. Finally, we model the synchrotron spectrum arising from a hotspot to show they can persist for 5–10 Myr at 5.5 GHz after the jets switch of—radio emission arising from such hotspots can therefore be expected in an appreciable fraction of genuine remnants.

Supermassive black holes at the centres of galaxies can cycle through periods of activity and quiescence. Characterising the duty cycle of active galactic nuclei is crucial for understanding the impact of the energy they release on the host galaxy. For radio AGN, this can be done by identifying dying (remnant) and restarted radio galaxies from their radio spectral properties. Using the combination of images at 1400 MHz produced by Apertif, the new phased-array feed receiver installed on the Westerbork Synthesis Radio Telescope, and images at 150 MHz provided by LOFAR, we have derived resolved spectral index images (at a resolution of ~15 arcsec) for all the sources within ~6 deg^2 area of the Lockman Hole region. We were able to select 15 extended radio sources with emission (partly or entirely) characterised by extremely steep spectral indices (steeper than 1.2). These objects represent radio sources in the remnant or the restarted phases of their life cycle. Our findings suggest this cycle to be relatively fast. They also show a variety of properties relevant for modelling the evolution of radio galaxies. For example, the restarted activity can occur while the remnant structure from a previous phase of activity is still visible. This provides constraints on the duration of the 'off' (dying) phase. In extended remnants with ultra-steep spectra at low frequencies, the activity likely stopped a few hundred megayears ago, and they correspond to the older tail of the age distribution of radio galaxies, in agreement with simulations of radio source evolution. We find remnant radio sources with a variety of structures (from double-lobed to amorphous), suggesting different types of progenitors. The present work sets the stage for exploiting low-frequency spectral index studies of extended sources by taking advantage of the large areas common to the LOFAR and the Apertif surveys., Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields. 17 pages, 10 figures

We present the data and initial results from the first Pilot Survey of the Evolutionary Map of the Universe (EMU), observed at 944 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The survey covers 270 \sqdeg of an area covered by the Dark Energy Survey, reaching a depth of 25--30 \ujybm\ rms at a spatial resolution of $\sim$ 11--18 arcsec, resulting in a catalogue of $\sim$ 220,000 sources, of which $\sim$ 180,000 are single-component sources. Here we present the catalogue of single-component sources, together with (where available) optical and infrared cross-identifications, classifications, and redshifts. This survey explores a new region of parameter space compared to previous surveys. Specifically, the EMU Pilot Survey has a high density of sources, and also a high sensitivity to low surface-brightness emission. These properties result in the detection of types of sources that were rarely seen in or absent from previous surveys. We present some of these new results here., Accepted by PASA

In this paper, we present the identification of five previously unknown giant radio galaxies (GRGs) using Data Release 1 of the Radio Galaxy Zoo citizen science project and a selection method appropriate to the training and validation of deep learning algorithms for new radio surveys. We associate one of these new GRGs with the brightest cluster galaxy (BCG) in the galaxy cluster GMBCG J251.67741+36.45295 and use literature data to identify a further 13 previously known GRGs as BCG candidates, increasing the number of known BCG GRGs by $\gt 60$ per cent. By examining local galaxy number densities for the number of all known BCG GRGs, we suggest that the existence of this growing number implies that GRGs are able to reside in the centres of rich (∼1014 M⊙) galaxy clusters and challenges the hypothesis that GRGs grow to such sizes only in locally underdense environments.

Radio galaxies are known to go through cycles of activity, where phases of apparent quiescence can be followed by repeated activity of the central supermassive black hole. A better understanding of this cycle is crucial for ascertaining the energetic impact that the jets have on the host galaxy, but little is known about it. We used deep LOFAR images at 150 MHz of the Lockman Hole extragalactic field to select a sample of 158 radio sources with sizes $> 60^{\prime\prime}$ in different phases of their jet life cycle. Using a variety of criteria (e.g. core prominence combined with low-surface brightness of the extended emission and steep spectrum of the central region) we selected a subsample of candidate restarted radio galaxies representing between 13% and 15% of the 158 sources of the main sample. We compare their properties to the rest of the sample, which consists of remnant candidates and active radio galaxies. Optical identifications and characterisations of the host galaxies indicate similar properties for candidate restarted, remnant, and active radio galaxies, suggesting that they all come from the same parent population. The fraction of restarted radio galaxies is slightly higher with respect to remnants, suggesting that the restarted phase can often follow after a relatively short remnant phase (the duration of the remnant phase being a few times 10$^{7}$ years). This confirms that the remnant and restarted phases are integral parts of the life cycle of massive elliptical galaxies. A preliminary investigation does not suggest a strong dependence of this cycle on the environment surrounding any given galaxy., 16 pages, 9 figures, 4 tables, A&A accepted, companion paper to Shabala et al., also in today arXiv posting

We present a detailed analysis of the radio galaxy PKS 2250-351, a giant of 1.2 Mpc projected size, its host galaxy, and its environment. We use radio data from the Murchison Widefield Array, the upgraded Giant Metre-wavelength Radio Telescope, the Australian Square Kilometre Array Pathfinder, and the Australia Telescope Compact Array to model the jet power and age. Optical and infra-red data come from the Galaxy And Mass Assembly (GAMA) survey and provide information on the host galaxy and environment. GAMA spectroscopy confirms that PKS 2250-351 lies at z=0.2115 in the irregular, and likely unrelaxed, cluster Abell 3936. We find its host is a massive, `red and dead' elliptical galaxy with negligible star formation but with a highly obscured active galactic nucleus dominating the mid-infrared emission. Assuming it lies on the local M-sigma relation it has an Eddington accretion rate of lambda_EDD~0.014. We find that the lobe-derived jet power (a time-averaged measure) is an order of magnitude greater than the hotspot-derived jet power (an instantaneous measure). We propose that over the lifetime of the observed radio emission (~300 Myr) the accretion has switched from an inefficient advection dominated mode to a thin-disc efficient mode, consistent with the decrease in jet power. We also suggest that the asymmetric radio morphology is due to its environment, with the host of PKS 2250-351 lying to the west of the densest concentration of galaxies in Abell 3936., Comment: 16 pages and 8 figures, accepted for publication in PASA

Supermassive black hole binaries may be detectable by an upcoming suite of gravitational wave experiments. Their binary nature can also be revealed by radio jets via a short-period precession driven by the orbital motion as well as the geodetic precession at typically longer periods. We have investigated Karl G. Jansky Very Large Array (VLA) and MERLIN radio maps of powerful jet sources for morphological evidence of geodetic precession. For perhaps the best studied source, Cygnus A, we find strong evidence for geodetic precession. Projection effects can enhance precession features, for which we find indications in strongly projected sources. For a complete sample of 33 3CR radio sources we find strong evidence for jet precession in 24 cases (73 per cent). The morphology of the radio maps suggests that the precession periods are of the order of 10^6 - 10^7 yr. We consider different explanations for the morphological features and conclude that geodetic precession is the best explanation. The frequently observed gradual jet angle changes in samples of powerful blazars can be explained by orbital motion. Both observations can be explained simultaneously by postulating that a high fraction of powerful radio sources have sub-parsec supermassive black hole binaries. We consider complementary evidence and discuss if any jetted supermassive black hole with some indication of precession could be detected as individual gravitational wave source in the near future. This appears unlikely, with the possible exception of M87., Comment: 24 pages, 36 figures, accepted by MNRAS

We have carried out numerical hydrodynamic simulations of radio jets from active galactic nuclei using the PLUTO simulation code, with the aim of investigating the effect of different environments and intermittency of energy injection on the resulting dynamics and observable properties of the jet-inflated lobes. Initially conical jets are simulated in poor group and cluster environments. We show that the environment into which a radio jet is propagating plays a large role in the resulting morphology, dynamics and observable properties of the radio source. The same jet collimates much later in a poor group compared to a cluster, which leads to pronounced differences in radio morphology. The intermittency of the jet also affects the observable properties of the radio source, and multiple hotspots are present for multiple outburst jets in the cluster environment. We quantify the detectability of active and quiescent phases, and find this to be strongly environment-dependent. We conclude that the dynamics and observational properties of jets depend strongly on the details of energy injection and environment., Comment: 22 pages, 24 figures. Accepted for publication in MNRAS

Sensitive continuum surveys with next-generation interferometers will characterise large samples of radio sources at epochs during which cosmological models predict feedback from radio jets to play an important role in galaxy evolution. Dynamical models of radio sources provide a framework for deriving from observations the radio jet duty cycles and energetics, and hence the energy budget available for feedback. Environment plays a crucial role in determining observable radio source properties, and I briefly summarise recent efforts to combine galaxy formation and jet models in a self-consistent framework. Galaxy clustering estimates from deep optical and NIR observations will provide environment measures needed to interpret the observed radio populations.

The optical Fundamental Plane of black hole activity relates radio continuum luminosity of Active Galactic Nuclei to [O III] luminosity and black hole mass. We examine the environments of low redshift ($z, 7 pages, 3 figures. MNRAS in press

This paper demonstrates a novel and efficient unsupervised clustering method with the combination of a Self-Organising Map (SOM) and a convolutional autoencoder. The rapidly increasing volume of radio-astronomical data has increased demand for machine learning methods as solutions to classification and outlier detection. Major astronomical discoveries are unplanned and found in the unexpected, making unsupervised machine learning highly desirable by operating without assumptions and labelled training data. Our approach shows SOM training time is drastically reduced and high-level features can be clustered by training on auto-encoded feature vectors instead of raw images. Our results demonstrate this method is capable of accurately separating outliers on a SOM with neighbourhood similarity and K-means clustering of radio-astronomical features complexity. We present this method as a powerful new approach to data exploration by providing a detailed understanding of the morphology and relationships of Radio Galaxy Zoo (RGZ) dataset image features which can be applied to new radio survey data., Accepted in Publications of the Astronomical Society of the Pacific, special issue on Machine Intelligence in Astronomy and Astrophysics. 23 pages, 8 full-page colour figures

Gaseous halos play a key role for understanding inflow, feedback and the overall baryon budget in galaxies. Literature models predict transitions of the state of the gaseous halo between cold and hot accretion, winds, fountains and hydrostatic halos at certain galaxy masses. Since luminosities of radio AGN are sensitive to halo densities, any significant transition would be expected to show up in the radio luminosities of large samples of galaxies. The Low Frequency Array (LOFAR) Two Metre Sky Survey (LoTSS) has indeed identified a galaxy stellar mass scale, $10^{11} M_\odot$ , above which the radio luminosities increase disproportionately. Here, we investigate, if radio luminosities of galaxies, especially the marked rise at galaxy masses around $10^{11} M_\odot$, can be explained with standard assumptions on jet powers, scaling between black hole-mass and galaxy mass and gaseous halos. We developed models for the radio luminosity of radio AGN in halos under infall, galactic wind and hydrostatic conditions based on observational data and theoretical constraints, and compared it to LoTSS data for a large sample of galaxies in the mass rangebetween $10^{8.5} M_\odot$ and $10^{12} M_\odot$. Assuming the same characteristic upper limit to jet powers as is known from high galaxy masses to hold at all masses, we find that the maximum radio luminosities for the hydrostatic gas halos fit well with the upper envelope of the distribution of the LOFAR data. The marked rise in radio luminosity at $10^{11} M_\odot$ is matched in our model, and is related to significant change in halo gas density around this galaxy mass, which is a consequence of the lower cooling rates at higher virial temperature. Wind and infall models overpredict the radio luminosities at small galaxy masses and have no particular steepening of the run of the radio luminosities predicted at any galaxy mass. [...], 6 pages, 3 figures, accepted by A&A, comments welcome

We study the impact of cluster environment on the morphology of a sample of 4304 extended radio galaxies from Radio Galaxy Zoo. A total of 87% of the sample lies within a projected 15 Mpc of an optically identified cluster. Brightest cluster galaxies (BCGs) are more likely than other cluster members to be radio sources, and are also moderately bent. The surface density as a function of separation from cluster center of non-BCG radio galaxies follows a power law with index $-1.10\pm 0.03$ out to $10~r_{500}$ ($\sim 7~$Mpc), which is steeper than the corresponding distribution for optically selected galaxies. Non-BCG radio galaxies are statistically more bent the closer they are to the cluster center. Within the inner $1.5~r_{500}$ ($\sim 1~$Mpc) of a cluster, non-BCG radio galaxies are statistically more bent in high-mass clusters than in low-mass clusters. Together, we find that non-BCG sources are statistically more bent in environments that exert greater ram pressure. We use the orientation of bent radio galaxies as an indicator of galaxy orbits and find that they are preferentially in radial orbits. Away from clusters, there is a large population of bent radio galaxies, limiting their use as cluster locators; however, they are still located within statistically overdense regions. We investigate the asymmetry in the tail length of sources that have their tails aligned along the radius vector from the cluster center, and find that the length of the inward-pointing tail is weakly suppressed for sources close to the center of the cluster., Comment: 23 pages, 17 figures, 2 tables. Supplemental data files available in The Astronomical Journal or contact author

We describe a physical model of the outflows produced as a result of gas accretion onto a black hole, and the resultant changes to star formation rates and efficiencies in galaxies, using the Radio-SAGE semi-analytic galaxy formation model. We show that the ratio of outflow rate to SFR of galaxies is mainly driven by black hole mass and virial halo mass, and show that the SFR is higher than the outflow rate at low black hole masses. The model consistently reproduces the observed evolution of star formation rate density from z = 6 to z = 0, as well as the trend of the stellar mass - halo mass relations. We show the characteristic growth of massive galaxies influenced by AGN feedback at different redshifts. We find feedback to be prevalent in the most massive galaxy halos, inhibiting the cooling catastrophe., 15 pages, 8 figures. Accepted for publication in MNRAS

We examine the processes triggering star formation and Active Galactic Nucleus (AGN) activity in a sample of 25 low redshift ($z10^7$ K) brightness temperature required for an mJIVE-20 detection allows us to unambiguously identify the radio AGN in our sample. We find three such objects. Our VLBI AGN identifications are classified as Seyferts or LINERs in narrow line optical diagnostic plots; mid-infrared colours of our targets and the comparison of H$��$ star formation rates with integrated radio luminosity are also consistent with the VLBI identifications. We reconstruct star formation histories in our galaxies using optical and UV photometry, and find that these radio AGN are not triggered promptly in the merger process, consistent with previous findings for non-VLBI samples of radio AGN. This delay can significantly limit the efficiency of feedback by radio AGN triggered in galaxy mergers. We find that radio AGN hosts have lower star formation rates than non-AGN radio-selected galaxies at the same starburst age. Conventional and VLBI radio imaging shows these AGN to be compact on arcsecond scales. Our modeling suggests that the actual sizes of AGN-inflated radio lobes may be much larger than this, but these are too faint to be detected in existing observations. Deep radio imaging is required to map out the true extent of the AGN, and to determine whether the low star formation rates in radio AGN hosts are a result of the special conditions required for radio jet triggering, or the effect of AGN feedback., 16 pages, 9 figures. Accepted for publication in MNRAS

The kinetic power of radio jets is a quantity of fundamental importance to studies of the AGN feedback process and radio galaxy physics. A widely used proxy for jet power is the extended radio luminosity. A number of empirical methods have been used to calibrate a scaling relationship between jet power (Q) and radio luminosity (L) of the form log(Q) = beta_L * log(L) + C. The regression slope has typically been found to be beta_L ~ 0.7 - 0.8. Here we show that the previously reported scaling relations are strongly affected by the confounding variable, distance. We find that in a sample of FRI X-ray cavity systems, after accounting for the mutual distance dependence, the jet power and radio luminosity are only weakly correlated, with slope beta_L ~ 0.3: significantly flatter than previously reported. We also find that in previously used samples of high-power sources, no evidence for an intrinsic correlation is present when the effect of distance is accounted for. Using a simple model we show that beta_L is expected to be significantly lower in samples of FRI radio galaxies than it is for FRIIs, due to the differing dynamics for these two classes of radio source. For FRI X-ray cavity systems the model predicts beta_L (FRI) ~ 0.5 in contrast to FRII radio galaxies, for which beta_L(FRII) ~ 0.8. We discuss the implications of our finding for studies of radio mode feedback, and radio galaxy physics., Comment: 14 pages, 6 figures, 2 tables. Accepted for publication in MNRAS

We investigate the role of environment on radio galaxy properties by constructing a sample of large ($\gtrsim100$~kpc), nearby ($z, 18 pages, 31 figures, 1 table; accepted in MNRAS

The upcoming next-generation large area radio continuum surveys can expect tens of millions of radio sources, rendering the traditional method for radio morphology classification through visual inspection unfeasible. We present ClaRAN - Classifying Radio sources Automatically with Neural networks - a proof-of-concept radio source morphology classifier based upon the Faster Region-based Convolutional Neutral Networks (Faster R-CNN) method. Specifically, we train and test ClaRAN on the FIRST and WISE images from the Radio Galaxy Zoo Data Release 1 catalogue. ClaRAN provides end users with automated identification of radio source morphology classifications from a simple input of a radio image and a counterpart infrared image of the same region. ClaRAN is the first open-source, end-to-end radio source morphology classifier that is capable of locating and associating discrete and extended components of radio sources in a fast (< 200 milliseconds per image) and accurate (>= 90 %) fashion. Future work will improve ClaRAN's relatively lower success rates in dealing with multi-source fields and will enable ClaRAN to identify sources on much larger fields without loss in classification accuracy., 22 pages, 16 figures, Accepted in Monthly Notices of the Royal Astronomical Society

We consider the problem of determining the host galaxies of radio sources by cross-identification. This has traditionally been done manually, which will be intractable for wide-area radio surveys like the Evolutionary Map of the Universe (EMU). Automated cross-identification will be critical for these future surveys, and machine learning may provide the tools to develop such methods. We apply a standard approach from computer vision to cross-identification, introducing one possible way of automating this problem, and explore the pros and cons of this approach. We apply our method to the 1.4 GHz Australian Telescope Large Area Survey (ATLAS) observations of the Chandra Deep Field South (CDFS) and the ESO Large Area ISO Survey South 1 (ELAIS-S1) fields by cross-identifying them with the Spitzer Wide-area Infrared Extragalactic (SWIRE) survey. We train our method with two sets of data: expert cross-identifications of CDFS from the initial ATLAS data release and crowdsourced cross-identifications of CDFS from Radio Galaxy Zoo. We found that a simple strategy of cross-identifying a radio component with the nearest galaxy performs comparably to our more complex methods, though our estimated best-case performance is near 100 per cent. ATLAS contains 87 complex radio sources that have been cross-identified by experts, so there are not enough complex examples to learn how to cross-identify them accurately. Much larger datasets are therefore required for training methods like ours. We also show that training our method on Radio Galaxy Zoo cross-identifications gives comparable results to training on expert cross-identifications, demonstrating the value of crowdsourced training data.

Gas-rich minor mergers contribute significantly to the gas reservoir of early-type galaxies (ETGs) at low redshift, yet the star formation efficiency (SFE; the star formation rate divided by the molecular gas mass) appears to be strongly suppressed following some of these events, in contrast to the more well-known merger-driven starbursts. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of six ETGs, which have each recently undergone a gas-rich minor merger, as evidenced by their disturbed stellar morphologies. These galaxies were selected because they exhibit extremely low SFEs. We use the resolving power of ALMA to study the morphology and kinematics of the molecular gas. The majority of our galaxies exhibit spatial and kinematical irregularities, such as detached gas clouds, warps, and other asymmetries. These asymmetries support the interpretation that the suppression of the SFE is caused by dynamical effects stabilizing the gas against gravitational collapse. Through kinematic modelling we derive high velocity dispersions and Toomre Q stability parameters for the gas, but caution that such measurements in edge-on galaxies suffer from degeneracies. We estimate merger ages to be about 100~Myr based on the observed disturbances in the gas distribution. Furthermore, we determine that these galaxies lie, on average, two orders of magnitude below the Kennicutt-Schmidt relation for star-forming galaxies as well as below the relation for relaxed ETGs. We discuss potential dynamical processes responsible for this strong suppression of star formation surface density at fixed molecular gas surface density., Comment: Accepted for publication in MNRAS

Thespatialstructureofsourcesmakingupthecelestialreferenceframe(CRF)atradiofrequencies is a systematic error source in very long baseline interferometry (VLBI) measurements. Using simulations, we investigate the effects of source structure on the CRF, determined by the actual observational programme of the International VLBI Service for Geodesy and Astrometry. This is done using the source structure simulator of the Vienna VLBI Software. Applying various mock two-component source models, systematic displacements of 10–80 μas in median source position offsets are found. These offsets are predominantly aligned with the direction of the jet. The simulations further show that slight changes in the source model can significantly change the estimated positions. We finally present a new parametrization of sourcepositionsintheanalysis,namelyalongthejetdirectionandperpendiculartoit,allowing us to significantly mitigate the effects of source structure on an estimated CRF.

Inhomogeneous networks and reference frames are an important issue for Very Long Baseline Interferometry (VLBI). In this work we examine the performance of southern stations and baselines in routine VLBI experiments. A positive impact on baseline length repeatabilities of the increased observing effort by the Australian AuScope VLBI antennas is found by analysing three years of global rapid-turnaround VLBI sessions: while worse results are found for southern baselines compared to northern baselines for the first half of the investigated sessions, the northern and southern baseline length repeatabilities are about the same in the second half of the period. In simulations, the actual observing plan with a significantly lower number of observations for southern stations is identified as a major reason for the worse length WRMS for southern baselines, though other factors seem to influence the results as well. Simulating radio source position uncertainties, effects of up to 10 mm are found on baseline length WRMS for long southern baselines. Improving all source position uncertainties through more frequent observations to better than 50 μ as could reduce this effect by up to 30 % .

We investigate the effects of quasar structure on geodetic very long baseline interferometry (VLBI) measurements. We create catalogues of simulated and real quasars with a range of structure indices, and use these to generate synthetic CONT11 observations with the Vienna VLBI Software simulator tool. We systematically investigate the effects of quasars with different amounts of source structure, and find that source structure can affect station positions at the one-millimetre level. This effect is stronger for isolated stations. Overall, source structure is found to contribute to about 10 % of the troposphere and clock effects. Our simulations confirm analytical predictions that source structure mitigation strategies must be developed in order to achieve millimetre-level VLBI position accuracy.

Monthly Notices of the Royal Astronomical Society, 453 (3), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Kinetic jet power estimates based exclusively on observed monochromatic radio luminosities are highly uncertain due to confounding variables and a lack of knowledge about some aspects of the physics of active galactic nuclei (AGNs). We propose a new methodology to calculate the jet powers of the largest, most powerful radio sources based on combinations of their size, lobe luminosity and shape of their radio spectrum; this approach avoids the uncertainties encountered by previous relationships. The outputs of our model are calibrated using hydrodynamical simulations and tested against independent X-ray inverse-Compton measurements. The jet powers and lobe magnetic field strengths of radio sources are found to be recovered using solely the lobe luminosity and spectral curvature, enabling the intrinsic properties of unresolved high-redshift sources to be inferred. By contrast, the radio source ages cannot be estimated without knowledge of the lobe volumes. The monochromatic lobe luminosity alone is incapable of accurately estimating the jet power or source age without knowledge of the lobe magnetic field strength and size respectively. We find that, on average, the lobes of the 3C radio sources have magnetic field strengths approximately a factor three lower than the equipartition value, inconsistent with equal energy in the particles and the fields at the 5$\sigma$ level. The particle content of 3C radio lobes is discussed in the context of complementary observations; we do not find evidence favouring an energetically-dominant proton population., Comment: 21 pages, 16 figures, 3 tables; accepted in MNRAS

We study the mutual alignment of radio sources within two surveys, Faint Images of the Radio Sky at Twenty-centimetres (FIRST) and TIFR GMRT Sky Survey (TGSS). This is done by producing two position angle catalogues containing the preferential directions of respectively 30 059 and 11 674 extended sources distributed over more than 7000 and 17 000 deg2. The identification of the sources in the FIRST sample was performed in advance by volunteers of the Radio Galaxy Zoo (RGZ) project, while for the TGSS sample it is the result of an automated process presented here. After taking into account systematic effects, marginal evidence of a local alignment on scales smaller than 2.5 deg is found in the FIRST sample. The probability of this happening by chance is found to be less than 2 per cent. Further study suggests that on scales up to 1.5 deg the alignment is maximal. For one third of the sources, the RGZ volunteers identified an optical counterpart. Assuming a flat Λ cold dark matter cosmology with Ωm=0.31,ΩΛ=0.69 Ωm=0.31,ΩΛ=0.69 , we convert the maximum angular scale on which alignment is seen into a physical scale in the range [19, 38] Mpc h−170 h70−1 . This result supports recent evidence reported by Taylor and Jagannathan of radio jet alignment in the 1.4 deg2 ELAIS N1 field observed with the Giant Metrewave Radio Telescope. The TGSS sample is found to be too sparsely populated to manifest a similar signal.

We model the triggering of Active Galactic Nuclei (AGN) in galaxy clusters using the semi- analytic galaxy formation model SAGE (?). We prescribe triggering methods based on the ram pressure galaxies experience as they move throughout the intracluster medium, which is hypothesized to trigger star formation and AGN activity. The clustercentric radius and velocity distribution of the simulated active galaxies produced by these models are compared with that of AGN and galaxies with intense star formation from a sample of low-redshift, relaxed clusters from the Sloan Digital Sky Survey. The ram pressure triggering model that best explains the clustercentric radius and velocity distribution of these observed galaxies has AGN and star formation triggered if $2.5\times10^{-14} < P_{ram} < 2.5\times10^{-13}$ Pa and $P_{ram} > 2P_{internal}$; this is consistent with expectations from hydrodynamical simulations of ram-pressure induced star formation. Our results show that ram pressure is likely to be an important mechanism for triggering star formation and AGN activity in clusters., 15 pages, 12 figures. Accepted for publication in MNRAS

We describe new efforts to model radio active galactic nuclei (AGN) in a cosmological context using the SAGE semi-analytic galaxy model. Our new method tracks the physical properties of radio jets in massive galaxies, including the evolution of radio lobes and their impact on the surrounding gas. This model also self consistently follows the gas cooling-heating cycle that significantly shapes star formation and the life and death of many galaxy types. Adding jet physics to SAGE adds new physical properties to the model output, which in turn allows us to make more detailed predictions for the radio AGN population. After calibrating the model to a set of core observations we analyse predictions for jet power, radio cocoon size, radio luminosity, and stellar mass. We find that the model is able to match the stellar mass--radio luminosity relation at $z\sim0$, and the radio luminosity function out to $z\sim1$. This updated model will make possible the construction of customised AGN-focused mock survey catalogues to be used for large-scale observing programs., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

With the advent of large scale-surveys the manual analysis and classification of individual radio source morphologies is rendered impossible as existing approaches do not scale. The analysis of complex morphological features in the spatial domain is a particularly important task. Here, we discuss the challenges of transferring crowdsourced labels obtained from the Radio Galaxy Zoo project and introduce a proper transfer mechanism via quantile random forest regression. By using parallelized rotation and flipping invariant Kohonen-maps, image cubes of Radio Galaxy Zoo selected galaxies formed from the Faint Images of the Radio Sky at Twenty-cm (FIRST) radio continuum and the Wide-field Infrared Survey Explorer (WISE) infrared all-sky surveys are first projected down to a two-dimensional embedding in an unsupervised way. This embedding can be seen as a discretized space of shapes with the coordinates reflecting morphological features as expressed by the automatically derived prototypes. We find that these prototypes have reconstructed physically meaningful processes across two channel images at radio and infrared wavelengths in an unsupervised manner. In the second step, images are compared with those prototypes to create a heat map, which is the morphological fingerprint of each object and the basis for transferring the user generated labels. These heat maps have reduced the feature space by a factor of 248, and are able to be used as the basis for subsequent machine-learning (ML) methods. Using an ensemble of decision trees we achieve upwards of 85.7% and 80.7% accuracy when predicting the number of components and peaks in an image, respectively, using these heat maps. We also question the currently used discrete classification schema and introduce a continuous scale that better reflects the uncertainty in transition between two classes, caused by sensitivity and resolution limits.

We examine the relationship between astrometric stability and astrophysical properties in flat-spectrum radio-loud quasars making up the celestial reference frame. We use position determinations from geodetic very long baseline interferometry measurements and develop a new metric for source position stability. We then compare this quantity to two measures of source compactness: structure index, which probes structure on milliarcsecond scales and the presence of interstellar scintillation associated with the quasar, which probes scales of tens of μas. We find that persistent scintillators have greater position stability than episodic scintillators, which are in turn more stable than non-scintillators. Scintillators are also more likely to be compact on milliarcsecond scales, as measured by the structure index. Persistent scintillators are therefore excellent candidates for inclusion in any future realization of the celestial reference frame. A list of these sources is presented in Appendix A. We find that slow (characteristic time-scale >3 d) scintillators have more stable positions than rapid scintillators, suggesting that they are more compact. High-cadence, long-term monitoring is therefore required to identify other members of this population of compact, high brightness temperature quasars.

The Australia Telescope Large Area Survey (ATLAS) has surveyed seven square degrees of sky around the Chandra Deep Field South (CDFS) and the European Large Area ISO Survey - South 1 (ELAIS-S1) fields at 1.4 GHz. ATLAS aims to reach a uniform sensitivity of $10 \mu$Jy beam$^{-1}$ rms over the entire region with data release 1 currently reaching $\sim30 \mu$Jy beam$^{-1}$ rms. Here we present 466 new spectroscopic redshifts for radio sources in ATLAS as part of our optical follow-up program. Of the 466 radio sources with new spectroscopic redshifts, 142 have star-forming optical spectra, 282 show evidence for AGN in their optical spectra, 10 have stellar spectra and 32 have spectra revealing redshifts, but with insufficient features to classify. We compare our spectroscopic classifications with two mid-infrared diagnostics and find them to be in broad agreement. We also construct the radio luminosity function for star-forming galaxies to z $= 0.5$ and for AGN to z $= 0.8$. The radio luminosity function for star-forming galaxies appears to be in good agreement with previous studies. The radio luminosity function for AGN appears higher than previous studies of the local AGN radio luminosity function. We explore the possibility of evolution, cosmic variance and classification techniques affecting the AGN radio luminosity function. ATLAS is a pathfinder for the forthcoming EMU survey and the data presented in this paper will be used to guide EMU's survey design and early science papers.