The MWA long baseline Epoch of reionisation survey—I. Improved source catalogue for the EoR 0 field

One of the principal systematic constraints on the Epoch of Reionisation (EoR) experiment is the accuracy of the foreground calibration model. Recent results have shown that highly accurate models of extended foreground sources, and including models for sources in both the primary beam and its sidelobes, are necessary for reducing foreground power. To improve the accuracy of the source models for the EoR fields observed by the Murchison Widefield Array (MWA), we conducted the MWA Long Baseline Epoch of Reionisation Survey (LoBES). This survey consists of multi-frequency observations of the main MWA EoR fields and their eight neighbouring fields using the MWA Phase II extended array. We present the results of the first half of this survey centred on the MWA EoR0 observing field (centred at RA (J2000) $0^\mathrm{h}$ , Dec (J2000) $-27^{\circ}$ ). This half of the survey covers an area of 3 069 degrees $^2$ , with an average rms of 2.1 mJy beam–1. The resulting catalogue contains a total of 80 824 sources, with 16 separate spectral measurements between 100 and 230 MHz, and spectral modelling for 78 $\%$ of these sources. Over this region we estimate that the catalogue is 90 $\%$ complete at 32 mJy, and 70 $\%$ complete at 10.5 mJy. The overall normalised source counts are found to be in good agreement with previous low-frequency surveys at similar sensitivities. Testing the performance of the new source models we measure lower residual rms values for peeled sources, particularly for extended sources, in a set of MWA Phase I data. The 2-dimensional power spectrum of these data residuals also show improvement on small angular scales—consistent with the better angular resolution of the LoBES catalogue. It is clear that the LoBES sky models improve upon the current sky model used by the Australian MWA EoR group for the EoR0 field.

Updated studies on exomoons in the HD 23079 system

We re-evaluate the outer edge of orbital stability for possible exomoons orbiting the radial velocity planet discovered in the HD 23079 system. In this system, a solar-type star hosts a Jupiter-mass planet in a nearly circular orbit in the outer stellar habitable zone. The outer stability limit of exomoons is deduced using N-body and tidal migration simulations considering a large range of initial conditions, encompassing both prograde and retrograde orbits. In particular, we extend previous works by evaluating many values in the satellite mean anomaly to identify and exclude regions of quasi-stability. Future observations of this system can make use of our results through a scale factor relative to the currently measured minimum mass. Using a constant time lag tidal model (Hut 1981), we find that plausible tidal interactions within the system are insufficient to induce significant outward migration toward the theoretical stability limit. While current technologies are incapable of detecting exomoons in this system, we comment on the detectability of putative moons through Doppler monitoring within direct imaging observations in view of future research capacities.

FIGARO simulation: FIlaments & GAlactic RadiO simulation

We produce the first low to mid-frequency radio simulation that incorporates both traditional extragalactic radio sources as well as synchrotron cosmic web emission. The FIlaments & GAlactic RadiO (FIGARO) simulation includes 10 unique 4° × 4° fields, incorporating active galactic nucleii (AGNs), star-forming galaxies (SFGs), and synchrotron cosmic web emission out to a redshift of z = 0.8 and over the frequency range 100–1 400 MHz. To do this, the simulation brings together a recent 1003 Mpc3 magnetohydrodynamic simulation (Vazza et al. 2019, A&A, 627, A5), calibrated to match observed radio relic population statistics, alongside updated ‘T-RECS’ code for simulating extragalactic radio sources (Bonaldi et al. 2019, MNRAS, 482, 2). Uniquely, the AGNs and SFGs are populated and positioned in accordance with the underlying matter density of the cosmological simulation. In this way, the simulation provides an accurate understanding of the apparent morphology, angular scales, and brightness of the cosmic web as well as—crucially—the clustering properties of the cosmic web with respect to the embedded extragalactic radio population. We find that the synchrotron cosmic web does not closely trace the underlying mass distribution of the cosmic web, but is instead dominated by shocked shells of emission surrounding dark matter halos and resembles a large, undetected population of radio relics. We also show that, with accurate kernels, the cosmic web radio emission is clearly detectable by cross-correlation techniques and this signal is separable from the embedded extragalactic radio population. We offer the simulation as a public resource towards the development of techniques for detecting and measuring the synchrotron cosmic web.

The GLEAM 200-MHz local radio luminosity function for AGN and star-forming galaxies

The GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) is a radio continuum survey at 76–227 MHz of the entire southern sky (Declination $<\!{+}30^{\circ}$ ) with an angular resolution of ${\approx}2$ arcmin. In this paper, we combine GLEAM data with optical spectroscopy from the 6dF Galaxy Survey to construct a sample of 1 590 local (median $z \approx 0.064$ ) radio sources with $S_{200\,\mathrm{MHz}} > 55$ mJy across an area of ${\approx}16\,700\,\mathrm{deg}^{2}$ . From the optical spectra, we identify the dominant physical process responsible for the radio emission from each galaxy: 73% are fuelled by an active galactic nucleus (AGN) and 27% by star formation. We present the local radio luminosity function for AGN and star-forming (SF) galaxies at 200 MHz and characterise the typical radio spectra of these two populations between 76 MHz and ${\sim}1$ GHz. For the AGN, the median spectral index between 200 MHz and ${\sim}1$ GHz, $\alpha_{\mathrm{high}}$ , is $-0.600 \pm 0.010$ (where $S \propto \nu^{\alpha}$ ) and the median spectral index within the GLEAM band, $\alpha_{\mathrm{low}}$ , is $-0.704 \pm 0.011$ . For the SF galaxies, the median value of $\alpha_{\mathrm{high}}$ is $-0.650 \pm 0.010$ and the median value of $\alpha_{\mathrm{low}}$ is $-0.596 \pm 0.015$ . Among the AGN population, flat-spectrum sources are more common at lower radio luminosity, suggesting the existence of a significant population of weak radio AGN that remain core-dominated even at low frequencies. However, around 4% of local radio AGN have ultra-steep radio spectra at low frequencies ( $\alpha_{\mathrm{low}} < -1.2$ ). These ultra-steep-spectrum sources span a wide range in radio luminosity, and further work is needed to clarify their nature.

Complex optical/UV and X-ray variability in Seyfert 1 galaxy Mrk 509

We performed a detailed spectral and timing analysis of a Seyfert 1 galaxy Mrk 509 using data from the Neil Gehrels Swift observatory that spanned over $\sim$ 13 years between 2006 and 2019. To study the variability properties from the optical/UV to X-ray emission, we used a total of 275 pointed observations in this work. The average spectrum over the entire duration exhibits a strong soft X-ray excess above the power law continuum. The soft X-ray excess is well described by two thermal components with temperatures of kT $_{\rm BB1}\sim$ 120 eV and kT $_{\rm BB2}\sim$ 460 eV. The warm thermal component is likely due to the presence of an optically thick and warm Comptonizing plasma in the inner accretion disk. The fractional variability amplitude is found to be decreasing with increasing wavelength, i.e., from the soft X-ray to UV/optical emission. However, the hard X-ray (2–8 keV) emission shows very low variability. The strength of the correlation within the UV and the optical bands (0.95–0.99) is found to be stronger than the correlation between the UV/optical and X-ray bands (0.40–0.53). These results clearly suggest that the emitting regions of the X-ray and UV/optical emission are likely distinct or partly interacting. Having removed the slow variations in the light curves, we find that the lag spectrum is well described by the 4/3 rule for the standard Shakura–Sunyaev accretion disk when we omit X-ray lags. All these results suggest that the real disk is complex, and the UV emission is likely reprocessed in the accretion disk to give X-ray and optical emission.

A comparison of Galactic electron density models using PyGEDM

Galactic electron density distribution models are crucial tools for estimating the impact of the ionised interstellar medium on the impulsive signals from radio pulsars and fast radio bursts. The two prevailing Galactic electron density models (GEDMs) are YMW16 (Yao et al. 2017, ApJ, 835, 29) and NE2001 (Cordes & Lazio 2002, arXiv e-prints, pp astro–ph/0207156). Here, we introduce a software package PyGEDM which provides a unified application programming interface for these models and the YT20 (Yamasaki & Totani 2020, ApJ, 888, 105) model of the Galactic halo. We use PyGEDM to compute all-sky maps of Galactic dispersion measure (DM) for YMW16 and NE2001 and compare the large-scale differences between the two. In general, YMW16 predicts higher DM values towards the Galactic anticentre. YMW16 predicts higher DMs at low Galactic latitudes, but NE2001 predicts higher DMs in most other directions. We identify lines of sight for which the models are most discrepant, using pulsars with independent distance measurements. YMW16 performs better on average than NE2001, but both models show significant outliers. We suggest that future campaigns to determine pulsar distances should focus on targets where the models show large discrepancies, so future models can use those measurements to better estimate distances along those line of sight. We also suggest that the Galactic halo should be considered as a component in future GEDMs, to avoid overestimating the Galactic DM contribution for extragalactic sources such as FRBs.

The nature and likely redshift of GLEAM J0917–0012

We previously reported a putative detection of a radio galaxy at $z=10.15$ , selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The redshift of this source, GLEAM J0917–0012, was based on three weakly detected molecular emission lines observed with the Atacama Large Millimetre Array (ALMA). In order to confirm this result, we conducted deep spectroscopic follow-up observations with ALMA and the Karl Jansky Very Large Array (VLA). The ALMA observations targeted the same CO lines previously reported in Band 3 (84–115 GHz) and the VLA targeted the CO(4-3) and [CI(1-0)] lines for an independent confirmation in Q-band (41 and 44 GHz). Neither observation detected any emission lines, removing support for our original interpretation. Adding publicly available optical data from the Hyper Suprime-Cam survey, Widefield Infrared Survey Explorer (WISE), and Herschel Space Observatory in the infrared, as well as $<$ 10 GHz polarisation and 162 MHz inter-planetary scintillation observations, we model the physical and observational characteristics of GLEAM J0917–0012 as a function of redshift. Comparing these predictions and observational relations to the data, we are able to constrain its nature and distance. We argue that if GLEAM J0917–0012 is at $z<3,$ then it has an extremely unusual nature, and that the more likely solution is that the source lies above $z=7$ .

