Search

We present the radio afterglow of short gamma-ray burst (GRB) 230217A, which was detected less than 1 day after the gamma-ray prompt emission with the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array (VLA). The ATCA rapid-response system automatically triggered an observation of GRB 230217A following its detection by the Neil Gehrels Swift Observatory and began observing the event just 32 minutes post-burst at 5.5 and 9 GHz for 7 hours. Dividing the 7-hour observation into three time-binned images allowed us to obtain radio detections with logarithmic central times of 1, 2.8 and 5.2 hours post-burst, the first of which represents the earliest radio detection of any GRB to date. The decline of the light curve is consistent with reverse shock emission if the observing bands are below the spectral peak and not affected by synchrotron self-absorption. This makes GRB 230217A the fifth short GRB with radio detections attributed to a reverse shock at early times ($<1$ day post-burst). Following brightness temperature arguments, we have used our early radio detections to place the highest minimum Lorentz factor (${\Gamma}_{min} > 50$ at $\sim1$ hour) constraints on a GRB in the radio band. Our results demonstrate the importance of rapid radio follow-up observations with long integrations and good sensitivity for detecting the fast-evolving radio emission from short GRBs and probing their reverse shocks., Comment: 11 pages, 3 figures, submitted to ApJL

Recent observations have revealed rare, previously unknown flashes of cosmic radio waves lasting from milliseconds to minutes, and with periodicity of minutes to an hour [1-4]. These transient radio signals must originate from sources in the Milky Way, and from coherent emission processes in astrophysical plasma. They are theorised to be produced in the extreme and highly magnetised environments around white dwarfs or neutron stars [5-8]. However, the astrophysical origin of these signals remains contested, and multiple progenitor models may be needed to explain their diverse properties. Here we present the discovery of a transient radio source, ILT J1101+5521, whose roughly minute-long pulses arrive with a periodicity of 125.5 minutes. We find that ILT J1101+5521 is an M dwarf - white dwarf binary system with an orbital period that matches the period of the radio pulses, which are observed when the two stars are in conjunction. The binary nature of ILT J1101+5521 establishes that some long-period radio transients originate from orbital motion modulating the observed emission, as opposed to an isolated rotating star. We conclude that ILT J1101+5521 is likely a polar system where magnetic interaction has synchronised the rotational and orbital periods of the white dwarf [9]. Magnetic interaction and plasma exchange between the two stars may generate the sporadic radio emission. Such mechanisms have been previously theorised [10-13], but not observationally established.

Short gamma-ray bursts (GRBs) are explosive transients caused by binary mergers of compact objects containing at least one neutron star. Multi-wavelength afterglow observations provide constraints on the physical parameters of the jet, its surrounding medium, and the microphysics of the enhanced magnetic fields and accelerated electrons in the blast wave at the front of the jet. The synchrotron radio emission can be tracked for much longer than in other spectral regimes, and it can pin down the evolution of the spectral peak. We present the results of a systematic observing campaign of eight short GRBs with the MeerKAT radio telescope. Additionally, we present observations of four of these short GRBs using the ATCA radio telescope and two of these short GRBs with the e-MERLIN radio telescope. Using these results we report one possible detection of a short GRB afterglow from GRB 230217A and deep upper limits for the rest of our short GRB observations. We use these observations to place constraints on some of the physical parameters, in particular those related to electron acceleration, the circumburst density, and gamma-ray energy efficiency. We discuss how deeper observations with new and upgraded telescopes should be able to determine if the gamma-ray efficiency differs between long and short GRBs. We also report detections of the likely host galaxies for four of the eight GRBs and upper limits for another GRB, increasing the number of detected host galaxies in the radio with implications for the star formation rate in these galaxies., Comment: to be published in the Monthly Notices of the Royal Astronomical Society

The main science aim of the BlackGEM array is to detect optical counterparts to gravitational wave mergers. Additionally, the array will perform a set of synoptic surveys to detect Local Universe transients and short time-scale variability in stars and binaries, as well as a six-filter all-sky survey down to ~22nd mag. The BlackGEM Phase-I array consists of three optical wide-field unit telescopes. Each unit uses an f/5.5 modified Dall-Kirkham (Harmer-Wynne) design with a triplet corrector lens, and a 65cm primary mirror, coupled with a 110Mpix CCD detector, that provides an instantaneous field-of-view of 2.7~square degrees, sampled at 0.564\arcsec/pixel. The total field-of-view for the array is 8.2 square degrees. Each telescope is equipped with a six-slot filter wheel containing an optimised Sloan set (BG-u, BG-g, BG-r, BG-i, BG-z) and a wider-band 440-720 nm (BG-q) filter. Each unit telescope is independent from the others. Cloud-based data processing is done in real time, and includes a transient-detection routine as well as a full-source optimal-photometry module. BlackGEM has been installed at the ESO La Silla observatory as of October 2019. After a prolonged COVID-19 hiatus, science operations started on April 1, 2023 and will run for five years. Aside from its core scientific program, BlackGEM will give rise to a multitude of additional science cases in multi-colour time-domain astronomy, to the benefit of a variety of topics in astrophysics, such as infant supernovae, luminous red novae, asteroseismology of post-main-sequence objects, (ultracompact) binary stars, and the relation between gravitational wave counterparts and other classes of transients, Comment: Published in PASP

In this paper, we present rapid follow-up observations of the short GRB 201006A, consistent with being a compact binary merger, using the LOw Frequency ARray (LOFAR). We have detected a candidate 5.6$\sigma$, short, coherent radio flash at 144 MHz at 76.6 mins post-GRB with a 3$\sigma$ duration of 38 seconds. This radio flash is 27 arcsec offset from the GRB location, which has a probability of occurring by chance of $\sim$0.05% (3.8$\sigma$) when accounting for measurement uncertainties. Despite the offset, we show that the probability of finding an unrelated transient within 40 arcsec of the GRB location is $<10^{-6}$ and conclude that this is a candidate radio counterpart to GRB 201006A. We performed image plane dedispersion and the radio flash is tentatively (2.4$\sigma$) shown to be highly dispersed, allowing a distance estimate, corresponding to a redshift of $0.58\pm0.06$. The corresponding luminosity of the event at this distance is $6.7^{+6.6}_{-4.4} \times 10^{32}$ erg s$^{-1}$ Hz$^{-1}$. If associated with GRB 201006A, this emission would indicate prolonged activity from the central engine that is consistent with being a newborn, supramassive, likely highly magnetised, millisecond spin neutron star (a magnetar)., Comment: MNRAS Accepted

