-
Like the magnetised planets in our Solar System, magnetised exoplanets should
emit strongly at radio wavelengths. Radio emission directly traces the
planetary magnetic fields and radio detections can place constraints on the
physical parameters of these features. Large comparative studies of predicted
radio emission characteristics for the known population of exoplanets help to
identify what physical parameters could be key for producing bright, observable
radio emission. Since the last comparative study, many thousands of exoplanets
have been discovered. We report new estimates for the radio flux densities and
maximum emission frequencies for the current population of known exoplanets
orbiting pre-main sequence and main-sequence stars with spectral types F-M. The
set of exoplanets predicted to produce observable radio emission are Hot
Jupiters orbiting young stars. The youth of these system predicts strong
stellar magnetic fields and/or dense winds, which are key for producing bright,
observable radio emission. We use a new all-sky circular polarisation Murchison
Widefield Array survey to place sensitive limits on 200 MHz emission from
exoplanets, with $3\sigma$ values ranging from 4.0 - 45.0 mJy. Using a targeted
Giant Metre Wave Radio Telescope observing campaign, we also report a $3\sigma$
upper limit of 4.5 mJy on the radio emission from V830 Tau b, the first Hot
Jupiter to be discovered orbiting a pre-main sequence star. Our limit is the
first to be reported for the low-frequency radio emission from this source.
-
Recently, low frequency, broadband radio emission has been observed
accompanying bright meteors by the Long Wavelength Array (LWA). The broadband
spectra between 20 and 60 MHz were captured for several events, while the
spectral index (dependence of flux density on frequency, with $S_\nu \propto
\nu^\alpha$) was estimated to be $-4\pm1$ during the peak of meteor afterglows.
Here we present a survey of meteor emission and other transient events using
the Murchison Widefield Array (MWA) at 72-103 MHz. In our 322-hour survey, down
to a $5\sigma$ detection threshold of 3.5 Jy/beam, no transient candidates were
identified as intrinsic emission from meteors. We derived an upper limit of
-3.7 (95% confidence limit) on the spectral index in our frequency range. We
also report detections of other transient events, like reflected FM broadcast
signals from small satellites, conclusively demonstrating the ability of the
MWA to detect and track space debris on scales as small as 0.1 m in low Earth
orbits.
-
The Murchison Widefield Array (MWA), and its recently-developed Voltage
Capture System (VCS), facilitates extending the low-frequency range of pulsar
observations at high-time and -frequency resolution in the Southern Hemisphere,
providing further information about pulsars and the ISM. We present the results
of an initial time-resolved census of known pulsars using the MWA. To
significantly reduce the processing load, we incoherently sum the detected
powers from the 128 MWA tiles, which yields ~10% of the attainable sensitivity
of the coherent sum. This preserves the large field-of-view (~450 deg2 at 185
MHz), allowing multiple pulsars to be observed simultaneously. We developed a
WIde-field Pulsar Pipeline (WIPP) that processes the data from each observation
and automatically folds every known pulsar located within the beam. We have
detected 50 pulsars to date, 6 of which are millisecond pulsars. This is
consistent with our expectation, given the telescope sensitivity and the sky
coverage of the processed data (~17,000 deg2). For ten pulsars, we present the
lowest-frequency detections published. For a subset of the pulsars, we present
multi-frequency pulse profiles by combining our data with published profiles
from other telescopes. Since the MWA is a low-frequency precursor to the Square
Kilometre Array (SKA), we use our census results to forecast that a survey
using Phase 1 of SKA-Low (SKA1-Low) can potentially detect around 9400 pulsars.
-
We search for an isotropic stochastic gravitational-wave background (GWB) in
the newly released $11$-year dataset from the North American Nanohertz
Observatory for Gravitational Waves (NANOGrav). While we find no significant
evidence for a GWB, we place constraints on a GWB from a population of
supermassive black-hole binaries, from cosmic strings, and from a primordial
GWB. For the first time, we find that the GWB upper limits and detection
statistics are sensitive to the Solar System ephemeris (SSE) model used, and
that SSE errors can mimic a GWB signal. To mitigate this effect, we developed
and implemented a novel approach that bridges systematic SSE differences,
producing the first PTA constraints that are robust against SSE uncertainties.
We place a $95\%$ upper limit on the GW strain amplitude of
$A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a
fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations
modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year
limit, calculated using the same procedure. Previous PTA upper limits on the
GWB will need revision in light of SSE systematic uncertainties. We also
characterize the combined influence of the mass-density of stars in galactic
cores, the eccentricity of binaries at formation, and the relation between the
mass of the central supermassive black hole and the galactic bulge (the
$M_\mathrm{BH}-M_\mathrm{bulge}$ relation). We constrain cosmic-string tension
on the basis of recent simulations, yielding an SSE-marginalized 95\% upper
limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of
$\sim 2$ better than the NANOGrav $9$-year constraints. We then use our new
Bayesian SSE model to limit the energy density of primordial GWBs,
corresponding to $\Omega_\mathrm{GWB}(f)h^2<3.4 \times 10^{-10}$ for a
radiation-dominated inflationary era. [ABRIDGED]
-
The Murchison Widefield Array (MWA), located in Western Australia, is one of
the low-frequency precursors of the international Square Kilometre Array (SKA)
project. In addition to pursuing its own ambitious science program, it is also
a testbed for wide range of future SKA activities ranging from hardware,
software to data analysis. The key science programs for the MWA and SKA require
very high dynamic ranges, which challenges calibration and imaging systems.
Correct calibration of the instrument and accurate measurements of source flux
densities and polarisations require precise characterisation of the telescope's
primary beam. Recent results from the MWA GaLactic Extragalactic All-sky MWA
(GLEAM) survey show that the previously implemented Average Embedded Element
(AEE) model still leaves residual polarisations errors of up to 10-20 % in
Stokes Q. We present a new simulation-based Full Embedded Element (FEE) model
which is the most rigorous realisation yet of the MWA's primary beam model. It
enables efficient calculation of the MWA beam response in arbitrary directions
without necessity of spatial interpolation. In the new model, every dipole in
the MWA tile (4 x 4 bow-tie dipoles) is simulated separately, taking into
account all mutual coupling, ground screen and soil effects, and therefore
accounts for the different properties of the individual dipoles within a tile.
We have applied the FEE beam model to GLEAM observations at 200 - 231 MHz and
used false Stokes parameter leakage as a metric to compare the models. We have
determined that the FEE model reduced the magnitude and declination-dependent
behaviour of false polarisation in Stokes Q and V while retaining low levels of
false polarisation in Stokes U.
-
The discovery of the first electromagnetic counterpart to a gravitational
wave signal has generated follow-up observations by over 50 facilities
world-wide, ushering in the new era of multi-messenger astronomy. In this
paper, we present follow-up observations of the gravitational wave event
GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label
AT2017gfo) by 14 Australian telescopes and partner observatories as part of
Australian-based and Australian-led research programs. We report early- to
late-time multi-wavelength observations, including optical imaging and
spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio
bursts. Our optical spectra reveal that the transient source afterglow cooled
from approximately 6400K to 2100K over a 7-day period and produced no
significant optical emission lines. The spectral profiles, cooling rate, and
photometric light curves are consistent with the expected outburst and
subsequent processes of a binary neutron star merger. Star formation in the
host galaxy probably ceased at least a Gyr ago, although there is evidence for
a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore
unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr
coalescence time could produce such a merger. The displacement (about 2.2 kpc)
of the binary star system from the centre of the main galaxy is not unusual for
stars in the host galaxy or stars originating in the merging galaxy, and
therefore any constraints on the kick velocity imparted to the progenitor are
poor.
-
Merging neutron stars offer an exquisite laboratory for simultaneously
studying strong-field gravity and matter in extreme environments. We establish
the physical association of an electromagnetic counterpart EM170817 to
gravitational waves (GW170817) detected from merging neutron stars. By
synthesizing a panchromatic dataset, we demonstrate that merging neutron stars
are a long-sought production site forging heavy elements by r-process
nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to
classical short gamma-ray bursts with ultra-relativistic jets. Instead, we
suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the
jet elegantly explains the low-luminosity gamma-rays, the high-luminosity
ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit
that all merging neutron stars may lead to a wide-angle cocoon breakout;
sometimes accompanied by a successful jet and sometimes a choked jet.
-
We present techniques developed to calibrate and correct Murchison Widefield
Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The
extremely wide field-of-view, excellent instantaneous (u, v)-coverage and
sensitivity to degree-scale structure that the MWA provides enable instrumental
calibration, removal of instrumental artefacts, and correction for ionospheric
Faraday rotation through imaging techniques. With the demonstrated polarimetric
capabilities of the MWA, we discuss future directions for polarimetric science
at low frequencies to answer outstanding questions relating to polarised source
counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the
nature of the interstellar and intergalactic media, and the solar environment.
-
We report on a search for Fast Radio Bursts (FRBs) with the Green Bank
Northern Celestial Cap (GBNCC) Pulsar Survey at 350 MHz. Pointings amounting to
a total on-sky time of 61 days were searched to a DM of 3000 pc cm$^{-3}$ while
the rest (23 days; 29% of the total time) were searched to a DM of 500 pc
cm$^{-3}$. No FRBs were detected in the pointings observed through May 2016. We
estimate a 95% confidence upper limit on the FRB rate of $3.6\times 10^3$ FRBs
sky$^{-1}$ day$^{-1}$ above a peak flux density of 0.63 Jy at 350 MHz for an
intrinsic pulse width of 5 ms. We place constraints on the spectral index
$\alpha$ by running simulations for different astrophysical scenarios and
cumulative flux density distributions. The non-detection with GBNCC is
consistent with the 1.4-GHz rate reported for the Parkes surveys for $\alpha >
+0.35 $ in the absence of scattering and free-free absorption and $\alpha >
-0.3$ in the presence of scattering, for a Euclidean flux distribution. The
constraints imply that FRBs exhibit either a flat spectrum or a spectral
turnover at frequencies above 400 MHz. These constraints also allow estimation
of the number of bursts that can be detected with current and upcoming surveys.
We predict that CHIME may detect anywhere from several to $\sim$50 FRBs a day
(depending on model assumptions), making it well suited for interesting
constraints on spectral index, the log $N$-log $S$ slope and pulse profile
evolution across its bandwidth (400-800 MHz).
-
We describe the design and performance of the Engineering Development Array
(EDA), which is a low frequency radio telescope comprising 256
dual-polarisation dipole antennas working as a phased-array. The EDA was
conceived of, developed, and deployed in just 18 months via re-use of Square
Kilometre Array (SKA) precursor technology and expertise, specifically from the
Murchison Widefield Array (MWA) radio telescope. Using drift scans and a model
for the sky brightness temperature at low frequencies, we have derived the
EDA's receiver temperature as a function of frequency. The EDA is shown to be
sky-noise limited over most of the frequency range measured between 60 and 240
MHz. By using the EDA in interferometric mode with the MWA, we used calibrated
visibilities to measure the absolute sensitivity of the array. The measured
array sensitivity matches very well with a model based on the array layout and
measured receiver temperature. The results demonstrate the practicality and
feasibility of using MWA-style precursor technology for SKA-scale stations. The
modular architecture of the EDA allows upgrades to the array to be rolled out
in a staged approach. Future improvements to the EDA include replacing the
second stage beamformer with a fully digital system, and to transition to using
RF-over-fibre for the signal output from first stage beamformers.
-
We have detected four flares from UV Ceti at 154 MHz using the Murchison
Widefield Array. The flares have flux densities between 10--65 mJy --- a factor
of 100 fainter than most flares in the literature at these frequencies --- and
are only detected in polarization. The circular polarized fractions are limited
to $>27$% at 3$\sigma$ confidence and two of the flares exhibit polarity
reversal. We suggest that these flares occur periodically on a time scale
consistent with the rotational period of UV Ceti. During the brightest observed
flare, we also detect significant linear polarization with polarization
fraction $>18$%. Averaging the data in 6-minute, 10 MHz frequency bins we find
that the flux density of these flares does not vary over the 30 MHz bandwidth
of the Murchison Widefield Array, however we cannot rule out finer
time-frequency structure. Using the measured flux densities for the flares, we
estimate brightness temperatures between $(10^{13}-10^{14})\,$K, indicative of
a coherent emission mechanism. The brightness temperature and polarization
characteristics point to the electron cyclotron maser mechanism. We also
calculate the flare rates given our four observed flares and compare them to
flare rates for the set of M dwarf stars with known 100--200 MHz flares. Our
measurement is the first for flares with intensities $<100$ mJy at 100-200 MHz.
-
We report the results of a 154 MHz survey to search for emission from
exoplanets located in the Upper Scorpius subgroup of the Sco Cen OB2
Association, the closest substantial region of recent star formation. This
survey is different from previous efforts in that it is the first to target
exoplanets orbiting Myr-old stars. Young exoplanet systems are expected to be
the best candidates for radio detections given the higher magnetic field
strengths predicted for young planets as well as the stronger and more dense
stellar wind expected for the host stars. The radio emission from exoplanets is
expected to be highly circularly polarised therefore we restricted our search
to the circular polarisation images rather than the total intensity images. We
carried out two different search methods using this data. The first method was
a targeted search for exoplanet emission using catalogues of known stars and
Hot Jupiters within the Upper Scorpius field. The second search method was a
blind search for highly circularly polarised sources in the field and for
sources identified only in our polarisation images. Both the blind and targeted
search resulted in non-detections with typical 3-sigma flux density limits of
4--235 mJy over timescales of 1.87--1000 minutes. In particular, we place the
first limits on low-frequency emission from the Hot Jupiter systems WASP-17 b
and K2-33 b. These are the first results from a larger program to
systematically search for low-frequency radio emission from planets orbiting
young stars.
-
Many astronomical sources produce transient phenomena at radio frequencies,
but the transient sky at low frequencies (<300 MHz) remains relatively
unexplored. Blind surveys with new widefield radio instruments are setting
increasingly stringent limits on the transient surface density on various
timescales. Although many of these instruments are limited by classical
confusion noise from an ensemble of faint, unresolved sources, one can in
principle detect transients below the classical confusion limit to the extent
that the classical confusion noise is independent of time. We develop a
technique for detecting radio transients that is based on temporal matched
filters applied directly to time series of images rather than relying on
source-finding algorithms applied to individual images. This technique has
well-defined statistical properties and is applicable to variable and transient
searches for both confusion-limited and non-confusion-limited instruments.
Using the Murchison Widefield Array as an example, we demonstrate that the
technique works well on real data despite the presence of classical confusion
noise, sidelobe confusion noise, and other systematic errors. We searched for
transients lasting between 2 minutes and 3 months. We found no transients and
set improved upper limits on the transient surface density at 182 MHz for flux
densities between ~20--200 mJy, providing the best limits to date for hour- and
month-long transients.
-
The Detection of redshifted 21 cm emission from the epoch of reionization
(EoR) is a challenging task owing to strong foregrounds that dominate the
signal. In this paper, we propose a general method, based on the delay spectrum
approach, to extract HI power spectra that is applicable to tracking
observations using an imaging radio interferometer (Delay Spectrum with Imaging
Arrays (DSIA)). Our method is based on modelling the HI signal taking into
account the impact of wide field effects such as the $w$-term which are then
used as appropriate weights in cross-correlating the measured visibilities. Our
method is applicable to any radio interferometer that tracks a phase center and
could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the
literature the delay spectrum approach has been implemented for near-redundant
baselines using drift scan observations. In this paper we explore the scheme
for non-redundant tracking arrays, and this is the first application of delay
spectrum methodology to such data to extract the HI signal. We analyze 3 hours
of MWA tracking data on the EoR1 field. We present both 2-dimensional
($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis.
Our results are in agreement with the findings of other pipelines developed to
analyse the MWA EoR data.
-
We present and evaluate several strategies to search for prompt,
low-frequency radio emission associated with gravitational wave transients
using the Murchison Widefield Array (MWA). As we are able to repoint the MWA on
timescales of tens of seconds, we can search for the dispersed radio signal
that has been predicted to originate along with or shortly after a neutron
star-neutron star merger. We find that given the large, 600 deg^2 instantaneous
field-of-view of the MWA we can cover a significant fraction of the predicted
gravitational wave error region, although due to the complicated geometry of
the latter we only cover >50% of the error region for approximately 5% of
events, and roughly 15% of events will be located <10 deg from the MWA pointing
center such that they will be covered in the radio images. For optimal
conditions our limiting flux density for a 10-s long transient would be 0.1 Jy,
increasing to about 1 Jy for a wider range of events. This corresponds to
luminosity limits of 1e38-1e39 erg/s based on expectations for the distances of
the gravitational wave transients, which should be sufficient to detect or
significantly constrain a range of models for prompt emission.
-
We report the discovery and timing results for five millisecond pulsars
(MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726,
J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted
with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3
yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low
eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1
$M_\odot$ companion and possible eclipses that last ~10% of the orbit. The
position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are
distant (>3.1 kpc) as determined from their dispersion measures, and none of
them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray
Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by
the PALFA survey to 26 and further demonstrate the power of this survey in
finding distant, highly dispersed MSPs deep in the Galactic plane.
-
The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch
of Reionization (EoR) data and now faces the challenge of overcoming foreground
and systematic contamination to reduce the data to a cosmological measurement.
We introduce several novel analysis techniques such as cable reflection
calibration, hyper-resolution gridding kernels, diffuse foreground model
subtraction, and quality control methods. Each change to the analysis pipeline
is tested against a two dimensional power spectrum figure of merit to
demonstrate improvement. We incorporate the new techniques into a deep
integration of 32 hours of MWA data. This data set is used to place a
systematic-limited upper limit on the cosmological power spectrum of $\Delta^2
\leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent
with other published limits, and a modest improvement (factor of 1.4) over
previous MWA results. From this deep analysis we have identified a list of
improvements to be made to our EoR data analysis strategies. These improvements
will be implemented in the future and detailed in upcoming publications.
-
We present deep polarimetric observations at 154 MHz with the Murchison
Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The
sensitivity available in our deep observations allows an in-band,
frequency-dependent analysis of polarized structure for the first time at long
wavelengths. Our analysis suggests that the polarized structures are dominated
by intrinsic emission but may also have a foreground Faraday screen component.
At these wavelengths, the compactness of the MWA baseline distribution provides
excellent snapshot sensitivity to large-scale structure. The observations are
sensitive to diffuse polarized emission at ~54' resolution with a sensitivity
of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a
sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The
sensitivity allows the effect of ionospheric Faraday rotation to be spatially
and temporally measured directly from the diffuse polarized background. Our
observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear
polarization clearly detectable in ~2 minute snapshots, which would remain
undetectable by interferometers with minimum baseline lengths >110 m at 154
MHz. The brightness temperature of these structures is on average 4 K in
polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that
the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with
a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc
to the polarized emission based on measurements of the in-field pulsar
J2330-2005. We detect four extragalactic linearly polarized point sources
within the field in our compact source survey. Based on the known polarized
source population at 1.4 GHz and non-detections at 154 MHz, we estimate an
upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz.
-
This Supplement provides supporting material for arXiv:1602.08492 . We
briefly summarize past electromagnetic (EM) follow-up efforts as well as the
organization and policy of the current EM follow-up program. We compare the
four probability sky maps produced for the gravitational-wave transient
GW150914, and provide additional details of the EM follow-up observations that
were performed in the different bands.
-
A gravitational-wave (GW) transient was identified in data recorded by the
Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors
on 2015 September 14. The event, initially designated G184098 and later given
the name GW150914, is described in detail elsewhere. By prior arrangement,
preliminary estimates of the time, significance, and sky location of the event
were shared with 63 teams of observers covering radio, optical, near-infrared,
X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In
this Letter we describe the low-latency analysis of the GW data and present the
sky localization of the first observed compact binary merger. We summarize the
follow-up observations reported by 25 teams via private Gamma-ray Coordinates
Network circulars, giving an overview of the participating facilities, the GW
sky localization coverage, the timeline and depth of the observations. As this
event turned out to be a binary black hole merger, there is little expectation
of a detectable electromagnetic (EM) signature. Nevertheless, this first
broadband campaign to search for a counterpart of an Advanced LIGO source
represents a milestone and highlights the broad capabilities of the transient
astronomy community and the observing strategies that have been developed to
pursue neutron star binary merger events. Detailed investigations of the EM
data and results of the EM follow-up campaign are being disseminated in papers
by the individual teams.
-
PSR J1024$-$0719 is a millisecond pulsar that was long thought to be
isolated. However, puzzling results concerning its velocity, distance, and low
rotational period derivative have led to reexamination of its properties. We
present updated radio timing observations along with new and archival optical
data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr)
binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z
\approx -0.9\,$ dex) main sequence star. Such a system can explain most of the
anomalous properties of this pulsar. We suggest that this system formed through
a dynamical exchange in a globular cluster that ejected it into a halo orbit,
consistent with the low observed metallicity for the stellar companion. Further
astrometric and radio timing observations such as measurement of the third
period derivative could strongly constrain the range of orbital parameters.
-
Detection of the Epoch of Reionization HI signal requires a precise
understanding of the intervening galaxies and AGN, both for instrumental
calibration and foreground removal. We present a catalogue of 7394
extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison
Widefield Array Epoch of Reionization observation programme. Motivated by
unprecedented requirements for precision and reliability we develop new methods
for source finding and selection. We apply machine learning methods to
self-consistently classify the relative reliability of 9490 source candidates.
A subset of 7466 are selected based on reliability class and signal-to-noise
ratio criteria. These are statistically cross-matched to four other radio
surveys using both position and flux density information. We find 7369 sources
to have confident matches, including 90 partially resolved sources that split
into a total of 192 sub-components. An additional 25 unmatched sources are
included as new radio detections. The catalogue sources have a median spectral
index of -0.85. Spectral flattening is seen toward lower frequencies with a
median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec.
compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers
approximately 1400 sq. deg. and is complete to approximately 80 mJy within half
beam power. This provides the most reliable discrete source sky model available
to date in the MWA EoR0 field for precision foreground subtraction.
-
Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active
galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by
their relative infrared faintness, resulting in enormous radio-to-infrared flux
density ratios of up to several thousand. We aim to test the hypothesis that
IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact
steep-spectrum (CSS) sources that have a low frequency turnover. We use the
rich radio data set available for the Australia Telescope Large Area Survey
fields, covering the frequency range between 150 MHz and 34 GHz with up to 19
wavebands from different telescopes, and build radio spectral energy
distributions (SEDs) for 34 IFRS. We then study the radio properties of this
class of object with respect to turnover, spectral index, and behaviour towards
higher frequencies. We also present the highest-frequency radio observations of
an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model
the multi-wavelength and radio-far-infrared SED of this source. We find IFRS
usually follow single power laws down to observed frequencies of around 150
MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In
particular, IFRS show statistically significantly steeper radio SEDs than the
broader RL AGN population. Our analysis reveals that the fractions of GPS and
CSS sources in the population of IFRS are consistent with the fractions in the
broader RL AGN population. We find that at least 18% of IFRS contain young AGN,
although the fraction might be significantly higher as suggested by the steep
SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED
modelling of one IFRS shows that it is different from ordinary AGN, although it
is consistent with a composite starburst-AGN model with a star formation rate
of 170 solar masses per year.
-
Model-independent distance constraints to binary millisecond pulsars (MSPs)
are of great value to both the timing observations of the radio pulsars, and
multiwavelength observations of their companion stars. Very Long Baseline
Interferometry (VLBI) astrometry can be employed to provide these
model-independent distances with very high precision via the detection of
annual geometric parallax. Using the Very Long Baseline Array, we have observed
two binary millisecond pulsars, PSR J1022+1001 and J2145-0750, over a two-year
period and measured their distances to be 700 +14 -10 pc and 613 +16 -14 pc
respectively. We use the well-calibrated distance in conjunction with revised
analysis of optical photometry to tightly constrain the nature of their massive
(M ~ 0.85 Msun) white dwarf companions. Finally, we show that several
measurements of their parallax and proper motion of PSR J1022+1001 and PSR
J2145-0750 obtained by pulsar timing array projects are incorrect, differing
from the more precise VLBI values by up to 5 sigma. We investigate possible
causes for the discrepancy, and find that imperfect modeling of the solar wind
is a likely candidate for the timing model errors given the low ecliptic
latitude of these two pulsars.
-
We present 154 MHz Murchison Widefield Array imaging observations and
variability information for a sample of pulsars. Over the declination range
$-80^{\circ} < {\delta} < 10^{\circ}$ we detect 17 known pulsars with mean flux
density greater than 0.3 Jy. We explore the variability properties of this
sample on timescales of minutes to years. For three of these pulsars, PSR
J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar
scintillation and variability on timescales of greater than 2 minutes. One
further pulsar, PSR J0034-0721, showed significant variability, the physical
origins of which are difficult to determine. The dynamic spectra for PSR
J0953+0755 and PSR J0437-4715 show discrete time and frequency structure
consistent with diffractive interstellar scintillation and we present the
scintillation bandwidth and timescales from these observations. The remaining
pulsars within our sample were statistically non-variable. We also explore the
spectral properties of this sample and find spectral curvature in pulsars PSR
J0835-4510, PSR J1752-2806 and PSR J0437-4715.
The Murchison Widefield Array (MWA), and its recently-developed Voltage Capture System (VCS), facilitates extending the low-frequency range of pulsar observations at high-time and -frequency resolution in the Southern Hemisphere, providing further information about pulsars and the ISM. We present the results of an initial time-resolved census of known pulsars using the MWA. To significantly reduce the processing load, we incoherently sum the detected powers from the 128 MWA tiles, which yields ~10% of the attainable sensitivity of the coherent sum. This preserves the large field-of-view (~450 deg2 at 185 MHz), allowing multiple pulsars to be observed simultaneously. We developed a WIde-field Pulsar Pipeline (WIPP) that processes the data from each observation and automatically folds every known pulsar located within the beam. We have detected 50 pulsars to date, 6 of which are millisecond pulsars. This is consistent with our expectation, given the telescope sensitivity and the sky coverage of the processed data (~17,000 deg2). For ten pulsars, we present the lowest-frequency detections published. For a subset of the pulsars, we present multi-frequency pulse profiles by combining our data with published profiles from other telescopes. Since the MWA is a low-frequency precursor to the Square Kilometre Array (SKA), we use our census results to forecast that a survey using Phase 1 of SKA-Low (SKA1-Low) can potentially detect around 9400 pulsars.
We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, from cosmic strings, and from a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. To mitigate this effect, we developed and implemented a novel approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We also characterize the combined influence of the mass-density of stars in galactic cores, the eccentricity of binaries at formation, and the relation between the mass of the central supermassive black hole and the galactic bulge (the $M_\mathrm{BH}-M_\mathrm{bulge}$ relation). We constrain cosmic-string tension on the basis of recent simulations, yielding an SSE-marginalized 95\% upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the NANOGrav $9$-year constraints. We then use our new Bayesian SSE model to limit the energy density of primordial GWBs, corresponding to $\Omega_\mathrm{GWB}(f)h^2<3.4 \times 10^{-10}$ for a radiation-dominated inflationary era. [ABRIDGED]
We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment.
We report on a search for Fast Radio Bursts (FRBs) with the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey at 350 MHz. Pointings amounting to a total on-sky time of 61 days were searched to a DM of 3000 pc cm$^{-3}$ while the rest (23 days; 29% of the total time) were searched to a DM of 500 pc cm$^{-3}$. No FRBs were detected in the pointings observed through May 2016. We estimate a 95% confidence upper limit on the FRB rate of $3.6\times 10^3$ FRBs sky$^{-1}$ day$^{-1}$ above a peak flux density of 0.63 Jy at 350 MHz for an intrinsic pulse width of 5 ms. We place constraints on the spectral index $\alpha$ by running simulations for different astrophysical scenarios and cumulative flux density distributions. The non-detection with GBNCC is consistent with the 1.4-GHz rate reported for the Parkes surveys for $\alpha > +0.35 $ in the absence of scattering and free-free absorption and $\alpha > -0.3$ in the presence of scattering, for a Euclidean flux distribution. The constraints imply that FRBs exhibit either a flat spectrum or a spectral turnover at frequencies above 400 MHz. These constraints also allow estimation of the number of bursts that can be detected with current and upcoming surveys. We predict that CHIME may detect anywhere from several to $\sim$50 FRBs a day (depending on model assumptions), making it well suited for interesting constraints on spectral index, the log $N$-log $S$ slope and pulse profile evolution across its bandwidth (400-800 MHz).
We describe the design and performance of the Engineering Development Array (EDA), which is a low frequency radio telescope comprising 256 dual-polarisation dipole antennas working as a phased-array. The EDA was conceived of, developed, and deployed in just 18 months via re-use of Square Kilometre Array (SKA) precursor technology and expertise, specifically from the Murchison Widefield Array (MWA) radio telescope. Using drift scans and a model for the sky brightness temperature at low frequencies, we have derived the EDA's receiver temperature as a function of frequency. The EDA is shown to be sky-noise limited over most of the frequency range measured between 60 and 240 MHz. By using the EDA in interferometric mode with the MWA, we used calibrated visibilities to measure the absolute sensitivity of the array. The measured array sensitivity matches very well with a model based on the array layout and measured receiver temperature. The results demonstrate the practicality and feasibility of using MWA-style precursor technology for SKA-scale stations. The modular architecture of the EDA allows upgrades to the array to be rolled out in a staged approach. Future improvements to the EDA include replacing the second stage beamformer with a fully digital system, and to transition to using RF-over-fibre for the signal output from first stage beamformers.
We have detected four flares from UV Ceti at 154 MHz using the Murchison Widefield Array. The flares have flux densities between 10--65 mJy --- a factor of 100 fainter than most flares in the literature at these frequencies --- and are only detected in polarization. The circular polarized fractions are limited to $>27$% at 3$\sigma$ confidence and two of the flares exhibit polarity reversal. We suggest that these flares occur periodically on a time scale consistent with the rotational period of UV Ceti. During the brightest observed flare, we also detect significant linear polarization with polarization fraction $>18$%. Averaging the data in 6-minute, 10 MHz frequency bins we find that the flux density of these flares does not vary over the 30 MHz bandwidth of the Murchison Widefield Array, however we cannot rule out finer time-frequency structure. Using the measured flux densities for the flares, we estimate brightness temperatures between $(10^{13}-10^{14})\,$K, indicative of a coherent emission mechanism. The brightness temperature and polarization characteristics point to the electron cyclotron maser mechanism. We also calculate the flare rates given our four observed flares and compare them to flare rates for the set of M dwarf stars with known 100--200 MHz flares. Our measurement is the first for flares with intensities $<100$ mJy at 100-200 MHz.
We report the results of a 154 MHz survey to search for emission from exoplanets located in the Upper Scorpius subgroup of the Sco Cen OB2 Association, the closest substantial region of recent star formation. This survey is different from previous efforts in that it is the first to target exoplanets orbiting Myr-old stars. Young exoplanet systems are expected to be the best candidates for radio detections given the higher magnetic field strengths predicted for young planets as well as the stronger and more dense stellar wind expected for the host stars. The radio emission from exoplanets is expected to be highly circularly polarised therefore we restricted our search to the circular polarisation images rather than the total intensity images. We carried out two different search methods using this data. The first method was a targeted search for exoplanet emission using catalogues of known stars and Hot Jupiters within the Upper Scorpius field. The second search method was a blind search for highly circularly polarised sources in the field and for sources identified only in our polarisation images. Both the blind and targeted search resulted in non-detections with typical 3-sigma flux density limits of 4--235 mJy over timescales of 1.87--1000 minutes. In particular, we place the first limits on low-frequency emission from the Hot Jupiter systems WASP-17 b and K2-33 b. These are the first results from a larger program to systematically search for low-frequency radio emission from planets orbiting young stars.
Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between ~20--200 mJy, providing the best limits to date for hour- and month-long transients.
The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data.
We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1 $M_\odot$ companion and possible eclipses that last ~10% of the orbit. The position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.
This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.
A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.
PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.
Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction.
Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year.
Model-independent distance constraints to binary millisecond pulsars (MSPs) are of great value to both the timing observations of the radio pulsars, and multiwavelength observations of their companion stars. Very Long Baseline Interferometry (VLBI) astrometry can be employed to provide these model-independent distances with very high precision via the detection of annual geometric parallax. Using the Very Long Baseline Array, we have observed two binary millisecond pulsars, PSR J1022+1001 and J2145-0750, over a two-year period and measured their distances to be 700 +14 -10 pc and 613 +16 -14 pc respectively. We use the well-calibrated distance in conjunction with revised analysis of optical photometry to tightly constrain the nature of their massive (M ~ 0.85 Msun) white dwarf companions. Finally, we show that several measurements of their parallax and proper motion of PSR J1022+1001 and PSR J2145-0750 obtained by pulsar timing array projects are incorrect, differing from the more precise VLBI values by up to 5 sigma. We investigate possible causes for the discrepancy, and find that imperfect modeling of the solar wind is a likely candidate for the timing model errors given the low ecliptic latitude of these two pulsars.
We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < {\delta} < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715.