Abstract We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes “Murriyang” radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with ∼ 3 years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release.

Abstract The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to compare and evaluate different SFs. Hydra, which currently includes the SFs Aegean, Caesar, ProFound, PyBDSF, and Selavy, provides for the addition of new SFs through containerisation and configuration files. The SF input RMS noise and island parameters are optimised to a 90% “percentage real detections” threshold (calculated from the difference between detections in the real and inverted images), to enable comparison between SFs. Hydra provides completeness and reliability diagnostics through observed-deep (D) and generated-shallow (S) images, as well as other statistics. In addition, it has a visual inspection tool for comparing residual images through various selection filters, such as S/N bins in completeness or reliability. The tool allows the user to easily compare and evaluate different SFs in order to choose their desired SF, or a combination thereof. This paper is part one of a two part series. In this paper we introduce the Hydra software suite and validate its D/S metrics using simulated data. The companion paper demonstrates the utility of Hydra by comparing the performance of SFs using both simulated and real images.

Abstract Active galactic nuclei (AGN) have been observed as far as redshift $z \sim 7$ . They are crucial in investigating the early Universe as well as the growth of supermassive black holes at their centres. Radio-loud AGN with their jets seen at a small viewing angle are called blazars and show relativistic boosting of their emission. Thus, their apparently brighter jets are easier to detect in the high-redshift Universe. DES J014132.4–542749.9 is a radio-luminous but X-ray weak blazar candidate at $z = 5$ . We conducted high-resolution radio interferometric observations of this source with the Australian Long Baseline Array at $1.7$ and $8.5$ GHz. A single, compact radio-emitting feature was detected at both frequencies with a flat radio spectrum. We derived the milliarcsecond-level accurate position of the object. The frequency dependence of its brightness temperature is similar to that of blazar sources observed at lower redshifts. Based on our observations, we can confirm its blazar nature. We compared its radio properties with those of two other similarly X-ray-weak and radio-bright AGN, and found that they show very different relativistic boosting characteristics.

Abstract Pulsar wind nebulae (PWN) are fascinating systems and archetypal sources for high-energy astrophysics in general. Due to their vicinity, brightness, to the fact that they shine at multi-wavelengths, and especially to their long-living emission at gamma rays, modelling their properties is particularly important for the correct interpretation of the visible Galaxy. A complication in this respect is the variety of properties and morphologies they show at different ages. Here, we discuss the differences among the evolutionary phases of PWN, how they have been modeled in the past and what progresses have been recently made. We approach the discussion from a phenomenological, theoretical (especially numerical) and observational point of view, with particular attention to the most recent results and open questions about the physics of such intriguing sources.

Abstract We present a systematic search for radio counterparts of novae using the Australian Square Kilometer Array Pathfinder (ASKAP). Our search used the Rapid ASKAP Continuum Survey, which covered the entire sky south of declination $+41^{\circ}$ ( $\sim$ $34000$ square degrees) at a central frequency of 887.5 MHz, the Variables and Slow Transients Pilot Survey, which covered $\sim$ $5000$ square degrees per epoch (887.5 MHz), and other ASKAP pilot surveys, which covered $\sim$ 200–2000 square degrees with 2–12 h integration times. We crossmatched radio sources found in these surveys over a two–year period, from 2019 April to 2021 August, with 440 previously identified optical novae, and found radio counterparts for four novae: V5668 Sgr, V1369 Cen, YZ Ret, and RR Tel. Follow-up observations with the Australian Telescope Compact Array confirm the ejecta thinning across all observed bands with spectral analysis indicative of synchrotron emission in V1369 Cen and YZ Ret. Our light-curve fit with the Hubble Flow model yields a value of $1.65\pm 0.17 \times 10^{-4} \rm \:M_\odot$ for the mass ejected in V1369 Cen. We also derive a peak surface brightness temperature of $250\pm80$ K for YZ Ret. Using Hubble Flow model simulated radio lightcurves for novae, we demonstrate that with a 5 $\sigma$ sensitivity limit of 1.5 mJy in 15-min survey observations, we can detect radio emission up to a distance of 4 kpc if ejecta mass is in the range $10^{-3}\rm \:M_\odot$ , and upto 1 kpc if ejecta mass is in the range $10^{-5}$ – $10^{-3}\rm \:M_\odot$ . Our study highlights ASKAP’s ability to contribute to future radio observations for novae within a distance of 1 kpc hosted on white dwarfs with masses $0.4$ – $1.25\:\rm M_\odot$ , and within a distance of 4 kpc hosted on white dwarfs with masses $0.4$ – $1.0\:\rm M_\odot$ .

