Search

Radio detections of stellar systems provide a window onto stellar magnetic activity and the space weather conditions of extrasolar planets, information that is difficult to attain at other wavelengths. There have been recent advances observing auroral emissions from radio-bright low-mass stars and exoplanets largely due to the maturation of low-frequency radio instruments and the plethora of wide-field radio surveys. To guide us in placing these recent results in context, we introduce the foremost local analogues for the field: Solar bursts and the aurorae found on Jupiter. We detail how radio bursts associated with stellar flares are foundational to the study of stellar coronae, and time-resolved radio dynamic spectra offers one of the best prospects of detecting and characterising coronal mass ejections from other stars. We highlight the prospects of directly detecting coherent radio emission from exoplanetary magnetospheres, and early tentative results. We bridge this discussion to the field of brown dwarf radio emission, in which their larger and stronger magnetospheres are amenable to detailed study with current instruments. Bright, coherent radio emission is also predicted from magnetic interactions between stars and close-in planets. We discuss the underlying physics of these interactions and implications of recent provisional detections for exoplanet characterisation. We conclude with an overview of outstanding questions in theory of stellar, star-planet interaction, and exoplanet radio emission, and the prospects of future facilities in answering them., Comment: Accepted to Nature Astronomy. The manuscript is designed to be a primer for new doctoral students and scholars to the field of radio stars and exoplanets. 36 pages, 3 figures

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.

Recent sensitive wide-field radio surveys, such as the LOFAR Two Meter Sky Survey (LoTSS), the LOFAR LBA Sky Survey (LoLSS), and the Very Large Array Sky Survey (VLASS), enable the selection of statistically large samples of peaked-spectrum (PS) sources. PS sources are radio sources that have a peak in their radio continuum spectrum and are observed to be compact. They are often considered to be the precursors to large radio galaxies. We present a sample of 8,032 gigahertz-peaked spectrum (GPS) sources with spectral turnovers near 1400 MHz, and a sample of 506 megahertz-peaked spectrum (MPS) sources with turnovers near 144 MHz. Our GPS sample is over five times larger than any previously known sample of PS sources. These large sample sizes allow us to make a robust comparison between GPS sources and MPS sources, such that we can investigate the differences between these types of sources, and study their lifetimes. The shape of the source counts of both samples match that of the general radio-loud active galactic nuclei (AGN) samples, scaled down by a factor 44 $\pm$ 2 for the MPS sample, and a factor 28 $\pm$ 1 for the GPS sample. Assuming no cosmological evolution, these offsets imply that both MPS and GPS sources have shorter duration than general radio-loud AGN, with MPS sources having an $\approx$1.6 times shorter lifespan than GPS sources. The shorter duration of MPS sources relative to GPS sources can be explained by the transition between GPS and MPS sources coinciding with the jet breakout phase of PS sources, such that GPS sources traverse through the surrounding medium at a lower speed than MPS sources. Such evolution has been observed in simulations of PS source evolution., Comment: Accepted for publication in Astronomy & Astrophysics. 16 pages, 9 figures

Merging galaxy clusters often host spectacular diffuse radio synchrotron sources. These sources can be explained by a non-thermal pool of relativistic electrons accelerated by shocks and turbulence in the intracluster medium. The origin of the pool and details of the cosmic ray transport and acceleration mechanisms in clusters are still open questions. Due to the often extremely steep spectral indices of diffuse radio emission, it is best studied at low frequencies. However, the lowest frequency window available to ground-based telescopes (10-30 MHz) has remained largely unexplored, as radio frequency interference and calibration problems related to the ionosphere become severe. Here, we present LOFAR observations from 16 to 168 MHz targeting the famous cluster Abell 2256. In the deepest-ever images at decametre wavelengths, we detect and resolve the radio halo, radio shock and various steep spectrum sources. We measure standard single power-law behaviour for the radio halo and radio shock spectra and find significant spectral index and curvature fluctuations across the radio halo, indicating an inhomogeneous emitting volume. In contrast to the straight power-law spectra of the large-scale diffuse sources, the various AGN-related sources often show extreme steepening towards higher frequencies and flattening towards low frequencies. We also discover a new fossil plasma source with a steep spectrum between 23 and 144 MHz, with $\alpha=-1.9\pm 0.1$. Finally, by comparing radio and gamma-ray observations, we rule out purely hadronic models for the radio halo origin in Abell 2256, unless the magnetic field strength in the cluster is exceptionally high, which is unsupportable by energetic arguments and inconsistent with the knowledge of other cluster magnetic fields., Comment: Accepted for publication in Astronomy & Astrophysics

The largely unexplored decameter radio band (10-30 MHz) provides a unique window for studying a range of astronomical topics, such as auroral emission from exoplanets, inefficient cosmic ray acceleration mechanisms, fossil radio plasma, and free-free absorption. The scarcity of low-frequency studies is mainly due to the severe perturbing effects of the ionosphere. Here we present a calibration strategy that can correct for the ionosphere in the decameter band. We apply this to an observation from the Low Frequency Array (LOFAR) between 16 to 30 MHz . The resulting image covers 330 square degrees of sky at a resolution of 45", reaching a sensitivity of 12 mJy/beam. Residual ionospheric effects cause additional blurring ranging between 60 to 100". This represents an order of magnitude improvement in terms of sensitivity and resolution compared to previous decameter band observations. In the region we surveyed, we have identified four fossil plasma sources. These rare sources are believed to contain old, possibly re-energised, radio plasma originating from previous outbursts of active galactic nuclei. At least three of them are situated near the center of low-mass galaxy clusters. Notably, two of these sources display the steepest radio spectral index among all the sources detected at 23 MHz. This indicates that fossil plasma sources constitute the primary population of steep-spectrum sources at these frequencies, emphasising the large discovery potential of ground-based decameter observations., Comment: 19 pages, 9 figures, 1 table. Submitted version, full version is published by Nature Astronomy

