We present deep near-infrared $K_\textrm{s}$-band imaging for 35 of the 53 sources from the high-redshift ($z \gt 2$) radio galaxy candidate sample defined in Broderick et al. (2022, PASA, 39, e061). These images were obtained using the High-Acuity Widefield K-band Imager (HAWK-I) on the Very Large Telescope. Host galaxies are detected for 27 of the sources, with $K_\textrm{s} \approx 21.6$–23.0 mag (2$^{\prime\prime}$ diameter apertures; AB). The remaining eight targets are not detected to a median $3\unicode{x03C3}$ depth of $K_\textrm{s} \approx 23.3$ mag (2$^{\prime\prime}$ diameter apertures). We examine the radio and near-infrared flux densities of the 35 sources, comparing them to the known $z \gt 3$ powerful radio galaxies with 500-MHz radio luminosities $L_{500\,\textrm{MHz}} \gt 10^{27}$ W Hz$^{-1}$. By plotting 150-MHz flux density versus $K_\textrm{s}$-band flux density, we find that, similar to the sources from the literature, these new targets have large radio to near-infrared flux density ratios, but extending the distribution to fainter flux densities. Five of the eight HAWK-I deep non-detections have a median $3\unicode{x03C3}$ lower limit of $K_\textrm{s} \gtrsim 23.8$ mag (1$.\!^{\prime\prime}$5 diameter apertures); these five targets, along with a further source from Broderick et al. (2022, PASA, 39, e061) with a deep non-detection ($K_\textrm{s} \gtrsim 23.7$ mag; $3\unicode{x03C3}$; 2$^{\prime\prime}$ diameter aperture) in the Southern H-ATLAS Regions $K_\textrm{s}$-band Survey, are considered candidates to be ultra-high-redshift ($z \gt 5$) radio galaxies. The extreme radio to near-infrared flux density ratios ($\gt 10^5$) for these six sources are comparable to TN J0924$-$2201, GLEAM J0856$+$0223 and TGSS J1530$+$1049, the three known powerful radio galaxies at $z \gt 5$. For a selection of galaxy templates with different stellar masses, we show that $z \gtrsim 4.2$ is a plausible scenario for our ultra-high-redshift candidates if the stellar mass $M_\textrm{*} \gtrsim 10^{10.5}$ M$_\odot$. In general, the 35 targets studied have properties consistent with the previously known class of infrared-faint radio sources. We also discuss the prospects for finding more UHzRG candidates from wide and deep near-infrared surveys.

This study presents an analysis of the optical variability of the blazar 1E 1458.8+2249 on diverse timescales using multi-band observations, including observations in the optical BVRI bands carried out with the T60 and T100 telescopes from 2020 to 2023 and ZTF gri data from 2018 to 2023. On seven nights, we searched for intraday variability using the power-enhanced F-test and the nested ANOVA test, but no significant variability was found. The long-term light curve shows a variability behaviour in the optical BVRI bands with amplitudes of $\sim$100% and in the gri bands with amplitudes of $\sim$120%, including short-term variability of up to $\sim$1.1 mag. Correlation analysis revealed a strong correlation between the optical multi-band emissions without any time lag. From 62 nightly spectral energy distributions, we obtained spectral indices between 0.826 and 1.360, with an average of $1.128\pm0.063$. The relationships of both spectral indices and colour with respect to brightness indicate a mild BWB trend throughout the observation period, both intraday and long-term. We also performed a periodicity search using the weighted wavelet Z-transform and Lomb–Scargle methods. A recurrent optical emission pattern with a quasi-periodicity of $\sim$340 days is detected in the combined V- and R-band light curves. The observational results indicate that the blazar 1E 1458.8+2249 has a complex variability, while emphasising the need for future observations to unravel its underlying mechanisms.