Measurement of the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn (CD) and Epoch of Reionisation (EoR) promises to unveil a wealth of information about the astrophysical processes during the first billion years of evolution of the universe. The AARTFAAC Cosmic Explorer (ACE) utilises the AARTFAAC wide-field imager of LOFAR to measure the power spectrum of the intensity fluctuations of the redshifted 21-cm signal from the CD at z~18. The RFI from various sources contaminates the observed data and it is crucial to exclude the RFI-affected data in the analysis for reliable detection. In this work, we investigate the impact of non-ground-based transient RFI using cross-power spectra and cross-coherence metrics to assess the correlation of RFI over time and investigate the level of impact of transient RFI on the ACE 21-cm power spectrum estimation. We detected moving sky-based transient RFI sources that cross the field of view within a few minutes and appear to be mainly from aeroplane communication beacons at the location of the LOFAR core in the 72-75 MHz band, by inspecting filtered images. This transient RFI is mostly uncorrelated over time and is only expected to dominate over the thermal noise for an extremely deep integration time of 3000 hours or more with a hypothetical instrument that is sky temperature dominated at 75 MHz. We find no visible correlation over different k-modes in Fourier space in the presence of noise for realistic thermal noise scenarios. We conclude that the sky-based transient RFI from aeroplanes, satellites and meteorites at present does not pose a significant concern for the ACE analyses at the current level of sensitivity and after integrating over the available 500 hours of observed data. However, it is crucial to mitigate or filter such transient RFI for more sensitive experiments aiming for significantly deeper integration., Comment: 13 pages, 9 figures, and 3 tables; accepted for publication in Astronomy and Astrophysics (A&A)

Context: GRB 211106A and GRB 211227A are recent gamma-ray bursts (GRBs) with initial X-ray positions suggesting associations with nearby galaxies (z < 0.7). Their prompt emission characteristics indicate GRB 211106A is a short-duration GRB and GRB 211227A is a short GRB with extended emission, likely originating from compact binary mergers. However, classifying solely based on prompt emission can be misleading. Aims: These short GRBs in the local Universe offer opportunities to search for associated kilonova (KN) emission and study host galaxy properties in detail. Methods: We conducted deep optical and NIR follow-up using ESO-VLT FORS2, HAWK-I, and MUSE for GRB 211106A, and ESO-VLT FORS2 and X-Shooter for GRB 211227A, starting shortly after the X-ray afterglow detection. We performed photometric analysis to look for afterglow and KN emissions associated with the bursts, along with host galaxy imaging and spectroscopy. Optical/NIR results were compared with Swift X-Ray Telescope (XRT) and other high-energy data. Results: For both GRBs we placed deep limits to the optical/NIR afterglow and KN emission. Host galaxies were identified: GRB 211106A at photometric z = 0.64 and GRB 211227A at spectroscopic z = 0.228. Host galaxy properties aligned with typical short GRB hosts. We also compared the properties of the bursts with the S-BAT4 sample to further examined the nature of these events. Conclusions: Study of prompt and afterglow phases, along with host galaxy analysis, confirms GRB 211106A as a short GRB and GRB 211227A as a short GRB with extended emission. The absence of optical/NIR counterparts is likely due to local extinction for GRB 211106A and a faint kilonova for GRB 211227A., Comment: Accepted to A&A on 08 August 2023, 21 pages, 24 figures

The composition of relativistic gamma-ray burst (GRB) jets and their emission mechanisms are still debated, and they could be matter or magnetically dominated. One way to distinguish these mechanisms arises because a Poynting flux dominated jet may produce low-frequency radio emission during the energetic prompt phase, through magnetic reconnection at the shock front. We present a search for radio emission coincident with three GRB X-ray flares with the LOw Frequency ARray (LOFAR), in a rapid response mode follow-up of long GRB 210112A (at z~2) with a 2 hour duration, where our observations began 511 seconds after the initial swift-BAT trigger. Using timesliced imaging at 120-168 MHz, we obtain upper limits at 3 sigma confidence of 42 mJy averaging over 320 second snapshot images, and 87 mJy averaging over 60 second snapshot images. LOFAR's fast response time means that all three potential radio counterparts to X-ray flares are observable after accounting for dispersion at the estimated source redshift. Furthermore, the radio pulse in the magnetic wind model was expected to be detectable at our observing frequency and flux density limits which allows us to disfavour a region of parameter space for this GRB. However, we note that stricter constraints on redshift and the fraction of energy in the magnetic field are required to further test jet characteristics across the GRB population., Comment: 13 pages, 7 figures, accepted for publication in MNRAS

The chemical enrichment of dust and metals in the interstellar medium (ISM) of galaxies throughout cosmic time is one of the key driving processes of galaxy evolution. Here we study the evolution of the gas-phase metallicities, dust-to-gas (DTG), and dust-to-metal (DTM) ratios of 36 star-forming galaxies at $1.7 < z < 6.3$ probed by gamma-ray bursts (GRBs). We compile all GRB-selected galaxies with intermediate (R=7000) to high (R>40,000) resolution spectroscopic data for which at least one refractory (e.g. Fe) and one volatile (e.g. S or Zn) element have been detected at S/N>3. This is to ensure that accurate abundances and dust depletion patterns can be obtained. We first derive the redshift evolution of the dust-corrected, absorption-line based gas-phase metallicity [M/H]$_{\rm tot}$ in these galaxies, for which we determine a linear relation with redshift ${\rm [M/H]_{tot}}(z) = (-0.21\pm 0.04)z -(0.47\pm 0.14)$. We then examine the DTG and DTM ratios as a function of redshift and through three orders of magnitude in metallicity, quantifying the relative dust abundance both through the direct line-of-sight visual extinction $A_V$ and the derived depletion level. We use a novel method to derive the DTG and DTM mass ratios for each GRB sightline, summing up the mass of all the depleted elements in the dust-phase. We find that the DTG and DTM mass ratios are both strongly correlated with the gas-phase metallicity and show a mild evolution with redshift as well. While these results are subject to a variety of caveats related to the physical environments and the narrow pencil-beam sightlines through the ISM probed by the GRBs, they provide strong implications for studies of dust masses to infer the gas and metal content of high-redshift galaxies, and particularly demonstrate the large offset from the average Galactic value in the low-metallicity, high-redshift regime., Comment: Accepted in A&A

Here, we present the results of our multi-wavelength campaign aimed at classifying \textit{Swift} J170800$-$402551.8 as part of the \textit{Swift} Deep Galactic Plane Survey (DGPS). We utilized Target of Opportunity (ToO) observations with \textit{Swift}, \textit{NICER}, \textit{XMM-Newton}, \textit{NuSTAR}, and the Southern African Large Telescope (SALT), as well as multi-wavelength archival observations from \textit{Gaia}, VPHAS, and VVV. The source displays a periodicity of 784 s in our \textit{XMM-Newton} observation. The X-ray spectrum (\textit{XMM-Newton} and \textit{NuSTAR}) can be described by thermal bremsstrahlung radiation with a temperature of $kT$\,$\approx$\,$30$ keV. The phase-folded X-ray lightcurve displays a double-peaked, energy-dependent pulse-profile. We used \textit{Chandra} to precisely localize the source, allowing us to identify and study the multi-wavelength counterpart. Spectroscopy with SALT identified a Balmer H$\alpha$ line, and potential HeI lines, from the optical counterpart. The faintness of the counterpart ($r$\,$\approx$\,$21$ AB mag) favors a low-mass donor star. Based on these criteria, we classify \textit{Swift} J170800$-$402551.8 as a candidate intermediate polar cataclysmic variable, where the spin period of the white dwarf is 784 s., Comment: Accepted in MNRAS

