No Low-frequency Emission from Extremely Bright Fast Radio Bursts

, , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2018 October 29 © 2018. The American Astronomical Society. All rights reserved.
, , Citation M. Sokolowski et al 2018 ApJL 867 L12 DOI 10.3847/2041-8213/aae58d

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

2041-8205/867/1/L12

Abstract

We present the results of a coordinated campaign conducted with the Murchison Widefield Array (MWA) to shadow fast radio bursts (FRBs) detected by the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz, which resulted in simultaneous MWA observations of seven ASKAP FRBs. We de-dispersed the 24 × 1.28 MHz MWA images across the 170–200 MHz band taken at 0.5 s time resolution at the known dispersion measures (DMs) and arrival times of the bursts and searched both within the ASKAP error regions (typically ∼10' × 10'), and beyond (4° × 4°). We identified no candidates exceeding a 5σ threshold at these DMs in the dynamic spectra. These limits are inconsistent with the mean fluence scaling of α = −1.8 ± 0.3 (${{ \mathcal F }}_{\nu }\propto {\nu }^{\alpha }$, where ν is the observing frequency) that is reported for ASKAP events, most notably for the three high-fluence (${{ \mathcal F }}_{1.4\mathrm{GHz}}\gtrsim 100$ Jy ms) FRBs 171020, 180110, and 180324. Our limits show that pulse broadening alone cannot explain our non-detections, and that there must be a spectral turnover at frequencies above 200 MHz. We discuss and constrain parameters of three remaining plausible spectral break mechanisms: free–free absorption, intrinsic spectral turnover of the radiative processes, and magnification of signals at ASKAP frequencies by caustics or scintillation. If free–free absorption were the cause of the spectral turnover, we constrain the thickness of the absorbing medium in terms of the electron temperature, T, to <0.03 (T/104 K)−1.35 pc for FRB 171020.

Export citation and abstract BibTeX RIS

1. Introduction

The origin of the bright, millisecond-timescale emission associated with fast radio bursts (FRBs) remains an open question. Many of the fundamental observational characteristics of these bursts remain stubbornly elusive, a fact exemplified by poor constraints on even the spectral extent of the radio emission. Until very recent FRB detections by the Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB; The CHIME/FRB Collaboration et al. 2018) down to 400 MHz (Boyle & Chime/Frb Collaboration 2018), the lowest-frequency FRB was observed at 700 MHz (Masui et al. 2015). Despite major efforts (including Coenen et al. 2014; Karastergiou et al. 2015; Tingay et al. 2015; Rowlinson et al. 2016; Keane et al. 2016; Amiri et al. 2017; Chawla et al. 2017; and Burke-Spolaor et al. 2016, to name a few), to date no FRB emission has been reported below 400 MHz or above 8 GHz (Gajjar et al. 2018). Currently the only published limit on the spectral index below 700 MHz resulting from simultaneous broadband observations is α < −3 (${{ \mathcal F }}_{\nu }\propto {\nu }^{\alpha }$, where ${{ \mathcal F }}_{\nu }$ is fluence at the observing frequency ν), for the Parkes FRB 150418 (Keane et al. 2016). Moreover, there have been no coincident detections of FRBs in any other waveband, despite extensive multi-wavelength simultaneous and follow-up observations with optical, infrared (IR), X-ray, and gamma-ray facilities (Petroff et al. 2015; Keane et al. 2016; Scholz et al. 2016; Law et al. 2017).

The dearth of FRB detections at low radio frequencies (≤400 MHz) presents a critical impediment to the analysis of the burst energetics. Burst energies up to E ∼ 1035 J are inferred by integrating the emission across the observing band and assuming that the emission is isotropic (Lorimer et al. 2007; Thornton et al. 2013; Bannister et al. 2017). However, reliable spectral measurements based on the 20 FRBs reported by the Commensal Realtime ASKAP Fast Transients (CRAFT) survey (Macquart et al. 2010) on the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz show that the average FRB fluence spectrum is steep, with a spectral index α = −1.8 ± 0.3 (Shannon et al. 2018). This indicates that the low-frequency cutoff likely dominates the energetics of the radio emission, notwithstanding the fact that the bursts detected by ASKAP often exhibit patchy spectral structure (Shannon et al. 2018), as does the repeating FRB 121102 (Spitler et al. 2016). High fluences and steep spectral indices of CRAFT FRBs make them ideal targets for low-frequency observations with the MWA.

