We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of g > 23.5 with a relative photometric calibration uncertainty of < 1 \%. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.

Jets from accreting black hole (BH) X-ray binaries (XRBs) are powerful outflows that release a large fraction of the accretion energy to the surrounding environment, providing a feedback mechanism that may alter the interstellar medium (ISM) properties. Studying accretion and feedback together enables estimates of matter and energy input/output around accreting BHs. We focus on the extended jet structures of the BH-XRB GRS1758-258. First seen in VLA data, these parsec-scale jets arise from jet-ISM interaction and show a Z-shaped morphology. Using the MeerKAT telescope we observed GRS1758-258 in L-band for a total exposure of 7 hr. Applying a calorimetry-based method developed for AGN and later used for XRBs, we estimated the properties of the jets and of the surrounding ISM. We detect a jet and counter-jet terminating in bow-shocks. Within the northern jet lobe we identify synchrotron and bremsstrahlung emission, while the southern lobe is dominated by thermal emission. We measure ISM densities between 10-40 cm-3 across both jets, slightly lower in the northern region. The estimated ages of the two lobes range from 6-51 kyr. The time-averaged jet power lies between 4.4x10^33 and 3.3x10^36 erg/s, with differences between north and south likely due to different local ISM conditions. Comparing new MeerKAT with archival VLA data, we measured a proper motion of 130 mas/yr in a portion of the northern jet. Jet-ISM interaction structures on both sides of GRS1758-258 reveal different ISM properties. The comparison between these structures and those from other XRBs suggests that the lobes in GRS1758-258 are younger and may result from different jet activity phases. The time-averaged energy transferred to the environment is slightly lower than in other XRBs, consistent with the younger age of the lobes in GRS1758-258 relative to those of other systems.

We revisit a class of coupled dark energy models where dark energy interacts with dark matter via phenomenological energy exchange terms. We include the perturbation of the Hubble rate in the interaction term, which was absent in previous works. We also consider more recent data sets such as cosmic microwave background (CMB) anisotropies from \text{Planck} 2018, type I-a supernovae (SNIa) measurements from Pantheon+ and data from baryon acoustic oscillations (BAO), and redshift space distortions (RSD). One of the models presents a strong incompatibility when different cosmological datasets are used. We analyze the influence of the SH0ES Cepheid host distances on the results and, although for one model the discrepancy of $H_0$ is reduced to $3.2\sigma$ when compared to the value obtained by the \textit{Planck} collaboration and $3.6\sigma$ when compared to the SH0ES team, joint analysis is incompatible. Including BAO with RSD shows incompatibility with SH0ES for all models considered here. We perform a model comparison and although there is a preference for interacting dark energy over $\Lambda$CDM for all the models for joint analysis CMB+BAO+RSD+SNIa, most of the 2D contours do not overlap. We conclude that the models of interactions in the dark sector considered in this paper are not flexible enough to fit all the cosmological data including values of $H_0$ from SH0ES in a statistically acceptable way. Therefore, the addition of one extra degree of freedom (i.e. the coupling between dark matter and dark matter) does not help enough to alleviate the already existing tension in the vanilla $\Lambda$CDM, suggesting that the models would need to be modified to include further flexibility of predictions to help elucidate this issue.

The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total intensity and spectral index mosaics of the region from 20 pointings (144 hours on-target in total), covering 6.5 square degrees with an angular resolution of 4$"$,at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT's high sensitivity, exceptional $u,v$-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new non-thermal filament complexes, and enhanced views of the Radio Arc Bubble, Sgr A and Sgr B regions. This project is a SARAO public legacy survey, and the image products are made available with this article.