We characterize the LHS 1610 system, a nearby ($d=9.7$ pc) M5 dwarf hosting a brown dwarf in a $10.6$ day, eccentric ($e \sim 0.37$) orbit. A joint fit of the available Gaia two-body solution, discovery radial velocities (RVs) from TRES, and new RVs obtained with the Habitable-zone Planet Finder, yields an orbital inclination of $117.2\pm0.9^\circ$ and a mass constraint of $50.9\pm0.9$ M$_J$. This gives LHS 1610 b the second most precise mass of brown dwarfs orbiting M stars within 25pc. We highlight a discrepancy between the Gaia two-body solution eccentricity ($e=0.52 \pm 0.03$) and that from the RVs ($e=0.3702\pm0.0003$), which requires the astrometric time-series release (Gaia DR4) for further diagnostics. With a flare rate of $0.28\pm 0.07$ flares/day from TESS photometry, and a rotation period of $84 \pm 8$ days, LHS 1610 joins other mid M stars -- including Proxima Centauri and YZ Ceti -- as nearby mid M dwarfs with flare rates on the higher end for their long rotation periods. These stars are promising candidates for searching for sub-Alfv\'enic star-companion interactions, raising the question whether LHS 1610 b could be driving the flares on its host star. However, the available TESS photometry is insufficient to confirm or rule out any orbital phase-dependence of the flares. We show that the LHS 1610 system, as a nearby mid M star with a large, short-period companion, is a promising target to look for evidence of star-companion interactions or aural emission from the brown dwarf at radio wavelengths., Comment: 24 pages, 7 figures, 3 tables. Submitted to AAS Journals on Oct 11, 2023

We present the low-frequency radio spectra of 9 high-redshift quasars at $5.6 \leq z \leq 6.6$ using the Giant Metre Radio Telescope band-3, -4, and -5 observations ($\sim$300-1200 MHz), archival Low Frequency Array (LOFAR; 144 MHz), and Very Large Array (VLA; 1.4 and 3 GHz) data. Five of the quasars in our sample have been discovered recently, representing some of the highest redshift radio bright quasars known at low-frequencies. We model their radio spectra to study their radio emission mechanism and age of the radio jets by constraining the spectral turnover caused by synchrotron self-absorption (SSA) or free-free absorption (FFA). Besides J0309+2717, a blazar at $z=6.1$, our quasars show no sign of a spectral flattening between 144 MHz and a few GHz, indicating there is no strong SSA or FFA absorption in the observed frequency range. However, we find a wide range of spectral indices between $-1.6$ and $0.05$, including the discovery of 3 potential ultra-steep spectrum quasars. Using further archival VLBA data, we confirm that the radio SED of the blazar J0309+2717 likely turns over at a rest-frame frequency of 0.6-2.3 GHz (90-330 MHz observed frame), with a high-frequency break indicative of radiative ageing of the electron population in the radio lobes. Ultra-low frequency data below 50 MHz are necessary to constrain the absorption mechanism for J0309+2717 and the turnover frequencies for the other high-$z$ quasars in our sample. A relation between linear radio jet size and turnover frequency has been established at low redshifts. If this relation were to hold at high redshifts, the limits on the turnover frequency of our sample suggest the radio jet sizes must be more extended than the typical sizes observed in other radio-bright quasars at similar redshift. To confirm this deep radio follow-up observations with high spatial resolution are required., Comment: 15 pages, 7 figures, accepted for publication in A&A, minor correction Tab. 1