Transient surveys are a vital tool in exploring the dynamic universe, with radio transients acting as beacons for explosive and highly energetic astrophysical phenomena. However, performing commensal transient surveys using radio imaging can require a significant amount of computing power, data storage and time. With the instrumentation available to us, and with new and exciting radio interferometers in development, it is essential that we develop efficient methods to probe the radio transient sky. In this paper, we present results from an commensal short duration transient survey, on time scales of 8 seconds, 128 seconds and 1 hour, using data from the MeerKAT radio telescope. The dataset used was obtained as part of a galaxy observing campaign, and we focus on the field of NGC 5068. We present a quick, wide field imaging strategy to enable fast imaging of large datasets, and develop methods to efficiently filter detected transient candidates. No transient candidates were identified on the time scales of 8 seconds, 128 seconds and 1 hour, leading to competitive limits on the transient surface densities of $6.7{\times}10^{-5}$ deg$^{-1}$, $1.1{\times}10^{-3}$ deg$^{-1}$, and $3.2{\times}10^{-2}$ deg$^{-1}$ at sensitivities of 56.4 mJy, 19.2 mJy, and 3.9 mJy for the respective time scales. We find one possible candidate that could be associated with a stellar flare, that was rejected due to strict image quality control. Further short time-scale radio observations of this candidate could give definite results to its origin., Comment: 11 pages (9 main, 2 appendix), 8 figures, 2 tables. Submitted to MNRAS

The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March 2017 and May 2021. The Survey is complete to depth $L_X$\,$>$\,$10^{34}$ erg s$^{-1}$ to the edge of the Galaxy. The main Survey goal is to produce a rich sample of new X-ray sources and transients, while also covering a broad discovery space. Here, we introduce the Survey strategy and present a catalog of sources detected during Phase-I observations. In total, we identify 928 X-ray sources, of which 348 are unique to our X-ray catalog. We report on the characteristics of sources in our catalog and highlight sources newly classified and published by the DGPS team., Comment: Accepted in ApJS. This is the final version

We present the results of our X-ray, ultraviolet, and optical follow-up campaigns of 1RXS J165424.6-433758, an X-ray source detected with the \textit{Swift} Deep Galactic Plane Survey (DGPS). The source X-ray spectrum (\textit{Swift} and \textit{NuSTAR}) is described by thermal bremsstrahlung radiation with a temperature of $kT=10.1\pm1.2$ keV, yielding an X-ray ($0.3-10$ keV) luminosity $L_X=(6.5\pm0.8)\times10^{31}$ erg s$^{-1}$ at a \textit{Gaia} distance of 460 pc. Spectroscopy with the Southern African Large Telescope (SALT) revealed a flat continuum dominated by emission features, demonstrating an inverse Balmer decrement, the $\lambda4640$ Bowen blend, almost a dozen HeI lines, and HeII $\lambda4541$, $\lambda4686$ and $\lambda 5411$. Our high-speed photometry demonstrates a preponderance of flickering and flaring episodes, and revealed the orbital period of the system, $P_\textrm{orb}=2.87$ hr, which fell well within the cataclysmic variable (CV) period gap between $2-3$ hr. These features classify 1RXS J165424.6-433758 as a nearby polar magnetic CV., Comment: Accepted for publication in ApJ

We present JWST and Hubble Space Telescope (HST) observations of the afterglow of GRB 221009A, the brightest gamma-ray burst (GRB) ever observed. This includes the first mid-IR spectra of any GRB, obtained with JWST/NIRSPEC (0.6-5.5 micron) and MIRI (5-12 micron), 12 days after the burst. Assuming that the intrinsic spectral slope is a single power-law, with $F_{\nu} \propto \nu^{-\beta}$, we obtain $\beta \approx 0.35$, modified by substantial dust extinction with $A_V = 4.9$. This suggests extinction above the notional Galactic value, possibly due to patchy extinction within the Milky Way or dust in the GRB host galaxy. It further implies that the X-ray and optical/IR regimes are not on the same segment of the synchrotron spectrum of the afterglow. If the cooling break lies between the X-ray and optical/IR, then the temporal decay rates would only match a post jet-break model, with electron index $p<2$, and with the jet expanding into a uniform ISM medium. The shape of the JWST spectrum is near-identical in the optical/nIR to X-shooter spectroscopy obtained at 0.5 days and to later time observations with HST. The lack of spectral evolution suggests that any accompanying supernova (SN) is either substantially fainter or bluer than SN 1998bw, the proto-type GRB-SN. Our HST observations also reveal a disc-like host galaxy, viewed close to edge-on, that further complicates the isolation of any supernova component. The host galaxy appears rather typical amongst long-GRB hosts and suggests that the extreme properties of GRB 221009A are not directly tied to its galaxy-scale environment., Comment: Accepted for publication to the Astrophysical Journal Letters for the GRB 221009A Special Issue. The results of this paper are under press embargo until March 28, 18 UT. 19 pages, 8 figures, 2 tables

We observed the rapid radio brightening of GRB 210702A with the Australian Telescope Compact Array (ATCA) just 11hr post-burst, tracking early-time radio variability over a 5hr period on ~15min timescales at 9.0, 16.7, and 21.2GHz. A broken power-law fit to the 9.0GHz light curve showed that the 5hr flare peaked at a flux density of 0.4+/-0.1mJy at ~13hr post-burst with a steep rise and decline. The observed temporal and spectral evolution are not expected in the standard internal-external shock model, where forward and reverse shock radio emission evolves on much longer timescales. The early-time (<1day) optical and X-ray light curves from the Neil Gehrels Swift Observatory demonstrated typical afterglow forward shock behaviour, allowing us to use blast wave physics to determine a likely homogeneous circumburst medium and an emitting electron population power-law index of p=2.9+/-0.1. We suggest the early-time radio flare is likely due to weak interstellar scintillation (ISS), which boosted the radio afterglow emission above the ATCA sensitivity limit on minute timescales. Using weak ISS relations, we were able to place an upper limit on the size of the blast wave of $\leq6 \times 10^{16}$cm in the plane of the sky, which is consistent with the theoretical forward shock size prediction of $8\times10^{16}$cm for GRB 210702A at ~13h post-burst. This represents the earliest ISS size constraint on a GRB blast wave to date, demonstrating the importance of rapid (<1day) radio follow-up of GRBs using several-hour integrations to capture the early afterglow evolution, and to track scintillation over a broad frequency range., Comment: Submitted to MNRAS 13 pages, 6 figures, 2 tables

We investigate pre-merger coherent radio emission from neutron star mergers arising due to the magnetospheric interaction between compact objects. We consider two plausible radiation mechanisms, and show that if one neutron star has a surface magnetic field $B_{\rm s} \ge 10^{12}$G, coherent millisecond radio bursts with characteristic temporal morphology and inclination angle dependence are observable to Gpc distances with next-generation radio facilities. We explore multi-messenger and multi-wavelength methods of identification of a NS merger origin of radio bursts, such as in fast radio burst surveys, triggered observations of gamma-ray bursts and gravitational wave events, and optical/radio follow-up of fast radio bursts in search of kilonova and radio afterglow emission. We present our findings for current and future observing facilities, and make recommendations for verifying or constraining the model., Comment: Re-upload to match journal accepted version. Clarifying corrections included to subsection 3.2