We present a novel application of cosmological rescaling, or "remapping," to generate 21 cm intensity mapping mocks for different cosmologies. The remapping method allows for computationally efficient generation of N-body catalogs by rescaling existing simulations. In this work, we employ the remapping method to construct dark matter halo catalogs, starting from the Horizon Run 4 simulation with WMAP5 cosmology, and apply it to different target cosmologies, including WMAP7, Planck18 and Chevallier-Polarski-Linder (CPL) models. These catalogs are then used to simulate 21 cm intensity maps. We use the halo occupation distribution (HOD) method to populate halos with neutral hydrogen (HI) and derive 21 cm brightness temperature maps. Our results demonstrate the effectiveness of the remapping approach in generating cosmological simulations for large-scale structure studies, offering an alternative for testing observational data pipelines and performing cosmological parameter forecasts without the need for computationally expensive full N-body simulations. We also analyze the precision and limitations of the remapping, in light of the rescaling parameters $s$ and $s_m$, as well as the effects of the halo mass and box size thresholds.

We investigate the capabilities of various stages of the SKA to perform world-leading weak gravitational lensing surveys. We outline a way forward to develop the tools needed for pursuing weak lensing in the radio band. We identify the key analysis challenges and the key pathfinder experiments that will allow us to address them in the run up to the SKA. We identify and summarize the unique and potentially very powerful aspects of radio weak lensing surveys, facilitated by the SKA, that can solve major challenges in the field of weak lensing. These include the use of polarization and rotational velocity information to control intrinsic alignments, and the new area of weak lensing using intensity mapping experiments. We show how the SKA lensing surveys will both complement and enhance corresponding efforts in the optical wavebands through cross-correlation techniques and by way of extending the reach of weak lensing to high redshift.

This work proposes a saliency-based attribution framework to evaluate and compare 10 state-of-the-art explainability methods for deep learning models in astronomy, focusing on the classification of radio galaxy images. While previous work has primarily emphasized classification accuracy, we prioritize model interpretability. Qualitative assessments reveal that Score-CAM, Grad-CAM, and Grad-CAM++ consistently produce meaningful attribution maps, highlighting the brightest regions of FRI and FRII galaxies in alignment with known astrophysical features. In contrast, other methods often emphasize irrelevant or noisy areas, reducing their effectiveness.

We discover what is in projection the largest known structure of galactic origin: a giant radio galaxy with a projected proper length of $4.99 \pm 0.04\ \mathrm{Mpc}$. The source, named Alcyoneus, was first identified in low-resolution LOFAR Two-metre Sky Survey images from which angularly compact sources had been removed. Being an extreme example in its class, Alcyoneus could shed light on the main mechanisms that drive radio galaxy growth. We find that - beyond geometry - Alcyoneus and its host galaxy appear suspiciously ordinary: the total low-frequency luminosity density, stellar mass and supermassive black hole mass are all lower than, though similar to, those of the medial giant radio galaxy (percentiles $45 \pm 3\%$, $25 \pm 9 \%$ and $23 \pm 11 \%$, respectively). The source resides in a filament of the Cosmic Web, with which it might have significant thermodynamic interaction. At $5 \cdot 10^{-16}\ \mathrm{Pa}$, the pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents one of the most promising radio galaxies yet to probe the warm-hot intergalactic medium.

As the scientific exploitation of the Square Kilometre Array (SKA) approaches, there is a need for new advanced data analysis and visualization tools capable of processing large high-dimensional datasets. In this study, we aim to generalize the YOLO-CIANNA deep learning source detection and characterization method for 3D hyperspectral HI emission cubes. We present the adaptations we made to the regression-based detection formalism and the construction of an end-to-end 3D convolutional neural network (CNN) backbone. We then describe a processing pipeline for applying the method to simulated 3D HI cubes from the SKA Observatory Science Data Challenge 2 (SDC2) dataset. The YOLO-CIANNA method was originally developed and used by the MINERVA team that won the official SDC2 competition. Despite the public release of the full SDC2 dataset, no published result has yet surpassed MINERVA's top score. In this paper, we present an updated version of our method that improves our challenge score by 9.5%. The resulting catalog exhibits a high detection purity of 92.3%, best-in-class characterization accuracy, and contains 45% more confirmed sources than concurrent classical detection tools. The method is also computationally efficient, processing the full ~1TB SDC2 data cube in 30 min on a single GPU. These state-of-the-art results highlight the effectiveness of 3D CNN-based detectors for processing large hyperspectral data cubes and represent a promising step toward applying YOLO-CIANNA to observational data from SKA and its precursors.

The Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) has been working for a decade to reduce the time and cost of designing, building and deploying new digital radio-astronomy instruments. Today, CASPER open-source technology powers over 45 scientific instruments worldwide, and is used by scientists and engineers at dozens of academic institutions. In this paper we catalog the current offerings of the CASPER collaboration, and instruments past and present built by CASPER users and developers. We describe the ongoing state of software development, as CASPER looks to support a broader range of programming environments and hardware and ensure compatibility with the latest vendor tools.

We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $\mu$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status.

Online deception and financial scams represent a pervasive threat in the digital age, yet a quantitative analysis and understanding of their propagation is lacking. This study introduces a novel model based on the framework of epidemiological models to describe the interaction between scammers and their victims. We propose a five-compartment deterministic model ($S-V-R-A_s-R_s$) calibrated using longitudinal data in fraud reports from the Canadian Anti-Fraud Centre. The model's theoretical properties are established, including the non-negativity of the state variables and the stability threshold defined by the basic reproduction number ($\mathcal{R}_0$). A non-standard finite difference scheme is developed for the numerical simulations to ensure dynamical consistency between the continuous deterministic model and its discrete equivalent. A key finding of the model sensitivity analysis indicates that the proliferation of scams is overwhelmingly driven by the lifecycle of scammers, their recruitment, attrition, and arrest, rather than the susceptibility of the victim population. The results of this study provide strong quantitative evidence that the most effective control strategies are those that directly disrupt the scammers' population. Overall, this study provides a crucial model for designing and evaluating evidence-based policies to combat the scourge of cybercrime.

We report on the lowest-frequency detection to date of three bursts from the fast radio burst FRB 180916, observed at 328 MHz with the Sardinia Radio Telescope (SRT). The SRT observed the periodic repeater FRB 180916 for five days from 2020 February 20 to 24 during a time interval of active radio bursting, and detected the three bursts during the first hour of observations; no more bursts were detected during the remaining ~ 30 hours. Simultaneous SRT observations at 1548 MHz did not detect any bursts. Burst fluences are in the range 37 to 13 Jy ms. No relevant scattering is observed for these bursts. We also present the results of the multi-wavelength campaign we performed on FRB 180916, during the five days of the active window. Simultaneously with the SRT observations, others with different time spans were performed with the Northern Cross at 408 MHz, with XMM-Newton, NICER, INTEGRAL, AGILE, and with the TNG and two optical telescopes in Asiago, which are equipped with fast photometers. XMM-Newton obtained data simultaneously with the three bursts detected by the SRT, and determined a luminosity upper limit in the 0.3-10 keV energy range of ~$10^{45}$ erg/s for the burst emission. AGILE obtained data simultaneously with the first burst and determined a fluence upper limit in the MeV range for millisecond timescales of $10^{-8}$ erg cm$^{-2}$.Our results show that absorption from the circumburst medium does not significantly affect the emission from FRB 180916, thus limiting the possible presence of a superluminous supernova around the source, and indicate that a cutoff for the bursting mechanism, if present, must be at lower frequencies. Our multi-wavelength campaign sensitively constrains the broadband emission from FRB 180916, and provides the best limits so far for the electromagnetic response to the radio bursting of this remarkable source of fast radio bursts.

