Search

We measured the harmonic-space power spectrum of galaxy clustering auto-correlation from the Evolutionary Map of the Universe Pilot Survey 1 data (EMU PS1) and its cross-correlation with the lensing convergence map of cosmic microwave background (CMB) from Planck Public Release 4 at the linear scale range from $\ell=2$ to 500. We applied two flux density cuts at $0.18$ and $0.4$mJy on the radio galaxies observed at 944MHz and considered two source detection algorithms. We found the auto-correlation measurements from the two algorithms at the 0.18mJy cut to deviate for $\ell\gtrsim250$ due to the different criteria assumed on the source detection and decided to ignore data above this scale. We report a cross-correlation detection of EMU PS1 with CMB lensing at $\sim$5.5$\sigma$, irrespective of flux density cut. In our theoretical modelling we considered two redshift distribution simulation models that yield consistent results, a linear and a non-linear matter power spectrum, and two linear galaxy bias models. That is a constant redshift-independent galaxy bias $b(z)=b_g$ and a constant amplitude galaxy bias $b(z)=b_g/D(z)$. By fixing a cosmology model and considering a non-linear matter power spectrum, we measured a constant galaxy bias at $0.18$mJy ($0.4$mJy) with $b_g=2.32^{+0.41}_{-0.33}$ ($2.18^{+0.17}_{-0.25}$) and a constant amplitude bias with $b_g=1.72^{+0.31}_{-0.21}$ ($1.78^{+0.22}_{-0.15}$). When $\sigma_8$ is a free parameter for the same models at $0.18$mJy ($0.4$mJy) with the constant model we found $\sigma_8=0.68^{+0.16}_{-0.14}$ ($0.82\pm0.10$), while with the constant amplitude model we measured $\sigma_8=0.61^{+0.18}_{-0.20}$ ($0.78^{+0.11}_{-0.09}$), respectively. Our results agree at $1\sigma$ with the measurements from Planck CMB and the weak lensing surveys and also show the potential of cosmology studies with future radio continuum survey data., Comment: To be submitted to PASA. 22 pages, 9 figures, 4 tables

We present the serendipitous discovery of a new radio-continuum ring-like object nicknamed Kyklos (J1802-3353), with MeerKAT UHF and L-band observations. The radio ring, which resembles the recently discovered odd radio circles (ORCs), has a diameter of 80 arcsec and is located just 6 deg from the Galactic plane. However, Kyklos exhibits an atypical thermal radio-continuum spectrum ({\alpha} = -0.1 +/- 0.3), which led us to explore different possible formation scenarios. We concluded that a circumstellar shell around an evolved massive star, possibly a Wolf-Rayet, is the most convincing explanation with the present data., Comment: 7 pages, 5 figures, accepted in A&A

Modern wide field radio surveys typically detect millions of objects. Techniques based on machine learning are proving to be useful for classifying large numbers of objects. The self-organizing map (SOM) is an unsupervised machine learning algorithm that projects a many-dimensional dataset onto a two- or three-dimensional lattice of neurons. This dimensionality reduction allows the user to visualize common features of the data better and develop algorithms for classifying objects that are not otherwise possible with large datasets. To this aim, we use the PINK implementation of a SOM. PINK incorporates rotation and flipping invariance so that the SOM algorithm may be applied to astronomical images. In this cookbook we provide instructions for working with PINK, including preprocessing the input images, training the model, and offering lessons learned through experimentation. The problem of imbalanced classes can be improved by careful selection of the training sample and increasing the number of neurons in the SOM (chosen by the user). Because PINK is not scale-invariant, structure can be smeared in the neurons. This can also be improved by increasing the number of neurons in the SOM. We also introduce pyink, a Python package used to read and write PINK binary files, assist in common preprocessing operations, perform standard analyses, visualize the SOM and preprocessed images, and create image-based annotations using a graphical interface. A tutorial is also provided to guide the user through the entire process. We present an application of PINK to VLA Sky Survey (VLASS) images. We demonstrate that the PINK is generally able to group VLASS sources with similar morphology together. We use the results of PINK to estimate the probability that a given source in the VLASS QuickLook Catalogue is actually due to sidelobe contamination., Comment: Accepted for publication in Astronomy and Computing

Generation of science-ready data from processed data products is one of the major challenges in next-generation radio continuum surveys with the Square Kilometre Array (SKA) and its precursors, due to the expected data volume and the need to achieve a high degree of automated processing. Source extraction, characterization, and classification are the major stages involved in this process. In this work we focus on the classification of compact radio sources in the Galactic plane using both radio and infrared images as inputs. To this aim, we produced a curated dataset of ~20,000 images of compact sources of different astronomical classes, obtained from past radio and infrared surveys, and novel radio data from pilot surveys carried out with the Australian SKA Pathfinder (ASKAP). Radio spectral index information was also obtained for a subset of the data. We then trained two different classifiers on the produced dataset. The first model uses gradient-boosted decision trees and is trained on a set of pre-computed features derived from the data, which include radio-infrared colour indices and the radio spectral index. The second model is trained directly on multi-channel images, employing convolutional neural networks. Using a completely supervised procedure, we obtained a high classification accuracy (F1-score>90%) for separating Galactic objects from the extragalactic background. Individual class discrimination performances, ranging from 60% to 75%, increased by 10% when adding far-infrared and spectral index information, with extragalactic objects, PNe and HII regions identified with higher accuracies. The implemented tools and trained models were publicly released, and made available to the radioastronomical community for future application on new radio data., Comment: 27 pages, 15 figures, 9 tables

