We present the first observational study of pulsars performed with the second-generation precursor stations to the low-frequency component of the Square Kilometre Array (SKA-Low): the Aperture Array Verification System 2 (AAVS2) and the Engineering Development Array 2 (EDA2). Using the SKA-Low stations, we have observed 100 southern-sky pulsars between 70–350 MHz, including follow-up observations at multiple frequencies for a selected sample of bright pulsars. These observations have yielded detections of 22 pulsars, including the lowest-frequency detections ever published for 6 pulsars, despite the modest sensitivity of initial system where the recording bandwidth is limited to ${\sim}{1}\,\mathrm{MHz}$ . By comparing simultaneous flux density measurements obtained with the SKA-Low stations and performing rigorous electromagnetic simulations, we verify the accuracy of the SKA-Low sensitivity simulation code presented in Sokolowski (2022, PASA, 39, e015). Furthermore, we perform model fits to the radio spectra of the detected pulsars using the method developed by Jankowski (2018, MNRAS, 473, 4436), including nine pulsars which were not fitted in the original work. We robustly classify the spectra into five morphological classes and find that all but one pulsar exhibit deviations from simple power-law behaviour. These findings suggest that pulsars with well-determined spectra are more likely to show spectral flattening or turnover than average. Our work demonstrates how SKA-Low stations can be meaningfully used for scientifically useful measurements and analysis of pulsar radio spectra, which are important inputs for informing pulsar surveys and science planned with the SKA-Low.

We present a study of scintillation induced by the mid-latitude ionosphere. By implementing methods currently used in Interplanetary Scintillation studies to measure amplitude scintillation at low frequencies, we have proven it is possible to use the Murchison Widefield Array to study ionospheric scintillation in the weak regime, which is sensitive to structures on scales ${\sim}300$ m at our observing frequency of 154 MHz, where the phase variance on this scale was $0.06\, \textrm{rad}^{2}$ in the most extreme case observed. Analysing over 1000 individual 2-min observations, we compared the ionospheric phase variance with that inferred with previous measurements of refractive shifts, which are most sensitive to scales almost an order of magnitude larger. The two measurements were found to be highly correlated (Pearson correlation coefficient 0.71). We observed that for an active ionosphere, the relationship between these two metrics is in line with what would be expected if the ionosphere’s structure is described by Kolmogorov turbulence between the relevant scales of 300 and 2000 m. In the most extreme ionospheric conditions, the refractive shifts were sometimes found to underestimate the small-scale variance by a factor of four or more, and it is these ionospheric conditions that could have significant effects on radio astronomy observations.

We describe a new low-frequency wideband radio survey of the southern sky. Observations covering 72–231 MHz and Declinations south of $+30^\circ$ have been performed with the Murchison Widefield Array “extended” Phase II configuration over 2018–2020 and will be processed to form data products including continuum and polarisation images and mosaics, multi-frequency catalogues, transient search data, and ionospheric measurements. From a pilot field described in this work, we publish an initial data release covering 1,447 $\mathrm{deg}^2$ over $4\,\mathrm{h}\leq \mathrm{RA}\leq 13\,\mathrm{h}$ , $-32.7^\circ \leq \mathrm{Dec} \leq -20.7^\circ$ . We process twenty frequency bands sampling 72–231 MHz, with a resolution of 2′–45′′, and produce a wideband source-finding image across 170–231 MHz with a root mean square noise of $1.27\pm0.15\,\mathrm{mJy\,beam}^{-1}$ . Source-finding yields 78,967 components, of which 71,320 are fitted spectrally. The catalogue has a completeness of 98% at ${{\sim}}50\,\mathrm{mJy}$ , and a reliability of 98.2% at $5\sigma$ rising to 99.7% at $7\sigma$ . A catalogue is available from Vizier; images are made available via the PASA datastore, AAO Data Central, and SkyView. This is the first in a series of data releases from the GLEAM-X survey.