Early Science from POSSUM: Shocks, turbulence, and a massive new reservoir of ionised gas in the Fornax cluster

We present the first Faraday rotation measure (RM) grid study of an individual low-mass cluster—the Fornax cluster—which is presently undergoing a series of mergers. Exploiting commissioning data for the POlarisation Sky Survey of the Universe’s Magnetism (POSSUM) covering a ${\sim}34$ square degree sky area using the Australian Square Kilometre Array Pathfinder (ASKAP), we achieve an RM grid density of ${\sim}25$ RMs per square degree from a 280-MHz band centred at 887 MHz, which is similar to expectations for forthcoming GHz-frequency ${\sim}3\pi$ -steradian sky surveys. These data allow us to probe the extended magnetoionic structure of the cluster and its surroundings in unprecedented detail. We find that the scatter in the Faraday RM of confirmed background sources is increased by $16.8\pm2.4$ rad m−2 within 1 $^\circ$ (360 kpc) projected distance to the cluster centre, which is 2–4 times larger than the spatial extent of the presently detectable X-ray-emitting intracluster medium (ICM). The mass of the Faraday-active plasma is larger than that of the X-ray-emitting ICM and exists in a density regime that broadly matches expectations for moderately dense components of the Warm-Hot Intergalactic Medium. We argue that forthcoming RM grids from both targeted and survey observations may be a singular probe of cosmic plasma in this regime. The morphology of the global Faraday depth enhancement is not uniform and isotropic but rather exhibits the classic morphology of an astrophysical bow shock on the southwest side of the main Fornax cluster, and an extended, swept-back wake on the northeastern side. Our favoured explanation for these phenomena is an ongoing merger between the main cluster and a subcluster to the southwest. The shock’s Mach angle and stand-off distance lead to a self-consistent transonic merger speed with Mach 1.06. The region hosting the Faraday depth enhancement also appears to show a decrement in both total and polarised radio emission compared to the broader field. We evaluate cosmic variance and free-free absorption by a pervasive cold dense gas surrounding NGC 1399 as possible causes but find both explanations unsatisfactory, warranting further observations. Generally, our study illustrates the scientific returns that can be expected from all-sky grids of discrete sources generated by forthcoming all-sky radio surveys.

Search for intermittent X-ray pulsations from neutron stars in low-mass X-ray binaries

We present the results of our extensive binary orbital motion corrected pulsation search for 13 low-mass X-ray binaries. These selected sources exhibit burst oscillations in X-rays with frequencies ranging from 45 to 1 122 Hz and have a binary orbital period varying from 2.1 to 18.9 h. We first determined episodes that contain weak pulsations around the burst oscillation frequency by searching all archival Rossi X-ray Timing Explorer data of these sources. Then, we applied Doppler corrections to these pulsation episodes to discard the smearing effect of the binary orbital motion and searched for recovered pulsations at the second stage. Here we report 75 pulsation episodes that contain weak but coherent pulsations around the burst oscillation frequency. Furthermore, we report eight new episodes that show relatively strong pulsations in the binary orbital motion corrected data.

Processing GOTO data with the Rubin Observatory LSST Science Pipelines I: Production of coadded frames

The past few decades have seen the burgeoning of wide-field, high-cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy; however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. One such example is the Gravitational-wave Optical Transient Observer (GOTO) whose primary science objective is the optical follow-up of gravitational wave events. The amount and rate of data production by GOTO and other wide-area, high-cadence surveys presents a significant challenge to data processing pipelines which need to operate in near-real time to fully exploit the time domain. In this study, we adapt the Rubin Observatory LSST Science Pipelines to process GOTO data, thereby exploring the feasibility of using this ‘off-the-shelf’ pipeline to process data from other wide-area, high-cadence surveys. In this paper, we describe how we use the LSST Science Pipelines to process raw GOTO frames to ultimately produce calibrated coadded images and photometric source catalogues. After comparing the measured astrometry and photometry to those of matched sources from PanSTARRS DR1, we find that measured source positions are typically accurate to subpixel levels, and that measured L-band photometries are accurate to $\sim50$ mmag at $m_L\sim16$ and $\sim200$ mmag at $m_L\sim18$ . These values compare favourably to those obtained using GOTO’s primary, in-house pipeline, gotophoto, in spite of both pipelines having undergone further development and improvement beyond the implementations used in this study. Finally, we release a generic ‘obs package’ that others can build upon, should they wish to use the LSST Science Pipelines to process data from other facilities.