Clusters of galaxies have been found to host Mpc-scale diffuse, non-thermal radio emission in the form of central radio halos and peripheral relics. Turbulence and shock-related processes in the intra-cluster medium are generally considered responsible for the emission, though details of these processes are still not clear. The low surface brightness makes detection of the emission a challenge, but with recent surveys with high-sensitivity radio telescopes we are beginning to build large samples of these sources. The Evolutionary Map of the Universe (EMU) is a Southern Sky survey being performed by the Australian SKA Pathfinder (ASKAP) over the next few years and is well-suited to detect and characterise such emission. To assess prospects of the full survey, we have performed a pilot search of diffuse sources in 71 clusters from the Planck Sunyaev-Zeldovich (SZ) cluster catalogue (PSZ2) found in archival ASKAP observations. After re-imaging the archival data and performing both (u,v)-plane and image-plane angular scale filtering, we detect 21 radio halos (12 for the first time, excluding an additional six candidates), 11 relics (in seven clusters, and six for the first time, excluding a further five candidate relics), along with 12 other, unclassified diffuse radio sources. From these detections, we predict the full EMU survey will uncover up to ~254 radio halos and ~85 radio relics in the 858 PSZ2 clusters that will be covered by EMU. The percentage of clusters found to host diffuse emission in this work is similar to the number reported in recent cluster surveys with the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (Botteon, et al. 2022a, A&A, 660, A78), suggesting EMU will complement similar searches being performed in the Northern Sky and provide us with statistically significant samples of halos and relics at the completion of the full survey., Comment: 46 pages, 116 individual figure files, accepted in PASA

We developed and trained a pipeline of three machine learning (ML) models than can predict which sources are more likely to be an AGN and to be detected in specific radio surveys. Also, it can estimate redshift values for predicted radio-detectable AGNs. These models, which combine predictions from tree-based and gradient-boosting algorithms, have been trained with multi-wavelength data from near-infrared-selected sources in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring field. Training, testing, calibration, and validation were carried out in the HETDEX field. Further validation was performed on near-infrared-selected sources in the Stripe 82 field. In the HETDEX validation subset, our pipeline recovers 96% of the initially labelled AGNs and, from AGNs candidates, we recover 50% of previously detected radio sources. For Stripe 82, these numbers are 94% and 55%. Compared to random selection, these rates are two and four times better for HETDEX, and 1.2 and 12 times better for Stripe 82. The pipeline can also recover the redshift distribution of these sources with $\sigma_{\mathrm{NMAD}}$ = 0.07 for HETDEX ($\sigma_{\mathrm{NMAD}}$ = 0.09 for Stripe 82) and an outlier fraction of 19% (25% for Stripe 82), compatible with previous results based on broad-band photometry. Feature importance analysis stresses the relevance of near- and mid-infrared colours to select AGNs and identify their radio and redshift nature. Combining different algorithms in ML models shows an improvement in the prediction power of our pipeline over a random selection of sources. Tree-based ML models (in contrast to deep learning techniques) facilitate the analysis of the impact that features have on the predictions. This prediction can give insight into the potential physical interplay between the properties of radio AGNs (e.g. mass of black hole and accretion rate)., Comment: Accepted for publication in A&A, 24 pages and 21 figures. Updated to include information on retrieval of data and models, which can be obtained from https://zenodo.org/doi/10.5281/zenodo.10220008

The 'Great Dimming' of the prototypical red supergiant Betelgeuse, which occurred between December 2019 and April 2020, gives us unprecedented insight into the processes occurring on the stellar surface and in the inner wind of this type of star. In particular it may bring further understanding of their dust nucleation and mass loss processes. Here, we present and analyse VLTI/MATISSE observations in the N-band (8 - 13 $\mu$m) taken near the brightness minimum in order to assess the status of the dusty circumstellar environment. We explore the compatibility of a dust clump obscuring the star with our mid-infrared interferometric observations using continuum 3D radiative transfer modelling, and probe the effect of adding multiple clumps close to the star on the observables. We also test the viability of a large cool spot on the stellar surface without dust present in the ambient medium. Using the visibility data, we derive a uniform disk diameter of 59.02 $\pm$ 0.64 mas in the spectral range 8 to 8.75 $\mu$m. We find that both the dust clump and the cool spot models are compatible with the data. Further to this, we note that the extinction and emission of our localised dust clump in the line of sight of the star, directly compensate each other making the clump undetectable in the spectral energy distribution and visibilities. The lack of infrared brightening during the 'Great Dimming' therefore does not exclude extinction due to a dust clump as one of the possible mechanisms. The visibilities can be reproduced by a spherical wind with dust condensing at 13 stellar radii and a dust mass-loss rate of (2.1 - 4.9) $\times$ 10$^{-10}$ $\mathit{M}_{\odot} {\rm yr}^{-1}$, however, in order to reproduce the complexity of the observed closure phases, additional surface features or dust clumps would be needed., Comment: 13 Pages, accepted for publication in A&A

We analyze the cooling and feedback properties of 48 galaxy clusters at redshifts $0.4 < z < 1.3$ selected from the South Pole Telescope (SPT) catalogs to evolve like the progenitors of massive and well-studied systems at $z{\sim}0$. We estimate the radio power at the brightest cluster galaxy (BCG) location of each cluster from an analysis of Australia Telescope Compact Array (ATCA) data. Assuming that the scaling relation between radio power and active galactic nucleus (AGN) cavity power $P_{\mathrm{cav}}$ observed at low redshift does not evolve with redshift, we use these measurements in order to estimate the expected AGN cavity power in the core of each system. We estimate the X-ray luminosity within the cooling radius $L_{\mathrm{cool}}$ of each cluster from a joint analysis of the available $Chandra$ X-ray and SPT Sunyaev-Zel'dovich (SZ) data. This allows us to characterize the redshift evolution of the $P_{\mathrm{cav}} / L_{\mathrm{cool}}$ ratio. When combined with low-redshift results, these constraints enable investigations of the properties of the feedback/cooling cycle across 9~Gyr of cluster growth. We model the redshift evolution of this ratio measured for cool core clusters by a log-normal distribution $\mathrm{Log}$-$\mathcal{N}(\alpha + \beta z, \sigma^2)$ and constrain the slope of the mean evolution $\beta = -0.05\pm 0.47$. This analysis improves the constraints on the slope of this relation by a factor of two. We find no evidence of redshift evolution of the feedback/cooling equilibrium in these clusters which suggests that the onset of radio-mode feedback took place at an early stage of cluster formation. High values of $P_{\mathrm{cav}} / L_{\mathrm{cool}}$ are found at the BCG location of non-cool core clusters which might suggest that the timescales of the AGN feedback cycle and the cool core / non-cool core transition are different., Comment: Submitted to The Astrophysical Journal

