We study in this paper bright late-type giants found in the First Byurakan Survey (FBS) data base. Phase dependent light-curves from large sky area variability data bases such as Catalina Sky Survey (CSS) and All-Sky Automated Survey for Supernovae (ASAS-SN), and the early installment of the third Gaia data release (Gaia EDR3) photometric and astrometric data have been used to characterize our sample of 1 100 M-type giants and 130 C-type stars found at high latitudes. Gaia radial velocities (RV) are available for 134 and luminosities for 158 stars out of 1 100. We show the behaviour of our sample stars in a Gaia color–absolute magnitude diagram (CaMD), the Gaia-2MASS-diagram from Lebzelter et al. with some alternative versions. In this way we explore the potential of these diagrams and their combination for the analysis and interpretation of datasets of LPVs. We show the possibility to classify stars into M- and C-types and to identify the mass of the bulk of the sample stars.

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$ .

The recent increase in well-localised fast radio bursts (FRBs) has facilitated in-depth studies of global FRB host properties, the source circumburst medium, and the potential impacts of these environments on the burst properties. The Australian Square Kilometre Array Pathfinder (ASKAP) has localised 11 FRBs with sub-arcsecond to arcsecond precision, leading to sub-galaxy localisation regions in some cases and those covering much of the host galaxy in others. The method used to astrometrically register the FRB image frame for ASKAP, in order to align it with images taken at other wavelengths, is currently limited by the brightness of continuum sources detected in the short-duration (‘snapshot’) voltage data captured by the Commensal Real-Time ASKAP Fast Transients (CRAFT) software correlator, which are used to correct for any frame offsets due to imperfect calibration solutions and estimate the accuracy of any required correction. In this paper, we use dedicated observations of bright, compact radio sources in the low- and mid-frequency bands observable by ASKAP to investigate the typical astrometric accuracy of the positions obtained using this so-called ‘snapshot’ technique. Having captured these data with both the CRAFT software and ASKAP hardware correlators, we also compare the offset distributions obtained from both data products to estimate a typical offset between the image frames resulting from the differing processing paths, laying the groundwork for future use of the longer duration, higher signal-to-noise ratio (S/N) data recorded by the hardware correlator. We find typical offsets between the two frames of ${\sim}0.6$ and ${\sim}0.3$ arcsec in the low- and mid-band data, respectively, for both RA and Dec. We also find reasonable agreement between our offset distributions and those of the published FRBs. We detect only a weak dependence in positional offset on the relative separation in time and elevation between target and calibrator scans, with the trends being more pronounced in the low-band data and in Dec. Conversely, the offsets show a clear dependence on frequency in the low band, which we compare to the frequency-dependent Dec. offsets found in FRB 200430. In addition, we present a refined methodology for estimating the overall astrometric accuracy of CRAFT FRBs.

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.

Light curves of photospheric radius expansion (PRE) bursts, a subset of type I X-ray bursts, have been used as standard candles to estimate the ‘nominal PRE distances’ for 63% of PRE bursters (bursters), assuming PRE burst emission is spherically symmetric. Model-independent geometric parallaxes of bursters provide a valuable chance to test models of PRE bursts (PRE models) and can be provided in some cases by Gaia astrometry of the donor stars in bursters. We searched for counterparts to 115 known bursters in the Gaia Early Data Release 3 and confirmed 4 bursters with Gaia counterparts that have detected ( $\gt\!3\,\sigma$ , prior to zero-point correction) parallaxes. We describe a generic approach to the Gaia parallax zero point as well as its uncertainty using an ensemble of Gaia quasars individually determined for each target. Assuming the spherically symmetric PRE model is correct, we refined the resultant nominal PRE distances of three bursters (i.e. Cen $\textrm{X}-4$ , Cyg $\textrm{X}-2$ , and $4\textrm{U}\,0919-54$ ) and put constraints on their compositions of the nuclear fuel powering the bursts. Finally, we describe a method for testing the correctness of the spherically symmetric PRE model using parallax measurements and provide preliminary results.