We present a catalogue of isolated field elliptical (IfE) galaxies drawn from the W1 field of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). 228 IfEs were identified from a flux-limited $(r<21.8)$ galaxy catalogue which corresponds to a density of 3 IfE/sq.deg. For comparison we consider a sample of elliptical galaxies living in dense environments, based on identification of the brightest cluster galaxies (BGCs) in the same survey. Using the same dataset for the comparison sample ensures a uniform selection, including in the redshift range as IfEs (i.e. $0.1<z<0.9$ ). A comparison of elliptical galaxies in different environments reveals that IfEs and BCGs have similar behaviours in their colours, star formation activities, and scaling relations of mass–size and size–luminosity. IfEs and BCGs have similar slopes in the scaling relations with respect to cluster ellipticals within the $-24 \leq M_{r} \leq -22$ magnitude and $10.2< \textrm{log}(M_{*}/ \textrm M_\odot)\leq12.0$ mass ranges. Three IfEs identified in this study can be associated with fossil groups found in the same survey area which gives clues for future studies.

Long-duration gamma-ray burst (GRB) afterglow observations offer cutting-edge opportunities to characterise the star formation history of the Universe back to the epoch of reionisation, and to measure the chemical composition of interstellar and intergalactic gas through absorption spectroscopy. The main barrier to progress is the low efficiency in rapidly and confidently identifying which bursts are high redshift ( $z > 5$ ) candidates before they fade, as this requires low-latency follow-up observations at near-infrared wavelengths (or longer) to determine a reliable photometric redshift estimate. Since no current or planned gamma-ray observatories carry near-infrared telescopes on-board, complementary facilities are needed. So far this task has been performed by instruments on the ground, but sky visibility and weather constraints limit the number of GRB targets that can be observed and the speed at which follow-up is possible. In this work we develop a Monte Carlo simulation framework to investigate an alternative approach based on the use of a rapid-response near-infrared nano-satellite, capable of simultaneous imaging in four bands from $0.8$ to $1.7\,\unicode{x03BC}$ m (a mission concept called SkyHopper). Using as reference a sample of 88 afterglows observed with the GROND instrument on the MPG/ESO telescope, we find that such a nano-satellite is capable of detecting in the H-band (1.6 $\unicode{x03BC}$ m) $72.5\% \pm 3.1\%$ of GRBs concurrently observable with the Swift satellite via its UVOT instrument (and $44.1\% \pm 12.3\%$ of high redshift ( $z>5$ ) GRBs) within 60 min of the GRB prompt emission. This corresponds to detecting ${\sim}55$ GRB afterglows per year, of which 1–3 have $z > 5$ . These rates represent a substantial contribution to the field of high-z GRB science, as only 23 $z > 5$ GRBs have been collectively discovered by the entire astronomical community over the last ${\sim}24$ yr. Future discoveries are critically needed to take advantage of next generation follow-up spectroscopic facilities such as 30m-class ground telescopes and the James Webb Space Telescope. Furthermore, a systematic space-based follow-up of afterglows in the near-infrared will offer new insight on the population of dusty (‘dark’) GRBs which are primarily found at cosmic noon ( $z\sim 1-3$ ). Additionally, we find that launching a mini-constellation of 3 near-infrared nano-satellites would increase the detection fraction of afterglows to ${\sim}83\%$ and substantially reduce the latency in the photometric redshift determination.

Recently Vernstrom et al. (2021, MNRAS) claimed the first definitive detection of the synchrotron cosmic web, obtained by ‘stacking’ hundreds of thousands of pairs of close-proximity clusters in low-frequency radio observations and looking for a residual excess signal spanning the intracluster bridge. A reproduction study by Hodgson et al. (2022, PASA, 39, e013), using both the original radio data as well as new observations with the Murchison Widefield Array, failed to confirm these findings. Whilst the detection remains unsure, we here turn to stacking a simulated radio sky to understand what kind of excess radio signal is predicted by our current best cosmological models of the synchrotron cosmic web. We use the FIlaments & GAlactic RadiO (FIGARO; Hodgson et al. 2021a, PASA, 38, e047) simulation, which models both the synchrotron cosmic web as well as various subtypes of active galactic nucleii and star-forming galaxies. Being a simulation, we have perfect knowledge of the location of clusters and galaxy groups which we use in our own stacking experiment. Whilst we do find an excess radio signature in our stacks that is attributable to the synchrotron cosmic web, its distribution is very different to that found by Vernstrom et al. (2021, MNRAS). Instead, we observe the appearance of excess emission on the immediate interiors of cluster pairs as a result of asymmetric, ‘radio relic’-like shocks surrounding cluster cores, whilst the excess emission spanning the intracluster region—attributable to filaments proper—is two orders of magnitude lower and undetectable in our experiment even under ideal conditions.