We present the discovery of highly-collimated radio jets spanning a total of 355 kpc around the nearby elliptical galaxy NGC 2663, and the possible first detection of recollimation on kiloparsec scales. The small distance to the galaxy (~28.5 Mpc) allows us to resolve portions of the jets to examine their structure. We combine multiwavelength data: radio observations by the Murchison Widefield Array (MWA), the Australian Square Kilometre Array Pathfinder (ASKAP) and the Australia Telescope Compact Array (ATCA), and X-ray data from Chandra, Swift and SRG/eROSITA. We present intensity, rotation measure, polarisation, spectral index and X-ray environment maps. Regions of the southern jet show simultaneous narrowing and brightening, which can be interpreted as a signature of the recollimation of the jet by external, environmental pressure, though it is also consistent with an intermittent Active Galactic Nuclei (AGN) or complex internal jet structure. X-ray data suggest that the environment is extremely poor; if the jet is indeed recollimating, the large recollimation scale (40 kpc) is consistent with a slow jet in a low-density environment., Comment: 16 pages, 14 figures, Accepted for publication in MNRAS

We present the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the Galaxy and Mass Assembly (GAMA)-23h field. The survey was carried out at 887.5 MHz and covers a 83 square degree field. We imaged the calibrated visibility data, taken as part of the Evolutionary Mapping of Universe (EMU) Early Science Programme, using the latest version of the ASKAPSoft pipeline. The final mosaic has an angular resolution of 10 arcsec and a central rms noise of around 38 $\mu$Jy beam$^{-1}$. The derived radio source catalogue has 39812 entries above a peak flux density threshold of 5$\sigma$. We searched for the radio source host galaxy counterparts using the GAMA spectroscopic (with an i-band magnitude limit of 19.2 mag) and multi-wavelength catalogues that are available as part of the collaboration. We identified hosts with GAMA spectroscopic redshifts for 5934 radio sources. We describe the data reduction, imaging, and source identification process, and present the source counts. Thanks to the wide area covered by our survey, we obtain very robust counts down to 0.2 mJy. ASKAP's exceptional survey speed, providing efficient, sensitive and high resolution mapping of large regions of the sky in conjunction with the multi-wavelength data available for the GAMA23 field, allowed us to discover 63 giant radio galaxies. The data presented here demonstrate the excellent capabilities of ASKAP in the pre-SKA era., Comment: 18 pages, 15 figures, accepted by MNRAS

Red supergiants are the most common final evolutionary stage of stars that have initial masses between 8 and 35 times that of the Sun. During this stage, which lasts roughly 100,000 years1, red supergiants experience substantial mass loss. However, the mechanism for this mass loss is unknown. Mass loss may affect the evolutionary path, collapse and future supernova light curve of a red supergiant, and its ultimate fate as either a neutron star or a black hole. From November 2019 to March 2020, Betelgeuse - the second-closest red supergiant to Earth (roughly 220 parsecs, or 724 light years, away) - experienced a historic dimming of its visible brightness. Usually having an apparent magnitude between 0.1 and 1.0, its visual brightness decreased to 1.614 +/- 0.008 magnitudes around 7-13 February 2020 - an event referred to as Betelgeuse's Great Dimming. Here we report high-angular-resolution observations showing that the southern hemisphere of Betelgeuse was ten times darker than usual in the visible spectrum during its Great Dimming. Observations and modelling support a scenario in which a dust clump formed recently in the vicinity of the star, owing to a local temperature decrease in a cool patch that appeared on the photosphere. The directly imaged brightness variations of Betelgeuse evolved on a timescale of weeks. Our findings suggest that a component of mass loss from red supergiants is inhomogeneous, linked to a very contrasted and rapidly changing photosphere

We report the discovery of J0624-6948, a low-surface brightness radio ring, lying between the Galactic Plane and the Large Magellanic Cloud (LMC). It was first detected at 888 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP), and with a diameter of ~196 arcsec. This source has phenomenological similarities to Odd Radio Circles (ORCs). Significant differences to the known ORCs - a flatter radio spectral index, the lack of a prominent central galaxy as a possible host, and larger apparent size - suggest that J0624-6948 may be a different type of object. We argue that the most plausible explanation for J0624-6948 is an intergalactic supernova remnant due to a star that resided in the LMC outskirts that had undergone a single-degenerate type Ia supernova, and we are seeing its remnant expand into a rarefied, intergalactic environment. We also examine if a massive star or a white dwarf binary ejected from either galaxy could be the supernova progenitor. Finally, we consider several other hypotheses for the nature of the object, including the jets of an active galactic nucleus (AGN) or the remnant of a nearby stellar super-flare., Comment: 20 pages accepted to MNRAS

The Shapley Supercluster ($\langle z \rangle\approx0.048$) contains several tens of gravitationally bound clusters and groups, making it it is an ideal subject for radio studies of cluster mergers. We used new high sensitivity radio observations to investigate the less energetic events of mass assembly in the Shapley Supercluster from supercluster down to galactic scales. We created total intensity images of the full region between A 3558 and A 3562, from $\sim 230$ to $\sim 1650$ MHz, using ASKAP, MeerKAT and the GMRT, with sensitivities ranging from $\sim 6$ to $\sim 100$ $\mu$Jy beam$^{-1}$. We performed a detailed morphological and spectral study of the extended emission features, complemented with ESO-VST optical imaging and X-ray data from XMM-Newton. We report the first GHz frequency detection of extremely low brightness intercluster diffuse emission on a $\sim 1$ Mpc scale connecting a cluster and a group, namely: A 3562 and the group SC 1329--313. It is morphologically similar to the X-ray emission in the region. We also found (1) a radio tail generated by ram pressure stripping in the galaxy SOS 61086 in SC 1329-313; (2) a head-tail radio galaxy, whose tail is broken and culminates in a misaligned bar; (3) ultrasteep diffuse emission at the centre of A 3558. Finally (4), we confirm the ultra-steep spectrum nature of the radio halo in A 3562. Our study strongly supports the scenario of a flyby of SC 1329-313 north of A 3562 into the supercluster core. [abridged...], Comment: 23 pages, 17 figures, A&A accepted