We present a collection of double-lobed sources towards a $20\,\mathrm{deg}^2$ area of the Cygnus region at the northern sky, observed at 325 and 610 MHz with the Giant Metrewave Radio Telescope. The 10 $^{\prime\prime}$ resolution achieved at 325 MHz is 5.5 times better than previous studies, while at 610 MHz, these are the first results ever of such a large area, mapped with 6 $^{\prime\prime}$ angular resolution. After a thorough visual inspection of the images at the two bands, we found 43 double-lobed source candidates, proposed as such due to the presence of 2 bright peaks, within a few arcminutes apart, joined by a bridge or a central nucleus. All but two are presented here as a double-lobed candidates for the first time. Thirty nine of the candidates were covered at both bands, and we provide the spectral index information for them. We have searched for positional coincidences between the detected sources/components and other objects from the literature, along the electromagnetic spectrum. Twenty-three candidates possess radio counterpart(s), 12 present infrared counterparts, and 1 showed an overlapping X-ray source. We analysed each candidate considering morphology, counterparts, and spectral indices. Out of the 43 candidates, 37 show characteristics compatible with an extragalactic nature, 2 of probably Galactic origin, 3 remain as dubious cases, though with feature(s) compatible with an extragalactic nature, and the remaining one, evidence of physically unrelated components. The median spectral index of the 40 putative extragalactic sources is $-1.0$ . Their celestial surface density at 610 MHz resulted in $1.9\,\mathrm{per\ deg}^2$ , across a region lying at the Galactic plane.

We present the data and initial results from the first pilot survey of the Evolutionary Map of the Universe (EMU), observed at 944 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The survey covers $270 \,\mathrm{deg}^2$ of an area covered by the Dark Energy Survey, reaching a depth of 25–30 $\mu\mathrm{Jy\ beam}^{-1}$ rms at a spatial resolution of $\sim$ 11–18 arcsec, resulting in a catalogue of $\sim$ 220 000 sources, of which $\sim$ 180 000 are single-component sources. Here we present the catalogue of single-component sources, together with (where available) optical and infrared cross-identifications, classifications, and redshifts. This survey explores a new region of parameter space compared to previous surveys. Specifically, the EMU Pilot Survey has a high density of sources, and also a high sensitivity to low surface brightness emission. These properties result in the detection of types of sources that were rarely seen in or absent from previous surveys. We present some of these new results here.

We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535–571, which underwent a prolonged outburst beginning on 2017 September 2. We monitored MAXI J1535–571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 2017 September 20 to 2017 October 14. The source was quasi-simultaneously observed over the frequency range 0.84–19 GHz by UTMOST (the Upgraded Molonglo Observatory Synthesis Telescope) the Australian Square Kilometre Array Pathfinder (ASKAP), the Australia Telescope Compact Array (ATCA), and the Australian Long Baseline Array (LBA). Using the LBA observations from 2017 September 23, we measured the source size to be $34\pm1$ mas. During the brightest radio flare on 2017 September 21, the source was detected down to 119 MHz by the MWA, and the radio spectrum indicates a turnover between 250 and 500 MHz, which is most likely due to synchrotron self-absorption (SSA). By fitting the radio spectrum with a SSA model and using the LBA size measurement, we determined various physical parameters of the jet knot (identified in ATCA data), including the jet opening angle ( $\phi_{\rm op} = 4.5\pm1.2^{\circ}$ ) and the magnetic field strength ( $B_{\rm s} = 104^{+80}_{-78}$ mG). Our fitted magnetic field strength agrees reasonably well with that inferred from the standard equipartition approach, suggesting the jet knot to be close to equipartition. Our study highlights the capabilities of the Australian suite of radio telescopes to jointly probe radio jets in black hole X-ray binaries via simultaneous observations over a broad frequency range, and with differing angular resolutions. This suite allows us to determine the physical properties of X-ray binary jets. Finally, our study emphasises the potential contributions that can be made by the low-frequency part of the Square Kilometre Array (SKA-Low) in the study of black hole X-ray binaries.

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.