We present the first unbiased survey of neutral hydrogen absorption in the Small Magellanic Cloud. The survey utilises pilot neutral hydrogen observations with the Australian Square Kilometre Array Pathfinder telescope as part of the Galactic Australian Square Kilometre Array Pathfinder neutral hydrogen project whose dataset has been processed with the Galactic Australian Square Kilometre Array Pathfinder-HI absorption pipeline, also described here. This dataset provides absorption spectra towards 229 continuum sources, a 275% increase in the number of continuum sources previously published in the Small Magellanic Cloud region, as well as an improvement in the quality of absorption spectra over previous surveys of the Small Magellanic Cloud. Our unbiased view, combined with the closely matched beam size between emission and absorption, reveals a lower cold gas faction (11%) than the 2019 ATCA survey of the Small Magellanic Cloud and is more representative of the Small Magellanic Cloud as a whole. We also find that the optical depth varies greatly between the Small Magellanic Cloud’s bar and wing regions. In the bar we find that the optical depth is generally low (correction factor to the optically thin column density assumption of $\mathcal{R}_{\mathrm{HI}} \sim 1.04$ ) but increases linearly with column density. In the wing however, there is a wide scatter in optical depth despite a tighter range of column densities.

Understanding the physical and evolutionary properties of Hot Stellar Systems (HSS) is a major challenge in astronomy. We studied the dataset on 13 456 HSS of Misgeld & Hilker (2011, MNRAS, 414, 3 699) that includes 12 763 candidate globular clusters using stellar mass ( $M_s$ ), effective radius ( $R_e$ ) and mass-to-luminosity ratio ( $M_s/L_\nu$ ), and found multi-layered homogeneous grouping among these stellar systems. Our methods elicited eight homogeneous ellipsoidal groups at the finest sub-group level. Some of these groups have high overlap and were merged through a multi-phased syncytial algorithm motivated from Almodóvar-Rivera & Maitra (2020, JMLR, 21, 1). Five groups were merged in the first phase, resulting in three complex-structured groups. Our algorithm determined further complex structure and permitted another merging phase, revealing two complex-structured groups at the highest level. A nonparametric bootstrap procedure was also used to estimate the confidence of each of our group assignments. These assignments generally had high confidence in classification, indicating great degree of certainty of the HSS assignments into our complex-structured groups. The physical and kinematic properties of the two groups were assessed in terms of $M_s$ , $R_e$ , surface density and $M_s/L_\nu$ . The first group consisted of older, smaller and less bright HSS while the second group consisted of brighter and younger HSS. Our analysis provides novel insight into the physical and evolutionary properties of HSS and also helps understand physical and evolutionary properties of candidate globular clusters. Further, the candidate globular clusters (GCs) are seen to have very high chance of really being GCs rather than dwarfs or dwarf ellipticals that are also indicated to be quite distinct from each other.

MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsar (MSPs) to high precision ( ${<} 1 \unicode{x03BC} \mathrm{s}$ ) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing Array (IPTA). In order to plan for the remainder of the programme and to use the allocated time most efficiently, we have conducted an initial census with the MeerKAT ‘L-band’ receiver of 189 MSPs visible to MeerKAT and here present their dispersion measures, polarisation profiles, polarisation fractions, rotation measures, flux density measurements, spectral indices, and timing potential. As all of these observations are taken with the same instrument (which uses coherent dedispersion, interferometric polarisation calibration techniques, and a uniform flux scale), they present an excellent resource for population studies. We used wideband pulse portraits as timing standards for each MSP and demonstrated that the MeerTime Pulsar Timing Array (MPTA) can already contribute significantly to the IPTA as it currently achieves better than $1\,\unicode{x03BC}\mathrm{s}$ timing accuracy on 89 MSPs (observed with fortnightly cadence). By the conclusion of the initial five-year MeerTime programme in 2024 July, the MPTA will be extremely significant in global efforts to detect the gravitational wave background with a contribution to the detection statistic comparable to other long-standing timing programmes.