Low-frequency measurements also provide diagnostics of the plasma along the line of sight to the FRB. In particular, a number of propagation effects potentially influence the spectral characteristics at low frequencies, and their identification in FRB data would place constraints on the burst environment and the properties of the plasma encountered along the ray path. The two most obvious effects are scattering due to inhomogeneities in the plasma (most likely distant from the source, i.e., ≫1 pc), which causes temporal smearing of the signal proportional to ∼ν−4 (Bhat et al. 2004), and free–free absorption by dense, circumburst plasma, whose optical depth scales as ${\tau }_{\mathrm{ff}}\sim {\nu }^{-2.1}$.

The characteristics of the low-frequency emission bear heavily on FRB detection rates with future survey facilities such as the low-frequency component of the Square Kilometre Array (SKA-Low; Fender et al. 2015). The wide fields of view (FOVs) accessible by this telescope render it capable of detecting transient phenomena at extremely high rates Fender et al. 2015). Moreover, the recent CHIME/FRB detections show that FRBs can be observed at least down to 400 MHz, triggering even more interest in lower frequencies and spectral extent of FRB emission.

In this Letter we present the results of an observing campaign undertaken with the Murchison Widefield Array (MWA; Tingay et al. 2013) to detect low-frequency radio emission coincident with the bright FRBs detected by the CRAFT survey at 1.4 GHz (Macquart et al. 2010; Bannister et al. 2017; Shannon et al. 2018). The MWA shadowed (co-tracked) the pointing positions of the ASKAP antennas, so that the precise dispersion measure (DM), time of arrival, and the approximate position (typically a 10' × 10' region) of each burst are all known, greatly reducing the searched volume of parameter space relative to a blind survey.

2. Observations and Data Reduction

2.1. Observing Strategy

The observations were made with the MWA while the ASKAP antennas observed in the fly's-eye mode12 distributed along a certain Galactic latitude (see Bannister et al. 2017). The FOV of the MWA in the frequency band 170–200 MHz is approximately 450 deg2 (FWHM ∼ 21°), enabling nearly full coverage of the ASKAP fly's-eye FOV with sensitivity ≳50% of the primary beam. Moreover, the selected frequency band minimizes the effects of pulse broadening (due to scattering) and radio-frequency interference. Figure 1 shows the observing setup during FRB 180324.

Figure 1.

Figure 1. Illustration of the MWA shadowing strategy for FRB 180324 (on top of the image by Haslam et al. 1982 at 408 MHz). The large (≈7°) red circles represent coverage of ASKAP antenna beams distributed near Galactic latitude −20° and small white circles (≈1fdg2) are the individual 36 beams of ASKAP antenna 25, which detected FRB 180324 (yellow star). The green (at 0.5 and 0.1) and white-dashed (at 0.001 with angular size ≈42° × 42°) contours show the MWA beam (Sokolowski et al. 2017) normalized to maximum response at zenith.

Standard image High-resolution image

Table 1 summarizes the ASKAP FRBs (Macquart et al. 2018; Shannon et al. 2018, submitted) detected while MWA shadowed. The strategy succeeded for the first time with FRB 171020, and the MWA collected data before, during, and after the FRB detection by ASKAP (Table 1). The MWA data were collected during the transition to the extended array (Wayth et al. 2018), implying reduced sensitivity (only ∼70% of antennas). Since then, six more ASKAP FRB positions were observed by the MWA in a similar mode. In all cases MWA data were recorded in 10 kHz frequency and 0.5 s temporal resolutions.