Weak gravitational lensing measurements are traditionally made at optical wavelengths where many highly resolved galaxy images are readily available. However, the Square Kilometre Array (SKA) holds great promise for this type of measurement at radio wavelengths owing to its greatly increased sensitivity and resolution over typical radio surveys. The key to successful weak lensing experiments is in measuring the shapes of detected sources to high accuracy. In this document we describe a simulation pipeline designed to simulate radio images of the quality required for weak lensing, and will be typical of SKA observations. We provide as input, images with realistic galaxy shapes which are then simulated to produce images as they would have been observed with a given radio interferometer. We exploit this pipeline to investigate various stages of a weak lensing experiment in order to better understand the effects that may impact shape measurement. We first show how the proposed SKA1-Mid array configurations perform when we compare the (known) input and output ellipticities. We then investigate how making small changes to these array configurations impact on this input-outut ellipticity comparison. We also demonstrate how alternative configurations for SKA1-Mid that are smaller in extent, and with a faster survey speeds produce similar performance to those originally proposed. We then show how a notional SKA configuration performs in the same shape measurement challenge. Finally, we describe ongoing efforts to utilise our simulation pipeline to address questions relating to how applicable current (mostly originating from optical data) shape measurement techniques are to future radio surveys. As an alternative to such image plane techniques, we lastly discuss a shape measurement technique based on the shapelets formalism that reconstructs the source shapes directly from the visibility data.

Radio weak lensing, while a highly promising complementary probe to optical weak lensing, will require incredible precision in the measurement of galaxy shape parameters. In this paper, we extend the Bayesian Inference for Radio Observations model fitting approach to measure galaxy shapes directly from visibility data of radio continuum surveys, instead of from image data. We apply a Hamiltonian Monte Carlo (HMC) technique for sampling the posterior, which is more efficient than the standard Monte Carlo Markov Chain method when dealing with a large dimensional parameter space. Adopting the exponential profile for galaxy model fitting allows us to analytically calculate the likelihood gradient required by HMC, allowing a faster and more accurate sampling. The method is tested on SKA1-MID simulated observations at 1.4 GHz of a field containing up to 1000 star-forming galaxies. It is also applied to a simulated observation of the weak lensing precursor survey SuperCLASS. In both cases we obtain reliable measurements of the galaxies' ellipticity and size for all sources with SNR $\ge 10$, and we also find relationships between the convergence properties of the HMC technique and some source parameters. Direct shape measurement in the visibility domain achieves high accuracy at the expected source number densities of the current and next SKA precursor continuum surveys. The proposed method can be easily extended for the fitting of other galaxy and scientific parameters, as well as simultaneously marginalising over systematic and instrumental effects.

Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,\nu}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results.

Radio astronomical imaging arrays comprising large numbers of antennas, O(10^2-10^3) have posed a signal processing challenge because of the required O(N^2) cross correlation of signals from each antenna and requisite signal routing. This motivated the implementation of a Packetized Correlator architecture that applies Field Programmable Gate Arrays (FPGAs) to the O(N) "F-stage" transforming time domain to frequency domain data, and Graphics Processing Units (GPUs) to the O(N^2) "X-stage" performing an outer product among spectra for each antenna. The design is readily scalable to at least O(10^3) antennas. Fringes, visibility amplitudes and sky image results obtained during field testing are presented.

Despite the evidence that supermassive black holes (SMBHs) co-evolve with their host galaxy, and that most of the growth of these SMBHs occurs via merger-free processes, the underlying mechanisms which drive this secular co-evolution are poorly understood. We investigate the role that both strong and weak large-scale galactic bars play in mediating this relationship. Using 72,940 disc galaxies in a volume-limited sample from Galaxy Zoo DESI, we analyse the active galactic nucleus (AGN) fraction in strongly barred, weakly barred, and unbarred galaxies up to z = 0.1 over a range of stellar masses and colours. After controlling for stellar mass and colour, we find that the optically selected AGN fraction is 31.6 +/- 0.9 per cent in strongly barred galaxies, 23.3 +/- 0.8 per cent in weakly barred galaxies, and 14.2 +/- 0.6 per cent in unbarred disc galaxies. These are highly statistically robust results, strengthening the tantalising results in earlier works. Strongly barred galaxies have a higher fraction of AGNs than weakly barred galaxies, which in turn have a higher fraction than unbarred galaxies. Thus, while bars are not required in order to grow a SMBH in a disc galaxy, large-scale galactic bars appear to facilitate AGN fuelling, and the presence of a strong bar makes a disc galaxy more than twice as likely to host an AGN than an unbarred galaxy at all galaxy stellar masses and colours.