Bayesian inference is a powerful tool in gravitational-wave astronomy. It enables us to deduce the properties of merging compact-object binaries and to determine how these mergers are distributed as a population according to mass, spin, and redshift. As key results are increasingly derived using Bayesian inference, there is increasing scrutiny on Bayesian methods. In this review, we discuss the phenomenon of model misspecification, in which results obtained with Bayesian inference are misleading because of deficiencies in the assumed model(s). Such deficiencies can impede our inferences of the true parameters describing physical systems. They can also reduce our ability to distinguish the ‘best fitting’ model: it can be misleading to say that Model A is preferred over Model B if both models are manifestly poor descriptions of reality. Broadly speaking, there are two ways in which models fail. Firstly, models that fail to adequately describe the data (either the signal or the noise) have misspecified likelihoods. Secondly, population models—designed, for example, to describe the distribution of black hole masses—may fail to adequately describe the true population due to a misspecified prior. We recommend tests and checks that are useful for spotting misspecified models using examples inspired by gravitational-wave astronomy. We include companion python notebooks to illustrate essential concepts.

A plausible formation scenario for the Galactic globular clusters 47 Tucanae (47 Tuc) and Omega Centauri $(\omega$ Cen) is that they are tidally stripped remnants of dwarf galaxies, in which case they are likely to have retained a fraction of their dark matter cores. In this study, we have used the ultra-wide band receiver on the Parkes telescope (Murriyang) to place upper limits on the annihilation rate of exotic Light Dark Matter particles $(\chi)$ via the $\chi\chi\rightarrow e^+e^-$ channel using measurements of the recombination rate of positronium (Ps). This is an extension of a technique previously used to search for Ps in the Galactic Centre. However, by stacking of spectral data at multiple line frequencies, we have been able to improve sensitivity. Our measurements have resulted in $3-\sigma$ flux density (recombination rate) upper limits of 1.7 mJy $\left(1.4\times 10^{43}\, \mathrm{s}^{-1}\right)$ and 0.8 mJy $\left(1.1 \times 10^{43} \mathrm{s}^{-1}\right)$ for 47 Tuc and $\omega$ Cen, respectively. Within the Parkes beam at the cluster distances, which varies from 10–23 pc depending on the frequency of the recombination line, and for an assumed annihilation cross-section $\langle\sigma v\rangle = 3\times 10^{-29} \mathrm{cm}^3\, \mathrm{s}^{-1}$ , we calculate upper limits to the dark matter mass and rms dark matter density of ${\lesssim} 1.2-1.3\times 10^5 f_n^{-0.5}$ $\left(m_\chi/\mathrm{MeV\, c}^{-2}\right)$ $\mathrm{M}_{\odot}$ and ${\lesssim} 48-54 f_n^{-0.5}$ $\left(m_\chi/\mathrm{MeV\, c}^{-2}\right)$ $\mathrm{M}_{\odot} \mathrm{pc}^{-3}$ for the clusters, where $f_n=R_n/R_p$ is the ratio of Ps recombination transitions to annihilations, estimated to be ${\sim}0.01$ . The radio limits for $\omega$ Cen suggest that, for a fiducial dark/luminous mass ratio of ${\sim}0.05$ , any contribution from Light Dark Matter is small unless $\langle\sigma v\rangle < 7.9\times 10^{-28}\ \left(m_\chi/\mathrm{MeV\, c}^{-2}\right)^2 \mathrm{cm}^3 \mathrm{s}^{-1}$ . Owing to the compactness and proximity of the clusters, archival 511-keV measurements suggest even tighter limits than permitted by CMB anisotropies, $\langle\sigma v\rangle < 8.6\times 10^{-31}\ (m_\chi/\mathrm{MeV\, c}^{-2})^2 \mathrm{cm}^3 \mathrm{s}^{-1}$ . Due to the very low synchrotron radiation background, our recombination rate limits substantially improve on previous radio limits for the Milky Way.

We report the results of a sensitive search for water maser emission in the Local Group Galaxy NGC 6822 with the Karl G. Jansky Very Large Array. The observations provide tentative single-epoch detections of four candidates, associated with two infrared-bright star formation regions (Hubble I/III and Hubble IV). The candidate maser detections are all offset from the velocity range where strong emission from Hi neutral gas is observed towards NGC 6822, with the closest offset by $\sim\!40\, \mathrm{kms}^{-1}$ . Our observations include the location of NL1K, a previous tentative water maser detection in NGC 6822. We do not detect any emission from this location with a sensitivity limit approximately a factor of 5 better than the original Sardina Radio Telescope observations.