Table 1.  Details of the ASKAP FRBs Shadowed by the MWA

                ${{ \mathcal F }}_{185\mathrm{MHz}}$ (Jy ms) Expectedf  
FRB UTC ${\mathrm{DM}}_{\mathrm{tot}}$ a ${\mathrm{DM}}_{\mathrm{mw}}$ a tarrb tsweepc ${\tau }_{\mathrm{scat}}$ d ${{ \mathcal F }}_{1.4\mathrm{GHz}}$ e $\alpha =-1$ α = −2 α = −1.8g ${{ \mathcal U }}_{5\sigma }$ h
  Detection $(\mathrm{pc}/{\mathrm{cm}}^{3})$ (s) (s) (ms) (Jy ms )       (Jy ms)
171020 10:27:59.00 114.1 38.4 11.7 4.5 1.7 ${200}_{-100}^{+500}$ ${1500}_{-800}^{+4000}$ ${11400}_{-6000}^{+30000}$ ${7600}_{-4000}^{+19000}$ 2200
180110 07:34:34.95 715.7 38.8 73.0 28.0 4.5 ${420}_{-20}^{+20}$ ${3200}_{-150}^{+150}$ ${23900}_{-1100}^{+1100}$ ${16000}_{-800}^{+800}$ 3350i
                      or
                      6500j
180128.0 00:59:37.97 441.4 31.5 45.0 17.3 2.9 ${51}_{-2}^{+2}$ ${380}_{-15}^{+15}$ ${2900}_{-110}^{+110}$ ${1940}_{-80}^{+80}$ GLk
180128.2 04:53:26.80 495.9 41.0 50.6 19.40 2.3 ${66}_{-4}^{+4}$ ${500}_{-30}^{+30}$ ${3800}_{-230}^{+230}$ ${2500}_{-150}^{+150}$ SLk
180130 04:55:29.99 343.5 39.0 34.90 13.35 6.0 ${95}_{-3}^{+3}$ ${720}_{-20}^{+20}$ ${5400}_{-170}^{+170}$ ${3600}_{-110}^{+110}$ SLk
180315 05:05:30.99 479.0 101.7 48.66 18.63 2.4 ${56}_{-4}^{+4}$ ${420}_{-30}^{+30}$ ${3200}_{-230}^{+230}$ ${2100}_{-150}^{+150}$ SLk
180324 09:31:46.70 431.0 64.0 43.79 16.75 4.3 ${71}_{-3}^{+3}$ ${540}_{-20}^{+20}$ ${4000}_{-170}^{+170}$ ${2700}_{-110}^{+110}$ 450i

Notes.

a ${\mathrm{DM}}_{\mathrm{tot}}$: the total DM measured by ASKAP; ${\mathrm{DM}}_{\mathrm{mw}}$ is the contribution of the Milky Way (from NE2001; Cordes & Lazio 2002). btarr: the time delay between ASKAP detection at 1297 MHz, and the expected arrival at 200 MHz. ctsweep: the sweep time over the MWA observing band (170–200 MHz). d ${\tau }_{\mathrm{scat}}$: scattering time at 1.4 GHz fitted to ASKAP data for 180110 and 180130 and pulse width for other FRBs. eThe errors represent 90% confidence limits. fThe fluences are extrapolated to 185.6 MHz assuming power-law scaling ${{ \mathcal F }}_{\nu }\propto {\nu }^{\alpha }$. gThe mean spectral index of ASKAP FRBs. h ${{ \mathcal U }}_{5\sigma }$: the MWA 5σ upper limit on ${{ \mathcal F }}_{185\mathrm{MHz}}$. iThese limits are higher by a factor ≈1.2–1.5 because of de-dispersion in 1.28 MHz channels. jThis FRB was significantly temporally broadened (see Section 2.3); therefore, we also present the limit from de-dispersed images averaged over 10 s (Table 2). kThe data quality was too low to obtain meaningful upper limits due to the presence of the Sun in sidelobe (SL) or grating-lobe (GL).

Download table as:  ASCIITypeset image