We present the serendipitous discovery of 8 distant ($>$ 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M dwarf nature of our sources has been confirmed with optical spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5-1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities of the parent catalogues, we suggest that it is statistically probable the M dwarfs are associated with the radio emission. However, it remains plausible that for some of the sources the radio emission originates from an optically faint and red galaxy hiding behind the M dwarf. The isotropic radio luminosities ($\sim10^{17-18}$ erg s$^{-1}$ Hz$^{-1}$) of the M dwarfs suggest that if the association is real, the radio emission is likely driven by a coherent emission process produced via plasma or electron-cyclotron maser instability processes, which is potentially caused by binary interaction. Long term monitoring in the radio and high-resolution radio follow-up observations are necessary to search for any variability and pinpoint the radio emission to determine whether our tentative conclusion that these ultracool dwarfs are radio emitting is correct. If the low-frequency radio emission is conclusively associated with the M dwarfs, this would reveal a new population of optically faint and distant ($>$ 50 pc) radio emitting M dwarfs., Comment: 10 pages, 5 figures, accepted for publication in A&A

The Low Frequency Array (LOFAR) is the only existing radio interferometer able to observe at ultra-low frequencies (<100 MHz) with high resolution (<15") and high sensitivity (<1 mJy/beam). To exploit these capabilities, the LOFAR Surveys Key Science Project is using the LOFAR Low Band Antenna (LBA) to carry out a sensitive wide-area survey at 41-66 MHz named the LOFAR LBA Sky Survey (LoLSS). LoLSS is covering the whole northern sky above declination 24 deg with a resolution of 15" and a sensitivity of 1-2 mJy/beam (1 sigma) depending on declination, field properties, and observing conditions. Here we present the first data release. An automated pipeline was used to reduce the 95 fields included in this data release. The data reduction procedures developed for this project have general application and are currently being used to process LOFAR LBA interferometric observations. Compared to the preliminary release, direction-dependent errors have been corrected for during the calibration process. This results in a typical sensitivity of 1.55 mJy/beam at the target resolution of 15". The first data release of the LOFAR LBA Sky Survey covers 650 sqdeg in the HETDEX spring field. The resultant data products released to the community include mosaic images (I and V Stokes) of the region, and a catalogue of 42463 detected sources and related Gaussian components used to describe sources' morphologies. Separate catalogues for 6 in-band frequencies are also released. The first data release of LoLSS shows that, despite the influences of the ionosphere, LOFAR can conduct large-scale surveys in the frequency window 42-66 MHz with unprecedentedly high sensitivity and resolution. The data can be used to derive unique information on the low-frequency spectral properties of many thousands of sources with a wide range of applications in extragalactic and galactic astronomy., Comment: 20 pages, 22 figures, images and catalogues available at https://www.lofar-surveys.org/lolss.html

Brown dwarfs display Jupiter-like auroral phenomena such as magnetospheric H$\alpha$ emission and coherent radio emission. Coherent radio emission is a probe of magnetospheric acceleration mechanisms and provides a direct measurement of the magnetic field strength at the emitter's location, both of which are difficult to access by other means. Observations of the coldest brown dwarfs (spectral types T and Y) are particularly interesting as their magnetospheric phenomena may be very similar to those in gas-giant exoplanets. Here we present 144 MHz radio and infrared adaptive optics observations of the brown dwarf WISEP J101905.63+652954.2 made using the LOFAR and Keck telescopes respectively. The radio data shows pulsed highly circularly polarised emission which yields a rotation rate of $0.32\pm0.03$ hr$^{-1}$. The infrared imaging reveals the source to be a binary with a projected separation of $423.0\pm1.6$ mas between components of spectral type T5.$5\pm0.5$ and T7.$0\pm0.5$. With a simple "toy model" we show that the radio emission can in principle be powered by the interaction between the two dwarfs with a mass-loss rates of at least $25$ times the Jovian value. WISEP J101905.63+652954.2 is interesting because it is the first pulsed methane dwarf detected in a low radio-frequency search. Unlike previous gigahertz-frequency searches that were only sensitive to objects with kiloGauss fields, our low-frequency search is sensitive to surface magnetic fields of $\approx 50$ Gauss and above which might reveal the coldest radio-loud objects down to planetary mass-scales., Comment: Accepted for publication in A&A

Spectral variability offers a new technique to identify small scale structures from scintillation, as well as determining the absorption mechanism for peaked-spectrum (PS) radio sources. In this paper, we present very long baseline interferometry (VLBI) imaging using the Long Baseline Array (LBA) of two PS sources, MRC0225-065 and PMNJ0322-4820, identified as spectrally variable from observations with the Murchison Widefield Array (MWA). We compare expected milliarcsecond structures based on the detected spectral variability with direct LBA imaging. We find MRC0225-065 is resolved into three components, a bright core and two fainter lobes, roughly 430pc projected separation. A comprehensive analysis of the magnetic field, host galaxy properties, and spectral analysis implies that MRC0225-065 is a young radio source with recent jet activity over the last 10^2-10^3years. We find PMNJ0322-4820 is unresolved on milliarcsecond scales. We conclude PMNJ0322-4820 is a blazar with flaring activity detected in 2014 with the MWA. We use spectral variability to predict morphology and find these predictions consistent with the structures revealed by our LBA images., Comment: Accepted for publication in PASA. 11 pages, 4 figures

We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $\mu$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polarisation is measured to be $\approx$0.06%, and our survey is complete at flux densities $\geq1$ mJy. A detection is considered reliable when the circularly-polarised fraction exceeds 1%. We find the population of circularly-polarised sources is composed of four distinct classes: stellar systems, pulsars, active galactic nuclei, and sources unidentified in the literature. The stellar systems can be further separated into chromospherically-active stars, M dwarfs, and brown dwarfs. Based on the circularly-polarised fraction and lack of an optical counterpart, we show it is possible to infer whether the unidentified sources are likely unknown pulsars or brown dwarfs. By the completion of this survey of the northern sky, we expect to detect 300$\pm$100 circularly-polarised sources., Comment: Accepted for publication in A&A, 15 pages, 8 figures, 1 table. The catalogue will be publicly available at http://lofar-surveys.org/ and via Vizier shortly

We present a sample of 373 peaked-spectrum (PS) sources with spectral peaks around 150MHz, selected using a subset of two LOFAR all-sky surveys, the LOFAR Two Meter Sky Survey and the LOFAR LBA Sky Survey. These surveys are the most sensitive low-frequency widefield surveys to date, allowing us to select low-luminosity PS sources. Our sample increases the number of known PS sources in our survey area by a factor 50. The 5GHz luminosity distribution of our PS sample shows we sample the lowest luminosity PS sources to-date by nearly an order of magnitude. Since high-frequency PS sources and compact steep-spectrum sources are hypothesised to be the precursors to large radio galaxies, we investigate whether this is also the case for our sample of low-frequency PS sources. Using optical line emission criteria, we find that our PS sources are predominately high-excitation radio galaxies instead of low-excitation radio galaxies, corresponding to a quickly evolving population. We compute the radio source counts of our PS sample, and find they are scaled down by a factor of $\sim$40 compared to a general sample of radio-loud active galactic nuclei (AGN). This implies that the lifetimes of PS sources are 40 times shorter than large scale radio galaxies, if their luminosity functions are identical. To investigate this, we compute the first radio luminosity function for a homogeneously-selected PS sample. We find that for 144MHz luminosities $\gtrsim 10^{25}$W Hz$^{-1}$, the PS luminosity function has the same shape as an unresolved radio-loud AGN population but shifted down by a factor of $\sim$10. We interpret this as strong evidence that these high-luminosity PS sources evolve into large-scale radio-loud AGN. For local, low-luminosity PS sources, there is a surplus of PS sources, which we hypothesise to be the addition of frustrated PS sources that do not evolve into large-scale AGN., Comment: 21 pages, accepted in Astronomy & Astrophysics

Characterising spectral variability of radio sources is a technique that offers the ability to determine the astrophysics of the intervening media, source structure, emission and absorption processes. We present broadband (0.072--10 GHz) spectral variability of 15 peaked-spectrum (PS) sources with the Australia Telescope Compact Array (ATCA) and the Murchison Widefield Array (MWA). These 15 PS sources were observed quasi-contemporaneously with ATCA and the MWA four to six times during 2020 with approximately a monthly cadence. Variability was not detected at 1--10GHz frequencies but 13 of the 15 targets show significant variability with the MWA at megahertz frequencies. We conclude the majority of variability seen at megahertz frequencies is due to refractive interstellar scintillation of a compact component ~25 mas across. We also identify four PS sources that show a change in their spectral shape at megahertz frequencies. Three of these sources are consistent with a variable optical depth from an inhomogeneous free-free absorbing cloud around the source. One PS source with a variable spectral shape at megahertz frequencies is consistent with an ejection travelling along the jet. We present spectral variability as a method for determining the physical origins of observed variability and for providing further evidence to support absorption models for PS sources where spectral modelling alone is insufficient., Comment: Accepted for publication in MNRAS. 17 pages, 10 figures

Discovering radio pulsars, particularly millisecond pulsars (MSPs), is important for a range of astrophysical applications, such as testing theories of gravity or probing the magneto-ionic interstellar medium. We aim to discover pulsars that may have been missed in previous pulsar searches by leveraging known pulsar observables (primarily polarisation) in the sensitive, low-frequency radio images from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS), and have commenced the Targeted search, using LoTSS images, for polarised pulsars (TULIPP) survey. For this survey, we identified linearly and circularly polarised point sources with flux densities brighter than 2 mJy in LoTSS images at a centre frequency of 144 MHz with a 48 MHz bandwidth. Over 40 known pulsars, half of which are MSPs, were detected as polarised sources in the LoTSS images and excluded from the survey. We have obtained beam-formed LOFAR observations of 30 candidates, which were searched for pulsations using coherent de-dispersion. Here, we present the results of the first year of the TULIPP survey. We discovered two pulsars, PSRs J1049+5822 and J1602+3901, with rotational periods of P=0.73 s and 3.7 ms, respectively. We also detected a further five known pulsars (two slowly-rotating pulsars and three MSPs) for which accurate sky positions were not available to allow a unique cross-match with LoTSS sources. This targeted survey presents a relatively efficient method by which pulsars, particularly MSPs, may be discovered using the flexible observing modes of sensitive radio telescopes such as the Square Kilometre Array and its pathfinders/precursors, particularly since wide-area all-sky surveys using coherent de-dispersion are currently computationally infeasible., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in A&A

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$\mu$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$\mu$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data., Comment: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

The empirical relationship between the non-thermal 5GHz radio luminosity and the soft X-ray luminosity of active stellar coronae, canonically called the G\"udel-Benz relationship (G\"udel & Benz 1993), has been a cornerstone of stellar radio astronomy as it explicitly ties the radio emission to the coronal heating mechanisms. The relationship extends from microflares on the Sun to the coronae of the most active stars suggesting that active coronae are heated by a flare-like process (Benz & G\"udel 1994). The relationship is thought to originate from a consistent partition of the available flare energy into relativistic charges, that emit in the radio-band via the incoherent gyrosynchrotron mechanism, and heating of the bulk coronal plasma, that emits in the X-ray band via the Bremsstrahlung mechanism. Consequently, coherent emission from stellar and sub-stellar objects is not expected to adhere to this empirical relationship, as is observed in ultracool dwarf stars and brown dwarfs. Here we report a population of radio-detected chromospherically active stars that surprisingly follows the G\"udel-Benz relationship despite their radio emission being classified as coherent emission by virtue of its high circularly polarised fraction and high brightness temperature. Our results prompt a re-examination of the physics behind the G\"udel-Benz relationship, its implication for the mechanism of coronal heating and particle acceleration in active stars and the phenomenological connection between solar and stellar flares., Comment: Under review, ApJL

In this work, we study the optical properties of compact radio sources selected from the literature in order to determine the impact of the radio-jet in their circumnuclear environment. Our sample includes 58 Compact Steep Spectrum (CSS) and GigaHertz Peaked Spectrum (GPS) and 14 Megahertz-Peaked spectrum (MPS) radio sources located at $z\leq 1$. The radio luminosity ($L_R$) of the sample varies between Log\,L$_R\sim$ 23.2 and 27.7 W\,Hz$^{-1}$. We obtained optical spectra for all sources from SDSS-DR12 and performed a stellar population synthesis using the {\sc starlight} code. We derived stellar masses (M$_\star$), ages $\langle t_\star \rangle$, star formation rates (SFR), metallicities $\langle Z_\star \rangle$ and internal reddening A$_V$ for all young AGNs of our sample. A visual inspection of the SDSS images was made to assign a morphological class for each source. Our results indicate that the sample is dominated by intermediate to old stellar populations and there is no strong correlation between optical and radio properties of these sources. Also, we found that young AGNs can be hosted by elliptical, spiral and interacting galaxies, confirming recent findings. When comparing the optical properties of CSS/GPS and MPS sources, we do not find any significant difference. Finally, the Mid-Infrared WISE colours analysis suggest that compact radio sources defined as powerful AGNs are, in general, gas-rich systems., Comment: 44 pages, 12 figures, 5 tables, accepted for publication in MNRAS

Apep is the brightest and most luminous non-thermal colliding-wind binary by over an order of magnitude. It has been suggested from infrared observations that one of the Wolf-Rayet stars in Apep is launching an anisotropic wind. Here we present radio observations of Apep from 0.2 to 20 GHz taken over 33 years. The spectrum reveals an extremely steep turnover in the flux density at low frequencies, where the flux density decreases by two orders of magnitude over only 325 MHz of bandwidth. This exponential decline is best described by free-free absorption, with a turnover frequency at 0.54 $\pm$ 0.01 GHz. Above the turnover, the spectrum is well described by a power-law and a high-frequency cut-off likely caused by inverse-Compton cooling. The lightcurve of Apep shows significant variation over the observing period, with Apep brightening by over 50 mJy in a span of 25 years at 1.4 GHz. Models that assume spherical winds do not replicate all of the structure evident in the radio lightcurve. We derived a model that allows one of the winds in the system to be anisotropic. This anisotropic model recovers most of the structure of the lightcurve and is a significantly better statistical fit to the data than the spherical wind model. We suggest such a result is independent support that one of the Wolf-Rayet stars in Apep is launching an anisotropic wind. If the anisotropic wind model is correct, we predict a ~25% decrease of the 1.4 GHz flux density of Apep over the next five years., Comment: 14 pages, 6 figures, accepted for publication in MNRAS

Coherent low-frequency ($\lesssim 200$ MHz) radio emission from stars encodes the conditions of the outer corona, mass-ejection events, and space weather. Previous low-frequency searches for radio emitting stellar systems have lacked the sensitivity to detect the general population, instead largely focusing on targeted studies of anomalously active stars. Here we present 19 detections of coherent radio emission associated with known M~dwarfs from a blind flux-limited low-frequency survey. Our detections show that coherent radio emission is ubiquitous across the M~dwarf main sequence, and that the radio luminosity is independent of known coronal and chromospheric activity indicators. While plasma emission can generate the low-frequency emission from the most chromospherically active stars of our sample, the origin of the radio emission from the most quiescent sources is yet to be ascertained. Large-scale analogues of the magnetospheric processes seen in gas-giant planets likely drive the radio emission associated with these quiescent stars. The slowest-rotating stars of this sample are candidate systems to search for star-planet interaction signatures., Comment: Accepted to Nature Astronomy; 31 pages, 9 figures, and 2 tables

[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

Relatively little information is available about the Universe at ultra-low radio frequencies, i.e. below 50 MHz (ULF), although the ULF spectral window contains a wealth of unique diagnostics for studying galactic and extragalactic phenomena. Sub-arcsecond resolution imaging at these frequencies is extremely difficult, due to the long baselines (>1000 km) required and large ionospheric perturbations. We have conducted a pilot project to investigate the ULF performance and potential of the International LOFAR Telescope (ILT), a trans-European interferometric array with baselines up to ~2000 km and observing frequencies down to 10 MHz. We have successfully produced images with sub-arcsecond resolution for 6 radio sources at frequencies down to 30 MHz. This is more than an order of magnitude better resolution than pre-ILT observations at similar frequencies. The six targets that we have imaged (3C 196, 3C 225, 3C 273, 3C 295, 3C 298 and 3C 380) are bright radio sources with compact structures. By comparing our data of 3C 196 and 3C 273 with observations at higher frequencies, we investigate their spatially resolved radio spectral properties. Our success shows that at frequencies down to 30 MHz, sub-arcsecond imaging with the ILT is possible. Further analysis is needed to determine the feasibility of observations of fainter sources or sources with less compact emission., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 11 pages, 8 figures, 3 tables

The prominent radio source Hercules A features complex structures in its radio lobes. Although it is one of the most comprehensively studied sources in the radio sky, the origin of the ring structures in the Hercules A radio lobes remains an open question. We present the first sub-arcsecond angular resolution images at low frequencies (<300 MHz) of Hercules A, made with the International LOFAR Telescope. With the addition of data from the Karl G. Jansky Very Large Array, we mapped the structure of the lobes from 144 MHz to 7 GHz. We explore the origin of the rings within the lobes of Hercules A, and test whether their properties are best described by a shock model, where shock waves are produced by the jet propagating in the radio lobe, or by an inner-lobe model, where the rings are formed by decelerated jetted plasma. From spectral index mapping our large frequency coverage reveals that the curvature of the different ring spectra increases with distance away from the central active galactic nucleus. We demonstrate that the spectral shape of the rings is consistent with synchrotron aging, which speaks in favor of an inner-lobe model where the rings are formed from the deposition of material from past periods of intermittent core activity., Comment: Accepted for a special issue of A&A on sub-arcsecond imaging with LOFAR

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

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

Coherent radio emission from stars can be used to constrain fundamental coronal plasma parameters, such as plasma density and magnetic field strength. It is also a probe of chromospheric and magnetospheric acceleration mechanisms. Close stellar binaries, such as RS Canum Venaticorum (RS CVn) systems, are particularly interesting as their heightened level of chromospheric activity and possible direct magnetic interaction make them a unique laboratory to study coronal and magnetospheric acceleration mechanisms. RS CVn binaries are known to be radio-bright but coherent radio emission has only conclusively been detected previously in one system. Here, we present a population of 14 coherent radio emitting RS CVn systems. We identified the population in the ongoing LOFAR Two Metre Sky Survey as circularly polarised sources at 144MHz that are astrometrically associated with known RS CVn binaries. We show that the observed emission is powered by the electron cyclotron maser instability. We use numerical calculations of the maser's beaming geometry to argue that the commonly invoked 'loss-cone' maser cannot generate the necessary brightness temperature in some of our detections and that a more efficient instability, such as the shell or horseshoe maser, must be invoked. Such maser configurations are known to operate in planetary magnetospheres. We also outline two acceleration mechanisms that could produce coherent radio emission, one where the acceleration occurs in the chromosphere and one where the acceleration is due to an electrodynamic interaction between the stars. We propose radio and optical monitoring observations that can differentiate between these two mechanisms., Comment: This article consists of 14 pages, 10 normal figures and 4 gifs. The gifs themselves each consist of 32 figures. This article is to be published in Astronomy and Astrophysics, section 8. 'Stellar Atmospheres'

In recent years, thanks to new facilities such as LOFAR capable of sensitive observations, much work has been done on the detection of stellar radio emission at low frequencies. Such emission has commonly been shown to be coherent emission, generally attributed to electron-cyclotron maser emission, and has usually been detected from main-sequence M dwarfs. Here we report the first detection of coherent emission at low frequencies from T Tauri stars, which are known to be associated with high levels of stellar activity. Using LOFAR, we have detected several bright radio bursts at 150 MHz from two weak-line T Tauri stars: KPNO-Tau 14 and LkCa 4. All of the bursts have high brightness temperatures ($10^{13} - 10^{14}\ \mathrm{K}$) and high circular polarization fractions (60 - 90 \%), indicating that they must be due to a coherent emission mechanism. This could be either plasma emission or electron-cyclotron maser (ECM) emission. Due to the exceptionally high brightness temperatures seen in at least one of the bursts ($\geq 10^{14}\ \mathrm{K}$), and the high circular polarization levels, it seems unlikely that plasma emission could be the source and so ECM is favoured as the most likely emission mechanism. Assuming this is the case, the required magnetic field in the emission regions would be 40 - 70 G. We determine that the most likely method of generating ECM emission is plasma co-rotation breakdown in the stellar magnetosphere. There remains the possibility, however, it could be due to an interaction with an orbiting exoplanet., Comment: 16 pages, 9 figures, accepted for publication in A&A

The space-weather conditions that result from stellar winds significantly impact the habitability of exoplanets. The conditions can be calculated from first principles if the necessary boundary conditions -- namely on the plasma density in the outer corona and the radial distance at which the plasma forces the closed magnetic field into an open geometry -- are specified. Low frequency radio observations ($\nu \lesssim 200$ MHz) of plasma and cyclotron emission from stars probe these magneto-ionic conditions. Here we report the detection of low-frequency ($120-167\,{\rm MHz}$) radio emission associated with the dMe6 star WX UMa. If the emission originates in WX UMa's corona, we show that the closed field regions extends to at least $\approx 10$ stellar radii, that is about a factor of a few larger than the solar value, and possibly to $\gtrsim 20$ stellar radii. Our results suggest that the magnetic-field structure of M dwarfs is in between Sun-like and planet-like configurations, where compact over-dense coronal loops with X-ray emitting plasma co-exist with a large-scale magnetosphere with lower plasma density and closed magnetic geometry., Comment: Accepted for publication in A&A, 10 pages, 7 figures, 2 tables

LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire northern sky with 3170 pointings in the frequency range 42-66 MHz, at a resolution of 15 arcsec and at a sensitivity of 1 mJy/b. Here we outline the survey strategy, the observational status, the current calibration techniques, and briefly describe several scientific motivations. We also describe the preliminary public data release. The preliminary images were produced using a fully automated pipeline that aims to correct all direction-independent effects in the data. Whilst the direction-dependent effects, such as those from the ionosphere, are not yet corrected, the images presented in this work are still 10 times more sensitive than previous surveys available at these low frequencies. The preliminary data release covers 740 sqdeg around the HETDEX spring field region at a resolution of 47" with a median noise level of 5 mJy/b. The images and the catalogue with 25,247 sources are publicly released. We demonstrate that the system is capable of reaching an rms noise of 1 mJy/b and the resolution of 15" once direction-dependent effects are corrected for. LoLSS will provide the ultra-low-frequency information for hundreds of thousands of radio sources, providing critical spectral information and producing a unique dataset that can be used for a wide range of science topics such as: the search for high redshift galaxies and quasars, the study of the magnetosphere of exoplanets, and the detection of the oldest populations of cosmic-rays in galaxies, clusters of galaxies, and from AGN activity., Comment: 19 pages, 14 figures, Accepted A&A, catalogue and images on www.lofar-surveys.org

Recently detected coherent low-frequency radio emission from M dwarf systems shares phenomenological similarities with emission produced by magnetospheric processes from the gas giant planets of our Solar System. Such beamed electron-cyclotron maser emission can be driven by a star-planet interaction or a breakdown in co-rotation between a rotating plasma disk and a stellar magnetosphere. Both models suggest that the radio emission could be periodic. Here we present the longest low-frequency interferometric monitoring campaign of an M dwarf system, composed of twenty-one $\approx$8 hour epochs taken in two series of observing blocks separated by a year. We achieved a total on-source time of 6.5 days. We show that the M dwarf binary CR Draconis has a low-frequency 3$\sigma$ detection rate of 90$^{+5}_{-8}$% when a noise floor of $\approx$0.1 mJy is reached, with a median flux density of 0.92 mJy, consistent circularly polarised handedness, and a median circularly polarised fraction of 66%. We resolve three bright radio bursts in dynamic spectra, revealing the brightest is elliptically polarised, confined to 4 MHz of bandwidth centred on 170 MHz, and reaches a flux density of 205 mJy. The burst structure is mottled, indicating it consists of unresolved sub-bursts. Such a structure shares a striking resemblance with the low-frequency emission from Jupiter. We suggest the near-constant detection of high brightness temperature, highly-circularly-polarised radiation that has a consistent circular polarisation handedness implies the emission is produced via the electron-cyclotron maser instability. Optical photometric data reveal the system has a rotation period of 1.984$\pm$0.003 days. We observe no periodicity in the radio data, but the sampling of our radio observations produces a window function that would hide the near two-day signal., Comment: Accepted for publication in A&A, 16 pages, 7 figures, 2 tables

The recently discovered colliding-wind binary (CWB) Apep has been shown to emit luminously from radio to X-rays, with the emission driven by a binary composed of two Wolf-Rayet (WR) stars of one carbon-sequence (WC8) and one nitrogen-sequence (WN4-6b). Mid-infrared imaging revealed a giant spiral dust plume that is reminiscent of a pinwheel nebula but with additional features that suggest Apep is a unique system. We have conducted observations with the Australian Long Baseline Array to resolve Apep's radio emission on milliarcsecond scales, allowing us to relate the geometry of the wind-collision region to that of the spiral plume. The observed radio emission shows a bow-shaped structure, confirming its origin as a wind-collision region. The shape and orientation of this region is consistent with being originated by the two stars and with being likely dominated by the stronger wind of the WN4-6b star. This shape allowed us to provide a rough estimation of the opening angle of $\sim 150^\circ$ assuming ideal conditions. The orientation and opening angle of the emission also confirms it as the basis for the spiral dust plume. We also provide estimations for the two stars in the system to milliarcsecond precision. The observed radio emission, one order of magnitude brighter and more luminous than any other known non-thermal radio-emitting CWB, confirms it is produced by an extremely powerful wind collision. Such a powerful wind-collision region is consistent with Apep being a binary composed of two WR stars, so far the first unambiguously confirmed system of its kind., Comment: 9 pages, 3 figures, accepted for publication in MNRAS

Spectral variability of radio sources encodes information about the conditions of intervening media, source structure, and emission processes. With new low-frequency radio interferometers observing over wide fractional bandwidths, studies of spectral variability for a large population of extragalactic radio sources are now possible. Using two epochs of observations from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey that were taken one year apart, we search for spectral variability across 100--230 MHz for 21,558 sources. We present methodologies for detecting variability in the spectrum between epochs and for classifying the type of variability: either as a change in spectral shape or as a uniform change in flux density across the bandwidth. We identify 323 sources with significant spectral variability over a year-long timescale. Of the 323 variable sources, we classify 51 of these as showing a significant change in spectral shape. Variability is more prevalent in peaked-spectrum sources, analogous to gigahertz-peaked spectrum and compact steep-spectrum sources, compared to typical radio galaxies. We discuss the viability of several potential explanations of the observed spectral variability, such as interstellar scintillation and jet evolution. Our results suggest that the radio sky in the megahertz regime is more dynamic than previously suggested., Comment: Accepted for publication in MNRAS

The LOFAR Two-metre Sky Survey (LoTSS) will cover the full northern sky and, additionally, aims to observe the LoTSS deep fields to a noise level of ~10 microJy/bm over several tens of square degrees in areas that have the most extensive ancillary data. This paper presents the ELAIS-N1 deep field, the deepest of the LoTSS deep fields to date. With an effective observing time of 163.7 hours, it reaches a root mean square (RMS) noise level below 20 microJy/bm in the central region (and below 30 microJy/bm over 10 square degrees). The resolution is 6 arcsecs and 84862 radio sources were detected in the full area (68 sq. deg.) with 74127 sources in the highest quality area at less than 3 degrees from the pointing centre. The observation reaches a sky density of more than 5000 sources per sq. deg. in the central ~5 sq. deg. region. We present the calibration procedure, which addresses the special configuration of some observations and the extended bandwidth covered (115 to 177 MHz; central frequency 146.2 MHz) compared to standard LoTSS. We also describe the methods used to calibrate the flux density scale using cross-matching with sources detected by other radio surveys in the literature. We find the flux density uncertainty related to the flux density scale to be ~6.5%. By studying the variations of the flux density measurements between different epochs, we show that relative flux density calibration is reliable out to about a 3 degree radius, but that additional flux density uncertainty is present for all sources at about the 3 per cent level; this is likely to be associated with residual calibration errors, and is shown to be more significant in datasets with poorer ionosphere conditions. We also provide intra-band spectral indices, which can be useful to detect sources with unusual spectral properties. The final uncertainty in the flux densities is estimated to be ~10% for ELAIS-N1., Comment: Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication. 21 pages, 17 figures

Magnetospheric processes seen in gas-giants such as aurorae and circularly-polarized cyclotron maser radio emission have been detected from some brown dwarfs. However, previous radio observations targeted known brown dwarfs discovered via their infrared emission. Here we report the discovery of BDR J1750+3809, a circularly polarized radio source detected around 144 MHz with the LOFAR telescope. Follow-up near-infrared photometry and spectroscopy show that BDR J1750+3809 is a cold methane dwarf of spectral type T$6.5\pm 1$ at a distance of $65^{+9}_{-8}\,{\rm pc}$. The quasi-quiescent radio spectral luminosity of BDR J1750+3809 is $\approx 5\times 10^{15}\,{\rm erg}\,{\rm s}^{-1}\,{\rm Hz}^{-1}$ which is over two orders of magnitude larger than that of the known population of comparable spectral type. This could be due to a preferential geometric alignment or an electrodynamic interaction with a close companion. In addition, as the emission is expected to occur close to the electron gyro-frequency, the magnetic field strength at the emitter site in BDR J1750+3809 is $B\gtrsim 25\,{\rm G}$, which is comparable to planetary-scale magnetic fields. Our discovery suggests that low-frequency radio surveys can be employed to discover sub-stellar objects that are too cold to be detected in infrared surveys., Comment: Accepted for publication in ApJL

The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark discrepancy was found in the observed outflow speed of the gas as measured spectroscopically in the line-of-sight direction compared to the proper motion expansion of the dust in the sky plane. This enigmatic behaviour arises at the wind base within the central Wolf-Rayet binary: a system that has so far remained spatially unresolved. Here we present an updated proper motion study deriving the expansion speed of Apep's dust plume over a two-year baseline that is four times slower than the spectroscopic wind speed, confirming and strengthening the previous finding. We also present the results from high-angular-resolution near-infrared imaging studies of the heart of the system, revealing a close binary with properties matching a Wolf-Rayet colliding-wind system. Based on these new observational constraints, an improved geometric model is presented yielding a close match to the data, constraining the orbital parameters of the Wolf-Rayet binary and lending further support to the anisotropic wind model., Comment: This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society. 17 pages, 9 figures, 6 tables

Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star., Comment: Accepted for publication in MNRAS, 10 pages, 9 figures, 1 table

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc., Comment: 19 pages, 17 figures; accepted by MNRAS

Low frequency ($\nu\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active flare-star (Lynch et al. 2017) that is not representative of the wider stellar population. Here we report the detection of low-frequency radio emission from a quiescent star, GJ 1151--- a member of the most common stellar type (red dwarf or spectral class M) in the Galaxy. The characteristics of the emission are similar to those of planetary auroral emissions (Zarka 1998) (e.g. Jupiter's decametric emission), suggesting a coronal structure dominated by a global magnetosphere with low plasma density. Our results show that large-scale currents that power radio aurorae operate over a vast range of mass and atmospheric composition, ranging from terrestrial planets to main-sequence stars. The Poynting flux required to produce the observed radio emission cannot be generated by GJ 1151's slow rotation, but can originate in a sub-Alfv\'{e}nic interaction of its magnetospheric plasma with a short-period exoplanet. The emission properties are consistent with theoretical expectations (Zarka 2007; Lanza 2009; Saur et al. 2013; Turnpenney et al. 2018) for interaction with an Earth-size planet in a $\sim 1-5$ day-long orbit., Comment: Accepted version, Nat Astron (2020)

Serendipitously, we have rediscovered radio emission from the binary system Ross 867 (M4.5V) and Ross 868 (M3.5V) while inspecting archival Giant Metrewave Radio Telescope (GMRT) observations. The binary system consists of two M-dwarf stars that share common characteristics such as spectral type, astrometric parameters, age and emission at infrared, optical and X-rays frequencies. The GMRT data at 610 MHz taken on July 2011 shows that the radio emission from Ross 867 is polarized and highly variable on hour time scales with a peak flux of 10.4 $\pm$ 0.7 mJy/beam. Additionally, after reviewing archival data from several observatories (VLA, GMRT, JVLA and LOFAR), we confirm that although both stars are likely coeval, only Ross 867 has been detected, while Ross 868 remains undetected at radio wavelengths. As they have a a large orbital separation, this binary stellar system provides a coeval laboratory to examine and constrain the stellar properties linked to radio activity in M dwarfs. We speculate that the observed difference in radio activity between the dwarfs could be due to vastly different magnetic field topologies or that Ross 867 has an intrinsically different dynamo., Comment: Accepted for publication in A&A. Language edition included

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories., Comment: 38 pages, 12 figures, accepted for publication in PASA

Obscuration of quasars by accreted gas and dust, or dusty intervening galaxies, can cause active galactic nuclei (AGN) to be missed in optically-selected surveys. Radio observations can overcome this dust bias. In particular, radio surveys searching for HI absorption inform us on how the AGN can impact on the cold neutral gas medium within the host galaxy, or the population of intervening galaxies through the observed line of sight gas kinematics. We present the results of a HI absorption line survey at $0.4 < z < 1$ towards 34 obscured quasars with the Australian SKA Pathfinder (ASKAP) commissioning array. We detect three HI absorption lines, with one of these systems previously unknown. Through optical follow-up for two sources, we find that in all detections the HI gas is associated with the AGN, and hence that these AGN are obscured by material within their host galaxies. Most of our sample are compact, and in addition, are either gigahertz peaked spectrum (GPS), or steep spectrum (CSS) sources, both thought to represent young or recently re-triggered radio AGN. The radio spectral energy distribution classifications for our sample agree with galaxy evolution models in which the obscured AGN has only recently become active. Our associated HI detection rate for GPS and compact SS sources matches those of other surveys towards such sources. We also find shallow and asymmetric HI absorption features, which agrees with previous findings that the cold neutral medium in compact radio galaxies is typically kinematically disturbed by the AGN., Comment: Accepted for publication in MNRAS. 19 pages, 8 figures

We have examined the giant radio galaxy 3C~236 using LOFAR at 143 MHz down to an angular resolution of 7", in combination with observations at higher frequencies. We have used the low frequency data to derive spectral index maps with the highest resolution yet at these low frequencies. We confirm a previous detection of an inner hotspot in the north-west lobe and for the first time observe that the south-east lobe hotspot is in fact a triple hotspot, which may point to an intermittent source activity. Also, the spectral index map of 3C 236 shows that the spectral steepening at the inner region of the northern lobe is prominent at low frequencies. The outer regions of both lobes show spectral flattening, in contrast with previous high frequency studies. We derive spectral age estimates for the lobes, as well as particle densities of the IGM at various locations. We propose that the morphological differences between the lobes are driven by variations in the ambient medium density as well as the source activity history., Comment: 10 pages, 6 6figures

Extragalactic peaked-spectrum radio sources are thought to be the progenitors of larger, radio-loud active galactic nuclei (AGN). Synchrotron self-absorption (SSA) has often been identified as the cause of their spectral peak. The identification of new megahertz-peaked spectrum sources from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey provides an opportunity to test how radio sources with spectral peaks below 1 GHz fit within this evolutionary picture. We observed six peaked-spectrum sources selected from the GLEAM survey, three that have spectral characteristics which violate SSA and three that have spectral peaks below 230 MHz, with the Very Long Baseline Array at 1.55 and 4.96 GHz. We present milliarcsecond resolution images of each source and constrain their morphology, linear size, luminosity, and magnetic field strength. Of the sources that are resolved by our study, the sources that violate SSA appear to be compact doubles, while the sources with peak frequencies below 230 MHz have core-jet features. We find that all of our sources are smaller than expected from SSA by factors of >20. We also find that component magnetic field strengths calculated from SSA are likely inaccurate, differing by factors of >5 from equipartition estimates. The calculated equipartition magnetic field strengths more closely resemble estimates from previously studied gigahertz-peaked spectrum sources. Exploring a model of the interaction between jets and the interstellar medium, we demonstrate that free-free absorption (FFA) can accurately describe the linear sizes and peak frequencies of our sources. Our findings support the theory that there is a fraction of peaked-spectrum sources whose spectral peaks are best modelled by FFA, implying our understanding of the early stages of radio AGN is incomplete., Comment: Accepted for publication in Astronomy & Astrophysics (A&A) on 16 July 2019. 13 pages, 6 figures

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$., Comment: Accepted for publication in A&A

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS., Comment: 16 figures, 1 table and 22 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

The LOFAR Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey of the Northern sky. The LoTSS First Data Release (DR1) presents 424 square degrees of radio continuum observations over the HETDEX Spring Field (10h45m00s $<$ right ascension $<$ 15h30m00s and 45$^\circ$00$'$00$'$ $<$ declination $<$ 57$^\circ$00$'$00$''$) with a median sensitivity of 71$\mu$Jy/beam and a resolution of 6$''$. In this paper we present photometric redshifts (photo-$z$) for 94.4% of optical sources over this region that are detected in the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 3$\pi$ steradian survey. Combining the Pan-STARRS optical data with mid-infrared photometry from the Wide-field Infrared Survey Explorer, we estimate photo-$z$s using a novel hybrid photometric redshift methodology optimised to produce the best possible performance for the diverse sample of radio continuum selected sources. For the radio-continuum detected population, we find an overall scatter in the photo-$z$ of 3.9% and an outlier fraction ($\left | z_{\rm{phot}} - z_{\rm{spec}} \right | / (1+z_{\rm{spec}}) > 0.15$) of 7.9%. We also find that, at a given redshift, there is no strong trend in photo-$z$ quality as a function of radio luminosity. However there are strong trends as a function of redshift for a given radio luminosity, a result of selection effects in the spectroscopic sample and/or intrinsic evolution within the radio source population. Additionally, for the sample of sources in the LoTSS First Data Release with optical counterparts, we present rest-frame optical and mid-infrared magnitudes based on template fits to the consensus photometric (or spectroscopic when available) redshift., Comment: 16 figures, 1 table and 18 pages. This paper is part of the LOFAR surveys Data Release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, Volume 622. The catalogues and images from the data release will be publicly available at http://www.lofar-surveys.org upon publication of the volume

Historically, the blazar population has been poorly understood at low frequencies because survey sensitivity and angular resolution limitations have made it difficult to identify megahertz counterparts. We used the LOFAR Two-Metre Sky Survey (LoTSS) first data release value-added catalogue (LDR1) to study blazars in the low-frequency regime with unprecedented sensitivity and resolution. We identified radio counterparts to all $98$ known sources from the Third \textit{Fermi}-LAT Point Source Catalogue (3FGL) or Roma-BZCAT Multi-frequency Catalogue of Blazars ($5^{\mathrm{th}}$ edition) that fall within the LDR1 footprint. Only the 3FGL unidentified $\gamma$-ray sources (UGS) could not be firmly associated with an LDR1 source; this was due to source confusion. We examined the redshift and radio luminosity distributions of our sample, finding flat-spectrum radio quasars (FSRQs) to be more distant and more luminous than BL Lacertae objects (BL Lacs) on average. Blazars are known to have flat spectra in the gigahertz regime but we found this to extend down to $144$ MHz, where the radio spectral index, $\alpha$, of our sample is $-0.17 \pm 0.14$. For BL Lacs, $\alpha = -0.13 \pm 0.16$ and for FSRQs, $\alpha = -0.15 \pm 0.17$. We also investigated the radio-to-$\gamma$-ray connection for the $30$ $\gamma$-ray-detected sources in our sample. We find Pearson's correlation coefficient is $0.45$ ($p = 0.069$). This tentative correlation and the flatness of the spectral index suggest that the beamed core emission contributes to the low-frequency flux density. We compare our sample distribution with that of the full LDR1 on colour-colour diagrams, and we use this information to identify possible radio counterparts to two of the four UGS within the LDR1 field. We will refine our results as LoTSS continues., Comment: 7 figures, 6 tables and 11 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

We constructed a sample of 23,344 radio-loud active galactic nuclei (RLAGN) from the catalogue derived from the LOFAR Two-Metre Sky Survey (LoTSS) survey of the HETDEX Spring field. Although separating AGN from star-forming galaxies remains challenging, the combination of spectroscopic and photometric techniques we used gives us one of the largest available samples of candidate RLAGN. We used the sample, combined with recently developed analytical models, to investigate the lifetime distribution of RLAGN. We show that large or giant powerful RLAGN are probably the old tail of the general RLAGN population, but that the low-luminosity RLAGN candidates in our sample, many of which have sizes $<100$ kpc, either require a very different lifetime distribution or have different jet physics from the more powerful objects. We then used analytical models to develop a method of estimating jet kinetic powers for our candidate objects and constructed a jet kinetic luminosity function based on these estimates. These values can be compared to observational quantities, such as the integrated radiative luminosity of groups and clusters, and to the predictions from models of RLAGN feedback in galaxy formation and evolution. In particular, we show that RLAGN in the local Universe are able to supply all the energy required per comoving unit volume to counterbalance X-ray radiative losses from groups and clusters and thus prevent the hot gas from cooling. Our computation of the kinetic luminosity density of local RLAGN is in good agreement with other recent observational estimates and with models of galaxy formation., Comment: 21 pages, 11 colour figures. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available at https://lofar-surveys.org/ upon publication of the journal

Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections., Comment: 7 figures, 2 tables and 10 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release., Comment: 8 figures, 5 tables and 23 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available at https://lofar-surveys.org/ upon publication in the journal