We present LOFAR imaging observations from the April/May 2020 active episode of magnetar SGR 1935+2154. We place the earliest radio limits on persistent emission following the low-luminosity fast radio burst FRB 200428 from the magnetar. We also perform an image-plane search for transient emission and find no radio flares during our observations. We examine post-FRB radio upper limits in the literature and find that all are consistent with the multi-wavelength afterglow predicted by the synchrotron maser shock model interpretation of FRB 200428. However, early optical observations appear to rule out the simple versions of the afterglow model with constant-density circumburst media. We show that these constraints may be mitigated by adapting the model for a wind-like environment, but only for a limited parameter range. In addition, we suggest that late-time non-thermal particle acceleration occurs within the afterglow model when the shock is no longer relativistic, which may prove vital for detecting afterglows from other Galactic FRBs. We also discuss future observing strategies for verifying either magnetospheric or maser shock FRB models via rapid radio observations of Galactic magnetars and nearby FRBs., Comment: Accepted for publication in MNRAS

Many transient and variable sources detected at multiple wavelengths are also observed to vary at radio frequencies. However, these samples are typically biased towards sources that are initially detected in wide-field optical, X-ray or gamma-ray surveys. Many sources that are insufficiently bright at higher frequencies are therefore missed, leading to potential gaps in our knowledge of these sources and missing populations that are not detectable in optical, X-rays or gamma-rays. Taking advantage of new state-of-the-art radio facilities that provide high quality wide-field images with fast survey speeds, we can now conduct unbiased surveys for transient and variable sources at radio frequencies. In this paper, we present an unbiased survey using observations obtained by MeerKAT, a mid-frequency ($\sim$1.4 GHz) radio array in South Africa's Karoo Desert. The observations used were obtained as part of a weekly monitoring campaign for X-ray binaries (XRBs) and we focus on the field of MAXI J1820+070. We develop methods to optimally filter transient and variable candidates that can be directly applied to other datasets. In addition to MAXI J1820+070, we identify four likely active galactic nuclei, one source that could be a Galactic source (pulsar or quiescent X-ray binary) or an AGN, and one variable pulsar. No transient sources, defined as being undetected in deep images, were identified leading to a transient surface density of $<3.7\times10^{-2}$ deg$^{-2}$ at a sensitivity of 1 mJy on timescales of one week at 1.4 GHz., Comment: MNRAS Accepted, 18 pages, 17 figures

We present 21 new long-term variable radio sources found commensally in two years of weekly MeerKAT monitoring of the low-mass X-ray binary GX 339-4. The new sources vary on time scales of weeks to months and have a variety of light curve shapes and spectral index properties. Three of the new variable sources are coincident with multi-wavelength counterparts; and one of these is coincident with an optical source in deep MeerLICHT images. For most sources, we cannot eliminate refractive scintillation of active galactic nuclei as the cause of the variability. These new variable sources represent $2.2\pm0.5$ per cent of the unresolved sources in the field, which is consistent with the 1-2 per cent variability found in past radio variability surveys. However, we expect to find short-term variable sources in the field as well as these 21 new long-term variable sources. We present the radio light curves and spectral index variability of the new variable sources, as well as the absolute astrometry and matches to coincident sources at other wavelengths., Comment: 26 pages, 17 figures, 7 tables. Accepted for publication in MNRAS 2022 March 14. Received 2022 February 17; in original form 2021 July 11

Aims: Contamination from bright diffuse Galactic thermal and non-thermal radio emission poses crucial challenges in experiments aiming to measure the 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization. If not included in calibration, this diffuse emission can severely impact the analysis and signal extraction in 21-cm experiments. We examine large-scale diffuse Galactic emission at 122~MHz, around the North Celestial Pole, using the AARTFAAC-HBA system. Methods: In this pilot project, we present the first-ever wide-field image produced with a single sub-band of the data recorded with the AARTFAAC-HBA. We demonstrate two methods: multiscale CLEAN and shapelet decomposition, to model the diffuse emission revealed in the image. We use angular power spectrum metrics to quantify different components of the emission and compare the performance of the two diffuse structure modelling approaches. Results: We observe that the point sources dominate the angular power spectrum ($\ell(\ell+1)C_{\ell}/2\pi\equiv \Delta^2(\ell)$) of the emission in the field on scales $\ell\gtrsim 60$ ($\lesssim 3$~degree). The angular power spectrum after subtraction of compact sources is flat within $20\lesssim \ell \lesssim200$ range, suggesting that the residual power is dominated by the diffuse emission on scales $\ell\lesssim200$. The residual diffuse emission has a brightness temperature variance of $\Delta^2_{\ell=180} = (145.64 \pm 13.61)~{\rm K}^2$ at 122~MHz on angular scales of 1~degree, and is consistent with a power-law following $C_{\ell}\propto \ell^{-2.0}$ in $20\lesssim \ell \lesssim200$ range. We also find that, in the current setup, the multiscale CLEAN is suitable to model the compact and diffuse structures on a wide range of angular scales, whereas the shapelet decomposition method better models the large scales, which are of the order of a few degrees and wider., Comment: 11 pages, 9 figures, accepted (A&A)

We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore, by comparing the infrared counterpart's spectral energy distribution to known BeXRBs, we confirm this classification and set a distance of approximately $10-15$ kpc for the source. The source's broadband X-ray spectrum ($1.5-50$ keV) is described by an absorbed power-law with photon index $\Gamma$\,$\sim$\,$0.5$ and cutoff energy at $\sim$\,$13$ keV., Comment: Revised manuscript. Accepted to ApJ

We provide constraints on the nature of particle bunches that power fast radio bursts (FRBs) in the coherent curvature radiation model. It has been shown that current-induced perturbation to the motion of individual particles results in a high-energy, incoherent component of emission. We consider photo-magnetic interactions and show that the high-energy radiation can produce pairs which screen the accelerating electric field. We find that to avoid catastrophic cascades that quench emission, bunches capable of producing FRBs must have a modest density $n_e \approx 10^{13-14}\, {\rm cm^{-3}}$, and likely propagate along field lines with large curvature radii, $\rho > 10^8 \, {\rm cm}$. This rules out rapidly rotating magnetars as FRB sources within the coherent curvature radiation model., Comment: 9 pages, 9 figures. Submitting to MNRAS, comments welcome

High brightness temperature radio transients such as pulsars and fast radio bursts require the coherent radiation of particles. The antenna class of coherent radiation models require a large number of charged particles radiating in phase, therefore the particles must be spatially confined and have well-aligned velocities. Given these necessary conditions, we look at the magnetic field induced by the currents associated with coherently emitting accelerated particles and consider the interaction between the radiating particles and the induced magnetic field. We find a maximum luminosity of coherent curvature radiation that depends on source parameters such as surface magnetic field and neutron star spin period. We find that coherent radio emission across all luminosities can be explained by coherent curvature radiation and suggest it could be universally responsible for both FRBs and extreme galactic sources. Using the Crab Pulsar as an example, we constrain the emission parameters and origin of the most extreme nanoshots to within 60km of the pulsar surface assuming coherent curvature radiation. In agreement with recent observations, we also predict simultaneous X-ray emission from small-scale particle gyration due to the induced field., Comment: Accepted to MNRAS Letters

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