The shadowing program was performed mostly during daytime. Hence, only FRBs 171020 and 180324 were detected by ASKAP after sunset, while the FRB 180110 field was observed with the Sun close to a "null" of the MWA's primary beam (a direction with very low sensitivity). In the other cases the data quality was too low to derive meaningful limits due to the presence of the Sun in the sidelobe of the primary beam.

2.2. Calibration

The MWA data were calibrated with the Common Astronomy Software Applications (CASA; McMullin et al. 2007), using observations of a calibrator source (3C444 or Pictor A). We applied calibration solutions to FRB field observations and created dirty images in 0.5 s temporal (shortest possible) and 1.28 MHz frequency resolutions using wsclean (Offringa et al. 2014) and natural weighting.

We also calibrated and imaged MWA observations (292 s) collected before the FRB observations in order to create reference images of the FRB fields (typical standard deviation of noise σ ∼ 20–40 $\mathrm{mJy}\,{\mathrm{beam}}^{-1}$). The flux density scale was calibrated using sources from the GaLactic Extragalactic All-sky MWA (GLEAM) survey (Hurley-Walker et al. 2017) identified in the reference images.

2.3. FRB Detection Algorithm

The transient detection algorithm used the reference images (Section 2.2) to create a list of reference sources above certain threshold (5σ). The 24 × 1.28 MHz, 0.5 s images across the observing band were time-aligned according to the ASKAP DM and summed to provide a 0.5 s resolution de-dispersed image (example images are shown in Figure 2). We also examined the de-dispersed images before and after the expected burst arrival time (we inspected all 0.5 s images within the analyzed MWA observation). The details of candidates exceeding a 5σ threshold, which were not present in the list of reference sources, were saved for further visual inspection. We did not observe any >5σ event within the ASKAP error boxes.13

Figure 2.

Figure 2. Example images of the FRB 180324 position (image centers). Left panel: reference image of the field (the mean of six observations without the FRB), with noise in the center σ ≈ 17 $\mathrm{mJy}\,{\mathrm{beam}}^{-1}$; the green circles are NRAO Very Large Array Sky Survey (NVSS; Condon et al. 1998) sources brighter than 10 mJy at 1.4 GHz. The image scale spans −0.2–0.2 $\mathrm{Jy}\,{\mathrm{beam}}^{-1}$. Center panel: an example 0.5 s dirty image obtained by stacking the appropriately time-shifted 24 × 1.28 MHz spectral channels (i.e., de-dispersing) to DM = 431.0 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$; the noise in the image center is σ ≈ 180 $\mathrm{mJy}\,{\mathrm{beam}}^{-1}$ and the green circles are NVSS sources brighter than 100 mJy at 1.4 GHz. Right panel: an example single 0.5 s dirty image in a single 1.28 MHz channel; noise in the image center is σ ≈ 760 $\mathrm{mJy}\,{\mathrm{beam}}^{-1}$ and the green circles are NVSS sources brighter than 200 mJy at 1.4 GHz. The scales in the center and right images are as shown in the colorbar (in $\mathrm{Jy}\,{\mathrm{beam}}^{-1}$).

Standard image High-resolution image

The 5σ transient candidates identified within a larger (4° × 4°) field were visually inspected on de-dispersed images and dynamic spectra of the candidate pixels, and none of them showed any signs of dispersion sweep in the dynamic spectra. As a final check, we visually inspected all 0.5 s, 1.28 MHz images.

FRB 180110 was significantly scattered (Table 1), with an expected pulse width at the MWA's frequency, ${\tau }_{\mathrm{MWA}}\approx 5.4$ s (Section 3). Therefore, in the second part of the algorithm we averaged over multiple de-dispersed (and non-de-dispersed) images on timescales of 5, 10, 15, 20, 25, and 30 s. Yet, we did not identify any transient events exceeding the 5σ threshold with any signs of dispersion sweep in the dynamic spectra of the candidate pixels. These time-averaged images were also visually inspected. The algorithm was executed on Stokes I and V images because noise in 0.5 s V images was slightly lower (∼20%–30%). The above procedure resulted in upper limits on the flux densities and consequently fluences of low-frequency counterparts of three ASKAP FRBs, which are summarized in Tables 1 and 2.