We present source detection and catalogue construction pipelines to build the first catalogue of radio galaxies from the 270 $\rm deg^2$ pilot survey of the Evolutionary Map of the Universe (EMU-PS) conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The detection pipeline uses Gal-DINO computer vision networks (Gupta et al. 2024, PASA, 41, e001) to predict the categories of radio morphology and bounding boxes for radio sources, as well as their potential infrared host positions. The Gal-DINO network is trained and evaluated on approximately 5 000 visually inspected radio galaxies and their infrared hosts, encompassing both compact and extended radio morphologies. We find that the Intersection over Union (IoU) for the predicted and ground-truth bounding boxes is larger than 0.5 for 99% of the radio sources, and 98% of predicted host positions are within $3^{\prime \prime}$ of the ground-truth infrared host in the evaluation set. The catalogue construction pipeline uses the predictions of the trained network on the radio and infrared image cutouts based on the catalogue of radio components identified using the Selavy source finder algorithm. Confidence scores of the predictions are then used to prioritise Selavy components with higher scores and incorporate them first into the catalogue. This results in identifications for a total of 211 625 radio sources, with 201 211 classified as compact and unresolved. The remaining 10 414 are categorised as extended radio morphologies, including 582 FR-I, 5 602 FR-II, 1 494 FR-x (uncertain whether FR-I or FR-II), 2 375 R (single-peak resolved) radio galaxies, and 361 with peculiar and other rare morphologies. Each source in the catalogue includes a confidence score. We cross-match the radio sources in the catalogue with the infrared and optical catalogues, finding infrared cross-matches for 73% and photometric redshifts for 36% of the radio galaxies. The EMU-PS catalogue and the detection pipelines presented here will be used towards constructing catalogues for the main EMU survey covering the full southern sky.

Odd Radio Circles (ORCs) are a class of low surface brightness, circular objects approximately one arcminute in diameter. ORCs were recently discovered in the Australian Square Kilometre Array Pathfinder (ASKAP) data and subsequently confirmed with follow-up observations on other instruments, yet their origins remain uncertain. In this paper, we suggest that ORCs could be remnant lobes of powerful radio galaxies, re-energised by the passage of a shock. Using relativistic hydrodynamic simulations with synchrotron emission calculated in post-processing, we show that buoyant evolution of remnant radio lobes is alone too slow to produce the observed ORC morphology. However, the passage of a shock can produce both filled and edge-brightnened ORC-like morphologies for a wide variety of shock and observing orientations. Circular ORCs are predicted to have host galaxies near the geometric centre of the radio emission, consistent with observations of these objects. Significantly offset hosts are possible for elliptical ORCs, potentially causing challenges for accurate host galaxy identification. Observed ORC number counts are broadly consistent with a paradigm in which moderately powerful radio galaxies are their progenitors.

We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\mathrm{M}_*\leq10^{10}\mathrm{M}_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe (EMU) early science observations, and Wide-field Infrared Survey Explorer (WISE), we compare this technique with a selection of different AGN diagnostics to explore the similarities and differences in AGN classification. We find that diagnostics based on optical and near-infrared criteria (the standard BPT diagram, the WISE colour criterion, and the mass-excitation, or MEx diagram) tend to favour detection of AGN in high mass, high luminosity systems, while the “ProSpect” SED fitting tool can identify AGN efficiently in low mass systems. We investigate an explanation for this result in the context of proportionally lower mass black holes in lower mass galaxies compared to higher mass galaxies and differing proportions of emission from AGN and star formation dominating the light at optical and infrared wavelengths as a function of galaxy stellar mass. We conclude that SED-derived AGN classification is an efficient approach to identify low mass hosts with low radio luminosity AGN.

Cold, neutral interstellar gas, the reservoir for star formation, is traced through the absorption of the 21-cm continuum radiation by neutral hydrogen (H i). Although detected in one hundred cases in the host galaxies of distant radio sources, only recently have column densities approaching the maximum value observed in Lyman-$\alpha$ absorption systems ($N_{{\textrm{H}\,\scriptsize{\textrm{I}}}}\sim 10^{22}$ $\textrm{cm}^{-2}$) been found. Here, we explore the implications these have for the hypothesis that the detection rate of H i absorption is dominated by photo-ionisation from the active galactic nucleus (AGN). We find, with the addition all of the current searches for H i absorption at $z\geq0.1$, a strong correlation between the H i absorption strength and the ionising photon rate, with the maximum value at which H i is detected remaining close to the theoretical value in which all of the neutral gas would be ionised in a large spiral galaxy ($Q_{{\textrm{H}\,\scriptsize{\textrm{I}}}} = 2.9\times10^{56}$ ionising photons s$^{-1}$). We also rule out other effects (excitation by the radio continuum and changing gas properties) as the dominant cause for the decrease in the detection rate with redshift. Furthermore, from the maximum theoretical column density, we find that the five high column density systems have spin temperatures close to those of the Milky Way ($T_{\textrm{spin}} \lesssim 300$ K), whereas, from our model of a gaseous galactic disc, the H i detection at $Q_{{\textrm{H}\,\scriptsize{\textrm{I}}}} =2.9\times10^{56}$ s$^{-1}$ yields $T_{\textrm{spin}}\sim10\,000$ K, consistent with the gas being highly ionised.