The Murchison Widefield Array (MWA) is a low-frequency aperture array capable of high-time and frequency resolution astronomy applications such as pulsar studies. The large field-of-view of the MWA (hundreds of square degrees) can also be exploited to attain fast survey speeds for all-sky pulsar search applications, but to maximise sensitivity requires forming thousands of tied-array beams from each voltage-capture observation. The necessity of using calibration solutions that are separated from the target observation both temporally and spatially makes pulsar observations vulnerable to uncorrected, frequency-dependent positional offsets due to the ionosphere. These offsets may be large enough to move the source away from the centre of the tied-array beam, incurring sensitivity drops of ${\sim}30{-}50\%$ in Phase II extended array configuration. We analyse these offsets in pulsar observations and develop a method for mitigating them, improving both the source position accuracy and the sensitivity. This analysis prompted the development of a multi-pixel beamforming functionality that can generate dozens of tied-array beams simultaneously, which runs a factor of ten times faster compared to the original single-pixel version. This enhancement makes it feasible to observe multiple pulsars within the vast field of view of the MWA and supports the ongoing large-scale pulsar survey efforts with the MWA. We explore the extent to which ionospheric offset correction will be necessary for the MWA Phase III and the low-frequency square kilometre array (SKA-low).

One of the major priorities of international radio astronomy is to study the early universe through the detection of the 21 cm HI line from the epoch of reionisation (EoR). Due to the weak nature of the 21 cm signal, an important part in the detection of the EoR is removing contaminating foregrounds from our observations as they are multiple orders of magnitude brighter. In order to achieve this, sky maps spanning a wide range of frequencies and angular scales are required for calibration and foreground subtraction. Complementing the existing low-frequency sky maps, we have constructed a Southern Sky map through spherical harmonic transit interferometry utilising the Engineering Development Array 2 (EDA2), a Square Kilometre Array (SKA) low-frequency array prototype system. We use the m-mode formalism to create an all-sky map at 159 MHz with an angular resolution of 3 degrees, with data from the EDA2 providing information over +60 degrees to –90 degrees in declination. We also introduce a new method for visualising and quantifying how the baseline distribution of an interferometer maps to the spherical harmonics and discuss how prior information can be used to constrain spherical harmonic components that the interferometer is not sensitive to.

