A SEARCH FOR FAST RADIO BURSTS WITH THE GBNCC SURVEY PRAGYA CHAWLA - - PowerPoint PPT Presentation

a search for fast radio
SMART_READER_LITE
LIVE PREVIEW

A SEARCH FOR FAST RADIO BURSTS WITH THE GBNCC SURVEY PRAGYA CHAWLA - - PowerPoint PPT Presentation

A SEARCH FOR FAST RADIO BURSTS WITH THE GBNCC SURVEY PRAGYA CHAWLA McGill University (On Behalf of the GBNCC Collaboration) P. Chawla 1 , V. M. Kaspi 1 , A. Josephy 1 , K. M. Rajwade 2 , D. R. Lorimer 2,3 , A. M. Archibald 4 , M. E. DeCesar 5 ,


slide-1
SLIDE 1

A SEARCH FOR FAST RADIO BURSTS WITH THE GBNCC SURVEY

PRAGYA CHAWLA McGill University (On Behalf of the GBNCC Collaboration)

slide-2
SLIDE 2

2 Pragya Chawla, McGill University December 9, 2016

  • P. Chawla1, V. M. Kaspi1, A. Josephy1, K. M. Rajwade2, D. R.

Lorimer2,3, A. M. Archibald4, M. E. DeCesar5, J. W. T. Hessels4,6,

  • D. L. Kaplan7, C. Karako-Argaman1, V. I. Kondratiev4,8, L. Levin9,
  • R. S. Lynch3, M. A. McLaughlin2, S. M. Ransom10, M. S. E.

Roberts11, I. H. Stairs12, K. Stovall13, J. K. Swiggum7 and J. van Leeuwen4,6

1Department of Physics & McGill Space Institute, McGill University 2Department of Physics and Astronomy, West Virginia University 3National Radio Astronomy Observatory, Green Bank 4ASTRON, the Netherlands Institute for Radio Astronomy 5Department of Physics, Lafayette College 6Anton Pannekoek Institute for Astronomy, University of Amsterdam 7Department of Physics, University of Wisconsin-Milwaukee 8Astro Space Center, Lebedev Physical Institute, Russian Academy of Sciences 9Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, The University of Manchester 10National Radio Astronomy Observatory, Charlottesville 11New York University, Abu Dhabi 12Department of Physics and Astronomy, University of British Columbia 13National Radio Astronomy Observatory, Socorro

2017, ApJ, Submitted

slide-3
SLIDE 3

3 Pragya Chawla, McGill University December 9, 2016

FAST RADIO BURSTS

  • Millisecond-duration events
  • 1.4 DMMW < DM < 35 DMMW
  • Origin unknown, likely

extragalactic

  • 18 FRBs discovered till date in

the freq. range of 700 – 1500 MHz.

  • Scattering, spectral index and

free-free absorption could be possible reasons for non- detection at low frequencies.

DM = ׬

𝟏 𝒆 𝒐𝒇 ⅆ𝒎

𝐮𝐭𝐝𝐛𝐮𝐮 ∝ 𝛏−𝟓.𝟓

Thornton et al. 2013

slide-4
SLIDE 4

4 Pragya Chawla, McGill University December 9, 2016

FAST RADIO BURSTS

  • Millisecond-duration events
  • 1.4 DMMW < DM < 35 DMMW
  • Origin unknown, likely

extragalactic

  • 18 FRBs discovered till date in

the freq. range of 700 – 1500 MHz.

  • Scattering, spectral index and

free-free absorption could be possible reasons for non- detection at low frequencies.

DMMW Contours Based on NE2001 Model (Cordes & Lazio 2002)

slide-5
SLIDE 5

5 Pragya Chawla, McGill University December 9, 2016

FAST RADIO BURSTS

  • Millisecond-duration events
  • 1.4 DMMW < DM < 35 DMMW
  • Origin unknown, likely

extragalactic

  • 18 FRBs discovered till date in

the freq. range of 700 – 1500 MHz.

  • Scattering, spectral index and

free-free absorption could be possible reasons for non- detection at low frequencies.

slide-6
SLIDE 6

6 Pragya Chawla, McGill University December 9, 2016

THE GBNCC PULSAR SURVEY

  • Operates at 350 MHz with a

bandwidth of 100 MHz.

  • The analysis pipeline (Stovall et al.

2014) based on PRESTO (Ransom 2001) includes:

  • RFI Removal
  • Dedispersion
  • Single Pulse Search
  • Dispersion delays are corrected

for by dedispersing the data at a large number of trial DMs.

  • Single pulses are searched for by

convolving the dedispersed time series with boxcars of varying duration.

An Example Plot Generated by the Analysis Pipeline

slide-7
SLIDE 7

7 Pragya Chawla, McGill University December 9, 2016

THE GBNCC PULSAR SURVEY

  • Maximum DM for FRB search =

3000 pc cm-3 (for 45000 pointings).

  • 17000 pointings searched to a

DM of 500 pc cm-3 have maximum Galactic DM along LOS < 100 pc cm-3.

  • Total observing time for all

pointings = 84 days.

GBNCC Pointings Searched to DM = 3000 pc cm-3

slide-8
SLIDE 8

8 Pragya Chawla, McGill University December 9, 2016

THE GBNCC PULSAR SURVEY

  • Maximum DM for FRB search =

3000 pc cm-3 (for 45000 pointings).

  • 17000 pointings searched to a

DM of 500 pc cm-3 have maximum Galactic DM along LOS < 100 pc cm-3.

  • Total observing time for all

pointings = 84 days.

GBNCC Pointings Searched to DM = 500 pc cm-3

slide-9
SLIDE 9