Gamma-ray bursts (GRBs) are ultra-relativistic collimated outflows, which emit synchrotron radiation throughout the entire electromagnetic spectrum when they interact with their environment. This afterglow emission enables us to probe the dynamics of relativistic blast waves, the microphysics of shock acceleration, and environments of GRBs. We perform Bayesian inference on a sample of GRB afterglow data sets consisting of 22 long GRBs and 4 short GRBs, using the afterglow model "scalefit", which is based on 2D relativistic hydrodynamic simulations. We make use of Gaussian processes to account for systematic deviations in the data sets, which allows us to obtain robust estimates for the model parameters. We present the inferred parameters for the sample of GRBs, and make comparisons between short GRBs and long GRBs in constant-density and stellar-wind-like environments. We find that in almost all respects such as energy and opening angle, short and long GRBs are statistically the same. Short GRBs however have a markedly lower prompt $\gamma$-ray emission efficiency than long GRBs. We also find that for long GRBs in ISM-like ambient media there is a significant anti-correlation between the fraction of thermal energy in the magnetic fields, $\epsilon_B$, and the beaming corrected kinetic energy. Furthermore, we find no evidence that the mass-loss rates of the progenitor stars are lower than those of typical Wolf-Rayet stars., Comment: Submitted to MNRAS. Accepted 2022 January 17. Received 2021 December 29; in original form 2021 June 28

The electromagnetic counterparts to gravitational wave (GW) merger events are highly sought after, but difficult to find owing to large localization regions. In this study, we present a strategy to search for compact object merger radio counterparts in wide-field data collected by the Low-Frequency Array (LOFAR). In particular, we use multi-epoch LOFAR observations centred at 144 MHz spanning roughly 300 deg$^2$ at optimum sensitivity of a since retracted neutron star-black hole merger candidate detected during O2, the second Advanced Ligo-Virgo GW observing run. The minimum sensitivity of the entire (overlapping) 1809 deg$^2$ field searched is 50 mJy and the false negative rate is 0.1 per cent above 200 mJy. We do not find any transients and thus place an upper limit at 95 per cent confidence of 0.02 transients per square degree above 20 mJy on one, two and three month timescales, which are the most sensitive limits available to date. Finally, we discuss the prospects of observing GW events with LOFAR in the upcoming GW observing run and show that a single multi-beam LOFAR observation can probe the full projected median localization area of binary neutron star mergers down to a median sensitivity of at least 8 mJy., Comment: Accepted version (in MNRAS)

Dense aperture arrays provide key benefits in modern astrophysical research. They are flexible, employing cheap receivers, while relying on the ever more sophisticated compute back-end to deal with the complexities of signal processing required for optimal use. Their advantage is that they offer very large fields of view and are readily scalable to any size, all other things being equal. Since they represent "software telescopes", the science cases these arrays can be applied to are quite broad. Here, we describe the calibration and performance of the AARTFAAC-12 instrument, which is composed of the twelve centrally located stations of the LOFAR array. We go into the details of the data acquisition and pre-processing, we describe the newly developed calibration pipeline as well as the noise properties of the resulting images and present radio source counts at 41.7 MHz and 61 MHz. We find that AARTFAAC-12 is confusion limited at 0.9 Jy/PSF at 61 MHz with a PSF size of 17'x11' and that the normalized source counts agree with the scaled VLSSr and 6C survey counts. The median spectral index of the sources between the two frequencies we observed at is -0.78. Further, we have used the derived source counts to estimate any excess cosmic radio background, and we do not find evidence for it at our observing frequencies compared to published literature values., Comment: Accepted for publication in MNRAS

Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT optical telescope in the days and months after the April 28 event. We did not detect any significant single radio pulses down to fluence limits between 25 mJy ms and 18 Jy ms. Some observing epochs overlapped with times when X-ray bursts were detected. Radio images made on four days using the MeerKAT telescope revealed no point-like persistent or transient emission at the location of the magnetar. No transient or persistent optical emission was detected over seven days. Using the multi-colour MeerLICHT images combined with relations between DM, NH and reddening we constrain the distance to SGR J1935+2154, to be between 1.5 and 6.5 kpc. The upper limit is consistent with some other distance indicators and suggests that the April 28 burst is closer to two orders of magnitude less energetic than the least energetic FRBs. The lack of single-pulse radio detections shows that the single pulses detected over a range of fluences are either rare, or highly clustered, or both. It may also indicate that the magnetar lies somewhere between being radio-quiet and radio-loud in terms of its ability to produce radio emission efficiently., Comment: Accepted for publication in MNRAS. Corresponding author B. W. Stappers

We introduce the Australia Telescope Compact Array (ATCA) rapid-response mode by presenting the first successful trigger on the short-duration gamma-ray burst (GRB) 181123B. Early-time radio observations of short GRBs may provide vital insights into the radio afterglow properties of Advanced LIGO- and Virgo-detected gravitational wave events, which will in turn inform follow-up strategies to search for counterparts within their large positional uncertainties. The ATCA was on target within 12.6 hr post-burst, when the source had risen above the horizon. While no radio afterglow was detected during the 8.3 hr observation, we obtained force-fitted flux densities of $7 \pm 12$ and $15 \pm 11~\mu$Jy at 5.5 and 9 GHz, respectively. Afterglow modelling of GRB 181123B showed that the addition of the ATCA force-fitted radio flux densities to the Swift X-ray Telescope detections provided more stringent constraints on the fraction of thermal energy in the electrons (log$\epsilon_e = -0.75^{+0.39}_{-0.40}$ rather than log$\epsilon_e = -1.13^{+0.82}_{-1.2}$ derived without the inclusion of the ATCA values), which is consistent with the range of typical $\epsilon_e$ derived from GRB afterglow modelling. This allowed us to predict that the forward shock may have peaked in the radio band $\sim10$ days post-burst, producing detectable radio emission $\gtrsim3-4$ days post-burst. Overall, we demonstrate the potential for extremely rapid radio follow-up of transients and the importance of triggered radio observations for constraining GRB blast wave properties, regardless of whether there is a detection, via the inclusion of force-fitted radio flux densities in afterglow modelling efforts., Comment: 15 pages, 7 figures, accepted for publication in MNRAS