Equilibrium configurations of the internal magnetic field of a pulsar play a key role in modelling astrophysical phenomena from glitches to gravitational wave emission. In this paper, we present a numerical scheme for solving the Grad–Shafranov equation and calculating equilibrium configurations of pulsars, accounting for superconductivity in the core of the neutron star, and for the Hall effect in the crust of the star. Our numerical code uses a finite difference method in which the source term appearing in the Grad–Shafranov equation, which is used to model the magnetic equilibrium is non-linear. We obtain solutions by linearising the source and applying an under-relaxation scheme at each step of computation to improve the solver’s convergence. We have developed our code in both C++ and Python, and our numerical algorithm can further be adapted to solve any non-linear PDEs appearing in other areas of computational astrophysics. We produce mixed toroidal–poloidal field configurations, and extend the portion of parameter space that can be investigated with respect to previous studies. We find that in even in the more extreme cases, the magnetic energy in the toroidal component does not exceed approximately 5% of the total. We also find that if the core of the star is superconducting, the toroidal component is entirely confined to the crust of the star, which has important implications for pulsar glitch models which rely on the presence of a strong toroidal field region in the core of the star, where superfluid vortices pin to superconducting fluxtubes.

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.

The radio sky at lower frequencies, particularly below 20 MHz, is expected to be a combination of increasingly bright non-thermal emission and significant absorption from intervening thermal plasma. The sky maps at these frequencies cannot therefore be obtained by simple extrapolation of those at higher frequencies. However, due to severe constraints in ground-based observations, this spectral window still remains greatly unexplored. In this paper, we propose and study, through simulations, a novel minimal configuration for a space interferometer system which would enable imaging of the radio sky at frequencies well below 20 MHz with angular resolutions comparable to those achieved at higher radio frequencies in ground-based observations by using the aperture synthesis technique. The minimal configuration consists of three apertures aboard Low Earth Orbit (LEO) satellites orbiting the Earth in mutually orthogonal orbits. Orbital periods for the satellites are deliberately chosen to differ from each other so as to obtain maximum (u,v) coverage in short time spans with baselines greater than 15 000 km, thus, giving us angular resolutions finer than 10 arcsec even at these low frequencies. The sensitivity of the (u,v) coverage is assessed by varying the orbit and the initial phase of the satellites. We discuss the results obtained from these simulations and highlight the advantages of such a system.

One of the key open questions in extragalactic astronomy is what stops star formation in galaxies. While it is clear that the cold gas reservoir, which fuels the formation of new stars, must be affected first, how this happens and what are the dominant physical mechanisms involved is still a matter of debate. At least for satellite galaxies, it is generally accepted that internal processes alone cannot be responsible for fully quenching their star formation, but that environment should play an important, if not dominant, role. In nearby clusters, we see examples of cold gas being removed from the star-forming discs of galaxies moving through the intracluster medium, but whether active stripping is widespread and/or necessary to halt star formation in satellites, or quenching is just a consequence of the inability of these galaxies to replenish their cold gas reservoirs, remains unclear. In this work, we review the current status of environmental studies of cold gas in star-forming satellites in the local Universe from an observational perspective, focusing on the evidence for a physical link between cold gas stripping and quenching of the star formation. We find that stripping of cold gas is ubiquitous in satellite galaxies in both group and cluster environments. While hydrodynamical mechanisms such as ram pressure are important, the emerging picture across the full range of dark matter halos and stellar masses is a complex one, where different physical mechanisms may act simultaneously and cannot always be easily separated. Most importantly, we show that stripping does not always lead to full quenching, as only a fraction of the cold gas reservoir might be affected at the first pericentre passage. We argue that this is a key point to reconcile apparent tensions between statistical and detailed analyses of satellite galaxies, as well as disagreements between various estimates of quenching timescales. We conclude by highlighting several outstanding questions where we expect to see substantial progress in the coming decades, thanks to the advent of the Square Kilometre Array and its precursors, as well as the next-generation optical and millimeter facilities.

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.