Table 2.  MWA 1σ Noise of Flux Density ($\mathrm{mJy}\,{\mathrm{beam}}^{-1}$) in the Centers of Beam-corrected Images on Average of 24 De-dispersed 1.28 MHz Channels (${\sigma }_{{dd}}$), Single 0.5 s Images (${\sigma }_{0.5s}$) and Averaged over 5 (${\sigma }_{5s}$), 10 (${\sigma }_{10s}$), 20 (${\sigma }_{20s}$), and 30 s (${\sigma }_{30s}$)

FRB ${\sigma }_{{dd}}$ ${\sigma }_{0.5s}$ ${\sigma }_{5s}$ ${\sigma }_{10s}$ ${\sigma }_{15s}$ ${\sigma }_{20s}$
20171020 880 4570 400 350 330 310
20180110 1340 6640 700 650 620 600
20180324 180 840 120 120 110 110

Download table as:  ASCIITypeset image

2.4. Verification of the FRB Pipeline

The algorithm was verified using a relatively high-DM pulsar (147.29 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$), PSR J0837-4135, with a pulse period P ≈ 0.751 s (Manchester et al. 2005). The pipeline detected three individual bright pulses above 5σ in a 112 s (≈149 pulses) observation. Example images with and without a pulse, the dynamic spectrum, and the lightcurve observed in the de-dispersed images are shown in Figure 3.

Figure 3.

Figure 3. Example of de-dispersed images of PSR J0837-4135 (left panel) without a pulse, and (center panel) with a pulse detected by the algorithm. Right panel: the corresponding dynamic spectrum of the central pixel showing the frequency sweep due to DM = 147.29 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$ (right panel, upper area). The three sweeps (sweep time ≈5.8 s in the 170–200 MHz range) can be seen at times ranges 97–105, 160–166, and 183–191 s between the red arrows and the data points surrounded by the red diamonds, marking the three brightest points in the lightcurve from the de-dispersed images (right panel, lower area).

Standard image High-resolution image

The efficiency of the algorithm was verified on a 112 s observation of a lower-DM pulsar (34.425 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$), PSR J0630-2834, with P ≈ 1.244 s (Manchester et al. 2005). In the 112 s (90 pulses) observation, the algorithm identified 31 pulses above the 5σ threshold in the de-dispersed images (110 pulses with the 3σ threshold). Due to a lack of pulsars with a suitable combination of DM and scattering time, we have verified the algorithm on averaged images by "injecting" a simulated FRB signal into 0.5 s MWA images and confirming that it was detected in the averaged images. These tests confirmed that the algorithm is capable of detecting highly dispersed, FRB-like transients in single and averaged 0.5 s MWA images.

3. Discussion

The absence of low-frequency emission coincident with any of the shadowed CRAFT FRBs places constraints on the emission characteristics of these bursts and their environments. We provide a brief discussion here in the context of the four most obvious interpretations.

(i) Temporal smearing due to scattering through an inhomogeneous plasma potentially dominates over the intrinsic pulse width at low frequencies and decreases the burst detectability. An upper limit to the low-frequency pulse width is inferred by attributing the ASKAP-derived pulse width entirely to scattering, and assuming that pulse broadening scales as ν−3.5 (or ν−4).14 The measured pulse widths of FRBs 171020, 180110, and 180324 at 1400 MHz in turn imply limits on the pulse width at 185 MHz of 2.1 (5.6), 5.4 (14.8), and 5.1 (14) s, respectively. Because pulse fluence is unaltered by scattering, the fluence upper limits relative to the expected fluence (Tables 1 and 2) show that any pulse with a spectrum steeper than α ≈ −1 would have been detected. We conclude that pulse broadening alone is insufficient to explain the MWA non-detections because the spectral index of the emission would have to flatten substantially between 185 MHz and 1295 MHz from the mean value of α = −1.8 measured for the ASKAP FRB population (Shannon et al. 2018).