All very massive early-type galaxies contain supermassive blackholes, but are these blackholes all sufficiently active to produce detectable radio continuum sources? We have used the 887.5 MHz Rapid ASKAP Continuum Survey DR1 to measure the radio emission from morphological early-type galaxies brighter than $K_S=9.5$ selected from the 2MASS Redshift Survey, HyperLEDA, and RC3. In line with previous studies, we find median radio power increases with infrared luminosity, with $P_{1.4} \propto L_K^{2.2}$, although the scatter about this relation spans several orders of magnitude. All 40 of the $M_K<-25.7$ early-type galaxies in our sample have measured radio flux densities that are more than $2\sigma$ above the background noise, with $1.4\,{\rm GHz}$ radio powers spanning ${\sim} 3 \times 10^{20}$ to ${\sim} 3\times 10^{25}\,{\rm W/Hz^{-1}}$. Cross-matching our sample with integral field spectroscopy of early-type galaxies reveals that the most powerful radio sources preferentially reside in galaxies with relatively low angular momentum (i.e. slow rotators). While the infrared colours of most galaxies in our early-type sample are consistent with passive galaxies with negligible star formation and the radio emission produced by active galactic nuclei or AGN remnants, very low levels of star formation could power the weakest radio sources with little effect on many other star formation rate tracers.

Creating radio galaxy catalogues from next-generation deep surveys requires automated identification of associated components of extended sources and their corresponding infrared hosts. In this paper, we introduce RadioGalaxyNET, a multimodal dataset, and a suite of novel computer vision algorithms designed to automate the detection and localization of multi-component extended radio galaxies and their corresponding infrared hosts. The dataset comprises 4 155 instances of galaxies in 2 800 images with both radio and infrared channels. Each instance provides information about the extended radio galaxy class, its corresponding bounding box encompassing all components, the pixel-level segmentation mask, and the keypoint position of its corresponding infrared host galaxy. RadioGalaxyNET is the first dataset to include images from the highly sensitive Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope, corresponding infrared images, and instance-level annotations for galaxy detection. We benchmark several object detection algorithms on the dataset and propose a novel multimodal approach to simultaneously detect radio galaxies and the positions of infrared hosts.

We here present 0.02–0.04″ resolution ALMA observation of the compact obscured nucleus (CON) of IRAS 17578-0400. A dusty torus within the nucleus, approximately 4 pc in radius, has been uncovered, exhibiting a usually flat spectral index at ALMA band 3, likely due to the millimeter corona emission from the central supermassive black hole (SMBH). The dense gas disk, traced by 13CO(1-0), spans 7 pc in radius and suggests an outflow driven by a disk wind due to its asymmetrical structure along the minor axis. Collimated molecular outflows (CMO), traced by the low-velocity components of the HCN(3-2) and HCO+(3-2) lines, align with the minor axis gas disk. Examination of position-velocity plots of HCN(3-2) and HCO+(3-2) reveals a flared dense gas disk extended a radius of ∼60 pc, infalling and rotating at speeds of about 200 km s−1 and 300 km s−1 respectively. A centrifugal barrier, located around 4 pc from the dynamical center, implies an SMBH mass of approximately 108 Mȯ, consistent with millimeter corona emission estimates. The CMO maintains a steady rotation speed of 200 km s−1 over the 100 pc scale along the minor axis. The projected speed of the CMO is about 80 km s−1, corresponding to around ∼500 km s−1, assuming an inclination angle of 80°. Such a kinematics structure of disk-driven collimated rotating molecular outflow with gas supplies from a falling rotating disk indicates that the feedback of the compact obscured nucleus is likely regulated by the momentum transfer of the molecular gas that connects to both the feeding of the nuclear starburst and supermassive black hole.