Recent observations using several different telescopes and sky surveys showed patterns of asymmetry in the distribution of galaxies by their spin directions as observed from Earth. These studies were done with data imaged from the Northern hemisphere, showing excellent agreement between different telescopes and different analysis methods. Here, data from the DESI Legacy Survey was used. The initial dataset contains $\sim\!2.2\times 10^7$ galaxy images, reduced to $\sim\!8.1\times 10^5$ galaxies annotated by their spin direction using a symmetric algorithm. That makes it not just the first analysis of its kind in which the majority of the galaxies are in the Southern hemisphere, but also by far the largest dataset used for this purpose to date. The results show strong agreement between opposite parts of the sky, such that the asymmetry in one part of the sky is similar to the inverse asymmetry in the corresponding part of the sky in the opposite hemisphere. Fitting the distribution of galaxy spin directions to cosine dependence shows a dipole axis with probability of 4.66 $\sigma$ . Interestingly, the location of the most likely axis is within close proximity to the CMB Cold Spot. The profile of the distribution is nearly identical to the asymmetry profile of the distribution identified in Pan-STARRS, and it is within 1 $\sigma$ difference from the distribution profile in SDSS and HST. All four telescopes show similar large-scale profile of asymmetry.

In order to further develop and implement novel drift scan imaging experiments to undertake wide-field, high time resolution surveys for millisecond optical transients, an appropriate telescope drive system is required. This paper describes the development of a simple and inexpensive hardware and software system to monitor, characterise, and correct the primary category of telescope drive errors and periodic errors due to imperfections in the drive and gear chain. A model for the periodic errors is generated from direct measurements of the telescope drive shaft rotation, verified by comparison to astronomical measurements of the periodic errors. The predictive model is generated and applied in real time in the form of corrections to the drive rate. A demonstration of the system shows that inherent periodic errors of peak-to-peak amplitude ${\sim}{100}''$ are reduced to below the seeing limit of ${\sim}3''$ . This demonstration allowed an estimate of the uncertainties on the transient sensitivity timescales of the prototype survey of Tingay $\&$ Joubert (2021), with the nominal timescale sensitivity of 21 ms revised to be in the range of $20\!-\!22$ ms, which does not significantly affect the results of the experiment. The correction system will be adopted into the final version of high-cadence imaging experiment, which is currently under construction. The correction system is inexpensive ( $<\!{\$}$ A100) and composed of readily available hardware and is readily adaptable to other applications. Design details and codes are therefore made publicly available.

In 1978, Bracewell suggested the technique of nulling interferometry to directly image exoplanets which would enable characterisation of their surfaces, atmospheres, weather, and possibly determine their capacity to host life. The contrast needed to discriminate starlight reflected by a terrestrial-type planet from the glare of its host star lies at or beyond a forbidding $10^{-10}$ for an exo-Earth in the habitable zone around a Sun-like star at near-infrared wavelengths, necessitating instrumentation with extremely precise control of the light. Guided Light Interferometric Nulling Technology (GLINT) is a testbed for new photonic devices conceived to overcome the challenges posed by nulling interferometry. At its heart, GLINT employs a single-mode nulling photonic chip fabricated by direct-write technology to coherently combine starlight from an arbitrarily large telescope at 1 550 nm. It operates in combination with an actuated segmented mirror in a closed-loop control system, to produce and sustain a deep null throughout observations. The GLINT South prototype interfaces the 3.9-m Anglo-Australian Telescope and was tested on a sample of bright Mira variable stars. Successful and continuous starlight injection into the photonic chip was achieved. A statistical model of the data was constructed, enabling a data reduction algorithm to retrieve contrast ratios of about $10^{-3}$ . As a byproduct of this analysis, stellar angular diameters that were below the telescope diffraction limit ( $\sim$ 100 mas) were recovered with 1 $\sigma$ accuracy and shown to be in agreement with literature values despite working in the seeing-limited regime. GLINT South serves as a demonstration of the capability of direct-write photonic technology for achieving coherent, stable nulling of starlight, which will encourage further technological developments towards the goal of directly imaging exoplanets with future large ground based and space telescopes.