We present a comprehensive multi-frequency catalogue of radio sources behind the Large Magellanic Cloud between 0.2 and 20 GHz, gathered from a combination of new and legacy radio continuum surveys. This catalogue covers an area of $\sim$144~deg$^2$ at angular resolutions from 45 arcsec to $\sim$3 arcmin. We find 6434 discrete radio sources in total, of which 3789 are detected at two or more radio frequencies. We estimate the median spectral index ($\alpha$; where $S_{v}\sim\nu^\alpha$) of $\alpha = -0.89 $ and mean of $-0.88 \pm 0.48$ for 3636 sources detected exclusively at two frequencies (0.843 and 1.384 GHz) with similar resolution (FWHM $\sim$40-45 arcsec). The large frequency range of the surveys makes it an effective tool to investigate Gigahertz Peak Spectrum (GPS), Compact Steep Spectrum (CSS) and Infrared Faint Radio sources populations within our sample. We find 10 GPS candidates with peak frequencies near 5 GHz, from which we estimate their linear size. 1866 sources from our catalogue are (CSS) candidates with $\alpha <-0.8$. We found six candidates for High Frequency Peaker (HFP) sources, whose radio fluxes peak above 5 GHz and no sources with unconstrained peaks and $\alpha~>0.5$. We found optical counterparts for 343 of the radio continuum sources, of which 128have a redshift measurement. Finally, we investigate the population of 123 Infrared Faint Radio Sources (IFRSs) found in this study., Comment: Accepted for publication in MNRAS

We report a detailed CO(1-0) survey of a galaxy protocluster field at $z=2.16$, based on 475 hours of observations with the Australia Telescope Compact Array. We constructed a large mosaic of 13 individual pointings, covering an area of 21 arcmin$^2$ and $\pm6500$ km/s range in velocity. We obtain a robust sample of 46 CO(1-0) detections spanning $z=2.09-2.22$, constituting the largest sample of molecular gas measurements in protoclusters to date. The CO emitters show an overdensity at $z=2.12-2.21$, suggesting a galaxy super-protocluster or a protocluster connected to large-scale filaments with ~120 cMpc size. We find that 90% CO emitters have distances $>0'.5-4'$ to the center galaxy, indicating that small area surveys would miss the majority of gas reservoirs in similar structures. Half of the CO emitters have velocities larger than escape velocities, which appears gravitationally unbound to the cluster core. These unbound sources are barely found within the $R_{200}$ radius around the center, which is consistent with a picture in which the cluster core is collapsed while outer regions are still in formation. Compared to other protoclusters, this structure contains relatively more CO emitters with relatively narrow line width and high luminosity, indicating galaxy mergers. We use these CO emitters to place the first constraint on the CO luminosity function and molecular gas density in an overdense environment. The amplitude of the CO luminosity function is 1.6$\pm$0.5 orders of magnitudes higher than observed for field galaxy samples at $z\sim2$, and one order of magnitude higher than predictions for galaxy protoclusters from semi-analytical SHARK models. We derive a high molecular gas density of $0.6-1.3\times10^{9}$ $M_\odot$ cMpc$^{-3}$ for this structure, consistent with predictions for cold gas density of massive structures from hydro-dynamical DIANOGA simulations., Comment: 21 pages + Appendices, 9 figures, 6 tables. Accepted for publication in A&A

In this paper we describe the system design and capabilities of the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope at the conclusion of its construction project and commencement of science operations. ASKAP is one of the first radio telescopes to deploy phased array feed (PAF) technology on a large scale, giving it an instantaneous field of view that covers 31 square degrees at 800 MHz. As a two-dimensional array of 36x12m antennas, with baselines ranging from 22m to 6km, ASKAP also has excellent snapshot imaging capability and 10 arcsecond resolution. This, combined with 288 MHz of instantaneous bandwidth and a unique third axis of rotation on each antenna, gives ASKAP the capability to create high dynamic range images of large sky areas very quickly. It is an excellent telescope for surveys between 700 MHz and 1800 MHz and is expected to facilitate great advances in our understanding of galaxy formation, cosmology and radio transients while opening new parameter space for discovery of the unknown., Comment: 38 pages, 24 figures, accepted for publication in PASA

We present observations of a region of the Galactic plane taken during the Early Science Program of the Australian Square Kilometre Array Pathfinder (ASKAP). In this context, we observed the SCORPIO field at 912 MHz with an uncompleted array consisting of 15 commissioned antennas. The resulting map covers a square region of ~40 deg^2, centred on (l, b)=(343.5{\deg}, 0.75{\deg}), with a synthesized beam of 24"x21" and a background rms noise of 150-200 {\mu}Jy/beam, increasing to 500-600 {\mu}Jy/beam close to the Galactic plane. A total of 3963 radio sources were detected and characterized in the field using the CAESAR source finder. We obtained differential source counts in agreement with previously published data after correction for source extraction and characterization uncertainties, estimated from simulated data. The ASKAP positional and flux density scale accuracy were also investigated through comparison with previous surveys (MGPS, NVSS) and additional observations of the SCORPIO field, carried out with ATCA at 2.1 GHz and 10" spatial resolution. These allowed us to obtain a measurement of the spectral index for a subset of the catalogued sources and an estimated fraction of (at least) 8% of resolved sources in the reported catalogue. We cross-matched our catalogued sources with different astronomical databases to search for possible counterparts, finding ~150 associations to known Galactic objects. Finally, we explored a multiparametric approach for classifying previously unreported Galactic sources based on their radio-infrared colors., Comment: 31 pages, 15 figures, 15 tables

The pre-merging system of galaxy clusters Abell 3391-Abell 3395 located at a mean redshift of 0.053 has been observed at 1 GHz in an ASKAP/EMU Early Science observation as well as in X-rays with eROSITA. The projected separation of the X-ray peaks of the two clusters is $\sim$50$'$ or $\sim$ 3.1 Mpc. Here we present an inventory of interesting radio sources in this field around this cluster merger. While the eROSITA observations provide clear indications of a bridge of thermal gas between the clusters, neither ASKAP nor MWA observations show any diffuse radio emission coinciding with the X-ray bridge. We derive an upper limit on the radio emissivity in the bridge region of $\langle J \rangle_{1\,{\rm GHz}}< 1.2 \times 10^{-44} {\rm W}\, {\rm Hz}^{-1} {\rm m}^{-3}$. A non-detection of diffuse radio emission in the X-ray bridge between these two clusters has implications for particle-acceleration mechanisms in cosmological large-scale structure. We also report extended or otherwise noteworthy radio sources in the 30 deg$^2$ field around Abell 3391-Abell 3395. We identified 20 Giant Radio Galaxies, plus 7 candidates, with linear projected sizes greater than 1 Mpc. The sky density of field radio galaxies with largest linear sizes of $>0.7$ Mpc is $\approx 1.7$ deg$^{-2}$, three times higher than previously reported. We find no evidence for a cosmological evolution of the population of Giant Radio Galaxies. Moreover, we find seven candidates for cluster radio relics and radio halos., Comment: Astronomy & Astrophysics, in press

