Scintillation arcs in low-frequency observations of the timing-array millisecond pulsar PSR J04374715
Low-frequency observations of pulsars provide a powerful means for probing the microstructure in the turbulent interstellar medium (ISM). Here we report on high-resolution dynamic spectral analysis of our observations of the timing-array millisecond pulsar PSR J04374715 with the Murchison Widefield Array (MWA), enabled by our recently commissioned tied-array beam processing pipeline for voltage data recorded from the high time resolution mode of the MWA. A secondary spectral analysis reveals faint parabolic arcs, akin to those seen in high-frequency observations of pulsars with the Green Bank and Arecibo telescopes. Data from Parkes observations at a higher frequency of 732 MHz reveal a similar parabolic feature, with a curvature that scales approximately as the square of the observing wavelength () to the MWA’s frequency of 192 MHz. Our analysis suggests that scattering toward PSR J04374715 predominantly arises from a compact region about 115 pc from the Earth, which matches well with the expected location of the edge of the Local Bubble that envelopes the local Solar neighbourhood. As well as demonstrating new and improved pulsar science capabilities of the MWA, our analysis underscores the potential of low-frequency pulsar observations for gaining valuable insights into the local ISM and for characterising the ISM toward timing-array pulsars.
Subject headings:pulsars: general — pulsars: individual (PSR J0437-4715) — methods: observational — instrumentation: interferometers
Pulsar signals are subjected to a range of delays, distortions and amplitude modulations from dispersive and scattering effects due to the ionised interstellar medium (ISM).
These propagation effects scale steeply with the observing frequency, making low frequencies ( 400 MHz) less appealing for high-precision timing experiments such as pulsar timing arrays (PTAs; Manchester et al., 2013; Demorest et al., 2013; van Haasteren et al., 2011).
However, with PTA experiments approaching timing precisions of 0.1-0.8 s, there is renewed interest in understanding
the local ISM and its effects on high-precision timing, which may potentially limit achievable timing precision for most PTA pulsars (Arzoumanian et al., 2015; Lentati et al., 2015; Shannon et al., 2015).
While the effects of temporal variations in dispersion measure (DM) have been investigated to a certain extent
(You et al., 2007; Cordes & Shannon, 2010; Keith et al., 2013; Lee et al., 2014; Lam et al., 2015; Cordes et al., 2015), there exists only a limited understanding of the impact of scattering on timing precision.
PTA experiments currently rely on millisecond pulsars (MSPs) with low to moderate DMs ( 50 ) in order to minimise ISM effects on timing precision. Recent work that used the Parkes pulsar timing array (PPTA) data to place a limit on the strength of the stochastic gravitational-wave background (Shannon et al., 2015) advocates shorter-wavelength ( 3 GHz) observations to alleviate ISM effects. However, this is not currently feasible for the majority of PTA pulsars, for which 1-2 GHz remains the most practically viable choice due to sensitivity limitations of existing telescopes and instrumentation. The ISM effects, including multi-path scattering, may still be significant at those frequencies. Since scattering delays ( ) scale steeply with the frequency (, where is the observing frequency; Bhat et al., 2004), they are more readily measurable in observations with new low-frequency arrays such as the MWA (Tingay et al., 2013), the Long Wavelength Array (LWA; Taylor et al., 2012) and Low Frequency Array (LOFAR; van Haarlem et al., 2013). Early pulsar observations with these instruments already demonstrate their potential in this direction (Bhat et al., 2014; Dowell et al., 2013; Archibald et al., 2014).
Observations of “scintillation arcs” – faint, parabolic arc-like features seen in secondary spectral analysis of pulsar observations have provided new insights into both the micro-structure of the ISM and the interstellar scattering phenomenon (Cordes et al., 2006; Rickett, 2007; Stinebring, 2007). First recognised by Stinebring et al. (2001) in Arecibo data, detailed studies of these arcs have revealed a variety and richness in their observational manifestations; e.g. forward and reverse arcs, and a chain of arclets (Putney & Stinebring, 2006; Hill et al., 2005; Brisken et al., 2010), which also stimulated a great deal of theoretical and modelling work (Cordes et al., 2006; Walker et al., 2004; Brisken et al., 2010). In the context of PTAs, Hemberger & Stinebring (2008) measured scattering delays from their observations of scintillation arcs in PSR B1737+13 (DM=48.9 ); the delays varied from 0.2 to 2 s at a frequency of 1.3 GHz.
In this paper we present our observations of parabolic scintillation arcs in PSR J04374715, a high-priority target for PTAs. Details on processing and analysis are described in § 2 and § 3, while in §4 and §5 we describe the estimation of the arc curvature and the placement of the scattering screen. Our conclusions and future prospects are summarised in §6.
2. Observational data
The MWA data used in this paper were recorded with the voltage capture system (VCS) developed for the MWA. The VCS functionality allows recording up to 241.28 MHz from all 128 tiles (both polarisations) at the native 100-s, 10-kHz resolutions for up to 1.5 hr (Tremblay et al., 2015), and is the primary observing mode for observations of pulsars and fast transients. Only half the recording capability (121.28 MHz) was available in the early days of VCS commissioning when our observations were performed (MJD=56559). Further details are described in Bhat et al. (2014). These data have now been reprocessed using our new beamformer pipeline that coherently combines voltage signals from all 128 tiles. The Parkes data used are from observations made at a frequency of 732 MHz, at an epoch two weeks later (MJD=56573) than our MWA observations.
2.1. Tied-array beam processing of MWA observations
A tied-array beam is a coherent sum of voltage signals from individual tiles and is, theoretically, expected to yield a sensitivity improvement of over the incoherent addition of detected powers, where is the
number of tiles. For the MWA, this means potentially an order of magnitude boost in sensitivity, besides enabling high-time resolution polarimetric studies. It involves incorporating delay models to account for the geometric and cable lengths, as well as calibration for complex gains (amplitude and phase) of individual tiles, and proper accounting for the beam models that incorporate polarimetric response of individual tiles. The calibration and beam information are provided by an offline version of the real-time calibration and imaging system, RTS (Mitchell et al. in prep), which uses the visibilities generated from an off-line version of the MWA correlator (Ord et al., 2015). Calibration was performed using Pictor A, a bright in-beam source at 10 offset from PSR J04374715. The full processing pipeline runs on the Galaxy cluster of the Pawsey supercomputing facility
Fig. 1 shows the improvement in signal-to-noise from this tied-array beam processing for PSR J04374715 observations. The signal-to-noise ratio (S/N) of the integrated pulse profile has increased from 205 (incoherent detection) to 2100, i.e.
an improvement of a factor of 10, which is only 10% less than the theoretical expectation for a coherent sum from 126 tiles
2.2. High-resolution dynamic spectra
The beam-formed data are processed using the dspsr software package (van Straten & Bailes, 2011) to generate synchronously-folded pulse profiles over 10-second sub-integrations. With the sensitivity improvement provided by the tied-array beam processing, we were able to generate a dynamic spectrum at the native frequency resolution of 10 kHz for VCS – this is a dramatic improvement over our earlier analysis that used a spectral resolution of 640 kHz (Bhat et al., 2014). The resultant dynamic spectrum is shown in Fig. 2. It is dominated by a small number of bright scintles whose intensity maxima drift in the time-frequency plane – a consequence of refraction through the ISM. The increased sensitivity and spectral resolution results in higher sensitivity to subtle features, such as those caused by a combination of diffractive and refractive scattering effects.
Scintillation parameters including the characteristic scales in time and frequency (i.e. the scintillation bandwidth and the diffractive time scale ), can be obtained from a two-dimensional auto-correlation function analysis of dynamic spectrum. The results from such an analysis, along with a detailed comparison with the published measurements, is presented in our earlier paper (Bhat et al., 2014). For the data in Fig. 2, we obtain 1.7 MHz, 260 s, and a drift rate (in the time-frequency plane) 95 (with measurement uncertainties 25%). Our measured is discrepant with the majority of the published values, however it agrees with the larger scale of scintillation from Gwinn et al. (2006). Considering that all published measurements are from observations made at higher observing frequencies (300-600 MHz), it is possible many of them were underestimated, particularly when the observing bandwidth () is not large enough to allow reliable measurements (e.g., ).
3. Secondary spectral analysis
The dynamic spectrum is a record of the pulse intensity as a function of time and frequency, , and is the primary observable for scintillation analysis. Its two-dimensional power spectrum is the secondary spectrum, (where indicates two-dimensional Fourier transform). It is a powerful technique that captures interference patterns produced by different points in the image plane (e.g. Cordes et al., 2006; Stinebring et al., 2001). The scatter-broadened pulsar image is seen over a field of view , where is the effective distance to the screen and is the half-width angular size of the broadened pulsar image. If and are two arbitrary points in the image plane, the corresponding “fringe rates” in time and frequency, and , are given by,
where is the speed of light, is the observing wavelength, and is the fractional distance of the screen from the source; is a measure of the differential time delay between pairs of rays and is the temporal fringe frequency. Interference between the origin and pairs along an axis in the direction of net velocity vector produces parabolic scintillation arcs, represented by . In essence, parabolic arcs can be described as a natural consequence of small-angle forward scattering.
3.1. MWA observations at 192 MHz
Figure 3 shows the secondary spectrum from MWA observations at 192 MHz.
There is a clear, albeit faint, arc-like feature in the data, particularly on the left side of the zero Doppler frequency () axis. The feature is relatively more prominent in the lower one-third of the frequency band (170–175.12 MHz), however it is still visible, with somewhat reduced strength, in data over a larger, or the full frequency range.
Even though the full secondary spectrum spans fringe rates out to 50 mHz and 50 s for our dynamic spectral resolutions of
=10 kHz and =10 s in Fig. 2, the arc feature visible is largely restricted to a small region
( 5%) near the origin ( 5 mHz; 15s).
3.2. Parkes observations at 732 MHz
To confirm the parabolic arc seen in MWA data, we analysed archival Parkes data from observations at 732 MHz, the only PPTA observing frequency that is below the expected transition frequency (1 GHz) for this pulsar. Fig. 3 shows the dynamic spectrum from observations made at an epoch two weeks later than MWA observations. The data were recorded with the ATNF Parkes Digital Filterbank (DFB4), and pre-processed over 1-minute sub-integrations and have a spectral resolution of 125 kHz. A parabolic arc feature is clearly visible in the secondary spectrum of these data (Fig. 3). There are also hints of a “filled parabola,” as seen in some of the published data at lower frequencies (e.g. Stinebring, 2007).
4. Arc curvature and the placement of scattering screen
Theoretical treatments on parabolic arcs and the scattering geometry are discussed by Stinebring et al. (2001) and Cordes et al. (2006) (see also Brisken et al., 2010). The fringe frequency and delay parameters and can be related to the curvature of the arc (), the pulsar distance () and proper motion ( ), and the placement of scattering screen (). As discussed by Cordes et al. (2006), this relation depends on the number of arcs seen in observations and the scattering geometry, and in general, can be determined to within a pair of solutions that is symmetric about . For a screen located at a fractional distance from the pulsar, the curvature parameter is given by
where is the effective distance to the screen and is the angle between the net velocity vector and the orientation of the scattered image. The effective velocity is the velocity of the point in the screen intersected by a straight line from the pulsar to the observer, which is the weighted sum of the pulsar’s binary and proper motions, and the motion of the screen and the observer ( and respectively). Its transverse component is given by
where is the transverse pulsar motion (i.e. proper motion) and is the pulsar’s binary orbital motion (transverse component).
Thus, the measurement of can be used to determine the location of the scatterer, when all the contributing terms of Eq. 3 (and hence the net ) are precisely known.
The astrometric and binary orbital parameters are very well determined for PSR J04374715. Specifically, both the distance () and the proper motion ( ) are known at very high precisions from timing and interferometric observations (Verbiest et al., 2008; Deller et al., 2008); a parallax measurement of yields pc, which, when combined with the proper motion measurement of implies a transverse space motion . The three-dimensional sky geometry of the pulsar’s binary orbit is also well determined (van Straten et al., 2001; Verbiest et al., 2008), including the longitude of the ascending node . The screen velocity ( ) and the orientation angle () are generally unknown; ignoring these terms (i.e. assuming is small compared to all other terms in Eq. 3 and ), leads to a simplified form for Eq. 2, with the screen location being the sole unknown.
4.1. Estimation of the arc curvature
When observations are not limited by signal-to-noise (e.g. Arecibo data) and parabolic arcs are sharp and clearly visible, the curvature can be reliably estimated even as a best fit-by-eye (Stinebring et al., 2001; Hill et al., 2005). As the signal-to-noise of our detections are comparatively lower, we adopt a new technique that is more systematic and robust. It employs feature extraction via the application of the (one-dimensional) generalized Hough transform. We essentially parameterise the parabolic feature by the curvature parameter , and trial over a wide range of values within the constraints allowed by the data, each time summing the power along the parabolic segments outside the regions of low-frequency noise in order to minimise contamination from high spectral values near the origin and the zero axes. The mean power computed in this manner, (henceforth refered to as “arc strength”), is given by
where the summing procedure is performed along the points of arc outside the excluded low-frequency noise, and out to delays beyond which little power is detectable (i.e. =1…, where corresponds to =20 s for MWA data).
Fig. 4 shows a plot of the arc strength against from our analysis, where the point of maximum mean power corresponds to the best-fit curvature.
The measured values of and do not scale as per the theoretically-expected relation. The implied scaling index of (where ) is larger than the values ( to ) that were determined by Hill et al. (2003) from their observational data spanning a large frequency range (0.4 to 2.2 GHz). As we explain below, this departure from the scaling is due to the change in between the two observing epochs.
4.2. Placement of the scattering screen
In order to determine the location of the scatterer () from the measurements of , we need to compute the effective velocity at the two observing epochs (cf. Eq. 2 and 3). As the screen velocity and the angle
are generally unknown, we assume and in our analysis. For the observing epochs of MWA and Parkes data (MJD=56559.878 and MJD=56573.837, respectively), the corresponding true anomalies are and , respectively, i.e. a difference in the orbital phase of 0.43 cycles. The transverse pulsar and binary motions, when projected on to the plane of the sky (with the x axis defined along the line of nodes, positive toward the ascending node), are tabulated in Table 1.
There is a substantial change in between the two epochs due the pulsar’s binary motion alone; accounting for just this term, i.e.
, we obtain the fractional distance from the pulsar, , from MWA measurements and from Parkes measurements.
The contribution from the Earth’s orbital motion around the Sun ( ) can also be significant depending on the pulsar’s line of sight and the observing epoch, although it will be weighted down for a screen that is located closer to the pulsar (Eq. 3).
The x and y components of in the coordinate system that we have employed are tabulated in Table 1.
The change in between the two observing epochs is relatively small ( 5 ) in comparison to that from the pulsar’s binary motion. Accounting also for this term involves solving for in Eq. 2, which yields a pair of solutions, and from , and and for . However, in light of our previous, independent analysis based on scintillation measurements (Bhat et al., 2014), which yields a scintillation velocity = , i.e. , and therefore the solution suggesting a screen closer to the pulsar is favoured.
|Observation||Frequency||Arc curvature,||Screen location ()|
|(MHz)||( )||( )||( )||( )|
Note. – The x and y components when projected on to the plane of sky, with the x axis along the line of nodes (see the text for details).
Our results in terms of the screen placements derived from MWA and Parkes observations are in very good agreement, despite the fact that the observations were made at widely separated observing frequencies and not contemporaneous. The time separation of two weeks is significantly longer than the expected refractive time scale ( ), i.e. the characteristic time on which a new volume of scattering material is expected to move across the pulsar’s line of sight. It is given by , where is the frequency of observation, and are the scintillation bandwidth and diffractive time scale, respectively, both of which are measurable from dynamic spectra. For our MWA observations, 1.7 MHz and 4.5 minutes (Bhat et al., 2014), and hence 17 hr 1 day. This estimate is however not so reliable since it is based on single-epoch measurements, but even then it is unlikely may be longer than a few days at the MWA’s frequency. Our observational results may therefore suggest that the underlying scattering structure persists over multiple refractive cycles. The existence of such large scattering structures in the ISM was also suggested by past observations including those where the drift slopes and multiple imaging episodes were seen to persist over time scales of several months (e.g. Gupta et al., 1994; Rickett et al., 1997; Bhat et al., 1999).
Another subtlety pertains to the ISM volume sampled by multipath scattering, which is a strong function of the observing frequency. As discussed in § 3, the scatter broadened pulsar image has a characteristic size , and consequently the ISM sampled by MWA observations is more than two orders of magnitude larger than Parkes observations (since ). Although often ignored in observational interpretations, this can be potentially an important effect. Recent work of Cordes et al. (2015) explores this in great detail in the context of frequency-dependent (chromatic) DMs in timing-array observations. Nonetheless, it is not yet clear how this may influence scattering and scintillation observables and their scalings with the frequency. There is no compelling observational evidence in support of chromatic DMs, and a wealth of observational data on scintillation and scattering measurements are seen follow the expected frequency scaling over a large range. In particular, we note the work of Hill et al. (2003), who experimentally verified the scaling relation for scintillation arcs; their derived scaling indices range from to , and are consistent with the scaling, despite observational data spanning a large frequency range (from 0.4 to 2.2 GHz). Contemporaneous observations at multiple different frequencies, similar to those advocated by Lam et al. (2015) for improved DM corrections in PTA observations, will be useful for gaining further insights into
Aside from these subtleties, our observations of scintillations arcs are clear indications of scattering toward PSR J04374715 arising from a localized region (thin screen). The implied screen location of pc is, incidentally, consistent with the expected location of 100-120 pc to the edge of the Local Bubble (Snowden et al., 1990; Bhat et al., 1998; Cordes & Lazio, 2002; Spangler, 2009).
The possibility of the screen being closer to the pulsar was also hinted by our earlier, independent analysis, where our measured scintillation velocity ( = ) suggested a screen location of pc from the Earth (i.e. ) based on / 3 (Bhat et al., 2014). Scattering toward PSR J04374715 therefore most likely dominated by the material near the edge of the bubble.
Even as our observations of scintillation arcs suggesting a small fraction of the scattered radiation arriving at large delays, its impact on timing precision may be negligibly small for this pulsar at its timing frequencies of 1-3 GHz. This is because PSR J04374715 is a weakly-scattered pulsar, with the second lowest value for the measured strength of scattering (the wavenumber spectral coefficient, from our measurement of scintillation bandwidth). Based on MWA observations, a transition to weak scintillation can be expected near 1 GHz, and consequently scattering effects are no longer relevant at frequencies 1 GHz. However, this will not be the case for many other PTA pulsars. The DM range of PTA pulsars extends out to 300 , even though the majority of them are at DMs 50 . Since scattering delays ( ) are expected to scale as , timing perturbations 100 ns can be expected for PTA pulsars at the 1-2 GHz timing frequencies. DM variations may still be the dominant source of ISM noise in PTA data; however, scattering effects may also be important, particularly if DM corrections are to rely on observations at 1 GHz. Observations at the low frequencies with MWA, LWA and LOFAR, and eventually with SKA-LOW, can therefore prove to be very useful in assessing the importance of scattering delays and the nature of turbulent ISM toward PTA pulsars.
6. Conclusions and future work
A new processing pipeline for MWA high time resolution data enables forming a coherent combination of tile powers from recorded voltages, bringing an order of magnitude improvement in the sensitivity for pulsar observations. We have demonstrated one of its applications, through high-resolution dynamic spectral studies of PSR J04374715 from MWA observations at 192 MHz. A secondary spectral analysis reveals parabolic scintillation arcs, whose curvature scales as to Parkes observations at 732 MHz, once accounted for the change in the net effective velocity due to the pulsar’s binary orbital and the Earth’s motions. Our analysis suggests that scattering toward PSR J04374715 predominantly arises from a compact region located about 115 pc from the Earth, which is comparable to the distance to the edge of the Local Bubble (100-120 pc) that encapsulates the local Solar neighbourhood. Dedicated observational campaigns at the low frequencies of MWA and LOFAR, preferably contemporaneously with timing-array observations made at higher frequencies, can be potentially promising for a detailed characterisation of the ISM along the lines of sight, and for assessing the sources of ISM noise in timing-array data.
Acknowledgments: We thank an anonymous referee for several insightful comments which helped improve the content and presentation of this paper. We also thank J.-P. Macquart, R. M. Shannon, M. Bailes and H. Knight for several useful discussions. This scientific work makes use of the Murchison Radio-astronomy Observatory, operated by CSIRO. We acknowledge the Wajarri Yamatji people as the traditional owners of the Observatory site. NDRB is supported by a Curtin Research Fellowship. Support for the operation of the MWA is provided by the Australian Government Department of Industry and Science and Department of Education (National Collaborative Research Infrastructure Strategy: NCRIS), under a contract to Curtin University administered by Astronomy Australia Limited. We acknowledge the iVEC Petabyte Data Store and the Initiative in Innovative Computing and the CUDA Center for Excellence sponsored by NVIDIA at Harvard University, and support from the Centre for All-sky Astrophysics (CAASTRO) funded by grant CE110001020.
- Two tiles were excluded from tied-array beam-forming owing to poor calibration solutions.
- The computation of is restricted to where the arc feature is more prominent.
- =0.97 would imply 0.03 , which is not supported by observations.
- Archibald, A. M., Kondratiev, V. I., Hessels, J. W. T., & Stinebring, D. R. 2014, ApJ, 790, L22
- Arzoumanian, Z., Brazier, A., Burke-Spolaor, S., et al. 2015, arXiv:1508.03024
- Bhat, N. D. R., Cordes, J. M., Camilo, F., Nice, D. J., & Lorimer, D. R. 2004, ApJ, 605, 759
- Bhat, N. D. R., Gupta, Y., & Rao, A. P. 1998, ApJ, 500, 262
- Bhat, N. D. R., Rao, A. P., & Gupta, Y. 1999, ApJS, 121, 483
- Bhat, N. D. R., Ord, S. M., Tremblay, S. E., et al. 2014, ApJ, 791, L32
- Brisken, W. F., Macquart, J.-P., Gao, J. J., et al. 2010, ApJ, 708, 232
- Cordes, J. M., & Lazio, T. J. W. 2002, arXiv:astro-ph/0207156
- Cordes, J. M., Rickett, B. J., Stinebring, D. R., & Coles, W. A. 2006, ApJ, 637, 346
- Cordes, J. M., & Shannon, R. M. 2010, arXiv:1010.3785
- Cordes, J. M., Shannon, R. M., & Stinebring, D. R. 2015, arXiv:1503.08491
- Deller, A. T., Verbiest, J. P. W., Tingay, S. J., & Bailes, M. 2008, ApJ, 685, L67
- Demorest, P. B., Ferdman, R. D., Gonzalez, M. E., et al. 2013, ApJ, 762, 94
- Dowell, J., Ray, P. S., Taylor, G. B., et al. 2013, ApJ, 775, L28
- Gupta, Y., Rickett, B. J., & Lyne, A. G. 1994, MNRAS, 269, 1035
- Gwinn, C. R., Hirano, C., & Boldyrev, S. 2006, A&A, 453, 595
- Hemberger, D. A., & Stinebring, D. R. 2008, ApJ, 674, L37
- Hill, A. S., Stinebring, D. R., Barnor, H. A., Berwick, D. E., & Webber, A. B. 2003, ApJ, 599, 457
- Hill, A. S., Stinebring, D. R., Asplund, C. T., et al. 2005, ApJ, 619, L171
- Keith, M. J., Coles, W., Shannon, R. M., et al. 2013, MNRAS, 429, 2161
- Lam, M. T., Cordes, J. M., Chatterjee, S., & Dolch, T. 2015, ApJ, 801, 130
- Lee, K. J., Bassa, C. G., Janssen, G. H., et al. 2014, MNRAS, 441, 2831
- Lentati, L., Taylor, S. R., Mingarelli, C. M. F., et al. 2015, MNRAS, 453, 2576
- Manchester, R. N., Hobbs, G., Bailes, M., et al. 2013, PASA, 30, 17
- Ord, S. M., Crosse, B., Emrich, D., et al. 2015, PASA, 32, e006
- Putney, M. L., & Stinebring, D. R. 2006, Chinese Journal of Astronomy and Astrophysics Supplement, 6, 233
- Rickett, B. J. 2007, SINS - Small Ionized and Neutral Structures in the Diffuse Interstellar Medium, 365, 207
- Rickett, B. J., Lyne, A. G., & Gupta, Y. 1997, MNRAS, 287, 739
- Shannon, R. M., Ravi, V., Lentati, L. T., et al. 2015, arXiv:1509.07320
- Snowden, S. L., Cox, D. P., McCammon, D., & Sanders, W. T. 1990, ApJ, 354, 211
- Spangler, S. R. 2009, Space Sci. Rev., 143, 277
- Stinebring, D. R., McLaughlin, M. A., Cordes, J. M., et al. 2001, ApJ, 549, L97
- Stinebring, D. 2007, SINS - Small Ionized and Neutral Structures in the Diffuse Interstellar Medium, 365, 254
- Taylor, G. B., Ellingson, S. W., Kassim, N. E., et al. 2012, Journal of Astronomical Instrumentation , 1, 1250004
- Tingay, S. J., Goeke, R., Bowman, J. D., et al. 2013, PASA, 30, 7
- Tremblay, S. E., Ord, S. M., Bhat, N. D. R., et al. 2015, PASA, 32, e005
- van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2
- van Haasteren, R., Levin, Y., Janssen, G. H., et al. 2011, MNRAS, 414, 3117
- van Straten, W., & Bailes, M. 2011, PASA, 28, 1
- van Straten, W., Bailes, M., Britton, M., et al. 2001, Nature, 412, 158
- Verbiest, J. P. W., Bailes, M., van Straten, W., et al. 2008, ApJ, 679, 675
- Walker, M. A., Melrose, D. B., Stinebring, D. R., & Zhang, C. M. 2004, MNRAS, 354, 43
- You, X. P., Hobbs, G., Coles, W. A., et al. 2007, MNRAS, 378, 493