The 21-cm absorption feature reported by the EDGES collaboration is several times stronger than that predicted by traditional astrophysical models. If genuine, a deeper absorption may lead to stronger fluctuations on the 21-cm signal on degree scales (up to 1~Kelvin in rms), allowing these fluctuations to be detectable in nearly 50~times shorter integration times compared to previous predictions. We commenced the "AARTFAAC Cosmic Explorer" (ACE) program, that employs the AARTFAAC wide-field imager, to measure or set limits on the power spectrum of the 21-cm fluctuations in the redshift range $z = 17.9-18.6$ ($\Delta\nu = 72.36-75.09$~MHz) corresponding to the deep part of the EDGES absorption feature. Here, we present first results from two LST bins: 23.5-23.75h and 23.5-23.75h, each with 2~h of data, recorded in `semi drift-scan' mode. We demonstrate the application of the new ACE data-processing pipeline (adapted from the LOFAR-EoR pipeline) on the AARTFAAC data. We observe that noise estimates from the channel and time-differenced Stokes~$V$ visibilities agree with each other. After 2~h of integration and subtraction of bright foregrounds, we obtain $2\sigma$ upper limits on the 21-cm power spectrum of $\Delta_{21}^2 < (8139~\textrm{mK})^2$ and $\Delta_{21}^2 < (8549~\textrm{mK})^2$ at $k = 0.144~h\,\textrm{cMpc}^{-1}$ for the two LST bins. Incoherently averaging the noise bias-corrected power spectra for the two LST bins yields an upper limit of $\Delta_{21}^2 < (7388~\textrm{mK})^2$ at $k = 0.144~h\,\textrm{cMpc}^{-1}$. These are the deepest upper limits thus far at these redshifts., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

The mergers of two neutron stars are typically accompanied by broad-band electromagnetic emission from either a relativistic jet or a kilonova. It has also been long predicted that coherent radio emission will occur during the merger phase or from a newly formed neutron star remnant, however this emission has not been seen to date. This paper presents the deepest limits for this emission from a neutron star merger folowing triggered LOFAR observations of the short gamma-ray burst (SGRB) 181123B, starting 4.4 minutes after the GRB occurred. During the X-ray plateau phase, a signature of ongoing energy injection, we detect no radio emission to a 3$\sigma$ limit of 153 mJy at 144 MHz (image integration time of 136 seconds), which is significantly fainter than the predicted emission from a standard neutron star. At a redshift of 1.8, this corresponds to a luminosity of $2.5 \times 10^{44}$ erg s$^{-1}$. Snapshot images were made of the radio observation on a range of timescales, targeting short duration radio flashes similar to fast radio bursts (FRBs). No emission was detected in the snapshot images at the location of GRB 181123B enabling constraints to be placed on the prompt coherent radio emission model and emission predicted to occur when a neutron star collapses to form a black hole. At the putative host redshift of 1.8 for GRB 181123B, the non detection of the prompt radio emission is two orders of magnitude lower than expected for magnetic reconnection models for prompt GRB emission and no magnetar emission is expected., Comment: MNRAS Accepted, 11 pages, 6 figures

We present the results of MeerKAT radio observations of eleven nearby novalike cataclysmic variables. We have detected radio emission from IM Eri, RW Sex, V3885 Sgr and V603 Aql. While RW Sex, V3885 Sgr and V603 Aql had been previously detected, this is the first reported radio detection of IM Eri. Our observations have doubled the sample of non-magnetic CVs with sensitive radio data. We observe that at our radio detection limits, a specific optical luminosity $\gtrsim 2.2\times 10^{18}$ erg/s/Hz (corresponding to $M_V \lesssim 6.0$) is required to produce a radio detection. We also observe that the X-ray and radio luminosities of our detected novalikes are on an extension of the $L_{X} \propto L_{R}^{\sim 0.7}$ power law originally proposed for non-pulsating neutron star low-mass X-ray binaries. We find no other correlations between the radio emission and emission in other wavebands or any other system parameters for the existing sample of radio-detected non-magnetic CVs. We measure in-band (0.9-1.7 GHz) radio spectral indices that are consistent with reports from earlier work. Finally, we constructed broad spectral energy distributions for our sample from published multi-wavelength data, and use them to place constraints on the mass transfer rates of these systems., Comment: 16 pages, 9 figures, accepted to MNRAS

The prompt emission in long gamma-ray bursts arises from within relativistic outflows created during the collapse of massive stars, and the mechanism by which radiation is produced may be either magnetically- or matter-dominated. In this work we suggest an observational test of a magnetically-dominated Poynting flux model that predicts both gamma-ray and low-frequency radio pulses. A common feature among early light curves of long gamma-ray bursts are X-ray flares, which have been shown to arise from sites internal to the jet. Ascribing these events to the prompt emission, we take an established Swift XRT flare sample and apply a magnetically-dominated wind model to make predictions for the timing and flux density of corresponding radio pulses in the ~100-200 MHz band observable with radio facilities such as LOFAR. We find that 44 per cent of the X-ray flares studied would have had detectable radio emission under this model, for typical sensitivities reached using LOFAR's rapid response mode and assuming negligible absorption and scattering effects in the interstellar and intergalactic medium. We estimate the rate of Swift gamma-ray bursts displaying X-ray flares with detectable radio pulses, accessible to LOFAR, of order seven per year. We determine that LOFAR triggered observations can play a key role in establishing the long debated mechanism responsible for gamma-ray burst prompt emission., Comment: 7 pages, 4 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The afterglow emission from gamma-ray bursts (GRBs) is a valuable source of information to understand the physics of these energetic explosions. The blast wave model has become the standard to describe the evolution of the afterglow emission over time and frequency. Thanks to recent developments in the theory of afterglows and numerical simulations of relativistic outflows, we are able to model the afterglow emission with realistic dynamics and radiative processes. Although the models agree with observations remarkably well, the afterglow emission still contains additional physics, instrumental systematics, and propagation effects which make the modelling of these events challenging. In this work, we present a new approach to modelling GRB afterglows, using Gaussian processes (GPs) to take into account systematics in the afterglow data. We show that, using this new approach, it is possible to obtain more reliable estimates of the explosion and microphysical parameters of GRBs. We present fit results for 5 long GRBs and find a preliminary correlation between the isotropic energetics and opening angles of GRBs, which confirms the idea of a common energy reservoir for the kinetic energy of long GRBs., Comment: Accepted to MNRAS. 10 pages, 5 figures

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations., Comment: 9 pages, 2 figures, accepted for publication in MNRAS

Black holes in binary systems execute patterns of outburst activity where two characteristic X-ray states are associated with different behaviours observed at radio wavelengths. The hard state is associated with radio emission indicative of a continuously replenished, collimated, relativistic jet, whereas the soft state is rarely associated with radio emission, and never continuously, implying the absence of a quasi-steady jet. Here we report radio observations of the black hole transient MAXI J1820$+$070 during its 2018 outburst. As the black hole transitioned from the hard to soft state we observed an isolated radio flare, which, using high angular resolution radio observations, we connect with the launch of bi-polar relativistic ejecta. This flare occurs as the radio emission of the core jet is suppressed by a factor of over 800. We monitor the evolution of the ejecta over 200 days and to a maximum separation of 10$''$, during which period it remains detectable due to in-situ particle acceleration. Using simultaneous radio observations sensitive to different angular scales we calculate an accurate estimate of energy content of the approaching ejection. This energy estimate is far larger than that derived from state transition radio flare, suggesting a systematic underestimate of jet energetics., Comment: 57 pages, 8 figures, published in Nature Astronomy

The radio:X-ray correlation that characterises accreting black holes at all mass scales - from stellar mass black holes in binary systems to super-massive black holes powering Active Galactic Nuclei - is one of the most important pieces of observational evidence supporting the existence of a connection between the accretion process and the generation of collimated outflows - or jets - in accreting systems. Although recent studies suggest that the correlation extends down to low luminosities, only a handful of stellar mass black holes have been clearly detected, and in general only upper limits (especially at radio wavelengths) can be obtained during quiescence. We recently obtained detections of the black hole X-ray binary GX 339--4 in quiescence using the MeerKAT radio telescope and Swift X-ray Telescope instrument onboard the Neil Gehrels Swift Observatory, probing the lower end of the radio:X-ray correlation. We present the properties of accretion and of the connected generation of jets in the poorly studied low-accretion rate regime for this canonical black hole XRB system., Comment: 5 pages, 2 figure, accepted to MNRAS

We report the discovery of the first transient with MeerKAT, MKT J170456.2$-$482100, discovered in ThunderKAT images of the low mass X-ray binary GX339$-$4. MKT J170456.2$-$482100 is variable in the radio, reaching a maximum flux density of $0.71\pm0.11\,\mathrm{mJy}$ on 2019 Oct 12, and is undetected in 15 out of 48 ThunderKAT epochs. MKT J170456.2$-$482100 is coincident with the chromospherically active K-type sub-giant TYC 8332-2529-1, and $\sim18\,\mathrm{years}$ of archival optical photometry of the star shows that it varies with a period of $21.25\pm0.04\,\mathrm{days}$. The shape and phase of the optical light curve changes over time, and we detect both X-ray and UV emission at the position of MKT J170456.2$-$482100, which may indicate that TYC 8332-2529-1 has large star spots. Spectroscopic analysis shows that TYC 8332-2529-1 is in a binary, and has a line-of-sight radial velocity amplitude of $43\,\mathrm{km\,s^{-1}}$. We also observe a spectral feature in anti-phase with the K-type sub-giant, with a line-of-sight radial velocity amplitude of $\sim12\pm10\,\mathrm{km\,s^{-1}}$, whose origins cannot currently be explained. Further observations and investigation are required to determine the nature of the MKT J170456.2$-$482100 system., Comment: 16 pages, 12 figures, 7 tables

The nature of the central engines of gamma-ray bursts (GRBs) and the composition of their relativistic jets are still under debate. If the jets are Poynting flux dominated rather than baryon dominated, a coherent radio flare from magnetic re-connection events might be expected with the prompt gamma-ray emission. There are two competing models for the central engines of GRBs; a black hole or a newly formed milli-second magnetar. If the central engine is a magnetar it is predicted to produce coherent radio emission as persistent or flaring activity. In this paper, we present the deepest limits to date for this emission following LOFAR rapid response observations of GRB 180706A. No emission is detected to a 3$\sigma$ limit of 1.7 mJy beam$^{-1}$ at 144 MHz in a two-hour LOFAR observation starting 4.5 minutes after the gamma-ray trigger. A forced source extraction at the position of GRB 180706A provides a marginally positive (1 sigma) peak flux density of $1.1 \pm 0.9$ mJy. The data were time-sliced into different sets of snapshot durations to search for FRB like emission. No short duration emission was detected at the location of the GRB. We compare these results to theoretical models and discuss the implications of a non-detection., Comment: MNRAS accepted, 11 pages, 5 figures

We report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date., Comment: 17 pages, 6 figures, Version accepted by ApJ

Long gamma-ray bursts give us the chance to study both their extreme physics and the star-forming galaxies in which they form. GRB 100418A, at a z = 0.6239, had a bright optical and radio afterglow, and a luminous star-forming host galaxy. This allowed us to study the radiation of the explosion as well as the interstellar medium of the host both in absorption and emission. We collected photometric data from radio to X-ray wavelengths to study the evolution of the afterglow and the contribution of a possible supernova and three X-shooter spectra obtained during the first 60 hr. The light curve shows a very fast optical rebrightening, with an amplitude of 3 magnitudes, starting 2.4 hr after the GRB onset. This cannot be explained by a standard external shock model and requires other contributions, such as late central-engine activity. Two weeks after the burst we detect an excess in the light curve consistent with a SN with peak absolute magnitude M_V = -18.5 mag, among the faintest GRB-SNe detected to date. The host galaxy shows two components in emission, with velocities differing by 130 km s^-1, but otherwise having similar properties. While some absorption and emission components coincide, the absorbing gas spans much higher velocities, indicating the presence of gas beyond the star-forming regions. The host has a star-formation rate of 12.2 M_sol yr^-1, a metallicity of 12 + log(O/H) = 8.55 and a mass of 1.6x10^9 M_sol. GRB 100418A is a member of a class of afterglow light curves which show a steep rebrightening in the optical during the first day, which cannot be explained by traditional models. Its very faint associated SN shows that GRB-SNe can have a larger dispersion in luminosities than previously seen. Furthermore, we have obtained a complete view of the host of GRB 100418A owing to its spectrum, which contains a remarkable number of both emission and absorption lines., Comment: 23 pages, 14 figures, accepted for publication in A&A

The part played by stars in the ionization of the intergalactic medium remains an open question. A key issue is the proportion of the stellar ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst afterglows can be used to determine the neutral hydrogen column-density in their host galaxies and hence the opacity to extreme ultra-violet radiation along the lines-of-sight to the bursts. Thus, making the reasonable assumption that long-duration GRB locations are representative of the sites of massive stars that dominate EUV production, one can calculate an average escape fraction of ionizing radiation in a way that is independent of galaxy size, luminosity or underlying spectrum. Here we present a sample of NH measures for 138 GRBs in the range 1.6~0.005, with a 98% confidence upper limit of ~0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget of the IGM at z<5, where the bulk of the bursts lie. At higher redshifts, z>5, firm conclusions are limited by the small size of the GRB sample, but any decline in average HI column-density seems to be modest. We also find no indication of a significant correlation of NH with galaxy UV luminosity or host stellar mass, for the subset of events for which these are available. We discuss in some detail a number of selection effects and potential biases. Drawing on a range of evidence we argue that such effects, while not negligible, are unlikely to produce systematic errors of more than a factor ~2, and so would not affect the primary conclusions. Given that many GRB hosts are low metallicity, high specific star-formation rate, dwarf galaxies, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the universe., Comment: 31 pages

With their wide fields of view and often relatively long coverage of any position in the sky in imaging survey mode, modern radio telescopes provide a data stream that is naturally suited to searching for rare transients. However, Radio Frequency Interference (RFI) can show up in the data stream in similar ways to such transients, and thus the normal pre-treatment of filtering RFI (flagging) may also remove astrophysical transients from the data stream before imaging. In this paper we investigate how standard flagging affects the detectability of such transients by examining the case of transient detection in an observing mode used for Low Frequency Array (LOFAR; \citep{LOFAR}) surveys. We quantify the fluence range of transients that would be detected, and the reduction of their SNR due to partial flagging. We find that transients with a duration close to the integration sampling time, as well as bright transients with durations on the order of tens of seconds, are completely flagged. For longer transients on the order of several tens of seconds to minutes, the flagging effects are not as severe, although part of the signal is lost. For these transients, we present a modified flagging strategy which mitigates the effect of flagging on transient signals. We also present a script which uses the differences between the two strategies, and known differences between transient RFI and astrophysical transients, to notify the observer when a potential transient is in the data stream., Comment: 31 pages, 6 figures, accepted for publication in Astronomy and Computing

PySE is a Python software package for finding and measuring sources in radio telescope images. The software was designed to detect sources in the LOFAR telescope images, but can be used with images from other radio telescopes as well. We introduce the LOFAR Telescope, the context within which PySE was developed, the design of PySE, and describe how it is used. Detailed experiments on the validation and testing of PySE are then presented, along with results of performance testing. We discuss some of the current issues with the algorithms implemented in PySE and their inter- action with LOFAR images, concluding with the current status of PySE and its future development., Comment: 26 pages, 6 figures

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes., Comment: 19 pages, 8 figures, 4 tables, accepted for publication in MNRAS

The binary neutron star merger GW170817 was the first multi-messenger event observed in both gravitational and electromagnetic waves. The electromagnetic signal began approximately 2 seconds post-merger with a weak, short burst of gamma-rays, which was followed over the next hours and days by the ultraviolet, optical and near-infrared emission from a radioactively- powered kilonova. Later, non-thermal rising X-ray and radio emission was observed. The low luminosity of the gamma-rays and the rising non-thermal flux from the source at late times could indicate that we are outside the opening angle of the beamed relativistic jet. Alternatively, the emission could be arising from a cocoon of material formed from the interaction between a jet and the merger ejecta. Here we present late-time optical detections and deep near-infrared limits on the emission from GW170817 at 110 days post-merger. Our new observations are at odds with expectations of late-time emission from kilonova models, being too bright and blue. Instead, the emission arises from the interaction between the relativistic ejecta of GW170817 and the interstellar medium. We show that this emission matches the expectations of a Gaussian structured relativistic jet, which would have launched a high luminosity short GRB to an aligned observer. However, other jet structure or cocoon models can also match current data - the future evolution of the afterglow will directly distinguish the origin of the emission., Comment: Includes MCMC fitting

The merger of two dense stellar remnants including at least one neutron star (NS) is predicted to produce gravitational waves (GWs) and short duration gamma ray bursts (GRBs). In the process, neutron-rich material is ejected from the system and heavy elements are synthesized by r-process nucleosynthesis. The radioactive decay of these heavy elements produces additional transient radiation termed "kilonova" or "macronova". We report the detection of linear optical polarization P = (0.50 +/- 0.07)% at 1.46 days after detection of the GWs from GW170817, a double neutron star merger associated with an optical macronova counterpart and a short GRB. The optical emission from a macronova is expected to be characterized by a blue, rapidly decaying, component and a red, more slowly evolving, component due to material rich of heavy elements, the lanthanides. The polarization measurement was made when the macronova was still in its blue phase, during which there is an important contribution from a lanthanide-free outflow. The low degree of polarization is consistent with intrinsically unpolarized emission scattered by Galactic dust, suggesting a symmetric geometry of the emitting region and low inclination of the merger system. Stringent upper limits to the polarization degree from 2.45 - 9.48 days post-burst are consistent with the lanthanides-rich macronova interpretation., Comment: 18 pages, 1 figure, 2 tables, Nature Astronomy, in press

We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient light is consistent with predictions for the behaviour of a "kilonova/macronova", powered by the radioactive decay of massive neutron-rich nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In particular, evidence for this scenario is found from broad features seen in Hubble Space Telescope infrared spectroscopy, similar to those predicted for lanthanide dominated ejecta, and the much slower evolution in the near-infrared Ks-band compared to the optical. This indicates that the late-time light is dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to the 3rd r-process peak (atomic masses A~195). This discovery confirms that neutron-star mergers produce kilo-/macronovae and that they are at least a major - if not the dominant - site of rapid neutron capture nucleosynthesis in the universe.

We present Hubble Space Telescope and Chandra imaging, combined with Very Large Telescope MUSE integral field spectroscopy of the counterpart and host galaxy of the first binary neutron star merger detected via gravitational wave emission by LIGO & Virgo, GW170817. The host galaxy, NGC 4993, is an S0 galaxy at z=0.009783. There is evidence for large, face-on spiral shells in continuum imaging, and edge-on spiral features visible in nebular emission lines. This suggests that NGC 4993 has undergone a relatively recent (<1 Gyr) ``dry'' merger. This merger may provide the fuel for a weak active nucleus seen in Chandra imaging. At the location of the counterpart, HST imaging implies there is no globular or young stellar cluster, with a limit of a few thousand solar masses for any young system. The population in the vicinity is predominantly old with <1% of any light arising from a population with ages <500 Myr. Both the host galaxy properties and those of the transient location are consistent with the distributions seen for short-duration gamma-ray bursts, although the source position lies well within the effective radius (r_e ~ 3 kpc), providing an r_e-normalized offset that is closer than ~90% of short GRBs. For the long delay time implied by the stellar population, this suggests that the kick velocity was significantly less than the galaxy escape velocity. We do not see any narrow host galaxy interstellar medium features within the counterpart spectrum, implying low extinction, and that the binary may lie in front of the bulk of the host galaxy., Comment: ApJL in press, 13 pages

We present the Arcminute Microkelvin Imager (AMI) Large Array catalogue of 139 gamma-ray bursts (GRBs). AMI observes at a central frequency of 15.7 GHz and is equipped with a fully automated rapid-response mode, which enables the telescope to respond to high-energy transients detected by Swift. On receiving a transient alert, AMI can be on-target within two minutes, scheduling later start times if the source is below the horizon. Further AMI observations are manually scheduled for several days following the trigger. The AMI GRB programme probes the early-time (< 1 day) radio properties of GRBs, and has obtained some of the earliest radio detections (GRB 130427A at 0.36 and GRB 130907A at 0.51 days post-burst). As all Swift GRBs visible to AMI are observed, this catalogue provides the first representative sample of GRB radio properties, unbiased by multi-wavelength selection criteria. We report the detection of six GRB radio afterglows that were not previously detected by other radio telescopes, increasing the rate of radio detections by 50% over an 18-month period. The AMI catalogue implies a Swift GRB radio detection rate of >15%, down to ~0.2 mJy/beam. However, scaling this by the fraction of GRBs AMI would have detected in the Chandra & Frail sample (all radio-observed GRBs between 1997 - 2011), it is possible ~44 - 56% of Swift GRBs are radio-bright, down to ~0.1 - 0.15 mJy/beam. This increase from the Chandra & Frail rate (~30%) is likely due to the AMI rapid-response mode, which allows observations to begin while the reverse-shock is contributing to the radio afterglow., Comment: 47 pages (20 of which are catalogue tables), 10 figures, 4 tables, accepted for publication in MNRAS