We used dedicated SRG/eROSITA X-ray, ASKAP/EMU radio, and DECam optical observations of a 15 sq.deg region around the interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite. We trace the irregular morphology of warm-hot gas of the main clusters from their centers out to well beyond their characteristic radii, $r_{200}$. Between the two main cluster systems, we observe an emission bridge; thanks to eROSITA's unique soft response and large field of view, we discover tantalizing hints for warm gas. Several matter clumps physically surrounding the system are detected. For the "Northern Clump," we provide evidence that it is falling towards A3391 from the hot gas morphology and radio lobe structure of its central AGN. Many of the extended sources in the field detected by eROSITA are known clusters or new clusters in the background, including a known SZ cluster at redshift z=1. We discover an emission filament north of the virial radius, $r_{100}$, of A3391 connecting to the Northern Clump and extending south of A3395 towards another galaxy cluster. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA observation. The DECam galaxy density map shows galaxy overdensities in the same regions. The new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution compared to the Magneticum simulation. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure., Comment: 21 pages plus 16 figures in the main text and 13 pages plus 29 figures as appendix. Astronomy & Astrophysics, accepted on November 3, 2020. A press release, full resolution images plus additional images and movies are available at https://astro.uni-bonn.de/~reiprich/A3391_95/

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the ``DeepDrill'' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centered on 3.6 $\mu$m and 4.5 $\mu$m. These observations expand the area which was covered by an earlier set of observations in these three fields by the Spitzer Extragalactic Representative Volume Survey (SERVS). The combined DeepDrill and SERVS data cover the footprints of the LSST DDFs in the Extended Chandra Deep Field-South field (ECDFS), the ELAIS-S1 field (ES1), and the XMM Large-Scale Structure Survey field (XMM-LSS). The observations reach an approximate $5\sigma$ point-source depth of 2 $\mu$Jy (corresponding to an AB magnitude of 23.1; sufficient to detect a 10$^{11} M_{\odot}$ galaxy out to $z\approx 5$) in each of the two bands over a total area of $\approx 29\,$deg$^2$. The dual-band catalogues contain a total of 2.35 million sources. In this paper we describe the observations and data products from the survey, and an overview of the properties of galaxies in the survey. We compare the source counts to predictions from the SHARK semi-analytic model of galaxy formation. We also identify a population of sources with extremely red ([3.6]$-$[4.5] $>1.2$) colours which we show mostly consists of highly-obscured active galactic nuclei., Comment: 20 pages, 9 figures; MNRAS in press

We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole., Comment: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

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 >60%. 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 centers of rich ($\sim 10^{14}$ M$_{\odot}$) galaxy clusters and challenges the hypothesis that GRGs grow to such sizes only in locally under-dense environments.

We develop a new analysis approach towards identifying related radio components and their corresponding infrared host galaxy based on unsupervised machine learning methods. By exploiting PINK, a self-organising map algorithm, we are able to associate radio and infrared sources without the a priori requirement of training labels. We present an example of this method using $894,415$ images from the FIRST and WISE surveys centred towards positions described by the FIRST catalogue. We produce a set of catalogues that complement FIRST and describe 802,646 objects, including their radio components and their corresponding AllWISE infrared host galaxy. Using these data products we (i) demonstrate the ability to identify objects with rare and unique radio morphologies (e.g. 'X'-shaped galaxies, hybrid FR-I/FR-II morphologies), (ii) can identify the potentially resolved radio components that are associated with a single infrared host and (iii) introduce a "curliness" statistic to search for bent and disturbed radio morphologies, and (iv) extract a set of 17 giant radio galaxies between 700-1100 kpc. As we require no training labels, our method can be applied to any radio-continuum survey, provided a sufficiently representative SOM can be trained.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon., Comment: 19 pages, 9 figures, ApJ in press

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA)., Comment: Accepted for publication in ApSS (38 pages, 14 figures), see also https://www.atnf.csiro.au/research/WALLABY/ - Contact email: Baerbel.Koribalski@csiro.au

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

We present a new image of the 9.0 GHz radio emission from the extended Chandra Deep Field South. A total of 181 hours of integration with the Australia Telescope Compact Array has resulted in a 0.276 square degree image with a median sensitivity of $\sim$20 $\mu$Jy/beam rms, for a synthesised beam of 4.0 $\times$ 1.3 arcsec. We present a catalogue of the 9.0 GHz radio sources, identifying 70 source components and 55 individual radio galaxies. Source counts derived from this sample are consistent with those reported in the literature. The observed source counts are also generally consistent with the source counts from simulations of the faint radio population. Using the wealth of multiwavelength data available for this region, we classify the faint 9 GHz population and find that 91% are radio loud AGN, 7% are radio quiet AGN and 2% are star forming galaxies. The 9.0 GHz radio sources were matched to 5.5 and 1.4 GHz sources in the literature and we find a significant fraction of flat or inverted spectrum sources, with 36% of the 9 GHz sources having $\alpha_{5.5GHz}^{9.0GHz}$ $>$ -0.3 (for $S \propto \nu^\alpha$). This flat or inverted population is not well reproduced by current simulations of radio source populations., Comment: 17 pages, 12 figures. Accepted for publication in MNRAS