The Square Kilometre Array (SKA) will be the largest radio astronomy observatory ever built, providing unprecedented sensitivity over a very broad frequency band from 50 MHz to 15.3 GHz. The SKA’s low frequency component (SKA-Low), which will observe in the 50–350 MHz band, will be built at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. It will consist of 512 stations each composed of 256 dual-polarised antennas, and the sensitivity of an individual station is pivotal to the performance of the entire SKA-Low telescope. The answer to the question in the title is, it depends. The sensitivity of a low frequency array, such as an SKA-Low station, depends strongly on the pointing direction of the digitally formed station beam and the local sidereal time (LST), and is different for the two orthogonal polarisations of the antennas. The accurate prediction of the SKA-Low sensitivity in an arbitrary direction in the sky is crucial for future observation planning. Here, we present a sensitivity calculator for the SKA-Low radio telescope, using a database of pre-computed sensitivity values for two realisations of an SKA-Low station architecture. One realisation uses the log-periodic antennas selected for SKA-Low. The second uses a known benchmark, in the form of the bowtie dipoles of the Murchison Widefield Array. Prototype stations of both types were deployed at the MRO in 2019, and since then have been collecting commissioning and verification data. These data were used to measure the sensitivity of the stations at several frequencies and over at least 24 h intervals, and were compared to the predictions described in this paper. The sensitivity values stored in the SQLite database were pre-computed for the X, Y, and Stokes I polarisations in 10 MHz frequency steps, $\scriptsize{1/2}$ hour LST intervals, and $5^\circ$ resolution in pointing directions. The database allows users to quickly and easily estimate the sensitivity of SKA-Low for arbitrary observing parameters (your favourite object) using interactive web-based or command line interfaces. The sensitivity can be calculated using publicly available web interface (http://sensitivity.skalow.link) or a command line python package (https://github.com/marcinsokolowski/station_beam), which can also be used to calculate the sensitivity for arbitrary pointing directions, frequencies, and times without interpolations.

Closure phase is the phase of a closed-loop product of spatial coherences formed by a ${\ge}3$ -element interferometer array. Its invariance to phase corruption attributable to individual array elements acquired during the propagation and the measurement processes, subsequent calibration, and errors therein, makes it a valuable tool in interferometry applications that otherwise require high-accuracy phase calibration. However, its understanding has remained mainly mathematical and limited to the aperture plane (Fourier dual of the image plane). Here, we present a geometrical, image domain view of closure phase, which until now has been lacking. Using the principal triangle in a 3-element interference image formed by a triad of interferometer elements, we show that the properties of closure phase, particularly its invariance to multiplicative element-based corruption factors (even of a large magnitude) and to translation, are intricately related to the conserved properties of the triangle, namely, its shape, orientation, and size, which is referred herein as the ‘shape-orientation-size (SOS) conservation principle’. In the absence of a need for element-based amplitude calibration of the interferometer array (as is typical in optical interferometry), the principal triangle in any 3-element interference image formed from phase-uncalibrated spatial coherences is still a true and uncorrupted representation of the source object’s morphology, except for a possible shift. Based on this knowledge of the triangle SOS conservation principle, we present two geometric methods to measure the closure phase directly from a simple 3-element interference image (without requiring an aperture-plane view): (i) the closure phase is directly measurable from any one of the triangle’s heights, and (ii) the squared closure phase is proportional to the product of the areas enclosed by the triad of array elements and the principal triangle in the aperture and image planes, respectively. We validate the geometric understanding of closure phase in the image plane using observations with the Karl G. Jansky Very Large Array, and the Event Horizon Telescope. These results verify the SOS conservation principle across a wide range of radio interferometric conditions. This geometric insight can be potentially valuable to other interferometric applications, such as optical interferometry. We also generalise these geometric relationships to an N-element interferometer.

We follow up on the surprising recent announcement by Vernstrom et al. (2021, MNRAS) of the detection of the synchrotron cosmic web. We attempt to reproduce their detection with new observations with the Phase II, extended configuration of the Murchison Widefield Array at 118.5 MHz. We reproduce their detection methodology by stacking pairs of nearby luminous red galaxies (LRGs)—used as tracers for clusters and galaxy groups—contained in our low-frequency radio observations. We show that our observations are significantly more sensitive than those used in Vernstrom et al. and that our angular sensitivity is sufficient. And yet, we make no statistically significant detection of excess radio emission along the bridge spanning the LRG pairs. This non-detection is true both for the original LRG pair catalogue as used in Vernstrom et al., as well as for other larger catalogues with modified selection criteria. Finally, we return to the original data sets used in Vernstrom et al., and find that whilst we clearly reproduce the excess X-ray emission from ROSAT, we are not able to reproduce any kind of broad and extended excess intercluster filamentary emission using the original 118.5 MHz MWA survey data. In the interests of understanding this result, as part of this paper we release images of the 14 fields used in this study, the final stacked images, as well as key components of our stacking and modelling code.

Rest-frame mid- to far-infrared (IR) spectroscopy is a powerful tool to study how galaxies formed and evolved, because a major part of their evolution occurs in heavily dust enshrouded environments, especially at the so-called Cosmic Noon ( $1< z < 3$ ). Using the calibrations of IR lines and features, recently updated with Herschel and Spitzer spectroscopy, we predict their expected fluxes with the aim to measure the Star Formation (SF) and the Black Hole Accretion (BHA) rates in intermediate to high redshift galaxies. On the one hand, the recent launch of the James Webb Space Telescope (JWST) offers new mid-IR spectroscopic capabilities that will enable for the first time a detailed investigation of both the SF and the BHA obscured processes as a function of cosmic time. We make an assessment of the spectral lines and features that can be detected by JWST-MIRI in galaxies and active galactic nuclei up to redshift $z \sim 3$ . The fine structure lines of [MgIV]4.49 $\unicode{x03BC}\textrm{m}$ and [ArVI]4.53 $\unicode{x03BC}\textrm{m}$ can be used as BHA rate tracers for the $1 \lesssim z \lesssim 3$ range, and we propose the [NeVI]7.65 $\unicode{x03BC}\textrm{m}$ line as the best tracer for $z \lesssim 1.5$ . The [ArII]6.98 $\unicode{x03BC}\textrm{m}$ and [ArIII]8.99 $\unicode{x03BC}\textrm{m}$ lines can be used to measure the SF rate at $z \lesssim 3$ and $z \lesssim 2$ , respectively, while the stronger [NeII]12.8 $\unicode{x03BC}\textrm{m}$ line exits the JWST-MIRI spectral range above $z \gtrsim 1.2$ . At higher redshifts, the PAH features at 6.2 and 7.7 $\unicode{x03BC}\textrm{m}$ can be observed at $z \lesssim 3$ and $z \lesssim 2.7$ , respectively. On the other hand, rest-frame far-IR spectroscopic observations of high redshift galaxies ( $z \gtrsim 3$ ) have been collected with the Atacama Large Millimeter Array (ALMA) in the last few years. The observability of far-IR lines from ALMA depends on the observed frequency, due to the significant decrease of the atmospheric transmission at the highest frequencies ( $\gtrsim420\,\rm{Hz}$ ). The [CII]158 $\unicode{x03BC}\textrm{m}$ line is a reliable tracer of the SF rate and can in most cases ( $0.9 \lesssim z \lesssim 2$ and $2 \lesssim z \lesssim 9$ ) be observed. Additionally, we propose the use of the combination of [OIII]88 $\,\unicode{x03BC}$ m and [OI]145 $\,\unicode{x03BC}$ m lines as an alternative SF rate tracer, that can be detected above $z \gtrsim 3$ . Overall, we emphasize the importance of using multi-feature analysis to measure both BHA and SFR, since individual tracers can be strongly dependent on the local ISM conditions and vary from source to source. However, we conclude that the peak of the obscured SF and BHA activities at Cosmic Noon falls outside the wavelength coverage of facilities currently operating or under development. A new IR space telescope covering the full IR spectral range from about 10 to $300\,\unicode{x03BC}\textrm{m}$ and actively cooled to achieve high sensitivity, will be needed.

This paper is the fourth in a series of low-frequency searches for technosignatures. Using the Murchison Widefield Array over two nights, we integrate 7 h of data toward the Galactic Centre (centred on the position of Sagittarius $\mathrm{A}^{*}$ ) with a total field-of-view of $200\,\mathrm{deg}^{2}$ . We present a targeted search toward 144 exoplanetary systems, at our best yet angular resolution (75 arcsec). This is the first technosignature search at a central frequency of 155 MHz toward the Galactic Centre (our previous central frequencies have been lower). A blind search toward in excess of 3 million stars toward the Galactic Centre and Galactic bulge is also completed, placing an equivalent isotropic power limit $<\!1.1\times10^{19}\,\mathrm{W}$ at the distance to the Galactic Centre. No plausible technosignatures are detected.

Proposed next-generation networks of gravitational-wave observatories include dedicated kilohertz instruments that target neutron star science, such as the proposed Neutron Star Extreme Matter Observatory, NEMO. The original proposal for NEMO highlighted the need for it to exist in a network of gravitational-wave observatories to ensure detection confidence and sky localisation of sources. We show that NEMO-like observatories have significant utility on their own as coincident electromagnetic observations can provide the detection significance and sky localisation. We show that, with a single NEMO-like detector and expected electromagnetic observatories in the late 2020 s and early 2030 s such as the Vera C. Rubin observatory and SVOM, approximately 40% of all binary neutron star mergers detected with gravitational waves could be confidently identified as coincident multimessenger detections. We show that we expect $2^{+10}_{-1}{yr^{-1}}{}$ coincident observations of gravitational-wave mergers with gamma-ray burst prompt emission, $13^{+23}_{-10}{yr^{-1}}{}$ detections with kilonova observations, and $4^{+18}_{-3}{yr^{-1}}{}$ with broadband afterglows and kilonovae, where the uncertainties are 90% confidence intervals arising from uncertainty in current merger-rate estimates. Combined, this implies a coincident detection rate of $14^{+25}_{-11}{yr^{-1}}{}$ out to $300\,\mathrm{Mpc}$ . These numbers indicate significant science potential for a single kilohertz gravitational-wave detector operating without a global network of other gravitational-wave observatories.