Abstract Next-generation astronomical surveys naturally pose challenges for human-centred visualisation and analysis workflows that currently rely on the use of standard desktop display environments. While a significant fraction of the data preparation and analysis will be taken care of by automated pipelines, crucial steps of knowledge discovery can still only be achieved through various level of human interpretation. As the number of sources in a survey grows, there is need to both modify and simplify repetitive visualisation processes that need to be completed for each source. As tasks such as per-source quality control, candidate rejection, and morphological classification all share a single instruction, multiple data (SIMD) work pattern, they are amenable to a parallel solution. Selecting extragalactic neutral hydrogen (Hi) surveys as a representative example, we use system performance benchmarking and the visual data and reasoning methodology from the field of information visualisation to evaluate a bespoke comparative visualisation environment: the encube visual analytics framework deployed on the 83 Megapixel Swinburne Discovery Wall. Through benchmarking using spectral cube data from existing Hi surveys, we are able to perform interactive comparative visualisation via texture-based volume rendering of 180 three-dimensional (3D) data cubes at a time. The time to load a configuration of spectral cubes scale linearly with the number of voxels, with independent samples of 180 cubes (8.4 Gigavoxels or 34 Gigabytes) each loading in under 5 min. We show that parallel comparative inspection is a productive and time-saving technique which can reduce the time taken to complete SIMD-style visual tasks currently performed at the desktop by at least two orders of magnitude, potentially rendering some labour-intensive desktop-based workflows obsolete.

Abstract In order to gain more information on the 236 M dwarfs identified in the First Byurakan Survey (FBS) low-resolution (lr) spectroscopic database, Gaia EDR3 high-accuracy astrometric and photometric data and Transiting Exoplanet Survey Satellite (TESS) data are used to characterise these M dwarfs and their possible multiplicity. Among the sample of 236 relatively bright $(7.3 < K_S < 14.4)$ M dwarfs, 176 are new discoveries. The Gaia EDR3 G broadband magnitudes are in the range $11.3 < G < 17.1$ . New distance information based on the EDR3 parallaxes are used to estimate the G-band absolute magnitudes. Nine FBS M dwarfs out of 176 newly discovered lie within 25 pc of the Sun. The FBS 0909-082 is the most distant $(r=780$ pc) M dwarf of the analysed sample, with a G-band absolute magnitude $M(G) = 9.18$ , $M = 0.59$ M $_{\odot}$ , $L = 0.13597$ L $_{\odot}$ , and $T_{eff}$ = 3844 K; it can be classified as M1 - M2 subtype dwarf. The nearest is FBS 0250+167, a M7 subtype dwarf located at 3.83 pc from the Sun with a very high proper motion (5.13 arcsec yr $^{-1}$ ). The TESS estimated masses lie in the range 0.095 ( $\pm$ 0.02) M $_{\odot}\leq$ $M\leq$ 0.7 ( $\pm$ 0.1) M $_{\odot}$ and $T_{eff}$ in the range 4000 K < $T_{eff}$ < 2790 K. We analyse colour-colour and colour-absolute magnitude diagram (CaMD) diagrams for the M dwarfs. Results suggest that 27 FBS M dwarfs are double or multiple systems. The observed spectral energy distribution (SED) for some of the M dwarfs can be used to classify potential infrared excess. Using TESS light curves, flares are detected for some FBS M dwarfs. Finally, for early and late sub-classes of the M dwarfs, the detection range for survey is estimated for the first time.