We show that, contrary to simple predictions, most AGNs show at best only a small increase of lags with increasing wavelength in the J, H, K, and L bands. We suggest that a possible cause of this near simultaneity from the near-IR to the mid-IR is that the hot dust is in a hollow bi-conical outflow of which we preferentially see the near side. In the proposed model sublimation or re-creation of dust (with some delay relative luminosity variations) along our line of sight in the hollow cone as the flux varies could be a factor in explaining the AGN changing-look phenomenon (CL). Variations in the dust obscuration can help explain changes in relationship of Hβ time delay on Luv variability. The relative wavelength independence of IR lags simplifies the use of IR lags for estimating cosmological parameters.

Compact obscured nuclei (CONs) are relatively common in the centers of local (U)LIRGs, yet their nature remains unknown. Both AGN activity and extreme nuclear starbursts have been suggested as plausible nuclear power sources. The prevalence of outflows in these systems suggest that CONs represent a key phase in the nuclear feedback cycle, in which material is ejected from the central regions of the galaxy. Here, we present results from MUSE for the confirmed local CON galaxy NGC4418. For the first time we spatially map the spectral features and kinematics of the galaxy in the optical, revealing several previously unknown structures. In particular, we discover a bilateral outflow along the minor axis, an outflowing bubble, several knot structures and a receding outflow partially obscured by the galactic disk. Based on the properties of these features, we conclude that the CON in NGC4418 is most likely powered by an AGN.

Flows originating from black hole magnetospheres via Blandford-Znajek (BZ) process start highly relativistically, with very large Lorentz factors γ01, imprinted into the flow during pair production within the gaps. As a result, BZ-driven outflows would produce spine-brightened images, contrary to observations of the edge-brightened jet in M87. We conclude that M87 jet is not BZ-driven.

Ram-pressure stripping (RPS) is a process known to remove gas from satellite galaxies. Recent observational studies have found an increased ratio of active galactic nuclei (AGN) among the population of RPS galaxies compared to regular galaxies in the field. To test whether ram pressure (RP) can trigger an AGN, we perform a suite of hydrodynamical wind-tunnel simulations of a massive (Mstar = 1011 Mȯ) galaxy, with inclusion of star formation, stellar feedback and high resolution up to 39 pc. We find that RP increases the inflow of gas to the galaxy centre, which in turn can result in the enhanced BH accretion, as measured by the Bondi-Hoyle model. We also estimate pressure of outflows from our accretion rates and show that AGN feedback would play an important role on the early stages of stripping, while RP itself is not so strong.

The present work discusses the use of a weakly-supervised deep learning algorithm that reduces the cost of labelling pixel-level masks for complex radio galaxies with multiple components. The algorithm is trained on weak class-level labels of radio galaxies to get class activation maps (CAMs). The CAMs are further refined using an inter-pixel relations network (IRNet) to get instance segmentation masks over radio galaxies and the positions of their infrared hosts. We use data from the Australian Square Kilometre Array Pathfinder (ASKAP) telescope, specifically the Evolutionary Map of the Universe (EMU) Pilot Survey, which covered a sky area of 270 square degrees with an RMS sensitivity of 25–35 $\mu$Jy beam$^{-1}$. We demonstrate that weakly-supervised deep learning algorithms can achieve high accuracy in predicting pixel-level information, including masks for the extended radio emission encapsulating all galaxy components and the positions of the infrared host galaxies. We evaluate the performance of our method using mean Average Precision (mAP) across multiple classes at a standard intersection over union (IoU) threshold of 0.5. We show that the model achieves a mAP$_{50}$ of 67.5% and 76.8% for radio masks and infrared host positions, respectively. The network architecture can be found at the following link: https://github.com/Nikhel1/Gal-CAM

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.