We present a catalogue of a large sample of extended radio sources in the SCORPIO field, observed and resolved by the Australia Telescope Compact Array. SCORPIO, a pathfinder project for addressing the early operations of the Australia SKA Pathfinder, is a survey of ~5 square degrees between 1.4 and 3.1 GHz, centered at l=343.5{\deg}, b=0.75{\deg} and with an angular resolution of about 10 arcsec. It is aimed at understanding the scientific and technical challenges to be faced by future Galactic surveys. With a mean sensitivity around 100 $\mu$Jy/beam and the possibility to recover angular scales at least up to 4 arcmin, we extracted 99 extended sources, 35 of them detected for the first time. Among the 64 known sources 55 had at least a tentative classification in literature. Studying the radio morphology and comparing the radio emission with infrared we propose as candidates 6 new H II regions, 2 new planetary nebulae, 2 new luminous blue variable or Wolf--Rayet stars and 3 new supernova remnants. This study provides an overview of the potentiality of future radio surveys in terms of Galactic source extraction and characterization and a discussion on the difficulty to reduce and analyze interferometric data on the Galactic plane.

We report the discovery of a very young high-mass X-ray binary (HMXB) system associated with the supernova remnant (SNR) MCSNRJ0513-6724 in the Large Magellanic Cloud (LMC), using XMM-Newton X-ray observations. The HMXB is located at the geometrical centre of extended soft X-ray emission, which we confirm as an SNR. The HMXB spectrum is consistent with an absorbed power law with spectral index ~1.6 and a luminosity of 7x10^{33} ergs/s (0.2--12 keV). Tentative X-ray pulsations are observed with a periodicity of 4.4 s and the OGLE I-band light curve of the optical counterpart from more than 17.5 years reveals a period of 2.2324\pm0.0003 d, which we interpret as the orbital period of the binary system. The X-ray spectrum of the SNR is consistent with non-equilibrium shock models as expected for young/less evolved SNRs. From the derived ionisation time scale we estimate the age of the SNR to be <6 kyr. The association of the HMXB with the SNR makes it the youngest HMXB, in the earliest evolutionary stage known to date. A HMXB as young as this can switch on as an accreting pulsar only when the spin period has reached a critical value. Under this assumption, we obtain an upper limit to the magnetic field of < 5x10^{11} G. This implies several interesting possibilities including magnetic field burial, possibly by an episode of post-supernova hyper-critical accretion. Since these fields are expected to diffuse out on a timescale of 10^{3}-10^{4} years, the discovery of a very young HMXB can provide us the unique opportunity to observe the evolution of the observable magnetic field for the first time in X-ray binaries., Comment: 9 pages, 8 figures, accepted for publication in MNRAS

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories., Comment: 38 pages, 12 figures, accepted for publication in PASA

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies., Comment: 18 pages, 23 figues, 9 tables, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

A new era in radioastronomy will begin with the upcoming large-scale surveys planned at the Australian Square Kilometre Array Pathfinder (ASKAP). ASKAP started its Early Science program in October 2017 and several target fields were observed during the array commissioning phase. The SCORPIO field was the first observed in the Galactic Plane in Band 1 (792-1032 MHz) using 15 commissioned antennas. The achieved sensitivity and large field of view already allow to discover new sources and survey thousands of existing ones with improved precision with respect to previous surveys. Data analysis is currently ongoing to deliver the first source catalogue. Given the increased scale of the data, source extraction and characterization, even in this Early Science phase, have to be carried out in a mostly automated way. This process presents significant challenges due to the presence of extended objects and diffuse emission close to the Galactic Plane. In this context we have extended and optimized a novel source finding tool, named CAESAR , to allow extraction of both compact and extended sources from radio maps. A number of developments have been done driven by the analysis of the SCORPIO map and in view of the future ASKAP Galactic Plane survey. The main goals are the improvement of algorithm performances and scalability as well as of software maintainability and usability within the radio community. In this paper we present the current status of CAESAR and report a first systematic characterization of its performance for both compact and extended sources using simulated maps. Future prospects are discussed in light of the obtained results., Comment: 15 pages, 10 figures

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web., Comment: Published online in Science 27 June 2019

We have observed the G23 field of the Galaxy And Mass Assembly (GAMA) survey using the Australian Square Kilometre Array Pathfinder (ASKAP) in its commissioning phase, to validate the performance of the telescope and to characterize the detected galaxy populations. This observation covers $\sim$48 deg$^2$ with synthesized beam of 32.7$^{\prime\prime}$ by 17.8$^{\prime\prime}$ at 936 MHz, and $\sim$39 deg$^2$ with synthesized beam of 15.8$^{\prime\prime}$ by 12.0$^{\prime\prime}$ at 1320 MHz. At both frequencies, the r.m.s. (root-mean-square) noise is $\sim$0.1 mJy/beam. We combine these radio observations with the GAMA galaxy data, which includes spectroscopy of galaxies that are i-band selected with a magnitude limit of 19.2. Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry is used to determine which galaxies host an active galactic nucleus (AGN). In properties including source counts, mass distributions, and IR vs. radio luminosity relation, the ASKAP detected radio sources behave as expected. Radio galaxies have higher stellar mass and luminosity in IR, optical and UV than other galaxies. We apply optical and IR AGN diagnostics and find that they disagree for $\sim$30% of the galaxies in our sample. We suggest possible causes for the disagreement. Some cases can be explained by optical extinction of the AGN, but for more than half of the cases we do not find a clear explanation. Radio sources are more likely ($\sim$6%) to have an AGN than radio quiet galaxies ($\sim$1%), but the majority of AGN are not detected in radio at this sensitivity., Comment: 26 pages, 22 figures, 6 tables. Accepted for publication in Publications of the Astronomical Society of Australia

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 FIRST radio continuum and 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 discretised 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 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 report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age between 28-10 kyr. The radio spectral index for this SNR PWN pulsar system is flat (-0.29 $\pm$ 0.01) consistent with other similar objects. We infer that the putative pulsar has a radio spectral index of -1.8, which is typical for Galactic pulsars. We searched for dispersion measures (DMs) up to 1000 cm/pc^3 but found no convincing candidates with a S/N greater than 8. We produce a polarisation map for this PWN at 5500 MHz and find a mean fractional polarisation of P $\sim 23$ percent. The X-ray power-law spectrum (Gamma $\sim 2$) is indicative of non-thermal synchrotron emission as is expected from PWN-pulsar system. Finally, we detect DEM S5 in Infrared (IR) bands. Our IR photometric measurements strongly indicate the presence of shocked gas which is expected for SNRs. However, it is unusual to detect such IR emission in a SNR with a supersonic bow-shock PWN. We also find a low-velocity HI cloud of $\sim 107$ km/s which is possibly interacting with DEM S5. SNR DEM S5 is the first confirmed detection of a pulsar-powered bow shock nebula found outside the Galaxy., Comment: 20 pages, 15 figures, 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<0.3$) radio sources identified as part of the Radio Galaxy Zoo citizen science project. Our sample consists of 16 Fanaroff-Riley Type II (FR-II) sources, 6 FR-I sources, and one source with a hybrid morphology. FR-I sources appear to be hosted by more massive galaxies, consistent with previous studies. In the FR-II sample, we compare the degree of asymmetry in radio lobe properties to asymmetry in the radio source environment, quantified through optical galaxy clustering. We find that the length of radio lobes in FR-II sources is anti-correlated with both galaxy clustering and lobe luminosity. These results are in quantitative agreement with predictions from radio source dynamical models, and suggest that galaxy clustering provides a useful proxy for the ambient gas density distribution encountered by the radio lobes., Comment: 18 pages, 31 figures, 1 table; accepted in MNRAS

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.