Abstract Multi-messenger observations of the transient sky to detect cosmic explosions and counterparts of gravitational wave mergers critically rely on orbiting wide-FoV telescopes to cover the wide range of wavelengths where atmospheric absorption and emission limit the use of ground facilities. Thanks to continuing technological improvements, miniaturised space instruments operating as distributed-aperture constellations are offering new capabilities for the study of high-energy transients to complement ageing existing satellites. In this paper we characterise the performance of the upcoming joint SpIRIT and HERMES-TP/SP constellation for the localisation of high-energy transients through triangulation of signal arrival times. SpIRIT is an Australian technology and science demonstrator satellite designed to operate in a low-Earth Sun-synchronous Polar orbit that will augment the science operations for the equatorial HERMES-TP/SP constellation. In this work we simulate the improvement to the localisation capabilities of the HERMES-TP/SP constellation when SpIRIT is included in an orbital plane nearly perpendicular (inclination = 97.6°) to the HERMES-TP/SP orbits. For the fraction of GRBs detected by three of the HERMES satellites plus SpIRIT, we find that the combined constellation is capable of localising 60% of long GRBs to within ${\sim}30\,\textrm{deg}^{2}$ on the sky, and 60% of short GRBs within ${\sim}1850\,\textrm{deg}^{2}$ ( $1\sigma$ confidence regions), though it is beyond the scope of this work to characterise or rule out systematic uncertainty of the same order of magnitude. Based purely on statistical GRB localisation capabilities (i.e., excluding systematic uncertainties and sky coverage), these figures for long GRBs are comparable to those reported by the Fermi Gamma Burst Monitor instrument. These localisation statistics represents a reduction of the uncertainty for the burst localisation region for both long and short GRBs by a factor of ${\sim}5$ compared to the HERMES-TP/SP alone. Further improvements by an additional factor of 2 (or 4) can be achieved by launching an additional 4 (or 6) SpIRIT-like satellites into a Polar orbit, respectively, which would both increase the fraction of sky covered by multiple satellite elements, and also enable localisation of ${\geq} 60\%$ of long GRBs to within a radius of ${\sim}1.5^{\circ}$ (statistical uncertainty) on the sky, clearly demonstrating the value of a distributed all-sky high-energy transient monitor composed of nano-satellites.

Abstract Close galaxy flybys, interactions during which two galaxies inter-penetrate, are frequent and can significantly affect the evolution of individual galaxies. Equal-mass flybys are extremely rare and almost exclusively distant, while frequent flybys have mass ratios $q=0.1$ or lower, with a secondary galaxy penetrating deep into the primary. This can result in comparable strengths of interaction between the two classes of flybys and lead to essentially the same effects. To demonstrate this, emphasise and explore the role of the impact parameter further, we performed a series of N-body simulations of typical flybys with varying relative impact parameters $b/R_{\mathrm{vir},1}$ ranging from $0.114$ to $0.272$ of the virial radius of the primary galaxy. Two-armed spirals form during flybys, with radii of origin correlated with the impact parameter and strengths well approximated with an inverted S-curve. The impact parameter does not affect the shape of induced spirals, and the lifetimes of a distinguished spiral structure appear to be constant, $T_\mathrm{LF} \sim 2$ Gyr. Bars, with strengths anti-correlated with the impact parameter, form after the encounter is over in simulations with $b/R_{\mathrm{vir},1} \leq 0.178$ and interaction strengths $S\geq0.076$ , but they are short-lived except for the stronger interactions with $S\geq0.129$ . We showcase an occurrence of multiple structures (ring-like, double bar) that survives for an exceptionally long time in one of the simulations. Effects on the pre-existing bar instability, that develops much later, are diverse: from an acceleration of bar formation, little to no effect, to even bar suppression. There is no uniform correlation between these effects and the impact parameter, as they are secondary effects, happening later in a post-flyby stage. Classical bulges are resilient to flyby interactions, while dark matter halos can significantly spin up in the amount anti-correlated with the impact parameter. There is an offset angle between the angular momentum vector of the dark matter halo and that of a disc, and it correlates linearly with the impact parameter. Thus, flybys remain an important pathway for structural evolution within galaxies in the local Universe.

Abstract Asymmetric emission of gravitational waves during mergers of black holes (BHs) produces a recoil kick, which can set a newly formed BH on a bound orbit around the centre of its host galaxy, or even completely eject it. To study this population of recoiling BHs we extract properties of galaxies with merging BHs from Illustris TNG300 simulation and then employ both analytical and numerical techniques to model unresolved process of BH recoil. This comparative analysis between analytical and numerical models shows that, on cosmological scales, numerically modelled recoiling BHs have a higher escape probability and predict a greater number of offset active galactic nuclei (AGN). BH escaped probability $>$ 40% is expected in 25 $\%$ of merger remnants in numerical models, compared to 8 $\%$ in analytical models. At the same time, the predicted number of offset AGN at separations ${>}5$ kpc changes from 58 $\%$ for numerical models to 3 $\%$ for analytical models. Since BH ejections in major merger remnants occur in non-virialised systems, static analytical models cannot provide an accurate description. Thus we argue that numerical models should be used to estimate the expected number density of escaped BHs and offset AGN.