(ii) The viability of the hypothesis that free–free absorption is responsible for the absence of low-frequency emission may be cast in terms of a limit on the thickness of the absorbing material. An upper limit to the depth of the absorbing region is derived by assuming that all plasma not attributable to the Milky Way is confined to the absorbing medium (not to the intergalactic medium (IGM), thereby assuming the redshift term (1 + z) in Equation (1) is ≈1), whose emission measure is then $\mathrm{EM}=\int {N}_{e}^{2}{dl}={\mathrm{DM}}_{\mathrm{extra}}^{2}\,{\rm{\Delta }}{L}^{-1}\,$ pc cm−6, where ΔL is the region thickness in pc, Ne is electron density, and DMextra is the DM excess (in $\mathrm{pc}\,{\mathrm{cm}}^{-3}$). A burst whose predicted flux density is Spred but is undetected at a limit of Slim constrains the free–free optical depth to ${\tau }_{\mathrm{ff}}\gt \mathrm{ln}({S}_{\mathrm{pred}}/{S}_{\mathrm{lim}})$, and one has

Equation (1)

where T4 is the electron temperature in units of 104 K and ν185 is the observing frequency, normalized to 185 MHz. The redshift term is not included as the equation still leads to a valid upper limit. We assume the filling factor f = 1, and if f < 1 then the limit is ΔL ∝ f−1. The DM values in excess of the Milky Way contributions predicted by NE2001 (Cordes & Lazio 2002) are 75.7, 676.9, and 367.0 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$ for FRBs 171020, 180110, and 180324, respectively. We conservatively subtract a further 15 pc cm−3 for the Milky Way halo contribution (see the discussion in Shannon et al. 2018). For a burst spectral index of α = −1.8, and using 5σ limits for Slim derived from Table 2, this implies upper limits on ΔL of $3\times {10}^{-2}\,{T}_{4}^{-1.35}\,$ pc, $5\,{T}_{4}^{-1.35}$ pc (using ${{ \mathcal U }}_{5\sigma }=6500\,$ Jy ms) and $0.7\,{T}_{4}^{-1.35}$ pc for these three bursts, respectively. It is not possible to place significant constraint on τff for α > −1. The tightest constraint on τff from FRB 171020 is likely tighter because a fraction of the electron column density is associated with host galaxy interstellar medium (ISM), IGM, or both.

(iii) The flux densities of the bursts detected by ASKAP could be enhanced by plasma lensing effects due to caustics or scintillation (Cordes et al. 2017). It is improbable that any lensing effects, if present, extend over an order of magnitude in frequency, so it is unlikely that it plays a significant role both at ASKAP and low frequencies. The ratios of the measured ASKAP fluences to their 185 MHz upper limits yield a lower limit to the magnification at 1.4 GHz. Comparing against the 5σ upper limit at α = −1.8, we find magnifications μ > 3, 2, and 6 for FRBs 171020, 180110, and 180324, respectively.

Most ASKAP FRBs exhibit fine frequency structure, some with a factor >10 deviations from the mean, with a filling factor of ∼5% (Shannon et al. 2018). If this structure is due to diffractive scintillation, the decorrelation bandwidth at 185 MHz would be ∼10−3 smaller than that at 1.4 GHz (assuming ∼ν4 scaling), the spectrum would therefore appear spectrally smooth at low frequency at our resolution, and the low spectral occupancy observed at 1.4 GHz would not explain the MWA non-detections. However, even if the low spectral occupancy is an intrinsic feature of the emission, it is implausible to attribute the low-frequency non-detections for all three well-constrained FRBs to this effect alone, as the relative spectral bandwidths of ASKAP and the MWA are comparable at Δν/ν equal to 26% and 17%, respectively.