We present very long baseline interferometry (VLBI) observations of 179 radio sources in the COSMOS field with extremely high sensitivity using the Green Bank Telescope (GBT) together with the Very Long Baseline Array (VLBA) (VLBA+GBT) at 1.4 GHz, to explore the faint radio population in the flux density regime of tens of $\mu$Jy. Here, the identification of active galactic nuclei (AGN) is based on the VLBI detection of the source, i.e., it is independent of X-ray or infrared properties. The milli-arcsecond resolution provided by the VLBI technique implies that the detected sources must be compact and have large brightness temperatures, and therefore they are most likely AGN (when the host galaxy is located at z$\geq$0.1). On the other hand, this technique allows us to only positively identify when a radio-active AGN is present, i.e., we cannot affirm that there is no AGN when the source is not detected. For this reason, the number of identified AGN using VLBI should be always treated as a lower limit. We present a catalogue containing the 35 radio sources detected with the VLBA+GBT, 10 of which were not previously detected using only the VLBA. We have constructed the radio source counts at 1.4 GHz using the samples of the VLBA and VLBA+GBT detected sources of the COSMOS field to determine a lower limit for the AGN contribution to the faint radio source population. We found an AGN contribution of >40-75% at flux density levels between 150 $\mu$Jy and 1 mJy. This flux density range is characterised by the upturn of the Euclidean-normalised radio source counts, which implies a contribution of a new population. This result supports the idea that the sub-mJy radio population is composed of a significant fraction of radio-emitting AGN, rather than solely by star-forming galaxies, in agreement with previous studies., Comment: 14 pages, 8 figures, 4 tables. Accepted for publication in A&A

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts., Comment: MNRAS, accepted. 34 pages, 25 figures, and 8 tables. The data products are available at this http url: http://personal.psu.edu/wnb3/xmmservs/xmmservs.html

The circumgalactic medium (CGM) of the massive Spiderweb Galaxy, a conglomerate of merging proto-cluster galaxies at z=2.2, forms an enriched interface where feedback and recycling act on accreted gas. This is shown by observations of [CI], CO(1-0) and CO(4-3) performed with the Atacama Large Millimeter Array (ALMA) and Australia Telescope Compact Array (ATCA). [CI] and CO(4-3) are detected across ~50 kpc, following the distribution of previously detected low-surface-brightness CO(1-0) across the CGM. This confirms our previous results on the presence of a cold molecular halo. The central radio galaxy MRC1138-262 shows a very high global $L'_{\rm CO(4-3)}$/$L'_{\rm CO(1-0)}$ ~ 1, suggesting that mechanisms other than FUV-heating by star formation prevail at the heart of the Spiderweb Galaxy. Contrary, the CGM has $L'_{\rm CO(4-3)}$/$L'_{\rm CO(1-0)}$ and $L'_{\rm [CI]}$/$L'_{\rm CO(1-0)}$ similar to the ISM of five galaxies in the wider proto-cluster, and its carbon abundance, $X_{\rm [CI]}$/$X_{\rm H2}$, resembles that of the Milky Way and starforming galaxies. The molecular CGM is thus metal-rich and not diffuse, confirming a link between the cold gas and in-situ star formation. Thus, the Spiderweb Galaxy grows not directly through accretion of gas from the cosmic web, but from recycled gas in the GCM., Comment: Accepted for publication as a Letter in MNRAS (5 pages, 4 figures)

Context. Red supergiant stars are one of the latest stages in the evolution of massive stars. Their photospheric convection may play an important role in the launching mechanism of their mass loss; however, its characteristics and dynamics are still poorly constrained. Aims. By observing red supergiant stars with near infrared interferometry at different epochs, we expect to reveal the evolution of bright convective features on their stellar surface. Methods. We observed the M2Iab-Ib red supergiant star CE Tau with the VLTI/PIONIER instrument in the H band at two different epochs separated by one month. Results. We derive the angular diameter of the star and basic stellar parameters, and reconstruct two reliable images of its H-band photosphere. The contrast of the convective pattern of the reconstructed images is $5 \pm 1 \%$ and $6 \pm 1 \%$ for our two epochs of observation. Conclusions. The stellar photosphere shows few changes between the two epochs. The contrast of the convective pattern is below the average contrast variations obtained on 30 randomly chosen snapshots of the best matching 3D radiative hydrodynamics simulation: $23 \pm 1~\%$ for the original simulation images and $16 \pm 1~\%$ for the maps degraded to the reconstruction resolution. We offer two hypotheses to explain this observation. CE Tau may be experiencing a quiet convective activity episode or it could be a consequence of its warmer effective temperature (hence its smaller radius) compared to the simulation., Comment: 12 pages, 8 figures. Accepted for publication in Astronomy & Astrophysics