We examine the long-term stability (on decade-like timescales) of optical ‘high polarisation’ (HP) state with ${p_{opt}}$ ${> 3\%}$, which commonly occurs in flat-spectrum (i.e., beamed) radio quasars (FSRQs) and is a prominent marker of blazar state. Using this clue, roughly a quarter of the FSRQ population has been reported to undergo HP $\leftrightarrow$ non-HP state transition on year-like timescales. This work examines the extent to which HP (i.e., blazar) state can endure in a FSRQ, despite these ‘frequent’ state transitions. This is the first attempt to verify, using purely opto-polarimetric data for a much enlarged sample of blazars, the recent curious finding that blazar state in individual quasars persists for at least a few decades, despite its changing/swinging observed fairly commonly on year-like timescales. The present analysis is based on a well-defined sample of 83 radio quasars, extracted from the opto-polarimetric survey RoboPol (2013–2017), for which old opto-polarimetric data taken prior to 1990 could be found in the literature. By a source-wise comparison of these two datasets of the same observable ($p_{opt}$), we find that $\sim$90% of the 63 quasars found in blazar state in our RoboPol sample, were also observed to be in that state about three decades before. On the other hand, within the RoboPol survey itself, we find that roughly a quarter of the blazars in our sample migrated to the other polarisation state on year-like timescales, by crossing the customary $p_{opt}$ = 3% threshold. Evidently, these relatively frequent transitions (in either direction) do not curtail the propensity of a radio quasar to retain its blazar (i.e., HP) state for at least a few decades. The observed transitions/swings of polarisation state are probably manifestation of transient processes, like ejections of synchrotron plasma blobs (VLBI radio knots) from the active nucleus.

Spectral variability offers a new technique to identify small scale structures from scintillation, as well as determining the absorption mechanism for peaked-spectrum (PS) radio sources. In this paper, we present very long baseline interferometry (VLBI) imaging using the long baseline array (LBA) of two PS sources, MRC 0225–065 and PMN J0322–4820, identified as spectrally variable from observations with the Murchison Widefield Array (MWA). We compare expected milliarcsecond structures based on the detected spectral variability with direct LBA imaging. We find MRC 0225–065 is resolved into three components, a bright core and two fainter lobes, roughly 430 pc projected separation. A comprehensive analysis of the magnetic field, host galaxy properties, and spectral analysis implies that MRC 0225–065 is a young radio source with recent jet activity over the last $10^2$–$10^3$ yr. We find PMN J0322–4820 is unresolved on milliarcsecond scales. We conclude PMN J0322–4820 is a blazar with flaring activity detected in 2014 with the MWA. We use spectral variability to predict morphology and find these predictions consistent with the structures revealed by our LBA images.

We present dual-frequency VLBI observations of a nearby radio galaxy NGC 4261 at 22 and 43 GHz using the East Asia VLBI Network. In particular, the first sub-pc scale image of the 22 GHz water megamaser line in the circumnuclear region of NGC 4261 is shown. Our results suggest that the megamaser emission in NGC 4261 can be associated with the inner radius of the obscuring disk, as it is proposed for the nearest radio-loud megamaser source NGC 1052. An alternative hypothesis on the megamaser association is the shock region of the interaction between the jet and ambient molecular clouds.

If water megamaser disk activity is intimately related to the circumnuclear activity from accreting supermassive black holes, a thorough understanding of the co-evolution of galaxies with their central black holes should consider the degree to which the maser production correlates with traits of their host galaxies. This contribution presents an investigation of multiwavelength nuclear and host properties of galaxies with and without water megamasers, that reveals a rather narrow multi-dimensional parameter space associated with the megamaser emission. This “goldilocks” region embodies the availability of gas, the degree of dusty obscuration and reprocessing of the central emission, the black hole mass, and the accretion rate, suggesting that the disk megamaser emission in particular is linked to a short-lived phase in the intermediate-mass galaxy evolution, providing new tools for both 1) further constraining the growth process of the incumbent AGN and its host galaxy, and 2) significantly boosting the maser disk detection by efficiently confining the 22 GHz survey parameters.

High resolution (0."26 × 0."13 (70 × 35 pc)) L-band (18 cm) OH megamaser (OHM) e-Merlin observations of the LIRG Zw049.057 show that the emission is emerging from a low velocity outflowing structure - which is foreground to a fast, dense and collimated molecular outflow detected by ALMA. The extremely dusty compact obscured nucleus (CON) of Zw049.057 has no (or only little) OHM emission associated with it - possibly because of too high number densities that quench the OHM. In contrast we detect 6 cm H2CO emission primarily from the CON-region. We suggest that the OHM-region of Zw049.057 is not directly associated with star formation, but instead occurs in a wide-angle, slow outflow that surrounds the fast and dense outflow. The OHM is pumped by IR emission that likely stems from activities in the nucleus. We briefly discuss how OHM emission can be used as a probe of LIRG-CON galaxies.