(iv) A break in the intrinsic spectrum is a natural explanation for the non-detection and would be plausible on two grounds: (a) for a spectrum with α < −1, the finite burst energy is dominated by the low-energy cutoff; and (b) many pulsars, whose emission resembles the properties of FRBs in some respects, also exhibit spectral breaks in the region 100–300 MHz (see Figure 7 of Bilous et al. (2016) or Figure 6 of Murphy et al. (2017) for a recent summary). The limit on the spectral cutoff of ${\nu }_{\mathrm{lo}}\gt 200$ MHz for α = −1.8 implies that the total burst energy cannot exceed 18 times more than the pulse energy measured within the ASKAP band alone.

The presented simultaneous MWA observations of extremely bright ASKAP FRBs are a substantial advance over previous surveys conducted below 200 MHz. Therefore, the derived limits enable us to place tighter constraints on the physics. Indeed, no previous low-frequency survey has a sufficiently large exposure to have detected the low-frequency counterpart to the FRBs detected by CRAFT: the 37 ± 8 events day−1 sky−1 burst detection rate at ${{ \mathcal F }}_{1.4\mathrm{GHz}}\gt 26\,$ Jy ms is equivalent to one event every ${27}_{-5}^{+7}\times {10}^{3}$ ${\deg }^{2}\,\mathrm{hr}$  (Shannon et al. 2018). This substantially exceeds the 14 × 103 ${\deg }^{2}\,\mathrm{hr}$  exposure of the Low-Frequency Array (LOFAR) Pilot Pulsar Survey at 140 MHz (Coenen et al. 2014) with a fluence cutoff ${{ \mathcal F }}_{{\rm{c}}}=2.75\,\times {10}^{3}{\rm{\Delta }}{T}^{1/2}\,\mathrm{Jy}\,\mathrm{ms},$ for pulse widths ΔT ≤ 1.26 s in the range DM = 2–3000 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$. The ARTEMIS survey at 150 MHz (Karastergiou et al. 2015) examined DM ≤ 310 $\mathrm{pc}\,{\mathrm{cm}}^{-3}$, so would have missed six of the seven FRBs reported here. Moreover, it only searched for pulse durations below 21 ms, with a 10σ threshold ${{ \mathcal F }}_{{\rm{c}}}=4.47\times {10}^{3}{\rm{\Delta }}{T}^{1/2}\,\mathrm{Jy}\,\mathrm{ms}$, where ΔT is burst width in seconds.

The blind pilot survey of Tingay et al. (2015) had an exposure of only 4.2 × 103 ${\deg }^{2}\,\mathrm{hr}$  with a 7σ detection limit in a de-dispersed time series of 350 $\mathrm{mJy}\,{\mathrm{beam}}^{-1}$ on a timescale of 2 s (i.e., a limiting fluence of 700 Jy ms). The survey of Rowlinson et al. (2016) was sensitive to an event rate a factor of 8.5 lower than that of Tingay et al. (2015), but only resolved down to a time resolution of 28 s with a limiting fluence of 7980 Jy ms.

Our results exceed limits from previous MWA observations. The 3σ upper limit on the fluence of FRB 150418 at 185 MHz, based on shadowing of the Parkes radio telescope, is 1050 Jy ms (Keane et al. 2018). However, the Parkes burst fluence of ${2.0}_{-0.8}^{+1.2}\,$ Jy ms permits only a poor constraint on the spectral index between 1382 MHz and 185 MHz of α > −3.

4. Conclusions

We have used the MWA to attempt the co-detection of the low-frequency counterparts of the high-fluence FRBs detected by the CRAFT survey. None were found. The shadowing strategy employed here enabled us to significantly reduce the searched FRB parameter volume in space, DM, and time. We are able to derive strong upper limits based on the rare, very high ASKAP fluence (${{ \mathcal F }}_{1.4\mathrm{GHz}}\gtrsim 26$ Jy ms) FRBs, which would require "blind" surveys of enormous exposures (∼27000 deg2 hr per event) to obtain comparable limits.