9 Pragya Chawla, McGill University December 9, 2016

THE GBNCC PULSAR SURVEY

  • Selected pointings with at least
  • ne pulse with S/N > 10.
  • Single pulse events at

DM > 2 DMMW were processed with the grouping and rating algorithm, RRATtrap.

  • RRATtrap (Karako-Argaman et
  • al. 2015) was developed to

detect Rotating RAdio Transients (RRATs).

  • Enabled discovery of 10 new

RRATs in GBNCC data.

  • RRATtrap output was inspected

and no FRBs were detected.

All GBNCC Pointings Searched for FRBs

slide-10
SLIDE 10

10 Pragya Chawla, McGill University December 9, 2016

CALCULATION OF FRB RATE

  • Assuming Poisson statistics, field of view = 0.41 sq. deg for the GBT beam and

a threshold flux density of 0.63 Jy for a pulse of intrinsic width = 5 ms.

  • We place a 95% confidence upper limit on FRB rate = 3.6 x 103 FRBs sky-1 day-1

at 350 MHz.

slide-11
SLIDE 11

11 Pragya Chawla, McGill University December 9, 2016

CONSTRAINING SPECTRAL INDICES OF FRBs

  • Assuming a power-law flux density model for FRBs:
  • Varying the index γ of the log N–log S function of the FRB population

(γ = 0.8, 1.2 and 1.5) such that:

  • Performing Monte Carlo simulations of FRB flux density distribution consistent

with the 1.4-GHz rate estimate reported for the Parkes surveys by Crawford et

  • al. (2016) (3.3 x 103 FRBs sky-1 day-1)
slide-12
SLIDE 12

12 Pragya Chawla, McGill University December 9, 2016

ABSENCE OF SCATTERING AND FREE-FREE ABSORPTION

γ αlim

0.8 +0.19 1.2 +0.28 1.5 +0.35

For constraining spectral index, αlim, computed 350-MHz FRB rate = 95% confidence GBNCC upper limit

Generate 1.4-GHz flux density distribution (γ = 0.8, 1.2 or 1.5) Scale distribution to 350 MHz by sampling α from a normal distribution (σ = 0.5) Compute 350-MHz FRB rate above S = 0.63 Jy from the resulting distribution

slide-13
SLIDE 13

13 Generate 1.4-GHz flux density distribution Scale distribution to 350 MHz by sampling α from a normal distribution Assign tscatt drawn from a lognormal distribution to each detectable FRB (S > 0.65 Jy) Is tscatt for all FRBs > 100 ms? Increase mean of lognormal dist. FALSE Mean scattering time for spectral index α = mean of lognormal dist. TRUE

γ αlim 0.8

  • 0.9

1.2

  • 0.6

1.5

  • 0.3

Given the observed range of scattering times and FRB rate, R, at 1.4 GHz,

SCATTERING

slide-14
SLIDE 14

14 Pragya Chawla, McGill University December 9, 2016

PREDICTING FRB DETECTION RATES FOR OTHER SURVEYS

  • Based on flux

density distribution consistent with 1.4- GHz rate estimate for Parkes surveys.

  • Distribution scaled

to each survey’s center freq. by sampling α from a normal distribution.

  • Scattering time for

each FRB drawn from a lognormal distribution.

slide-15
SLIDE 15

15 Pragya Chawla, McGill University December 9, 2016

PREDICTING FRB DETECTION RATES FOR OTHER SURVEYS

  • Based on flux

density distribution consistent with 1.4- GHz rate estimate for Parkes surveys.

  • Distribution scaled

to each survey’s center freq. by sampling α from a normal distribution.

  • Scattering time for

each FRB drawn from a lognormal distribution.

slide-16
SLIDE 16

16 Pragya Chawla, McGill University December 9, 2016

PREDICTING FRB DETECTION RATES FOR OTHER SURVEYS

  • Limits on bursts per

hour arise because

  • f range of spectral

indices considered (αlim < α < +2).

  • No free-free

absorption assumed.

γ = 0.8

slide-17
SLIDE 17

17 Pragya Chawla, McGill University December 9, 2016

PREDICTING FRB DETECTION RATES FOR OTHER SURVEYS

  • Predictions

consistent with upper limits reported for FRB searches with MWA (Tingay et al. 2015, Rowlinson et al. 2016), LOFAR (Coenen et al. 2014, Karastergiou et al. 2016), AO327 (Deneva et al. 2016) and UTMOST (Caleb et al. 2016).

γ = 1.2

slide-18
SLIDE 18

18 Pragya Chawla, McGill University December 9, 2016

  • CHIME bandwidth

divided into 4 parts centered at 450, 550, 650 and 750 MHz.

  • For γ = 1.5,

simulations predict detection of 1-60 bursts a day for the freq. range

  • f 400-500 MHz.
  • For 700-800 MHz,

simulations predict 5-42 bursts a day.

PREDICTING FRB DETECTION RATES FOR CHIME

γ = 1.5

slide-19
SLIDE 19

19 Pragya Chawla, McGill University December 9, 2016

SUMMARY

  • No FRBs above a S/N of 10 were detected in GBNCC pointings

amounting to a total observing time of 84 days.

  • We place a 95% confidence upper limit on the FRB rate of 3.6 x 103

FRBs sky-1 day-1 above 0.63 Jy at 350 MHz.

  • Non-detection with GBNCC is consistent with 1.4-GHz rate estimate

for α > +0.35 in the absence of scattering and free-free absorption and α > -0.3 in the presence of scattering, for γ = 1.5.

  • We predict CHIME to detect 3-54 bursts per day assuming the

Crawford et al. rate estimate, for γ = 1.5.