Machine learning techniques have been increasingly useful in astronomical applications over the last few years, for example in the morphological classification of galaxies. Convolutional neural networks have proven to be highly effective in classifying objects in image data. The current work aims to establish when multiple components are present, in the astronomical context of synthesis imaging observations of radio sources. To this effect, we design a convolutional neural network to differentiate between different morphology classes using sources from the Radio Galaxy Zoo (RGZ) citizen science project. In this first step, we focus on exploring the factors that affect the performance of such neural networks, such as the amount of training data, number and nature of layers and the hyperparameters. We begin with a simple experiment in which we only differentiate between two extreme morphologies, using compact and multiple component extended sources. We found that a three convolutional layer architecture yielded very good results, achieving a classification accuracy of 97.4% on a test data set. The same architecture was then tested on a four-class problem where we let the network classify sources into compact and three classes of extended sources, achieving a test achieving a test accuracy of 93.5%. The best-performing convolutional neural network setup has been verified against RGZ Data Release 1 where a final test accuracy of 94.8% was obtained, using both original and augmented images. The use of sigma clipping does not offer a significant benefit overall, except in cases with a small number of training images., Comment: 17 pages, 11 figures. Accepted for publication in MNRAS

Hybrid morphology radio sources are a rare type of radio galaxy that display different Fanaroff-Riley classes on opposite sides of their nuclei. To enhance the statistical analysis of hybrid morphology radio sources, we embarked on a large-scale search of these sources within the international citizen science project, Radio Galaxy Zoo (RGZ). Here, we present 25 new candidate hybrid morphology radio galaxies. Our selected candidates are moderate power radio galaxies (L_median = 4.7x10^{24} W/(Hz sr) at redshifts 0.141Mpc) radio galaxies, one resides at a centre of a galaxy cluster, and one is hosted by a rare green bean galaxy. Although the origin of the hybrid morphology radio galaxies is still unclear, this type of radio source starts depicting itself as a rather diverse class. We discuss hybrid radio morphology formation in terms of the radio source environment (nurture) and intrinsically occurring phenomena (nature; activity cessation and amplification), showing that these peculiar radio galaxies can be formed by both mechanisms. While high angular resolution follow-up observations are still necessary to confirm our candidates, we demonstrate the efficacy of the Radio Galaxy Zoo in the pre-selection of these sources from all-sky radio surveys, and report the reliability of citizen scientists in identifying and classifying complex radio sources., Comment: 16 pages, published by Astronomical Journal

We have acquired radio continuum data between 70\,MHz and 48\,GHz for a sample of 19 southern starburst galaxies at moderate redshifts ($0.067 < z < 0.227$) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework we find the radio continuum is rarely characterised well by a single power law, instead often exhibiting low frequency turnovers below 500\,MHz, steepening at mid-to-high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 to 500\,MHz the radio-continuum at low frequency ($\nu < 200$\,MHz) could be overestimated by upwards of a factor of twelve if a simple power law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be $\alpha=-1.06$, which is steeper then the canonical value of $-0.8$ for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution.

Cross-identifying complex radio sources with optical or infra red (IR) counterparts in surveys such as the Australia Telescope Large Area Survey (ATLAS) has traditionally been performed manually. However, with new surveys from the Australian Square Kilometre Array Pathfinder (ASKAP) detecting many tens of million of radio sources such an approach is no longer feasible. This paper presents new software (LRPY - Likelihood Ratio in PYthon) to automate the process of cross-identifying radio sources with catalogues at other wavelengths. LRPY implements the Likelihood Ratio (LR) technique with a modification to account for two galaxies contributing to a sole measured radio component. We demonstrate LRPY by applying it to ATLAS DR3 and a {\it Spitzer}-based multi-wavelength fusion catalogue, identifying 3,848 matched sources via our LR-based selection criteria. A subset of 1987 sources have flux density values for all IRAC bands which allow us to use criteria to distinguish between active galactic nuclei (AGN) and star-forming galaxies (SFG). We find that 936 radio sources ($\approx47\,\%$) meet both of the Lacy and Stern AGN selection criteria. Of the matched sources, 295 have spectroscopic redshifts and we examine the radio to IR flux ratio vs redshift, proposing an AGN selection criterion below the Elvis radio-loud (RL) AGN limit for this dataset. Taking the union of all three AGN selection criteria we identify 956 as AGN ($\approx48\,\%$). From this dataset, we find a decreasing fraction of AGN with lower radio flux densities consistent with other results in the literature., Comment: 32 pages, 14 Figures, 9 Tables, accepted for publication in MNRAS

In the next few years the classification of radio sources observed by the large surveys will be a challenging problem, and spectral index is a powerful tool for addressing it. Here we present an algorithm to estimate the spectral index of sources from multiwavelength radio images. We have applied our algorithm to SCORPIO (Umana et al. 2015), a Galactic Plane survey centred around 2.1 GHz carried out with ATCA, and found we can measure reliable spectral indices only for sources stronger than 40 times the rms noise. Above a threshold of 1 mJy, the source density in SCORPIO is 20 percent greater than in a typical extra-galactic field, like ATLAS (Norris et al. 2006), because of the presence of Galactic sources. Among this excess population, 16 sources per square degree have a spectral index of about zero, suggesting optically thin thermal emission such as Hii regions and planetary nebulae, while 12 per square degree present a rising spectrum, suggesting optically thick thermal emission such as stars and UCHii regions., Comment: 12 pages, 11 figures, accepted by MNRAS

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume., Comment: Proceedings of Science, "MeerKAT Science: On the Pathway to the SKA", Stellenbosch, 25-27 May 2016

We study the faint radio population using wide-field very long baseline interferometry (VLBI) observations of 2865 known radio sources in the Cosmic Evolution Survey (COSMOS) field. The main objective of the project was to determine where active galactic nuclei (AGN) are present. The combination of number of sources, sensitivity, angular resolution and area covered by this project are unprecedented. We have detected 468 radio sources, expected to be AGNs, with the Very Long Baseline Array (VLBA) at 1.4 GHz. This is, to date, the largest sample assembled of VLBI detected sources in the sub-mJy regime. The input sample was taken from previous observations with the Very Large Array (VLA). We present the catalogue with additional multiwavelength information. We find a detection fraction of 20%, considering only those sources from the input catalogue which were in principle detectable with the VLBA (2361). As a function of redshift, we see no evolution of the detection fraction over the redshift range 0.51.5 we find that spiral galaxies become the most prevalent (48%). We demonstrate that wide-field VLBI observations, together with new calibration methods such as multi-source self-calibration and mosaicing, result in information which is difficult or impossible to obtain otherwise., Comment: 19 pages, 15 figures, 5 tables. Accepted for publication in A&A