No low-frequency counterparts of these FRBs were identified. The fluences of CRAFT FRBs are accurately measured, and the presented limits on broadband spectral indexes (α ≳ −1) supersede the only direct broadband limit existing in the literature (α > −3) obtained with the MWA and Parkes FRB 150418. These limits are strongly at variance with the mean α = −1.8 ± 0.3 spectral index for these FRBs measured within the ASKAP band.

We show that pulse broadening alone due to scattering is insufficient to explain the non-detections of the shadowed FRBs. We place limits on the thickness of a dense medium close to the progenitors under the assumption that free–free absorption is responsible for the spectral turnover; the tightest constraint is for FRB 171020, for which the medium thickness is $\lt 0.03\,{T}_{4}^{-1.35}$ pc. In this case the implied absorption region is so small that it casts considerable doubt that free–free absorption is responsible for the spectral turnover. The speed of the dense wind of any possible massive progenitor, or the shock associated with a young (older than 1 minute) supernova would easily encompass a larger region than that implied here. Moreover, (Mahony et al. 2018) examines the properties of the likely host galaxy of this FRB, and finds no obvious bright compact knots of continuum radio emission at frequencies up to 21 GHz down to 31.8 μJy.

If FRBs are enhanced by caustics in the ASKAP band, we constrain the magnification (between 1295 and 185 MHz) to be μ > 6, in the case of FRB 180324. The largely unconstrained properties of the host galaxy ISM and intervening IGM make it difficult to assess the viability of these implied lensing magnifications. However, the analysis in Macquart et al. (2018) shows that the spectral structure of the bursts at 1.4 GHz appears to be qualitatively inconsistent with the effect of plasma lensing, which disfavors lensing as the explanation of the MWA non-detections. A turnover in the intrinsic FRB spectrum, however, remains a plausible explanation.

For α ≤ −1 the low-frequency cutoff dominates the burst energy. Thus, for α = −1.8 our non-detections limit the total burst energy to <18 times the ∼1035 J values deduced for FRB detections in the ASKAP and Parkes bands alone.

To date, no FRB has been detected simultaneously at high and low frequencies. The ongoing enhancements to the MWA will improve the detection sensitivity for FRB searches by about an order of magnitude (owing to much higher temporal and spectral resolutions possible via the voltage buffer mode). These improvements will significantly increase the prospects of detecting low-frequency emission if it is much fainter than our currently set limits, which may warrant a review of FRB models and physics.

This scientific work makes use of the Murchison Radio-astronomy Observatory (MRO), operated by CSIRO. We acknowledge the Wajarri Yamatji people as the traditional owners of the Observatory site. Support for the operation of the MWA is provided by the Australian Government (NCRIS), under a contract to Curtin University administered by Astronomy Australia Limited. The ASKAP is part of the Australia Telescope National Facility which is managed by CSIRO. Operation of ASKAP is funded by the Australian Government with support from the National Collaborative Research Infrastructure Strategy. ASKAP uses the resources of the Pawsey Supercomputing Centre. Establishment of ASKAP, the MRO and the Pawsey Supercomputing Centre are initiatives of the Australian Government, with support from the Government of Western Australia and the Science and Industry Endowment Fund. J.P.M., R.M.S., and K.W.B. acknowledge Australian Research Council grant DP180100857. R.M.S. acknowledges salary support through Australian Research Council (ARC) grants FL150100148 and CE170100004. Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. Parts of this research were conducted by the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020.

Footnotes

  • 12 

    Individual antennas pointing in different directions and covering ≈300 deg2 (with 10 ASKAP antennas).

  • 13 

    The expected number of candidates exceeding 5σ due to Gaussian noise fluctuations is ≪1 for ∼100 × 100 pixels.

  • 14 

    Shannon et al. (2018) measured ${\tau }_{\mathrm{scat}}\propto {\nu }^{-3.5\pm 0.5}$ in the ASKAP band, so we examine the consequences of both a ν−3.5 and a ν−4 scaling.

Please wait… references are loading.
10.3847/2041-8213/aae58d