sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Shannon, R.

Normalized to: Shannon, R.

135 article(s) in total. 1954 co-authors, from 1 to 48 common article(s). Median position in authors list is 5,0.

[1]  oai:arXiv.org:2006.16502  [pdf] - 2129465
A population analysis of pulse broadening in ASKAP Fast Radio Bursts
Comments: 10 pages, 6 figures, Accepted to MNRAS
Submitted: 2020-06-29, last modified: 2020-07-07
The pulse morphology of fast radio bursts (FRBs) provides key information in both understanding progenitor physics and the plasma medium through which the burst propagates. We present a study of the profiles of 33 bright FRBs detected by the Australian Square Kilometre Array Pathfinder. We identify seven FRBs with measureable intrinsic pulse widths, including two FRBs that have been seen to repeat. In our modest sample we see no evidence for bimodality in the pulse width distribution. We also identify five FRBs with evidence of millisecond timescale pulse broadening caused by scattering in inhomogeneous plasma. We find no evidence for a relationship between pulse broadening and extragalactic dispersion measure. The scattering could be either caused by extreme turbulence in the host galaxy or chance propagation through foreground galaxies. With future high time resolution observations and detailed study of host galaxy properties we may be able to probe line-of-sight turbulence on gigaparsec scales.
[2]  oai:arXiv.org:2007.03128  [pdf] - 2129518
Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network
Ackley, K.; Adya, V. B.; Agrawal, P.; Altin, P.; Ashton, G.; Bailes, M.; Baltinas, E.; Barbuio, A.; Beniwal, D.; Blair, C.; Blair, D.; Bolingbroke, G. N.; Bossilkov, V.; Boublil, S. Shachar; Brown, D. D.; Burridge, B. J.; Bustillo, J. Calderon; Cameron, J.; Cao, H. Tuong; Carlin, J. B.; Casey, A.; Chang, S.; Charlton, P.; Chatterjee, C.; Chattopadhyay, D.; Chen, X.; Chi, J.; Chow, J.; Chu, Q.; Ciobanu, A.; Clarke, T.; Clearwater, P.; Cooke, J.; Coward, D.; Crisp, H.; Dattatri, R. J.; Deller, A. T.; Dobie, D. A.; Dunn, L.; Easter, P. J.; Eichholz, J.; Evans, R.; Flynn, C.; Foran, G.; Forsyth, P.; Gai, Y.; Galaudage, S.; Galloway, D. K.; Gendre, B.; Goncharov, B.; Goode, S.; Gozzard, D.; Grace, B.; Graham, A. W.; Heger, A.; Vivanco, F. Hernandez; Hirai, R.; Holland, N. A.; Holmes, Z. J.; Howard, E.; Howell, E.; Howitt, G.; Hübner, M. T.; Hurley, J.; Ingram, C.; Hamedan, V. Jaberian; Jenner, K.; Ju, L.; Kapasi, D. P.; Kaur, T.; Kijbunchoo, N.; Kovalam, M.; Choudhary, R. Kumar; Lasky, P. D.; Lau, M. Y. M.; Leung, J.; Liu, J.; Loh, K.; Mailvagan, A.; Mandel, I.; McCann, J. J.; McClelland, D. E.; McKenzie, K.; McManus, D.; McRae, T.; Melatos, A.; Meyers, P.; Middleton, H.; Miles, M. T.; Millhouse, M.; Mong, Y. Lun; Mueller, B.; Munch, J.; Musiov, J.; Muusse, S.; Nathan, R. S.; Naveh, Y.; Neijssel, C.; Neil, B.; Ng, S. W. S.; Oloworaran, V.; Ottaway, D. J.; Page, M.; Pan, J.; Pathak, M.; Payne, E.; Powell, J.; Pritchard, J.; Puckridge, E.; Raidani, A.; Rallabhandi, V.; Reardon, D.; Riley, J. A.; Roberts, L.; Romero-Shaw, I. M.; Roocke, T. J.; Rowell, G.; Sahu, N.; Sarin, N.; Sarre, L.; Sattari, H.; Schiworski, M.; Scott, S. M.; Sengar, R.; Shaddock, D.; Shannon, R.; SHI, J.; Sibley, P.; Slagmolen, B. J. J.; Slaven-Blair, T.; Smith, R. J. E.; Spollard, J.; Steed, L.; Strang, L.; Sun, H.; Sunderland, A.; Suvorova, S.; Talbot, C.; Thrane, E.; Töyrä, D.; Trahanas, P.; Vajpeyi, A.; van Heijningen, J. V.; Vargas, A. F.; Veitch, P. J.; Vigna-Gomez, A.; Wade, A.; Walker, K.; Wang, Z.; Ward, R. L.; Ward, K.; Webb, S.; Wen, L.; Wette, K.; Willcox, R.; Winterflood, J.; Wolf, C.; Wu, B.; Yap, M. Jet; You, Z.; Yu, H.; Zhang, J.; Zhang, J.; Zhao, C.; Zhu, X.
Comments: Submitted
Submitted: 2020-07-06
Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.
[3]  oai:arXiv.org:2006.14906  [pdf] - 2122755
A search for fast radio burst-like emission from Fermi gamma-ray bursts
Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal
Submitted: 2020-06-26
We report the results of the rapid follow-up observations of gamma-ray bursts (GRBs) detected by the Fermi satellite to search for associated fast radio bursts. The observations were conducted with the Australian Square Kilometre Array Pathfinder at frequencies from 1.2-1.4 GHz. A set of 20 bursts, of which four were short GRBs, were followed up with a typical latency of about one minute, for a duration of up to 11 hours after the burst. The data was searched using 4096 dispersion measure trials up to a maximum dispersion measure of 3763 pc cm$^{-3}$, and for pulse widths $w$ over a range of duration from 1.256 to 40.48 ms. No associated pulsed radio emission was observed above $26 {\rm Jy ms} (w/1 {\rm ms})^{-1/2}$ for any of the 20 GRBs.
[4]  oai:arXiv.org:2005.13159  [pdf] - 2121501
A search for supernova-like optical counterparts to ASKAP-localised Fast Radio Bursts
Comments: 7 pages, 3 figures. Accepted to Astronomy & Astrophysics on 03 June 2020
Submitted: 2020-05-27, last modified: 2020-06-24
Fast radio bursts (FRBs) are millisecond-scale radio pulses, which originate in distant galaxies and are produced by unknown sources. The mystery remains partially because of the typical difficulty in localising FRBs to host galaxies. Accurate localisations delivered by the Commensal Real-time ASKAP Fast Transients (CRAFT) survey now provide an opportunity to study the host galaxies and potential transient counterparts of FRBs at a large range of wavelengths. In this work, we investigate whether the first three FRBs accurately localised by CRAFT have supernova-like transient counterparts. We obtained two sets of imaging epochs with the Very Large Telescope for three host galaxies, one soon after the burst detection and one several months later. After subtracting these images no optical counterparts were identified in the associated FRB host galaxies, so we instead place limits on the brightness of any potential optical transients. A Monte Carlo approach, in which supernova light curves were modelled and their base properties randomised, was used to estimate the probability of a supernova associated with each FRB going undetected. We conclude that Type Ia and IIn supernovae are unlikely to accompany every apparently non-repeating FRB.
[5]  oai:arXiv.org:2002.12533  [pdf] - 2121464
First Constraints on Compact Dark Matter from Fast Radio Burst Microstructure
Comments: 9 pages, 5 Figures, 1 table, Submitted to ApJ Letters
Submitted: 2020-02-27, last modified: 2020-06-23
Despite existing constraints, it remains possible that up to $35\%$ of all dark matter is comprised of compact objects, such as the black holes in the 10-100\,M$_\odot$ range whose existence has been confirmed by LIGO. The strong gravitational lensing of transients such as FRBs and GRBs has been suggested as a more sensitive probe for compact dark matter than intensity fluctuations observed in microlensing experiments. Recently ASKAP has reported burst substructure down to $15\mu$s timescales in FRBs in the redshift range $0.3-0.5$. We investigate here the implications of this for the detectability of compact dark matter by FRBs. We find that a sample size of $\sim130$ FRBs would be required to constrain compact dark matter to less than the existing 35$\%$ limit with 95$\%$ confidence, if it were distributed along $\gtrsim 1\,$Gpc-long FRB sightlines through the cosmic web. Conversely, existing constraints on the fraction of compact dark matter permit as many as 1 in $\approx 40$ of all $z \lesssim 0.4$ FRBs to exhibit micro-lensed burst structure. Approximately $170$ FRBs intercepting halos within $\sim 50\,$kpc would be required to exclude the fraction of compact dark matter in each intercepted halo to a similar level. Furthermore, we consider the cumulative effects of lensing of the FRB signal by a macroscopic dark matter distribution. We conclude that lensing from a uniform distribution of compact objects is likely not observable, but suggest that FRBs may set meaningful limits on power-law distributions of dark matter.
[6]  oai:arXiv.org:2004.11522  [pdf] - 2117641
Spectropolarimetric properties of Swift J1818.0$-$1607: a 1.4 s radio magnetar
Comments: 10 pages, 5 figures and 1 table. Revised version accepted for publication in The Astrophysical Journal Letters
Submitted: 2020-04-23, last modified: 2020-06-01
The soft-gamma repeater Swift J1818.0$-$1607 is only the fifth magnetar found to exhibit pulsed radio emission. Using the Ultra-Wideband Low receiver system of the Parkes radio telescope, we conducted a 3 h observation of Swift J1818.0$-$1607. Folding the data at a rotation period of $P = 1.363$ s, we obtained wideband polarization profiles and flux density measurements covering radio frequencies between 704-4032 MHz. After measuring, and then correcting for the pulsar's rotation measure of $1442.0 \pm 0.2$ rad m$^{-2}$, we find the radio profile is between 80-100 per cent linearly polarized across the wide observing band, with a small amount of depolarization at low frequencies that we ascribe to scatter broadening. We also measure a steep spectral index of $\alpha = -2.26^{+0.02}_{-0.03}$ across our large frequency range, a significant deviation from the flat or inverted spectra often associated with radio-loud magnetars. The steep spectrum and temporal rise in flux density bears some resemblance to the behaviour of the magnetar-like, rotation-powered pulsar PSR J1119$-$6127. This leads us to speculate that Swift J1818.0$-$1607 may represent an additional link between rotation-powered pulsars and magnetars.
[7]  oai:arXiv.org:2005.14366  [pdf] - 2104094
The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime
Comments: 20 pages, 16 Figures, 4 Tables, accepted for publication in PASA
Submitted: 2020-05-28
We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.
[8]  oai:arXiv.org:2005.13161  [pdf] - 2103002
A census of baryons in the Universe from localized fast radio bursts
Comments: Published online in Nature 27 May, 2020
Submitted: 2020-05-27
More than three quarters of the baryonic content of the Universe resides in a highly diffuse state that is difficult to observe, with only a small fraction directly observed in galaxies and galaxy clusters. Censuses of the nearby Universe have used absorption line spectroscopy to observe these invisible baryons, but these measurements rely on large and uncertain corrections and are insensitive to the majority of the volume, and likely mass. Specifically, quasar spectroscopy is sensitive either to only the very trace amounts of Hydrogen that exists in the atomic state, or highly ionized and enriched gas in denser regions near galaxies. Sunyaev-Zel'dovich analyses provide evidence of some of the gas in filamentary structures and studies of X-ray emission are most sensitive to gas near galaxy clusters. Here we report the direct measurement of the baryon content of the Universe using the dispersion of a sample of localized fast radio bursts (FRBs), thus utilizing an effect that measures the electron column density along each sight line and accounts for every ionised baryon. We augment the sample of published arcsecond-localized FRBs with a further four new localizations to host galaxies which have measured redshifts of 0.291, 0.118, 0.378 and 0.522, completing a sample sufficiently large to account for dispersion variations along the line of sight and in the host galaxy environment to derive a cosmic baryon density of $\Omega_{b} = 0.051_{-0.025}^{+0.021} \, h_{70}^{-1}$ (95% confidence). This independent measurement is consistent with Cosmic Microwave Background and Big Bang Nucleosynthesis values.
[9]  oai:arXiv.org:2005.13158  [pdf] - 2102999
Dissecting the Local Environment of FRB 190608 in the Spiral Arm of its Host Galaxy
Comments: 16 pages, 8 figure files, 7 figures in paper, 3 tables. Submitted to ApJ
Submitted: 2020-05-27
We present a high-resolution analysis of the host galaxy of fast radio burst FRB 190608, an SBc galaxy at $z=0.11778$ (hereafter HG 190608), to dissect its local environment and its contributions to the FRB properties. Our Hubble Space Telescope WFC3/UVIS image reveals that the sub-arcsecond localization of FRB 190608 is coincident with a knot of star-formation ($\Sigma_{SFR} = 1.2 \times 10^{-2}~ M_{\odot} \, kpc^{-2}$) in one of the prominent spiral arms of HG 190608. This is confirmed by H$\beta$ emission present in our Keck/KCWI integral field spectrum of the galaxy with a surface brightness of $\mu_{H\beta} = (3.35\pm0.18)\times10^{-17}\;erg\;s^{-1}\;cm^{-2}\;arcsec^{-2}$. We infer an extinction-corrected H$\alpha$ surface brightness and compute a dispersion measure from the interstellar medium of HG 190608 of ${DM}_{Host,ISM} = 82 \pm 35~ pc \, cm^{-3}$. The galaxy rotates with a circular velocity $v_{circ} = 141 \pm 8~ km \, s^{-1}$ at an inclination $i_{gas} = 37 \pm 3^\circ$, giving a dynamical mass $M_{halo}^{dyn} \approx 10^{11.96 \pm 0.08}~ M_{\odot}$. A surface photometric analysis of the galaxy using FORS2 imaging suggests a stellar disk inclination of $i_{stellar} = 26 \pm 3^\circ$. The dynamical mass estimate implies a halo contribution to the dispersion measure of ${DM}_{Host,Halo} = 55 \pm 25\; pc \, cm^{-3}$ subject to assumptions on the density profile and fraction of baryons retained. The relatively high temporal broadening ($\tau = 3.3 \pm 0.2 \; ms$ at 1.28 GHz) and rotation measure ($ RM = 353 \pm 2\; rad \; m^{-2}$) (Day et al. 2020) of FRB 190608 may be attributable to both turbulent gas within the spiral arm and gas local to the FRB progenitor. In contrast to previous high-resolution studies of FRB progenitor environments, we find no evidence for disturbed morphology, emission, nor kinematics for FRB 190608.
[10]  oai:arXiv.org:2005.13157  [pdf] - 2102998
Disentangling the Cosmic Web Towards FRB 190608
Comments: 14 pages, 9 figures, 3 tables. The full version of Table 1 is available as a LaTeX file. Only the first 10 entries are present in the print version. Submitted to ApJ
Submitted: 2020-05-27
FRB 190608 was detected by ASKAP and localized to a spiral galaxy at $z_{host}=0.11778$ in the SDSS footprint. The burst has a large dispersion measure ($DM_{FRB}=339.8$ $pc/cm^3$) compared to the expected cosmic average at its redshift. It also has a large rotation measure ($RM_{FRB}=353$ $rad/m^2$) and scattering timescale ($\tau=3.3$ $ms$ at $1.28$ $GHz$). Chittidi et al (2020) perform a detailed analysis of the ultraviolet and optical emission of the host galaxy and estimate the host DM contribution to be $137\pm 43$ $pc/cm^3$. This work complements theirs and reports the analysis of the optical data of galaxies in the foreground of FRB 190608 to explore their contributions to the FRB signal. Together, the two manuscripts delineate an observationally driven, end-to-end study of matter distribution along an FRB sightline; the first study of its kind. Combining KCWI observations and public SDSS data, we estimate the expected cosmic dispersion measure $DM_{cosmic}$ along the sightline to FRB 190608. We first estimate the contribution of hot, ionized gas in intervening virialized halos ($DM_{halos} \approx 7-28$ $pc/cm^3$). Then, using the Monte Carlo Physarum Machine (MCPM) methodology, we produce a 3D map of ionized gas in cosmic web filaments and compute the DM contribution from matter outside halos ($DM_{IGM} \approx 91-126$ $pc/cm^3$). This implies a greater fraction of ionized gas along this sightline is extant outside virialized halos. We also investigate whether the intervening halos can account for the large FRB rotation measure and pulse width and conclude that it is implausible. Both the pulse broadening and the large Faraday rotation likely arise from the progenitor environment or the host galaxy.
[11]  oai:arXiv.org:2005.13160  [pdf] - 2110104
The host galaxies and progenitors of Fast Radio Bursts localized with the Australian Square Kilometre Array Pathfinder
Comments: 15 pages, 5 figures, 5 tables. Accepted for publication in ApJL
Submitted: 2020-05-27
The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11<z<0.48$. All four are massive galaxies (log( $M_{*}/ M_{\odot}$) $\sim 9.4 -10.4$) with modest star-formation rates of up to $2M_{\odot}$yr$^{-1}$ -- very different to the host galaxy of the first repeating FRB 121102, which is a dwarf galaxy with a high specific star-formation rate. The FRBs localized by ASKAP typically lie in the outskirts of their host galaxies, which appears to rule out FRB progenitor models that invoke active galactic nuclei (AGN) or free-floating cosmic strings. The stellar population seen in these host galaxies also disfavors models in which all FRBs arise from young magnetars produced by superluminous supernovae (SLSNe), as proposed for the progenitor of FRB 121102. A range of other progenitor models (including compact-object mergers and magnetars arising from normal core-collapse supernovae) remain plausible.
[12]  oai:arXiv.org:2005.13162  [pdf] - 2103003
High time resolution and polarisation properties of ASKAP-localised fast radio bursts
Comments: 18 pages, 4 figures, 5 tables, submitted to MNRAS
Submitted: 2020-05-27
Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescope has provided the first ensemble of bursts with this information. In this paper, we present the high time and spectral resolution, full polarisation observations of five localised FRBs to complement the results published for the previously studied ASKAP FRB~181112. We find that every FRB is highly polarised, with polarisation fractions ranging from 80 -- 100\%, and that they are generally dominated by linear polarisation. While some FRBs in our sample exhibit properties associated with an emerging archetype (i.e., repeating or apparently non-repeating), others exhibit characteristic features of both, implying the existence of a continuum of FRB properties. When examined at high time resolution, we find that all FRBs in our sample have evidence for multiple sub-components and for scattering at a level greater than expected from the Milky Way. We find no correlation between the diverse range of FRB properties (e.g., scattering time, intrinsic width, and rotation measure) and any global property of their host galaxy. The most heavily scattered bursts reside in the outskirts of their host galaxies, suggesting that the source-local environment rather than the host interstellar medium is likely the dominant origin of the scattering in our sample.
[13]  oai:arXiv.org:2005.05558  [pdf] - 2129331
Radio observations of two intermittent pulsars: PSRs J1832+0029 and J1841-0500
Comments: 7 pages, 7 figures, accepted for publication in ApJ
Submitted: 2020-05-12
We present long-term observations of two intermittent pulsars, PSRs~J1832+0029 and J1841$-$0500 using the Parkes 64\,m radio telescope. The radio emission for these pulsars switches "off" for year-long durations. Our new observations have enabled us to improve the determination of the on-off timescales and the spin down rates during those emission states. In general our results agree with previous studies of these pulsars, but we now have significantly longer data spans. We have identified two unexpected signatures in the data. Weak emission was detected in a single observation of PSR~J1832$+$0029 during an "off" emission state. For PSR~J1841$-$0500, we identified a quasi-periodic fluctuation in the intensities of the detectable single pulses, with a modulation period between 21 and 36 pulse periods.
[14]  oai:arXiv.org:2004.14668  [pdf] - 2119976
Very long baseline astrometry of PSR J1012+5307 and its implications on alternative theories of gravity
Comments: Accepted for publication in ApJ
Submitted: 2020-04-30
PSR J1012+5307, a millisecond pulsar in orbit with a helium white dwarf (WD), has been timed with high precision for about 25 years. One of the main objectives of this long-term timing is to use the large asymmetry in gravitational binding energy between the neutron star and the WD to test gravitational theories. Such tests, however, will be eventually limited by the accuracy of the distance to the pulsar. Here, we present VLBI (very long baseline interferometry) astrometry results spanning approximately 2.5 years for PSR J1012+5307, obtained with the Very Long Baseline Array as part of the MSPSRPI project. These provide the first proper motion and absolute position for PSR J1012+5307 measured in a quasi-inertial reference frame. From the VLBI results, we measure a distance of $0.83^{+0.06}_{-0.02}$kpc (all the estimates presented in the abstract are at 68% confidence) for PSR J1012+5307, which is the most precise obtained to date. Using the new distance, we improve the uncertainty of measurements of the unmodeled contributions to orbital period decay, which, combined with three other pulsars, places new constraints on the coupling constant for dipole gravitational radiation $\kappa_D=(-1.7\pm1.7)\times 10^{-4}$ and the fractional time derivative of Newton's gravitational constant $\dot{G}/G = -1.8^{\,+5.6}_{\,-4.7}\times 10^{-13}\,{\rm yr^{-1}}$ in the local universe. As the uncertainties of the observed decays of orbital period for the four leading pulsar-WD systems become negligible in $\approx10$ years, the uncertainties for $\dot{G}/G$ and $\kappa_D$ will be improved to $\leq1.5\times10^{-13}\,{\rm yr^{-1}}$ and $\leq1.0\times10^{-4}$, respectively, predominantly limited by the distance uncertainties.
[15]  oai:arXiv.org:2003.13303  [pdf] - 2077140
Timing of young radio pulsars II. Braking indices and their interpretation
Comments: 15 pages, 10 figures and 5 tables
Submitted: 2020-03-30
In Paper I of this series, we detected a significant value of the braking index ($n$) for 19 young, high-$\dot{E}$ radio pulsars using $\sim$ 10 years of timing observations from the 64-m Parkes radio telescope. Here we investigate this result in more detail using a Bayesian pulsar timing framework to model timing noise and to perform selection to distinguish between models containing exponential glitch recovery and braking index signatures. We show that consistent values of $n$ are maintained with the addition of substantial archival data, even in the presence of glitches. We provide strong arguments that our measurements are unlikely due to exponential recovery signals from unseen glitches even though glitches play a key role in the evolution of a pulsar's spin frequency. We conclude that, at least over decadal time scales, the value of $n$ can be significantly larger than the canonical 3 and discuss the implications for the evolution of pulsars.
[16]  oai:arXiv.org:1912.07847  [pdf] - 2101335
Which bright fast radio bursts repeat?
Comments: 12 pages, 7 figures, 5 tables. Updated after error discovered in original (v1) submission
Submitted: 2019-12-17, last modified: 2020-03-23
A handful of fast radio bursts (FRBs) are now known to repeat. However, the question remains --- do they all? We report on an extensive observational campaign with the Australian Square Kilometre Array Pathfinder (ASKAP), Parkes, and Robert C. Byrd Green Bank Telescope, searching for repeat bursts from FRBs detected by the Commensal Real-time ASKAP Fast Transients survey. In 383.2 hr of follow-up observations covering 27 FRBs initially detected as single bursts, only two repeat bursts from a single FRB, FRB 171019, were detected, which have been previously reported by Kumar et al. We use simulations of repeating FRBs that allow for clustering in burst arrival times to calculate new estimates for the repetition rate of FRB 171019, finding only slight evidence for incompatibility with the properties of FRB 121102. Our lack of repeat bursts from the remaining FRBs set limits on the model of all bursts being attributable to repeating FRBs. Assuming a reasonable range of repetition behaviour, at most 60% (90% C.L.) of these FRBs having an intrinsic burst distribution similar to FRB~121102. This result is shown to be robust against different assumptions on the nature of repeating FRB behaviour, and indicates that if indeed all FRBs repeat, the majority must do so very rarely.
[17]  oai:arXiv.org:2003.09780  [pdf] - 2110000
The Parkes Pulsar Timing Array Project: Second data release
Comments: 26 pages, 6 figures. Accepted for publication in PASA. The data release can be downloaded from https://doi.org/10.25919/5db90a8bdeb59
Submitted: 2020-03-21
We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 ${\mu}$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files.
[18]  oai:arXiv.org:2002.12481  [pdf] - 2059713
The UTMOST pulsar timing programme II: Timing noise across the pulsar population
Comments: Accepted by MNRAS. 28 pages, 8 figures, 8 tables
Submitted: 2020-02-27
While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesis Telescope (MOST). We reanalysed the spin and spin-down changes associated with nine previously reported pulsar glitches, report the discovery of three new glitches and four unusual glitch-like events in the rotational evolution of PSR J1825$-$0935. We develop a refined Bayesian framework for determining how red noise strength scales with pulsar spin frequency ($\nu$) and spin-down frequency ($\dot{\nu}$), which we apply to a sample of 280 non-recycled pulsars. With this new method and a simple power-law scaling relation, we show that red noise strength scales across the non-recycled pulsar population as $\nu^{a} |\dot{\nu}|^{b}$, where $a = -0.84^{+0.47}_{-0.49}$ and $b = 0.97^{+0.16}_{-0.19}$. This method can be easily adapted to utilise more complex, astrophysically motivated red noise models. Lastly, we highlight our timing of the double neutron star PSR J0737$-$3039, and the rediscovery of a bright radio pulsar originally found during the first Molonglo pulsar surveys with an incorrectly catalogued position.
[19]  oai:arXiv.org:2002.12539  [pdf] - 2065447
Spectropolarimetric analysis of FRB 181112 at microsecond resolution: Implications for Fast Radio Burst emission mechanism
Comments: Accepted to ApJ Letters; 13 pages, 5 figures
Submitted: 2020-02-27
We have developed a new coherent dedispersion mode to study the emission of Fast Radio Bursts that trigger the voltage capture capability of the Australian SKA Pathfinder (ASKAP) interferometer. In principle the mode can probe emission timescales down to 3 ns with full polarimetric information preserved. Enabled by the new capability, here we present a spectropolarimetric analysis of FRB 181112 detected by ASKAP, localized to a galaxy at redshift 0.47. At microsecond time resolution the burst is resolved into four narrow pulses with a rise time of just $15 \mu$s for the brightest. The pulses have a diversity of morphology, but do not show evidence for temporal broadening by turbulent plasma along the line of sight, nor is there any evidence for periodicity in their arrival times. The pulses are highly polarized (up to 95%), with the polarization position angle varying both between and within pulses. The pulses have apparent rotation measures that vary by $15\pm 2\, {\rm rad \,m^{-2}}$ and apparent dispersion measures that vary by $0.041\pm 0.004\,{\rm pc\,cm^{-3}}$. Conversion between linear and circular polarization is observed across the brightest pulse. We conclude that the FRB 181112 pulses are most consistent with being a direct manifestation of the emission process or the result of propagation through a relativistic plasma close to the source. This demonstrates that our method, which facilitates high-time-resolution polarimetric observations of FRBs, can be used to study not only burst emission processes, but also a diversity of propagation effects present on the gigaparsec paths they traverse.
[20]  oai:arXiv.org:2002.10250  [pdf] - 2057948
The Thousand-Pulsar-Array programme on MeerKAT I: Science objectives and first results
Comments: Accepted by MNRAS
Submitted: 2020-02-18
We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.
[21]  oai:arXiv.org:1912.10827  [pdf] - 2020696
Detection of a Glitch in PSR J0908$-$4913 by UTMOST
Comments: 2 pages, 1 table
Submitted: 2019-12-17
We report the first detection of a glitch in the radio pulsar PSR J0908$-$4913 (PSR B0906$-$49) during regular timing observations by the Molonglo Observatory Synthesis Telescope (MOST) as part of the UTMOST project.
[22]  oai:arXiv.org:1810.03227  [pdf] - 2014965
Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter
Comments: Fixed a typo in Fig. 5; Published in Phys. Rev. D as Editor's Suggestion article, 17 pages, 7 figures, 4 tables
Submitted: 2018-10-07, last modified: 2019-12-16
It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities.
[23]  oai:arXiv.org:1908.10026  [pdf] - 2017218
Faint Repetitions from a Bright Fast Radio Burst Source
Comments: 9 pages, 4 figures; accepted for publication in ApJL; updated plots, revised scattering analysis
Submitted: 2019-08-27, last modified: 2019-12-08
We report the detection of repeat bursts from the source of FRB 171019, one of the brightest fast radio bursts (FRBs) detected in the Australian Square Kilometre Array Pathfinder (ASKAP) fly's eye survey. Two bursts from the source were detected with the Green Bank Telescope in observations centered at 820 MHz. The repetitions are a factor of $\sim 590$ fainter than the ASKAP-discovered burst. All three bursts from this source show no evidence of scattering and have consistent pulse widths. The pulse spectra show modulation that could be evidence for either steep spectra or patchy emission. The two repetitions were the only ones found in an observing campaign for this FRB totaling 1000 hr, which also included ASKAP and the 64-m Parkes radio telescope, over a range of frequencies (720$-$2000 MHz) at epochs spanning two years. The inferred scaling of repetition rate with fluence of this source agrees with the other repeating source, FRB 121102. The detection of faint pulses from FRB 171019 shows that at least some FRBs selected from bright samples will repeat if follow-up observations are conducted with more sensitive telescopes.
[24]  oai:arXiv.org:1911.12024  [pdf] - 2042276
Probing the emission states of PSR J1107-5907
Comments: 7 pages, 6 figures, accepted by ApJ
Submitted: 2019-11-27
The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state.
[25]  oai:arXiv.org:1911.00656  [pdf] - 2076708
An ultra-wide bandwidth (704 to 4032 MHz) receiver for the Parkes radio telescope
Comments: submitted to PASA
Submitted: 2019-11-02
We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.
[26]  oai:arXiv.org:1910.13628  [pdf] - 2026196
A pulsar-based timescale from the International Pulsar Timing Array
Comments: accepted by MNRAS
Submitted: 2019-10-29
We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.
[27]  oai:arXiv.org:1909.11681  [pdf] - 2068935
The low density and magnetization of a massive galaxy halo exposed by a fast radio burst
Comments: Published in Science on 2019 September 26; Main (3 figures; 1 Table) + Supp (12 figures; 7 Tables)
Submitted: 2019-09-25
Present-day galaxies are surrounded by cool and enriched halo gas extending to hundreds of kiloparsecs. This halo gas is thought to be the dominant reservoir of material available to fuel future star formation, but direct constraints on its mass and physical properties have been difficult to obtain. We report the detection of a fast radio burst (FRB 181112) with arcsecond precision, which passes through the halo of a foreground galaxy. Analysis of the burst shows the halo gas has low net magnetization and turbulence. Our results imply predominantly diffuse gas in massive galactic halos, even those hosting active supermassive black holes, contrary to some previous results.
[28]  oai:arXiv.org:1909.05779  [pdf] - 1967091
A fast radio burst in the direction of the Virgo cluster
Comments: 9 pages, 6 figures, Accepted for publication in MNRAS
Submitted: 2019-09-12
The rate of fast radio bursts (FRBs) in the direction of nearby galaxy clusters is expected to be higher than the mean cosmological rate if intrinsically faint FRBs are numerous. In this paper, we describe a targeted search for faint FRBs near the core of the Virgo cluster using the Australian Square Kilometer Array Pathfinder telescope. During 300 hr of observations, we discovered one burst, FRB 180417, with dispersion measure DM $=474.8$ cm$^{-3}$pc. The FRB was promptly followed up by several radio telescopes for 27 h, but no repeat bursts were detected. An optical follow-up of FRB 180417 using the PROMPT5 telescope revealed no new sources down to an $R$-band magnitude of 20.1. We argue that FRB 180417 is likely behind the Virgo cluster as the Galactic and intracluster DM contribution are small compared to the DM of the FRB, and there are no galaxies in the line of sight. The non-detection of FRBs from Virgo constrains the faint-end slope, $\alpha<1.52$ (at 68\% confidence limit), and the minimum luminosity, $L_{\rm min}\gtrsim 2\times 10^{40}$ erg s$^{-1}$ (at 68\% confidence limit), of the FRB luminosity function assuming cosmic FRB rate of $10^4$ FRBs sky$^{-1}$ day$^{-1}$ with flux above 1 Jy located out to redshift of 1. Further FRB surveys of galaxy clusters with high-sensitivity instruments will tighten the constraints on the faint end of the luminosity function and, thus, are strongly encouraged.
[29]  oai:arXiv.org:1909.04534  [pdf] - 1983956
The International Pulsar Timing Array: Second data release
Comments: Submitted to MNRAS and in review, 23 pages, 5 figures
Submitted: 2019-09-10
In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.
[30]  oai:arXiv.org:1908.11709  [pdf] - 1967033
Timing of young radio pulsars I: Timing noise, periodic modulation and proper motion
Comments:
Submitted: 2019-08-29
The smooth spin-down of young pulsars is perturbed by two non-deterministic phenomenon, glitches and timing noise. Although the timing noise provides insights into nuclear and plasma physics at extreme densities, it acts as a barrier to high-precision pulsar timing experiments. An improved methodology based on Bayesian inference is developed to simultaneously model the stochastic and deterministic parameters for a sample of 85 high-$\dot{E}$ radio pulsars observed for $\sim$ 10 years with the 64-m Parkes radio telescope. Timing noise is known to be a red process and we develop a parametrization based on the red-noise amplitude ($A_{\rm red}$) and spectral index ($\beta$). We measure the median $A_{\rm red}$ to be $-10.4^{+1.8}_{-1.7}$ yr$^{3/2}$ and $\beta$ to be $-5.2^{+3.0}_{-3.8}$ and show that the strength of timing noise scales proportionally to $\nu^{1}|\dot{\nu}|^{-0.6\pm0.1}$, where $\nu$ is the spin frequency of the pulsar and $\dot{\nu}$ its spin-down rate. Finally, we measure significant braking indices for 19 pulsars, proper motions for two pulsars and discuss the presence of periodic modulation in the arrival times of five pulsars.
[31]  oai:arXiv.org:1908.02911  [pdf] - 1958248
The emission and scintillation properties of RRAT J2325-0530 at 154 MHz and 1.4 GHz
Comments: 18 pages, 8 figures, 5 tables, accepted for publication in PASA
Submitted: 2019-08-07
Rotating Radio Transients (RRATs) represent a relatively new class of pulsar, primarily characterised by their sporadic bursting emission of single pulses on time scales of minutes to hours. In addition to the difficulty involved in detecting these objects, low-frequency ($<$300 MHz) observations of RRATs are sparse, which makes understanding their broadband emission properties in the context of the normal pulsar population problematic. Here, we present the simultaneous detection of RRAT J2325-0530 using the Murchison Widefield Array (154 MHz) and Parkes radio telescope (1.4 GHz). On a single-pulse basis, we produce the first polarimetric profile of this pulsar, measure the spectral index ($\alpha=-2.2\pm 0.1$), pulse energy distributions, and present the pulse rates in the context of detections in previous epochs. We find that the distribution of time between subsequent pulses is consistent with a Poisson process and find no evidence of clustering over the $\sim$1.5 hr observations. Finally, we are able to quantify the scintillation properties of RRAT J2325-0530 at 1.4 GHz, where the single pulses are modulated substantially across the observing bandwidth, and show that this characterisation is feasible even with irregular time sampling as a consequence of the sporadic emission behaviour.
[32]  oai:arXiv.org:1907.08916  [pdf] - 1966872
A high time resolution study of the millisecond pulsar J2241-5236 at frequencies below 300 MHz
Comments: 09 pages, 4 figures, 2 tables
Submitted: 2019-07-21
One of the major challenges for pulsar timing array (PTA) experiments is the mitigation of the effects of the turbulent interstellar medium (ISM) from timing data. These can potentially lead to measurable delays and/or distortions in the pulse profiles and scale strongly with the inverse of the radio frequency. Low-frequency observations are therefore highly appealing for characterizing them. However, in order to achieve the necessary time resolution to resolve profile features of short-period millisecond pulsars, phase-coherent de-dispersion is essential, especially at frequencies below $300$ MHz. We present the lowest-frequency ($80$-$220$ MHz), coherently de-dispersed detections of one of the most promising pulsars for current and future PTAs, PSR J2241$-$5236, using our new beam-former software for the MWA's voltage capture system (VCS), which reconstructs the time series at a much higher time resolution of $\sim 1 \mu$s by re-synthesizing the recorded voltage data at $10$-kHz/$100$-$\mu$s native resolutions. Our data reveal a dual-precursor type feature in the pulse profile that is either faint or absent in high-frequency observations from Parkes. The resultant high-fidelity detections have enabled dispersion measure (DM) determinations with very high precision, of the order of $(2$-$6)\times10^{-6}$ $\rm pc\,cm^{-3}$, owing to the microsecond level timing achievable for this pulsar at the MWA's low frequencies. This underscores the usefulness of low-frequency observations for probing the ISM toward PTA pulsars and informing optimal observing strategies for PTA experiments.
[33]  oai:arXiv.org:1906.11476  [pdf] - 1907636
A single fast radio burst localized to a massive galaxy at cosmological distance
Comments: Published online in Science 27 June 2019
Submitted: 2019-06-27
Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.
[34]  oai:arXiv.org:1906.09793  [pdf] - 1912956
Commensal discovery of four Fast Radio Bursts during Parkes Pulsar Timing Array observations
Comments: 8 pages, accepted for publication in MNRAS
Submitted: 2019-06-24
The Parkes Pulsar Timing Array (PPTA) project monitors two dozen millisecond pulsars (MSPs) in order to undertake a variety of fundamental physics experiments using the Parkes 64m radio telescope. Since June 2017 we have been undertaking commensal searches for fast radio bursts (FRBs) during the MSP observations. Here, we report the discovery of four FRBs (171209, 180309, 180311 and 180714). The detected events include an FRB with the highest signal-to-noise ratio ever detected at the Parkes observatory, which exhibits unusual spectral properties. All four FRBs are highly polarized. We discuss the future of commensal searches for FRBs at Parkes.
[35]  oai:arXiv.org:1905.02293  [pdf] - 1912829
Five new real-time detections of Fast Radio Bursts with UTMOST
Comments: 13 pages, 11 figures, submitted to MNRAS
Submitted: 2019-05-06
We detail a new fast radio burst (FRB) survey with the Molonglo Radio Telescope, in which six FRBs were detected between June 2017 and December 2018. By using a real-time FRB detection system, we captured raw voltages for five of the six events, which allowed for coherent dedispersion and very high time resolution (10.24 $\mu$s) studies of the bursts. Five of the FRBs show temporal broadening consistent with interstellar and/or intergalactic scattering, with scattering timescales ranging from 0.16 to 29.1 ms. One burst, FRB181017, shows remarkable temporal structure, with 3 peaks each separated by 1 ms. We searched for phase-coherence between the leading and trailing peaks and found none, ruling out lensing scenarios. Based on this survey, we calculate an all-sky rate at 843 MHz of $98^{+59}_{-39}$ events sky$^{-1}$ day$^{-1}$ to a fluence limit of 8 Jy-ms: a factor of 7 below the rates estimated from the Parkes and ASKAP telescopes at 1.4 GHz assuming the ASKAP-derived spectral index $\alpha=-1.6$ ($F_{\nu}\propto\nu^{\alpha}$). Our results suggest that FRB spectra may turn over below 1 GHz. Optical, radio and X-ray followup has been made for most of the reported bursts, with no associated transients found. No repeat bursts were found in the survey.
[36]  oai:arXiv.org:1903.06525  [pdf] - 1855811
A Southern sky search for repeating Fast Radio Bursts using the Australian SKA Pathfinder
Comments: 7 pages, 8 figures, accepted in MNRAS
Submitted: 2019-03-13
We have conducted a search for bright repeating Fast Radio Bursts in our nearby Universe with the Australian Square Kilometer Array Pathfinder (ASKAP) in single-dish mode. We used eight ASKAP 12-m dishes, each equipped with a Chequerboard Phased Array Feed forming 36 beams on the sky, to survey $\sim$30,000 deg$^{2}$ of the southern sky ($-90^{\circ} < \delta < +30^{\circ}$) in 158 antenna days. The fluence limit of the survey is 22 Jyms. We report the detection of FRB 180515 in our survey. We found no repeating FRBs in a total mean observation of 3hrs per pointing divided into one-hour intervals, which were separated in time ranging between a day to a month. Using our non-detection, we exclude the presence of a repeating FRB similar to FRB 121102 closer than $z=0.004$ in the survey area --- a volume of at least $9.4 \times 10^4$Mpc$^3$ --- at 95% confidence.
[37]  oai:arXiv.org:1903.04694  [pdf] - 1851769
A Survey of the Galactic Plane for Dispersed Pulses with the Australian Square Kilometre Array Pathfinder
Comments: 9 pages, 7 figures. Accepted for publication in MNRAS
Submitted: 2019-03-11
We report the results from a survey of the Galactic plane for dispersed single pulses using the Australian SKA Pathfinder (ASKAP). We searched for rare bright dispersed radio pulses comprising 160 pointings covering 4800 deg$^2$ of the Galactic plane within |b| $< 7 \deg$, each pointing with an exposure time of 10 hours. We detected one fast radio burst, FRB 180430, and single pulses from 11 pulsars. No rotating radio transients were detected. We detected FRB 180430 in the Galactic plane in the anticentre direction with a fluence of 216$\pm5$Jy ms a dispersion measure (DM) of 264.1 pc cm-3. We estimate the extragalactic DM of the object to be less than 86.7 $ \text{pc} \ \text{cm}^{-3} $ depending on the electron density model. One model suggests that this FRB may be a giant pulse within our galaxy; we discuss how this may not correctly represent the line-of-sight DM. Based on the single detection of FRB 180430 in 3.47 $\times 10^{4}$ deg$^2$ h we derive a FRB event rate in the Galactic plane at the 20 Jy ms threshold to be in the range 2-140 per sky per day at 95% confidence. Despite the necessarily large uncertainties from this single detection, this is consistent with the current ASKAP all-sky detection rate.
[38]  oai:arXiv.org:1810.04356  [pdf] - 1842376
The Performance and Calibration of the CRAFT Fly's Eye Fast Radio Burst Survey
Comments: 21 pages, 15 figures, 2 tables, accepted for publication in PASA
Submitted: 2018-10-10, last modified: 2019-01-20
Since January 2017, the Commensal Real-time ASKAP Fast Transients survey (CRAFT) has been utilising commissioning antennas of the Australian SKA Pathfinder (ASKAP) to survey for fast radio bursts (FRBs) in fly's eye mode. This is the first extensive astronomical survey using phased array feeds (PAFs), and a total of 20 FRBs have been reported. Here we present a calculation of the sensitivity and total exposure of this survey, using the pulsars B1641-45 (J1644-4559) and B0833-45 (J0835-4510, i.e.\ Vela) as calibrators. The design of the survey allows us to benchmark effects due to PAF beamshape, antenna-dependent system noise, radio-frequency interference, and fluctuations during commissioning on timescales from one hour to a year. Observation time, solid-angle, and search efficiency are calculated as a function of FRB fluence threshold. Using this metric, effective survey exposures and sensitivities are calculated as a function of the source counts distribution. The implied FRB rate is significantly lower than the $37$\,sky$^{-1}$\,day$^{-1}$ calculated using nominal exposures and sensitivities for this same sample by \citet{craft_nature}. At the Euclidean power-law index of $-1.5$, the rate is $10.7_{-1.8}^{+2.7}\,{\rm (sys)} \, \pm \, 3\,{\rm (stat)}$\,sky$^{-1}$\,day$^{-1}$ above a threshold of $57\pm6\,{\rm (sys)}$\,Jy\,ms, while for the best-fit index for this sample of $-2.1$, it is $16.6_{-1.5}^{+1.9} \,{\rm (sys)}\, \pm 4.7\,{\rm (stat)}$\,sky$^{-1}$\,day$^{-1}$ above a threshold of $41.6\pm1.5\,{\rm (sys)}$\,Jy\,ms. This strongly suggests that these calculations be performed for other FRB-hunting experiments, allowing meaningful comparisons to be made between them.
[39]  oai:arXiv.org:1811.01142  [pdf] - 1802766
Hunting for radio emission from the intermittent pulsar J1107-5907 at low frequencies
Comments: 19 pages, 7 figures, 4 tables, accepted for publication in ApJ; minor updates to keep in-line with published version
Submitted: 2018-11-02, last modified: 2018-12-19
The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molonglo Observatory Synthesis Telescope (UTMOST) at 835 MHz, as part of an on-going observing campaign. During a 30-minute simultaneous observation, we detected the pulsar in its bright emission state for approximately 15 minutes, where 86 and 283 pulses were detected above a signal-to-noise threshold of 6 with the MWA and UTMOST, respectively. Of the detected pulses, 51 had counterparts at both frequencies and exhibited steep spectral indices for both the bright main pulse component and the precursor component. We find that the bright state pulse energy distribution is best parameterised by a log-normal distribution at both frequencies, contrary to previous results which suggested a power law distribution. Further low-frequency observations are required in order to explore in detail aspects such as pulse-to-pulse variability, intensity modulations and to better constrain the signal propagation effects due to the interstellar medium and intermittency characteristics at these frequencies. The spectral index, extended profile emission covering a large fraction of pulse longitude, and the broadband intermittency of PSR J1107-5907 suggests that future low-frequency pulsar searches, for instance those planned with SKA-Low, will be in an excellent position to find and investigate new pulsars of this type.
[40]  oai:arXiv.org:1812.00719  [pdf] - 1822991
Searching a Thousand Radio Pulsars for Gamma-ray Emission
Comments: Accepted by the Astrophysical Journal
Submitted: 2018-12-03
Identifying as many gamma-ray pulsars as possible in the Fermi Large Area Telescope (LAT) data helps test pulsar emission models by comparing predicted and observed properties for a large, varied sample with as little selection bias as possible. It also improves extrapolations from the observed population to estimate the contribution of unresolved pulsars to the diffuse gamma-ray emission. We use a recently developed method to determine the probability that a given gamma-ray photon comes from a known position in the sky, convolving the photon's energy with the LAT's energy-dependent point-spread-function (PSF), without the need for an accurate spatial and spectral model of the gamma-ray sky around the pulsar. The method is simple and fast and, importantly, provides probabilities, or _weights_, for gamma rays from pulsars too faint for phase-integrated detection. We applied the method to over a thousand pulsars for which we obtained rotation ephemerides from radio observations, and discovered gamma-ray pulsations from 16 pulsars, 12 young and 4 recycled. PSR J2208+4056 has spindown power Edot = 8e32 erg/s, about three times lower than the previous observed gamma-ray emission "deathline". PSRs J2208+4056 and J1816-0755 have radio interpulses, constraining their geometry and perhaps enhancing their gamma-ray luminosity. We discuss whether the deathline is an artifact of selection bias due to the pulsar distance.
[41]  oai:arXiv.org:1810.04355  [pdf] - 1778698
No low-frequency emission from extremely bright Fast Radio Bursts
Comments: 9 pages, 3 figures , 2 tables (accepted for publication in ApJL)
Submitted: 2018-10-10
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$ $\times$ $1.28$ MHz MWA images across the $170-200$ MHz band taken at 0.5 second time resolution at the known dispersion measures (DMs) and arrival times of the bursts and searched both within the ASKAP error regions (typically $\sim$ $10$ arcmin $\times$ $10$ arcmin), and beyond ($4$ deg $\times$ $4$ deg). We identified no candidates exceeding a $5\sigma$ threshold at these DMs in the dynamic spectra. These limits are inconsistent with the mean fluence scaling of $\alpha=-1.8 \pm 0.3$ (${\cal F}_\nu \propto \nu^\alpha$, where $\nu$ is the observing frequency) that is reported for ASKAP events, most notably for the three high fluence (${\cal F}_{1.4\,{\rm 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 turn-over 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/10^4 K)^{-1.35}$ pc for FRB 171020.
[42]  oai:arXiv.org:1810.04354  [pdf] - 1779703
A search for the host galaxy of FRB171020
Comments: Submitted to ApJL
Submitted: 2018-10-10
We report on a search for the host galaxy of FRB171020, the fast radio burst with the smallest recorded dispersion measure (DM $=114$ pc cm$^{-3}$) of our on-ongoing ASKAP survey. The low DM confines the burst location within a sufficiently small volume to rigorously constrain the identity of the host galaxy. We identify 16 candidate galaxies in the search volume and single out ESO 601-G036, a Sc galaxy at redshift $z=0.00867$, as the most likely host galaxy. UV and optical imaging and spectroscopy reveal this galaxy has a star-formation rate of approximately 0.1 M$_\odot$ yr$^{-1}$ and oxygen abundance $12 + \log({\rm O/H}) = 8.3 \pm 0.2$, properties remarkably consistent with the galaxy hosting the repeating FRB121102. However, in contrast to FRB121102, follow-up radio observations of ESO 601-G036 show no compact radio emission above a 5$\sigma$ limit of $L_{2.1{\rm GHz}}=3.6\times 10^{19}$ W Hz$^{-1}$. Using radio continuum observations of the field, combined with archival optical imaging data, we find no analog to the persistent radio source associated with FRB121102 within the localization region of FRB171020 out to $z=0.06$. These results suggest that FRBs are not necessarily associated with a luminous and compact radio continuum source.
[43]  oai:arXiv.org:1810.04353  [pdf] - 1838206
The spectral properties of the bright FRB population
Comments: Submitted to ApJ Letters
Submitted: 2018-10-10
We examine the spectra of 23 fast radio bursts detected in a fly's-eye survey with the Australian SKA Pathfinder, including those of three bursts not previously reported. The mean spectral index of $\alpha = -1.6_{-0.2}^{+0.3}$ ($F_\nu \propto \nu^\alpha$) is close to that of the Galactic pulsar population. The sample is dominated by bursts exhibiting a large degree of spectral modulation: 17 exhibit fine-scale spectral modulation with an rms exceeding 50% of the mean, with decorrelation bandwidths (half-maximum) ranging from $\approx$ to 49 MHz. Most decorrelation bandwidths are an order of magnitude lower than the $\gtrsim 30\,$MHz expected of Galactic interstellar scintillation at the Galactic latitude of the survey, $|b| = 50 \pm 5 \deg$. A test of the amplitude distribution of the spectral fluctuations reveals only 12 bursts consistent at better than a 5% confidence level with the prediction of 100%-modulated diffractive scintillation. Moreover, five of six FRBs with a signal-to-noise ratio exceeding 18 are consistent with this prediction at less than 1% confidence. Nonetheless, there is weak evidence (88-95% confidence) that the amplitude of the fine-scale spectral modulation is anti-correlated with dispersion measure (DM) that would suggest it originates as a propagation effect. This effect appears to be corroborated by the smoothness of the higher-DM Parkes FRBs, and could arise due to quenching of diffractive scintillation (e.g. in the interstellar medium of the host galaxy) by angular broadening in the intergalactic medium.
[44]  oai:arXiv.org:1810.04357  [pdf] - 1791139
The Slope of the Source-Count Distribution for Fast Radio Bursts
Comments: 13 pages, 4 figures, 3 tables
Submitted: 2018-10-10
The slope of the source-count distribution of fast radio burst (FRB) fluences, $\alpha$, has been estimated using a variety of methods. Hampering all attempts have been the low number of detected FRBs, and the difficulty of defining a completeness threshold for FRB surveys. In this work, we extend maximum-likelihood methods for estimating $\alpha$, using detected and threshold signal-to-noise ratios applied to all FRBs in a sample without regard to a completeness threshold. Using this method with FRBs detected by the Parkes radio telescope, we find $\alpha=-1.18 \pm 0.24$ (68\% confidence interval, C.I.), i.e.\ consistent with a non-evolving Euclidean distribution ($\alpha=-1.5$). Applying these methods to the Australian Square Kilometre Array Pathfinder (ASKAP) Commensal Real-time ASKAP Fast Transients (CRAFT) FRB survey finds $\alpha=-2.2 \pm 0.47$ (68\% C.I.). A full maximum-likelihood estimate finds an inconsistency with the Parkes rate with a p-value of 0.86\% ($2.6\, \sigma$). If not due to statistical fluctuations or biases in Parkes data, this is the first evidence for deviations from a pure power law in the integral source-count distribution of FRBs. It is consistent with a steepening of the integral source-count distribution in the fluence range 5--40\,Jy\,ms, for instance due to a cosmological population of FRB progenitors evolving more rapidly than the star-formation rate, and peaking in the redshift range 1--3.
[45]  oai:arXiv.org:1809.10744  [pdf] - 1767703
Studying the solar system with the International Pulsar Timing Array
Comments: Accepted for publication by the Monthly Notices of the Royal Astronomical Society
Submitted: 2018-09-27
Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.
[46]  oai:arXiv.org:1809.01309  [pdf] - 1775671
Polarimetry of the eclipsing pulsar PSR J1748$-$2446A
Comments: 9 papes, 3 figures, ApJ accepted
Submitted: 2018-09-04
Observations with the Parkes radio telescope of the eclipsing millisecond binary pulsar PSR J1748$-$2446A, which is in the globular cluster Terzan 5, are presented. These include the first observations of this pulsar in the 3\,GHz frequency band, along with simultaneous observations in the 700\,MHz band and new observations around 1400\,MHz. We show that the pulsar signal is not eclipsed in the 3\,GHz band and observe the known eclipses in the lower frequency bands. We find that the observed pulse signal becomes depolarized during particular orbital phases and postulate that this depolarization occurs because of rotation-measure fluctuations resulting from turbulence in the stellar wind responsible for the eclipses.
[47]  oai:arXiv.org:1808.02580  [pdf] - 1730308
Detection of a glitch in the pulsar J1709-4429
Comments: 3 pages, 1 figure
Submitted: 2018-08-07
We report the detection of a glitch event in the pulsar J1709$-$4429 (also known as B1706$-$44) during regular monitoring observations with the Molonglo Observatory Synthesis Telescope (UTMOST). The glitch was found during timing operations, in which we regularly observe over 400 pulsars with up to daily cadence, while commensally searching for Rotating Radio Transients, pulsars, and FRBs. With a fractional size of $\Delta\nu/\nu \approx 52.4 \times10^{-9}$, the glitch reported here is by far the smallest known for this pulsar, attesting to the efficacy of glitch searches with high cadence using UTMOST.
[48]  oai:arXiv.org:1807.01302  [pdf] - 1790804
First Ground-based Measurement of Sub-20 GeV to 100 GeV $\gamma$-rays from the Vela Pulsar with H.E.S.S. II
Collaboration, H. E. S. S.; Abdalla, H.; Aharonian, F.; Benkhali, F. Ait; Angüner, E. O.; Arakawa, M.; Arcaro, C.; Arm, C.; Arrieta, M.; Backes, M.; Barnard, M.; Becherini, Y.; Tjus, J. Becker; Berge, D.; Bernhard, S.; Bernlöhr, K.; Blackwell, R.; Böttcher, M.; Boisson, C.; Bolmont, J.; Bonnefoy, S.; Bordas, P.; Bregeon, J.; Brun, F.; Brun, P.; Bryan, M.; Büchele, M.; Bulik, T.; Bylund, T.; Capasso, M.; Caroff, S.; Carosi, A.; Casanova, S.; Cerruti, M.; Chakraborty, N.; Ch, S.; Chaves, R. C. G.; Chen, A.; Colafrancesco, S.; Condon, B.; Davids, I. D.; Deil, C.; Devin, J.; deWilt, P.; Dirson, L.; Djannati-Ataï, A.; Dmytriiev, A.; Donath, A.; Doroshenko, V.; Drury, L. O'C.; Dyks, J.; Egberts, K.; Emery, G.; Ernenwein, J. -P.; Eschbach, S.; Fegan, S.; Fiasson, A.; Fontaine, G.; Funk, S.; Füßling, M.; Gabici, S.; Gallant, Y. A.; Gaté, F.; Giavitto, G.; Glawion, D.; Glicenstein, J. F.; Gottschall, D.; Grondin, M. -H.; Hahn, J.; Haupt, M.; Heinzelmann, G.; Henri, G.; Hermann, G.; Hinton, J. A.; Hofmann, W.; Hoischen, C.; Holch, T. L.; Holler, M.; Horns, D.; Huber, D.; Iwasaki, H.; Jacholkowska, A.; Jamrozy, M.; Jankowsky, D.; Jankowsky, F.; Jouvin, L.; Jung-Richardt, I.; Kastendieck, M. A.; Katarzyński, K.; Katsuragawa, M.; Katz, U.; Kerszberg, D.; Khangulyan, D.; Khélifi, B.; King, J.; Klepser, S.; Kluźniak, W.; Komin, Nu.; Kosack, K.; Krakau, S.; Kraus, M.; Krüger, P. P.; Lamanna, G.; Lau, J.; Lefaucheur, J.; Lemière, A.; Lemoine-Goumard, M.; Lenain, J. -P.; Leser, E.; Lohse, T.; Lorentz, M.; López-Coto, R.; Lypova, I.; Malyshev, D.; Mar, V.; Marcowith, A.; Mariaud, C.; Martí-Devesa, G.; Marx, R.; Maurin, G.; Meintjes, P. J.; Mitchell, A. M. W.; Moderski, R.; Mohamed, M.; Mohrmann, L.; Moulin, E.; Murach, T.; Nakashima, S.; de Naurois, M.; Ndiyavala, H.; Niederwanger, F.; Niemiec, J.; Oakes, L.; O'Brien, P.; Odaka, H.; Ohm, S.; Ostrowski, M.; Oya, I.; Padovani, M.; Panter, M.; Parsons, R. D.; Perennes, C.; Petrucci, P. -O.; Peyaud, B.; Piel, Q.; Pita, S.; Poireau, V.; Noel, A. Priyana; Prokhorov, D. A.; Prokoph, H.; Pühlhofer, G.; Punch, M.; Quirrenbach, A.; Raab, S.; Rauth, R.; Reimer, A.; Reimer, O.; Renaud, M.; Rieger, F.; Rinchiuso, L.; Romoli, C.; Rowell, G.; Rudak, B.; Ruiz-Velasco, E.; Sahakian, V.; Saito, S.; Sanchez, D. A.; Santangelo, A.; Sasaki, M.; Schlickeiser, R.; Schüssler, F.; Schulz, A.; Schwanke, U.; Schwemmer, S.; Seglar-Arroyo, M.; Senniappan, M.; Seyffert, A. S.; Shafi, N.; Shilon, I.; Shiningayamwe, K.; Simoni, R.; Sinha, A.; Sol, H.; Spanier, F.; Specovius, A.; Spir-Jacob, M.; Stawarz, Ł.; Steenkamp, R.; Stegmann, C.; Steppa, C.; Takahashi, T.; Tavernet, J. -P.; Tavernier, T.; Taylor, A. M.; Terrier, R.; Tibaldo, L.; Tiziani, D.; Tluczykont, M.; Trichard, C.; Tsirou, M.; Tsuji, N.; Tuffs, R.; Uchiyama, Y.; van der Walt, D. J.; van Eldik, C.; van Rensburg, C.; van Soelen, B.; Vasileiadis, G.; Veh, J.; Venter, C.; Vincent, P.; Vink, J.; Voisin, F.; Völk, H. J.; Vuillaume, T.; Wadiasingh, Z.; Wagner, S. J.; Wagner, R. M.; White, R.; Wierzcholska, A.; Yang, R.; Zaborov, D.; Zacharias, M.; Zanin, R.; Zdziarski, A. A.; Zech, A.; Zefi, F.; Ziegler, A.; Zorn, J.; Żywucka, N.; Kerr, M.; Johnston, S.; Shannon, R. M.
Comments: Accepted for publication in A&A. 16 Pages, 9 figures (slightly language edited version)
Submitted: 2018-07-03, last modified: 2018-07-24
We report on the measurement and investigation of pulsed high-energy $\gamma$-ray emission from the Vela pulsar, PSR B0833$-$45, based on 40.3 hours of observations with the largest telescope of H.E.S.S., CT5, in monoscopic mode, and on 8 years of data obtained with the Fermi-LAT. A dedicated very-low-threshold event reconstruction and analysis pipeline was developed and, together with the CT5 telescope response model, was validated using the Fermi-LAT data as reference. A pulsed $\gamma$-ray signal at a significance level of more than $15\sigma$ is detected from the P2 peak of the Vela pulsar light curve. Of a total of 15835 events, more than 6000 lie at an energy below 20 GeV, implying a significant overlap between H.E.S.S. II-CT5 and the Fermi-LAT. While the investigation of the pulsar light curve with the LAT confirms characteristics previously known up to 20 GeV, in the tens of GeV energy range, CT5 data show a change in the pulse morphology of P2, i.e., an extreme sharpening of its trailing edge, together with the possible onset of a new component at 3.4$\sigma$ significance level. Assuming a power-law model for the P2 spectrum, an excellent agreement is found for the photon indices ($\Gamma \simeq$ 4.1) obtained with the two telescopes above 10 GeV and an upper bound of 8% is derived on the relative offset between their energy scales. Using both instruments data, it is however shown that the spectrum of P2 in the 10-100 GeV has a pronounced curvature, i.e. a confirmation of the sub-exponential cutoff form found at lower energies with the LAT. This is further supported by the weak evidence for an emission above 100 GeV obtained with CT5. In contrast, converging indications are found from both CT5 and LAT data for the emergence of a hard component above 50 GeV in the leading wing (LW2) of P2, which possibly extends beyond 100 GeV.
[49]  oai:arXiv.org:1807.06989  [pdf] - 1740075
Observations of Low-Frequency Radio Emission from Millisecond Pulsars and Multipath Propagation in the Interstellar Medium
Comments: 23 pages, 10 Figures, 4 Tables, Accepted to ApJS
Submitted: 2018-07-18
Studying the gravitational-wave sky with pulsar timing arrays (PTAs) is a key science goal for the Square Kilometre Array (SKA) and its pathfinder telescopes. With current PTAs reaching sub-microsecond timing precision, making accurate measurements of interstellar propagation effects and mitigating them effectively has become increasingly important to realise PTA goals. As these effects are much stronger at longer wavelengths, low-frequency observations are most appealing for characterizing the interstellar medium (ISM) along the sight lines toward PTA pulsars. The Murchison Widefield Array (MWA) and the Engineering Development Array (EDA), which utilizes MWA technologies, present promising opportunities for undertaking such studies, particularly for PTA pulsars located in the southern sky. Such pulsars are also the prime targets for PTA efforts planned with the South African MeerKAT, and eventually with the SKA. In this paper we report on observations of two bright southern millisecond pulsars PSRs J0437-4715 and J2145-0750 made with these facilities; MWA observations sampling multiple frequencies across the 80-250 MHz frequency range, while the EDA providing direct-sampled baseband data to yield a large instantaneous usable bandwidth of $\sim$200 MHz. Using these observations, we investigate various aspects relating to pulsar emission and ISM properties, such as spectral evolution of the mean pulse shape, scintillation as a function of frequency, chromaticity in interstellar dispersion, and flux density spectra at low frequencies. Systematic and regular monitoring observations will help ascertain the role of low-frequency measurements in PTA experiments, while simultaneously providing a detailed characterization of the ISM toward the pulsars, which will be useful in devising optimal observing strategies for future PTA experiments.
[50]  oai:arXiv.org:1804.08402  [pdf] - 1717105
The geometric distance and binary orbit of PSR B1259-63
Comments: Accepted by MNRAS. v2 is the journal-accepted version. Changes since v1 include a Bayesian approach to distance estimation, and discussion of Gaia DR2 results. 13 pages, 7 figures, 3 tables
Submitted: 2018-04-23, last modified: 2018-07-03
The pulsar/massive star binary system PSR B1259-63 / LS 2883 is one of the best-studied gamma-ray binaries, a class of systems whose bright gamma-ray flaring can provide important insights into high-energy physics. Using the Australian Long Baseline Array we have conducted very long baseline interferometric observations of PSR B1259-63 over 4.4 years, fully sampling the 3.4-year orbital period. From our measured parallax of $0.38\pm0.05$ mas we use a Bayesian approach to infer a distance of $2.6^{+0.4}_{-0.3}$ kpc. We find that the binary orbit is viewed at an angle of $154\pm3$ degrees to the line of sight, implying that the pulsar moves clockwise around its orbit as viewed on the sky. Taking our findings together with previous results from pulsar timing observations, all seven orbital elements for the system are now fully determined. We use our measurement of the inclination angle to constrain the mass of the stellar companion to lie in the range 15-31$M_{\odot}$. Our measured distance and proper motion are consistent with the system having originated in the Cen OB1 association and receiving a modest natal kick, causing it to have moved $\sim$8 pc from its birthplace over the past $\sim3\times10^5$ years. The orientation of the orbit on the plane of the sky matches the direction of motion of the X-ray synchrotron-emitting knot observed by the Chandra X-ray Observatory to be moving away from the system.
[51]  oai:arXiv.org:1803.05697  [pdf] - 1674948
FRB microstructure revealed by the real-time detection of FRB170827
Comments: 10 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2018-03-15, last modified: 2018-05-01
We report a new Fast Radio Burst (FRB) discovered in real-time as part of the UTMOST project at the Molonglo Observatory Synthesis Radio Telescope (MOST). FRB170827 is the first detected with our low-latency ($< 24$ s), machine-learning-based FRB detection system. The FRB discovery was accompanied by the capture of voltage data at the native time and frequency resolution of the observing system, enabling coherent dedispersion and detailed off-line analysis, which have unveiled fine temporal and frequency structure. The dispersion measure (DM) of 176.80 $\pm$ 0.04 pc cm$^{-3}$, is the lowest of the FRB population. The Milky Way contribution along the line of sight is $\sim$ 40 pc cm$^{-3}$, leaving an excess DM of $\sim$ 145 pc cm$^{-3}$. The FRB has a fluence $>$ 20 $\pm$ 7 Jy ms, and is narrow, with a width of $\sim$ 400 $\mu$s at 10$\%$ of its maximum amplitude. However, the burst shows three temporal components, the narrowest of which is $\sim$ 30 $\mu$s, and a scattering timescale of $4.1 \pm 2.7$ $\mu$s. The FRB shows spectral modulations on frequency scales of 1.5 MHz and 0.1 MHz. Both are prominent in the dynamic spectrum, which shows a very bright region of emission between 841 and 843 MHz, and weaker, patchy emission across the entire band. We show the fine spectral structure could arise in the FRB host galaxy, or its immediate vicinity.
[52]  oai:arXiv.org:1712.02646  [pdf] - 1714761
Spectral flattening of Crab Giant Pulses
Comments: 2 pages, 1 figure, to appear in the Proceedings of IAU Symposium 337, Pulsar Astrophysics - The Next 50 Years
Submitted: 2017-12-06
The frequency dependence of normal pulsar radio emission is typically observed to be a power law, with some indications of a flattening or turnover at low frequencies ($<100$ MHz). The spectrum of the Crab pulsar's giant pulse emission has not been examined as closely. We conducted simultaneous wideband observations of the Crab pulsar, with the Parkes radio telescope and the Murchison Widefield Array, to study the spectral behaviour of its giant pulses. Our analysis shows that the mean spectral index of Crab giant pulses flattens at low frequencies, from $-2.6\pm0.5$ between the Parkes bands, to $-0.7\pm1.4$ between the lowest MWA subbands.
[53]  oai:arXiv.org:1712.00426  [pdf] - 1608731
Extreme Scattering Events Towards Two Young Pulsars
Comments: MNRAS accepted
Submitted: 2017-12-01
We have measured the scintillation properties of 151 young, energetic pulsars with the Parkes radio telescope and have identified two extreme scattering events (ESEs). Towards PSR J1057-5226 we discovered a three-year span of strengthened scattering during which the variability in flux density and the scintillation bandwidth decreased markedly. The transverse size of the scattering region is ~23 au, and strong flux density enhancement before and after the ESE may arise from refractive focusing. Long observations reveal scintillation arcs characteristic of interference between rays scattered at large angles, and the clearest arcs appear during the ESE. The arcs suggest scattering by a screen 100-200 pc from the earth, perhaps ionized filamentary structure associated with the boundary of the local bubble(s). Towards PSR J1740-3015 we observed a "double dip" in the measured flux density similar to ESEs observed towards compact extragalactic radio sources. The observed shape is consistent with that produced by a many-au scale diverging plasma lens with electron density ~500cm$^{-3}$ . The continuing ESE is at least 1500 d long, making it the longest detected event to date. These detections, with materially different observational signatures, indicate that well-calibrated pulsar monitoring is a keen tool for ESE detection and ISM diagnostics. They illustrate the strong r\^ole au-scale non-Kolmogorov density fluctuations and the local ISM structure play in such events and are key to understanding both their intrinsic physics and their impact on other phenomena, particularly fast radio bursts.
[54]  oai:arXiv.org:1701.07709  [pdf] - 1931377
First search for gravitational waves from known pulsars with Advanced LIGO
The LIGO Scientific Collaboration; the Virgo Collaboration; Abbott, B. P.; Abbott, R.; Abbott, T. D.; Abernathy, M. R.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R. X.; Adya, V. B.; Affeldt, C.; Agathos, M.; Agatsuma, K.; Aggarwal, N.; Aguiar, O. D.; Aiello, L.; Ain, A.; Ajith, P.; Allen, B.; Allocca, A.; Altin, P. A.; Ananyeva, A.; Anderson, S. B.; Anderson, W. G.; Appert, S.; Arai, K.; Araya, M. C.; Areeda, J. S.; Arnaud, N.; Arun, K. G.; Ascenzi, S.; Ashton, G.; Ast, M.; Aston, S. M.; Astone, P.; Aufmuth, P.; Aulbert, C.; Avila-Alvarez, A.; Babak, S.; Bacon, P.; Bader, M. K. M.; Baker, P. T.; Baldaccini, F.; Ballardin, G.; Ballmer, S. W.; Barayoga, J. C.; Barclay, S. E.; Barish, B. C.; Barker, D.; Barone, F.; Barr, B.; Barsotti, L.; Barsuglia, M.; Barta, D.; Bartlett, J.; Bartos, I.; Bassiri, R.; Basti, A.; Batch, J. C.; Baune, C.; Bavigadda, V.; Bazzan, M.; Beer, C.; Bejger, M.; Belahcene, I.; Belgin, M.; Bell, A. S.; Berger, B. K.; Bergmann, G.; Berry, C. P. L.; Bersanetti, D.; Bertolini, A.; Betzwieser, J.; Bhagwat, S.; Bhandare, R.; Bilenko, I. A.; Billingsley, G.; Billman, C. R.; Birch, J.; Birney, R.; Birnholtz, O.; Biscans, S.; Bisht, A.; Bitossi, M.; Biwer, C.; Bizouard, M. A.; Blackburn, J. K.; Blackman, J.; Blair, C. D.; Blair, D. G.; Blair, R. M.; Bloemen, S.; Bock, O.; Boer, M.; Bogaert, G.; Bohe, A.; Bondu, F.; Bonnand, R.; Boom, B. A.; Bork, R.; Boschi, V.; Bose, S.; Bouffanais, Y.; Bozzi, A.; Bradaschia, C.; Brady, P. R.; Braginsky, V. B.; Branchesi, M.; Brau, J. E.; Briant, T.; Brillet, A.; Brinkmann, M.; Brisson, V.; Brockill, P.; Broida, J. E.; Brooks, A. F.; Brown, D. A.; Brown, D. D.; Brown, N. M.; Brunett, S.; Buchanan, C. C.; Buikema, A.; Bulik, T.; Bulten, H. J.; Buonanno, A.; Buskulic, D.; Buy, C.; Byer, R. L.; Cabero, M.; Cadonati, L.; Cagnoli, G.; Cahillane, C.; Bustillo, J. Calder'on; Callister, T. A.; Calloni, E.; Camp, J. B.; Canepa, M.; Cannon, K. C.; Cao, H.; Cao, J.; Capano, C. D.; Capocasa, E.; Carbognani, F.; Caride, S.; Diaz, J. Casanueva; Casentini, C.; Caudill, S.; Cavagli`a, M.; Cavalier, F.; Cavalieri, R.; Cella, G.; Cepeda, C. B.; Baiardi, L. Cerboni; Cerretani, G.; Cesarini, E.; Chamberlin, S. J.; Chan, M.; Chao, S.; Charlton, P.; Chassande-Mottin, E.; Cheeseboro, B. D.; Chen, H. Y.; Chen, Y.; Cheng, H. -P.; Chincarini, A.; Chiummo, A.; Chmiel, T.; Cho, H. S.; Cho, M.; Chow, J. H.; Christensen, N.; Chu, Q.; Chua, A. J. K.; Chua, S.; Chung, S.; Ciani, G.; Clara, F.; Clark, J. A.; Cleva, F.; Cocchieri, C.; Coccia, E.; Cohadon, P. -F.; Colla, A.; Collette, C. G.; Cominsky, L.; Constancio, M.; Conti, L.; Cooper, S. J.; Corbitt, T. R.; Cornish, N.; Corsi, A.; Cortese, S.; Costa, C. A.; Coughlin, M. W.; Coughlin, S. B.; Coulon, J. -P.; Countryman, S. T.; Couvares, P.; Covas, P. B.; Cowan, E. E.; Coward, D. M.; Cowart, M. J.; Coyne, D. C.; Coyne, R.; Creighton, J. D. E.; Creighton, T. D.; Cripe, J.; Crowder, S. G.; Cullen, T. J.; Cumming, A.; Cunningham, L.; Cuoco, E.; Canton, T. Dal; Danilishin, S. L.; D'Antonio, S.; Danzmann, K.; Dasgupta, A.; Costa, C. F. Da Silva; Dattilo, V.; Dave, I.; Davier, M.; Davies, G. S.; Davis, D.; Daw, E. J.; Day, B.; Day, R.; De, S.; DeBra, D.; Debreczeni, G.; Degallaix, J.; De Laurentis, M.; Del'eglise, S.; Del Pozzo, W.; Denker, T.; Dent, T.; Dergachev, V.; De Rosa, R.; DeRosa, R. T.; DeSalvo, R.; Devine, R. C.; Dhurandhar, S.; D'iaz, M. C.; Di Fiore, L.; Di Giovanni, M.; Di Girolamo, T.; Di Lieto, A.; Di Pace, S.; Di Palma, I.; Di Virgilio, A.; Doctor, Z.; Dolique, V.; Donovan, F.; Dooley, K. L.; Doravari, S.; Dorrington, I.; Douglas, R.; 'Alvarez, M. Dovale; Downes, T. P.; Drago, M.; Drever, R. W. P.; Driggers, J. C.; Du, Z.; Ducrot, M.; Dwyer, S. E.; Edo, T. B.; Edwards, M. C.; Effler, A.; Eggenstein, H. -B.; Ehrens, P.; Eichholz, J.; Eikenberry, S. S.; Eisenstein, R. A.; Essick, R. C.; Etienne, Z.; Etzel, T.; Evans, M.; Evans, T. M.; Everett, R.; Factourovich, M.; Fafone, V.; Fair, H.; Fairhurst, S.; Fan, X.; Farinon, S.; Farr, B.; Farr, W. M.; Fauchon-Jones, E. J.; Favata, M.; Fays, M.; Fehrmann, H.; Fejer, M. M.; Galiana, A. Fern'andez; Ferrante, I.; Ferreira, E. C.; Ferrini, F.; Fidecaro, F.; Fiori, I.; Fiorucci, D.; Fisher, R. P.; Flaminio, R.; Fletcher, M.; Fong, H.; Forsyth, S. S.; Fournier, J. -D.; Frasca, S.; Frasconi, F.; Frei, Z.; Freise, A.; Frey, R.; Frey, V.; Fries, E. M.; Fritschel, P.; Frolov, V. V.; Fulda, P.; Fyffe, M.; Gabbard, H.; Gadre, B. U.; Gaebel, S. M.; Gair, J. R.; Gammaitoni, L.; Gaonkar, S. G.; Garufi, F.; Gaur, G.; Gayathri, V.; Gehrels, N.; Gemme, G.; Genin, E.; Gennai, A.; George, J.; Gergely, L.; Germain, V.; Ghonge, S.; Ghosh, Abhirup; Ghosh, Archisman; Ghosh, S.; Giaime, J. A.; Giardina, K. D.; Giazotto, A.; Gill, K.; Glaefke, A.; Goetz, E.; Goetz, R.; Gondan, L.; Gonz'alez, G.; Castro, J. M. Gonzalez; Gopakumar, A.; Gorodetsky, M. L.; Gossan, S. E.; Gosselin, M.; Gouaty, R.; Grado, A.; Graef, C.; Granata, M.; Grant, A.; Gras, S.; Gray, C.; Greco, G.; Green, A. C.; Groot, P.; Grote, H.; Grunewald, S.; Guidi, G. M.; Guo, X.; Gupta, A.; Gupta, M. K.; Gushwa, K. E.; Gustafson, E. K.; Gustafson, R.; Hacker, J. J.; Hall, B. R.; Hall, E. D.; Hammond, G.; Haney, M.; Hanke, M. M.; Hanks, J.; Hanna, C.; Hanson, J.; Hardwick, T.; Harms, J.; Harry, G. M.; Harry, I. W.; Hart, M. J.; Hartman, M. T.; Haster, C. -J.; Haughian, K.; Healy, J.; Heidmann, A.; Heintze, M. C.; Heitmann, H.; Hello, P.; Hemming, G.; Hendry, M.; Heng, I. S.; Hennig, J.; Henry, J.; Heptonstall, A. W.; Heurs, M.; Hild, S.; Hoak, D.; Hofman, D.; Holt, K.; Holz, D. E.; Hopkins, P.; Hough, J.; Houston, E. A.; Howell, E. J.; Hu, Y. M.; Huerta, E. A.; Huet, D.; Hughey, B.; Husa, S.; Huttner, S. H.; Huynh-Dinh, T.; Indik, N.; Ingram, D. R.; Inta, R.; Isa, H. N.; Isac, J. -M.; Isi, M.; Isogai, T.; Iyer, B. R.; Izumi, K.; Jacqmin, T.; Jani, K.; Jaranowski, P.; Jawahar, S.; Jim'enez-Forteza, F.; Johnson, W. W.; Jones, D. I.; Jones, R.; Jonker, R. J. G.; Ju, L.; Junker, J.; Kalaghatgi, C. V.; Kalogera, V.; Kandhasamy, S.; Kang, G.; Kanner, J. B.; Karki, S.; Karvinen, K. S.; Kasprzack, M.; Katsavounidis, E.; Katzman, W.; Kaufer, S.; Kaur, T.; Kawabe, K.; K'ef'elian, F.; Keitel, D.; Kelley, D. B.; Kennedy, R.; Key, J. S.; Khalili, F. Y.; Khan, I.; Khan, S.; Khan, Z.; Khazanov, E. A.; Kijbunchoo, N.; Kim, Chunglee; Kim, J. C.; Kim, Whansun; Kim, W.; Kim, Y. -M.; Kimbrell, S. J.; King, E. J.; King, P. J.; Kirchhoff, R.; Kissel, J. S.; Klein, B.; Kleybolte, L.; Klimenko, S.; Koch, P.; Koehlenbeck, S. M.; Koley, S.; Kondrashov, V.; Kontos, A.; Korobko, M.; Korth, W. Z.; Kowalska, I.; Kozak, D. B.; Kr"amer, C.; Kringel, V.; Krishnan, B.; Kr'olak, A.; Kuehn, G.; Kumar, P.; Kumar, R.; Kuo, L.; Kutynia, A.; Lackey, B. D.; Landry, M.; Lang, R. N.; Lange, J.; Lantz, B.; Lanza, R. K.; Lartaux-Vollard, A.; Lasky, P. D.; Laxen, M.; Lazzarini, A.; Lazzaro, C.; Leaci, P.; Leavey, S.; Lebigot, E. O.; Lee, C. H.; Lee, H. K.; Lee, H. M.; Lee, K.; Lehmann, J.; Lenon, A.; Leonardi, M.; Leong, J. R.; Leroy, N.; Letendre, N.; Levin, Y.; Li, T. G. F.; Libson, A.; Littenberg, T. B.; Liu, J.; Lockerbie, N. A.; Lombardi, A. L.; London, L. T.; Lord, J. E.; Lorenzini, M.; Loriette, V.; Lormand, M.; Losurdo, G.; Lough, J. D.; Lousto, C. O.; Lovelace, G.; L"uck, H.; Lundgren, A. P.; Lynch, R.; Ma, Y.; Macfoy, S.; Machenschalk, B.; MacInnis, M.; Macleod, D. M.; Magana-Sandoval, F.; Majorana, E.; Maksimovic, I.; Malvezzi, V.; Man, N.; Mandic, V.; Mangano, V.; Mansell, G. L.; Manske, M.; Mantovani, M.; Marchesoni, F.; Marion, F.; M'arka, S.; M'arka, Z.; Markosyan, A. S.; Maros, E.; Martelli, F.; Martellini, L.; Martin, I. W.; Martynov, D. V.; Mason, K.; Masserot, A.; Massinger, T. J.; Masso-Reid, M.; Mastrogiovanni, S.; Matichard, F.; Matone, L.; Mavalvala, N.; Mazumder, N.; McCarthy, R.; McClelland, D. E.; McCormick, S.; McGrath, C.; McGuire, S. C.; McIntyre, G.; McIver, J.; McManus, D. J.; McRae, T.; McWilliams, S. T.; Meacher, D.; Meadors, G. D.; Meidam, J.; Melatos, A.; Mendell, G.; Mendoza-Gandara, D.; Mercer, R. A.; Merilh, E. L.; Merzougui, M.; Meshkov, S.; Messenger, C.; Messick, C.; Metzdorff, R.; Meyers, P. M.; Mezzani, F.; Miao, H.; Michel, C.; Middleton, H.; Mikhailov, E. E.; Milano, L.; Miller, A. L.; Miller, A.; Miller, B. B.; Miller, J.; Millhouse, M.; Minenkov, Y.; Ming, J.; Mirshekari, S.; Mishra, C.; Mitra, S.; Mitrofanov, V. P.; Mitselmakher, G.; Mittleman, R.; Moggi, A.; Mohan, M.; Mohapatra, S. R. P.; Montani, M.; Moore, B. C.; Moore, C. J.; Moraru, D.; Moreno, G.; Morriss, S. R.; Mours, B.; Mow-Lowry, C. M.; Mueller, G.; Muir, A. W.; Mukherjee, Arunava; Mukherjee, D.; Mukherjee, S.; Mukund, N.; Mullavey, A.; Munch, J.; Muniz, E. A. M.; Murray, P. G.; Mytidis, A.; Napier, K.; Nardecchia, I.; Naticchioni, L.; Nelemans, G.; Nelson, T. J. N.; Neri, M.; Nery, M.; Neunzert, A.; Newport, J. M.; Newton, G.; Nguyen, T. T.; Nielsen, A. B.; Nissanke, S.; Nitz, A.; Noack, A.; Nocera, F.; Nolting, D.; Normandin, M. E. N.; Nuttall, L. K.; Oberling, J.; Ochsner, E.; Oelker, E.; Ogin, G. H.; Oh, J. J.; Oh, S. H.; Ohme, F.; Oliver, M.; Oppermann, P.; Oram, Richard J.; O'Reilly, B.; O'Shaughnessy, R.; Ottaway, D. J.; Overmier, H.; Owen, B. J.; Pace, A. E.; Page, J.; Pai, A.; Pai, S. A.; Palamos, J. R.; Palashov, O.; Palomba, C.; Pal-Singh, A.; Pan, H.; Pankow, C.; Pannarale, F.; Pant, B. C.; Paoletti, F.; Paoli, A.; Papa, M. A.; Paris, H. R.; Parker, W.; Pascucci, D.; Pasqualetti, A.; Passaquieti, R.; Passuello, D.; Patricelli, B.; Pearlstone, B. L.; Pedraza, M.; Pedurand, R.; Pekowsky, L.; Pele, A.; Penn, S.; Perez, C. J.; Perreca, A.; Perri, L. M.; Pfeiffer, H. P.; Phelps, M.; Piccinni, O. J.; Pichot, M.; Piergiovanni, F.; Pierro, V.; Pillant, G.; Pinard, L.; Pinto, I. M.; Pitkin, M.; Poe, M.; Poggiani, R.; Popolizio, P.; Post, A.; Powell, J.; Prasad, J.; Pratt, J. W. W.; Predoi, V.; Prestegard, T.; Prijatelj, M.; Principe, M.; Privitera, S.; Prix, R.; Prodi, G. A.; Prokhorov, L. G.; Puncken, O.; Punturo, M.; Puppo, P.; P"urrer, M.; Qi, H.; Qin, J.; Qiu, S.; Quetschke, V.; Quintero, E. A.; Quitzow-James, R.; Raab, F. J.; Rabeling, D. S.; Radkins, H.; Raffai, P.; Raja, S.; Rajan, C.; Rakhmanov, M.; Rapagnani, P.; Raymond, V.; Razzano, M.; Re, V.; Read, J.; Regimbau, T.; Rei, L.; Reid, S.; Reitze, D. H.; Rew, H.; Reyes, S. D.; Rhoades, E.; Ricci, F.; Riles, K.; Rizzo, M.; Robertson, N. A.; Robie, R.; Robinet, F.; Rocchi, A.; Rolland, L.; Rollins, J. G.; Roma, V. J.; Romano, R.; Romie, J. H.; Rosi'nska, D.; Rowan, S.; R"udiger, A.; Ruggi, P.; Ryan, K.; Sachdev, S.; Sadecki, T.; Sadeghian, L.; Sakellariadou, M.; Salconi, L.; Saleem, M.; Salemi, F.; Samajdar, A.; Sammut, L.; Sampson, L. M.; Sanchez, E. J.; Sandberg, V.; Sanders, J. R.; Sassolas, B.; Sathyaprakash, B. S.; Saulson, P. R.; Sauter, O.; Savage, R. L.; Sawadsky, A.; Schale, P.; Scheuer, J.; Schmidt, E.; Schmidt, J.; Schmidt, P.; Schnabel, R.; Schofield, R. M. S.; Sch"onbeck, A.; Schreiber, E.; Schuette, D.; Schutz, B. F.; Schwalbe, S. G.; Scott, J.; Scott, S. M.; Sellers, D.; Sengupta, A. S.; Sentenac, D.; Sequino, V.; Sergeev, A.; Setyawati, Y.; Shaddock, D. A.; Shaffer, T. J.; Shahriar, M. S.; Shapiro, B.; Shawhan, P.; Sheperd, A.; Shoemaker, D. H.; Shoemaker, D. M.; Siellez, K.; Siemens, X.; Sieniawska, M.; Sigg, D.; Silva, A. D.; Singer, A.; Singer, L. P.; Singh, A.; Singh, R.; Singhal, A.; Sintes, A. M.; Slagmolen, B. J. J.; Smith, B.; Smith, J. R.; Smith, R. J. E.; Son, E. J.; Sorazu, B.; Sorrentino, F.; Souradeep, T.; Spencer, A. P.; Srivastava, A. K.; Staley, A.; Steinke, M.; Steinlechner, J.; Steinlechner, S.; Steinmeyer, D.; Stephens, B. C.; Stevenson, S. P.; Stone, R.; Strain, K. A.; Straniero, N.; Stratta, G.; Strigin, S. E.; Sturani, R.; Stuver, A. L.; Summerscales, T. Z.; Sun, L.; Sunil, S.; Sutton, P. J.; Swinkels, B. L.; Szczepa'nczyk, M. J.; Tacca, M.; Talukder, D.; Tanner, D. B.; T'apai, M.; Taracchini, A.; Taylor, R.; Theeg, T.; Thomas, E. G.; Thomas, M.; Thomas, P.; Thorne, K. A.; Thrane, E.; Tippens, T.; Tiwari, S.; Tiwari, V.; Tokmakov, K. V.; Toland, K.; Tomlinson, C.; Tonelli, M.; Tornasi, Z.; Torrie, C. I.; T"oyr"a, D.; Travasso, F.; Traylor, G.; Trifir`o, D.; Trinastic, J.; Tringali, M. C.; Trozzo, L.; Tse, M.; Tso, R.; Turconi, M.; Tuyenbayev, D.; Ugolini, D.; Unnikrishnan, C. S.; Urban, A. L.; Usman, S. A.; Vahlbruch, H.; Vajente, G.; Valdes, G.; van Bakel, N.; van Beuzekom, M.; Brand, J. F. J. van den; Broeck, C. Van Den; Vander-Hyde, D. C.; van der Schaaf, L.; van Heijningen, J. V.; van Veggel, A. A.; Vardaro, M.; Varma, V.; Vass, S.; Vas'uth, M.; Vecchio, A.; Vedovato, G.; Veitch, J.; Veitch, P. J.; Venkateswara, K.; Venugopalan, G.; Verkindt, D.; Vetrano, F.; Vicer'e, A.; Viets, A. D.; Vinciguerra, S.; Vine, D. J.; Vinet, J. -Y.; Vitale, S.; Vo, T.; Vocca, H.; Vorvick, C.; Voss, D. V.; Vousden, W. D.; Vyatchanin, S. P.; Wade, A. R.; Wade, L. E.; Wade, M.; Walker, M.; Wallace, L.; Walsh, S.; Wang, G.; Wang, H.; Wang, M.; Wang, Y.; Ward, R. L.; Warner, J.; Was, M.; Watchi, J.; Weaver, B.; Wei, L. -W.; Weinert, M.; Weinstein, A. J.; Weiss, R.; Wen, L.; Wessels, P.; Westphal, T.; Wette, K.; Whelan, J. T.; Whiting, B. F.; Whittle, C.; Williams, D.; Williams, R. D.; Williamson, A. R.; Willis, J. L.; Willke, B.; Wimmer, M. H.; Winkler, W.; Wipf, C. C.; Wittel, H.; Woan, G.; Woehler, J.; Worden, J.; Wright, J. L.; Wu, D. S.; Wu, G.; Yam, W.; Yamamoto, H.; Yancey, C. C.; Yap, M. J.; Yu, Hang; Yu, Haocun; Yvert, M.; Zadrożny, A.; Zangrando, L.; Zanolin, M.; Zendri, J. -P.; Zevin, M.; Zhang, L.; Zhang, M.; Zhang, T.; Zhang, Y.; Zhao, C.; Zhou, M.; Zhou, Z.; Zhu, S. J.; Zhu, X. J.; Zucker, M. E.; Zweizig, J.; Buchner, S.; Cognard, I.; Corongiu, A.; Freire, P. C. C.; Guillemot, L.; Hobbs, G. B.; Kerr, M.; Lyne, A. G.; Possenti, A.; Ridolfi, A.; Shannon, R. M.; Stappers, B. W.; Weltevrede, P.
Comments: 19 pages, 4 tables, 4 figures, accepted by The Astrophysical Journal. Fixes Equation 4. Science summary of results available at http://www.ligo.org/science/Publication-O1KnownPulsar/index.php
Submitted: 2017-01-26, last modified: 2017-11-13
We present the result of searches for gravitational waves from 200 pulsars using data from the first observing run of the Advanced LIGO detectors. We find no significant evidence for a gravitational-wave signal from any of these pulsars, but we are able to set the most constraining upper limits yet on their gravitational-wave amplitudes and ellipticities. For eight of these pulsars, our upper limits give bounds that are improvements over the indirect spin-down limit values. For another 32, we are within a factor of 10 of the spin-down limit, and it is likely that some of these will be reachable in future runs of the advanced detector. Taken as a whole, these new results improve on previous limits by more than a factor of two.
[55]  oai:arXiv.org:1710.05846  [pdf] - 1608506
Follow up of GW170817 and its electromagnetic counterpart by Australian-led observing programs
Andreoni, I.; Ackley, K.; Cooke, J.; Acharyya, A.; Allison, J. R.; Anderson, G. E.; Ashley, M. C. B.; Baade, D.; Bailes, M.; Bannister, K.; Beardsley, A.; Bessell, M. S.; Bian, F.; Bland, P. A.; Boer, M.; Booler, T.; Brandeker, A.; Brown, I. S.; Buckley, D.; Chang, S. -W.; Coward, D. M.; Crawford, S.; Crisp, H.; Crosse, B.; Cucchiara, A.; Cupák, M.; de Gois, J. S.; Deller, A.; Devillepoix, H. A. R.; Dobie, D.; Elmer, E.; Emrich, D.; Farah, W.; Farrell, T. J.; Franzen, T.; Gaensler, B. M.; Galloway, D. K.; Gendre, B.; Giblin, T.; Goobar, A.; Green, J.; Hancock, P. J.; Hartig, B. A. D.; Howell, E. J.; Horsley, L.; Hotan, A.; Howie, R. M.; Hu, L.; Hu, Y.; James, C. W.; Johnston, S.; Johnston-Hollitt, M.; Kaplan, D. L.; Kasliwal, M.; Keane, E. F.; Kenney, D.; Klotz, A.; Lau, R.; Laugier, R.; Lenc, E.; Li, X.; Liang, E.; Lidman, C.; Luvaul, L. C.; Lynch, C.; Ma, B.; Macpherson, D.; Mao, J.; McClelland, D. E.; McCully, C.; Möller, A.; Morales, M. F.; Morris, D.; Murphy, T.; Noysena, K.; Onken, C. A.; Orange, N. B.; Oslowski, S.; Pallot, D.; Paxman, J.; Potter, S. B.; Pritchard, T.; Raja, W.; Ridden-Harper, R.; Romero-Colmenero, E.; Sadler, E. M.; Sansom, E. K.; Scalzo, R. A.; Schmidt, B. P.; Scott, S. M.; Seghouani, N.; Shang, Z.; Shannon, R. M.; Shao, L.; Shara, M. M.; Sharp, R.; Sokolowski, M.; Sollerman, J.; Staff, J.; Steele, K.; Sun, T.; Suntzeff, N. B.; Tao, C.; Tingay, S.; Towner, M. C.; Thierry, P.; Trott, C.; Tucker, B. E.; Väisänen, P.; Krishnan, V. Venkatraman; Walker, M.; Wang, L.; Wang, X.; Wayth, R.; Whiting, M.; Williams, A.; Williams, T.; Wolf, C.; Wu, C.; Wu, X.; Yang, J.; Yuan, X.; Zhang, H.; Zhou, J.; Zovaro, H.
Comments: 26 pages, 9 figures, 15 tables
Submitted: 2017-10-16
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.
[56]  oai:arXiv.org:1709.03651  [pdf] - 1600291
Spectral Flattening at Low Frequencies in Crab Giant Pulses
Comments: 21 pages, 12 figures, 8 tables, accepted for publication in ApJ
Submitted: 2017-09-11
We report on simultaneous wideband observations of Crab giant pulses with the Parkes radio telescope and the Murchison Widefield Array (MWA). The observations were conducted simultaneously at 732 and 3100 MHz with Parkes, and at 120.96, 165.76 and 210.56 MHz with the MWA. Flux density calibration of the MWA data was accomplished using a novel technique based on tied-array beam simulations. We detected between 90-648 giant pulses in the 120.96-210.56 MHz MWA subbands above a $ 5.5\sigma $ threshold while in the Parkes subbands we detected 6344 and 231 giant pulses above a threshold of $ 6\sigma $ at 732 and 3100 MHz, respectively. We show, for the first time over a wide frequency range, that the average spectrum of Crab giant pulses exhibits a significant flattening at low frequencies. The spectral index, $ \alpha $, for giant pulses evolves from a steep, narrow distribution with a mean $ \alpha=-2.6 $ and width $\sigma_\alpha=0.5 $ between 732 and 3100 MHz, to a wide, flat distribution of spectral indices with a mean $ \alpha=-0.7 $ and width $ \sigma_\alpha=1.4 $ between 120.96 and 165.76 MHz. We also comment on the plausibility of giant pulse models for Fast Radio Bursts based on this spectral information.
[57]  oai:arXiv.org:1705.07581  [pdf] - 1583683
The detection of an extremely bright fast radio burst in a phased array feed survey
Comments: 8 pages, 4 figures, accepted in ApJ letters. Version 2: Fixed galactic coordinates in Table 2
Submitted: 2017-05-22, last modified: 2017-05-23
We report the detection of an ultra-bright fast radio burst (FRB) from a modest, 3.4-day pilot survey with the Australian Square Kilometre Array Pathfinder. The survey was conducted in a wide-field fly's-eye configuration using the phased-array-feed technology deployed on the array to instantaneously observe an effective area of $160$ deg$^2$, and achieve an exposure totaling $13200$ deg$^2$ hr. We constrain the position of FRB 170107 to a region $8'\times8'$ in size (90% containment) and its fluence to be $58\pm6$ Jy ms. The spectrum of the burst shows a sharp cutoff above $1400$ MHz, which could be either due to scintillation or an intrinsic feature of the burst. This confirms the existence of an ultra-bright ($>20$ Jy ms) population of FRBs.
[58]  oai:arXiv.org:1704.00460  [pdf] - 1561251
Chandra Phase-Resolved Spectroscopy of the High Magnetic Field Pulsar B1509-58
Comments: 10 pages, 10 figures, 3 tables, accepted for publication in ApJ
Submitted: 2017-04-03
We report on a timing and spectral analysis of the young, high magnetic field rotation-powered pulsar (RPP) B1509-58 using Chandra continuous-clocking mode observation. The pulsar's X-ray light curve can be fit by the two Gaussian components and the pulsed fraction shows moderate energy dependence over the Chandra band. The pulsed X-ray spectrum is well described by a power law with a photon index 1.16(4), which is harder than the values measured with RXTE PCA and NuSTAR. This result supports the log-parabolic model for the broadband X-ray spectrum. With the unprecedented angular resolution of Chandra, we clearly identified off-pulse X-ray emission from the pulsar, and its spectrum is best fit by a power law plus blackbody model. The latter component has a temperature of ~0.14 keV with a bolometric luminosity comparable to the luminosities of other young and high magnetic field RPPs, and it lies between the temperature of magnetars and typical RPPs. In addition, we found that the nonthermal X-ray emission of PSR B1509-58 is significantly softer in the off-pulse phase than in the pulsed phase, with the photon index varying between 1.0 and 1.8 and anticorrelated with the flux. This is similar to the behavior of three other young pulsars. We interpreted it as different contributions of pair-creation processes at different altitudes from the neutron star surface according to the outer-gap model.
[59]  oai:arXiv.org:1704.01011  [pdf] - 1582337
Comparison of Pulsar Positions from Timing and Very Long Baseline Astrometry
Comments: accepted by MNRAS,22 pages, 5 figures
Submitted: 2017-03-31
Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radio telescopes in order to relate the two reference frames. We find that the timing positions differ significantly from the VLBI positions (and also differ between different ephemerides). A statistically significant change in the obliquity of the ecliptic of $2.16\pm0.33$\,mas is found for the JPL ephemeris DE405, but no significant rotation is found in subsequent JPL ephemerides. The accuracy with which we can relate the two frames is limited by the current uncertainties in the VLBI reference source positions and in matching the pulsars to their reference source. Not only do the timing positions depend on the ephemeris used in computing them, but also different segments of the timing data lead to varying position estimates. These variations are mostly common to all ephemerides, but slight changes are seen at the 10$\mu$as level between ephemerides.
[60]  oai:arXiv.org:1703.02108  [pdf] - 1561204
Robust Estimation of Scattering in Pulsar Timing Analysis
Comments: 12 pages, 6 figures, Accepted to MNRAS
Submitted: 2017-03-06
We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsar's timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4\sigma$ level, for example, changing proper motion in right ascension by $0.50 \pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited.
[61]  oai:arXiv.org:1612.05258  [pdf] - 1533476
Wide-band Profile Domain Pulsar Timing Analysis
Comments: 22 pages, submitted to MNRAS
Submitted: 2016-12-15
We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $\sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band.
[62]  oai:arXiv.org:1610.01731  [pdf] - 1531913
The NANOGrav Nine-Year Data Set: Excess Noise in Millisecond Pulsar Arrival Times
Comments: 20 pages, 16 figures, to be published in ApJ
Submitted: 2016-10-06, last modified: 2016-12-02
Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon \& Cordes (2010). We briefly discuss our results in terms of detection of gravitational waves at nanohertz frequencies.
[63]  oai:arXiv.org:1611.05758  [pdf] - 1532687
The magnetic field and turbulence of the cosmic web measured using a brilliant fast radio burst
Comments: Published online in Science First Release. 41 pages including the mean text and supplementary material
Submitted: 2016-11-17
Fast radio bursts (FRBs) are millisecond-duration events thought to originate beyond the Milky Way galaxy. Uncertainty surrounding the burst sources, and their propagation through intervening plasma, has limited their use as cosmological probes. We report on a mildly dispersed (dispersion measure 266.5+-0.1 pc cm^-3), exceptionally intense (120+-30 Jy), linearly polarized, scintillating burst (FRB 150807) that we directly localize to 9 arcmin^2. Based on a low Faraday rotation (12.0+-0.7 rad m^-2), we infer negligible magnetization in the circum-burst plasma and constrain the net magnetization of the cosmic web along this sightline to <21 nG, parallel to the line-of-sight. The burst scintillation suggests weak turbulence in the ionized intergalactic medium.
[64]  oai:arXiv.org:1611.05580  [pdf] - 1546299
Radio-interferometric monitoring of FRB 131104: A coincident AGN flare, but no evidence for a cosmic fireball
Comments: 8 pages, 3 Figures, submitted to ApJ Letters
Submitted: 2016-11-17
The localization of fast radio bursts (FRBs) has been hindered by the poor angular resolution of the detection observations and inconclusive identification of transient or variable counterparts. Recently a $\gamma$-ray pulse of $380$ s duration has been associated with the fast radio burst FRB 131104. We report on radio-continuum imaging observations of the original localization region of the FRB, beginning three days after the event and comprising 25 epochs over 2.5 yr. Besides probabilistic arguments that suggest that the association between the $\gamma$-ray transient and the FRB is not compelling, we provide upper limits on a putative radio afterglow of this transient that are at odds with standard models for its progenitor. We further report the discovery of an unusual variable radio source spatially and temporally coincident with FRB 131104, but not spatially coincident with the $\gamma$-ray event. The radio variable flares by a factor of $3$ above its long term average within $10$ d of the FRB at 7.5 GHz, with a factor-of-$2$ increase at $5.5$ GHz. Within our observations, the variable has persisted with only modest modulation and never approached the flux density observed in the days after the FRB. We identify an optical counterpart to the variable. Optical and infrared photometry, and deep optical spectroscopy, suggest that the object is a narrow-line radio AGN.
[65]  oai:arXiv.org:1609.07573  [pdf] - 1513706
Modelling and mitigating refractive propagation effects in precision pulsar timing observations
Comments: 16 pages, 8 figures. Accepted for publication in MNRAS
Submitted: 2016-09-24
To obtain the most accurate pulse arrival times from radio pulsars, it is necessary to correct or mitigate the effects of the propagation of radio waves through the warm and ionised interstellar medium. We examine both the strength of propagation effects associated with large-scale electron-density variations and the methodology used to estimate infinite-frequency arrival times. Using simulations of two-dimensional phase-varying screens, we assess the strength and non-stationarity of timing perturbations associated with large-scale density variations. We identify additional contributions to arrival times that are stochastic in both radio frequency and time and therefore not amenable to correction solely using times of arrival. We attribute this to the frequency dependence of the trajectories of the propagating radio waves. We find that this limits the efficacy of low-frequency (metre-wavelength) observations. Incorporating low-frequency pulsar observations into precision timing campaigns is increasingly problematic for pulsars with larger dispersion measures.
[66]  oai:arXiv.org:1509.05446  [pdf] - 1941997
Single-Source Gravitational Wave Limits from the J1713+0747 24-hr Global Campaign
Comments: 4 pages, 2 figures, published in the proceedings of Amaldi 11
Submitted: 2015-09-17, last modified: 2016-08-20
Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.
[67]  oai:arXiv.org:1608.02163  [pdf] - 1470812
The disturbance of a millisecond pulsar magnetosphere
Comments: Accepted to ApJ Letters
Submitted: 2016-08-06
Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability.
[68]  oai:arXiv.org:1606.06795  [pdf] - 1498310
The Fluence and Distance Distributions of Fast Radio Bursts
Comments: In press (ApJ)
Submitted: 2016-06-21
Fast radio bursts (FRB) are millisecond-duration radio pulses with apparent extragalactic origins. All but two of the FRBs have been discovered using the Parkes dish which employs multiple beams formed by an array of feed horns on its focal plane. In this paper, we show that (i) the preponderance of multiple-beam detections, and (ii) the detection rates for varying dish diameters, can be used to infer the index $\alpha$ of the cumulative fluence distribution function (the log$N$-log$F$ function: $\alpha=1.5$ for a non-evolving population in a Euclidean universe). If all detected FRBs arise from a single progenitor population, multiple-beam FRB detection rates from the Parkes telescope yield the constraint $0.52<\alpha<1.0$ with $90$% confidence. Searches at other facilities with different dish sizes refine the constraint to $0.66<\alpha<0.96$. Our results favor FRB searches with smaller dishes, because for $\alpha<1$, the gain in field-of-view for a smaller dish is more important than the reduction in sensitivity. Further, our results suggest that (i) FRBs are not standard candles, and (ii) the distribution of distances to the detected FRBs is weighted towards larger distances. If FRBs are extragalactic, these results are consistent with a cosmological population, which would make FRBs excellent probes of the baryonic content and geometry of the Universe.
[69]  oai:arXiv.org:1604.05322  [pdf] - 1393032
Radio polarimetry of Galactic centre pulsars
Comments: 7 pages, 3 figures. Accepted for publication in MNRAS
Submitted: 2016-04-18
To study the strength and structure of the magnetic field in the Galactic centre (GC) we measured Faraday rotation of the radio emission of pulsars which are seen towards the GC. Three of these pulsars have the largest rotation measures (RMs) observed in any Galactic object with the exception of Sgr A*. Their large dispersion measures, RMs and the large RM variation between these pulsars and other known objects in the GC implies that the pulsars lie in the GC and are not merely seen in projection towards the GC. The large RMs of these pulsars indicate large line-of-sight magnetic field components between ~ 16-33 microgauss; combined with recent model predictions for the strength of the magnetic field in the GC this implies that the large-scale magnetic field has a very small inclination angle with respect to the plane of the sky (~ 12 degrees). Foreground objects like the Radio Arc or possibly an ablated, ionized halo around the molecular cloud G0.11-0.11 could contribute to the large RMs of two of the pulsars. If these pulsars lie behind the Radio Arc or G0.11-0.11 then this proves that low-scattering corridors with lengths >~ 100 pc must exist in the GC. This also suggests that future, sensitive observations will be able to detect additional pulsars in the GC. Finally, we show that the GC component in our most accurate electron density model oversimplifies structure in the GC.
[70]  oai:arXiv.org:1604.02661  [pdf] - 1397005
Characterising the rotational irregularities of the Vela pulsar from 21 yr of phase-coherent timing
Comments: Accepted by MNRAS. 10 pages, 1 Figure, 2 tables
Submitted: 2016-04-10
Pulsars show two classes of rotational irregularities that can be used to understand neutron-star interiors and magnetospheres: glitches and timing noise. Here we present an analysis of the Vela pulsar spanning nearly 21 yr of observation and including 8 glitches. We identify the relative pulse number of all of the observations between glitches, with the only pulse-number ambiguities existing over glitch events. We use the phase coherence of the timing solution to simultaneously model the timing noise and glitches in a Bayesian framework, allowing us to select preferred models for both. We find the glitches can be described using only permanent and transient changes in spin frequency, i.e., no step changes in frequency derivative. For all of the glitches, we only need two exponentially decaying changes in spin frequency to model the transient components. In contrast to previous studies, we find that the dominant transient components decay on a common $\approx$ 1300 d time scale, and that a larger fraction ( $\gtrsim 25\%$) of glitch amplitudes are associated with these transient components. We also detect shorter-duration transient components of $\approx$ 25 d, as previously observed, but are limited in sensitivity to events with shorter durations by the cadence of our observations. The timing noise is well described by a steep power-law process that is independent of the glitches and subdominant to the glitch recovery. The braking index is constrained to be $<$ 8 with 95% confidence. This methodology can be used to robustly measure the properties of glitches and timing noise in other pulsars.
[71]  oai:arXiv.org:1511.05994  [pdf] - 1385411
Gravitational-wave cosmology across 29 decades in frequency
Comments: accepted for publication in Physical Review X
Submitted: 2015-11-18, last modified: 2016-02-28
Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.
[72]  oai:arXiv.org:1602.05570  [pdf] - 1370407
From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release
Comments: 29 pages. 16 figures. Accepted for publication in MNRAS
Submitted: 2016-02-16
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.
[73]  oai:arXiv.org:1602.03640  [pdf] - 1362729
The International Pulsar Timing Array: First Data Release
Comments: 25 pages, 6 tables, 5 figures. Accepted for publication in MNRAS
Submitted: 2016-02-11
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.
[74]  oai:arXiv.org:1602.02872  [pdf] - 1359306
Gravitational wave astronomy: the current status
Comments: published (2015)
Submitted: 2016-02-09
In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry.
[75]  oai:arXiv.org:1512.02702  [pdf] - 1342986
A pilot ASKAP survey of radio transient events in the region around the intermittent pulsar PSR J1107-5907
Comments: Accepted by MNRAS
Submitted: 2015-12-08
We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.
[76]  oai:arXiv.org:1511.05481  [pdf] - 1347636
Emission-rotation correlation in pulsars: new discoveries with optimal techniques
Comments: 24 pages, 25 figures
Submitted: 2015-11-17, last modified: 2015-11-18
Pulsars are known to display short-term variability. Recently, examples of longer-term emission variability have emerged that are often correlated with changes in the rotational properties of the pulsar. To further illuminate this relationship, we have developed techniques to identify emission and rotation variability in pulsar data, and determine correlations between the two. Individual observations may be too noisy to identify subtle changes in the pulse profile. We use Gaussian process (GP) regression to model noisy observations and produce a continuous map of pulse profile variability. Generally, multiple observing epochs are required to obtain the pulsar spin frequency derivative ($\dot{\nu}$). GP regression is, therefore, also used to obtain $\dot{\nu}$, under the hypothesis that pulsar timing noise is primarily caused by unmodelled changes in $\dot{\nu}$. Our techniques distinguish between two types of variability: changes in the total flux density versus changes in the pulse shape. We have applied these techniques to 168 pulsars observed by the Parkes radio telescope, and see that although variations in flux density are ubiquitous, substantial changes in the shape of the pulse profile are rare. We reproduce previously published results and present examples of profile shape changing in seven pulsars; in particular, a clear new example of correlated changes in profile shape and rotation is found in PSR~J1602$-$5100. In the shape changing pulsars, a more complex picture than the previously proposed two state model emerges. We conclude that our simple assumption that all timing noise can be interpreted as $\dot{\nu}$ variability is insufficient to explain our dataset.
[77]  oai:arXiv.org:1510.08068  [pdf] - 1323551
Versatile Directional Searches for Gravitational Waves with Pulsar Timing Arrays
Comments: Accepted for publication by MNRAS
Submitted: 2015-10-27
By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "phased-up" to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise.
[78]  oai:arXiv.org:1510.06078  [pdf] - 1313354
Periodic modulation in pulse arrival times from young pulsars: a renewed case for neutron star precession
Comments: 11 pages, 8 figures, 3 tables, accepted in MNRAS
Submitted: 2015-10-20
In a search for periodic variation in the arrival times of pulses from 151 young, energetic pulsars, we have identified seven cases of modulation consistent with one or two harmonics of a single fundamental with time-scale 0.5-1.5 yr. We use simulations to show that these modulations are statistically significant and of high quality (sinusoidal) even when contaminated by the strong stochastic timing noise common to young pulsars. Although planetary companions could induce such modulation, the large implied masses and 2:1 mean motion resonances challenge such an explanation. Instead, the modulation is likely to be intrinsic to the pulsar, arising from quasi-periodic switching between stable magnetospheric states, and we propose that precession of the neutron star may regulate this switching.
[79]  oai:arXiv.org:1510.05099  [pdf] - 1319682
Timing Gamma-ray Pulsars with the Fermi Large Area Telescope: Timing Noise and Astrometry
Comments: 13 pages, 11 figures, 2 tables; accepted in ApJ
Submitted: 2015-10-17
We have constructed timing solutions for 81 gamma-ray pulsars covering more than five years of Fermi data. The sample includes 37 radio-quiet or radio-faint pulsars which cannot be timed with other telescopes. These timing solutions and the corresponding pulse times of arrival are prerequisites for further study, e.g. phase-resolved spectroscopy or searches for mode switches. Many gamma-ray pulsars are strongly affected by timing noise, and we present a new method for characterizing the noise process and mitigating its effects on other facets of the timing model. We present an analysis of timing noise over the population using a new metric for characterizing its strength and spectral shape, namely its time-domain correlation. The dependence of the strength on spin frequency and spin-down rate is in good agreement with previous studies. We find that noise process power spectra $S(f)$ for unrecycled pulsars are steep, with strong correlations over our entire data set and spectral indices $S(f)\propto f^{-\alpha}$ of 5 to 9. One possible explanation for these results is the occurrence of unmodelled, episodic 'microglitches'. Finally, we show that our treatment of timing noise results in robust parameter estimation, and in particular we measure a precise timing position for each pulsar. We extensively validate our results with multi-wavelength astrometry, and using our updated position, we firmly identify the X-ray counterpart of PSR J1418-6058.
[80]  oai:arXiv.org:1510.04434  [pdf] - 1294139
Timing analysis for 20 millisecond pulsars in the Parkes Pulsar Timing Array
Comments: 21 pages, 5 figures, 7 tables. Accepted for publication in MNRAS
Submitted: 2015-10-15
We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with $M_p=1.44\pm0.07$ $M_\odot$ and $M_p=1.47\pm0.03$ $M_\odot$ respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16% level of precision, $D=156.79\pm0.25$ pc, one of the most fractionally precise distance measurements of any star to date.
[81]  oai:arXiv.org:1510.02363  [pdf] - 1327503
A study of spatial correlations in pulsar timing array data
Comments: 13 pages, 7 figures. Accepted for publication in MNRAS
Submitted: 2015-10-08
Pulsar timing array experiments search for phenomena that produce angular correlations in the arrival times of signals from millisecond pulsars. The primary goal is to detect an isotropic and stochastic gravitational wave background. We use simulated data to show that this search can be affected by the presence of other spatially correlated noise, such as errors in the reference time standard, errors in the planetary ephemeris, the solar wind and instrumentation issues. All these effects can induce significant false detections of gravitational waves. We test mitigation routines to account for clock errors, ephemeris errors and the solar wind. We demonstrate that it is non-trivial to find an effective mitigation routine for the planetary ephemeris and emphasise that other spatially correlated signals may be present in the data.
[82]  oai:arXiv.org:1509.07320  [pdf] - 1283271
Gravitational waves from binary supermassive black holes missing in pulsar observations
Comments: Published in Science 25 September 2015
Submitted: 2015-09-24
Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.
[83]  oai:arXiv.org:1509.07505  [pdf] - 1283277
Profile Stochasticity in PSR J1909-3744
Comments: 16 pages, 13 figures, 2 tables
Submitted: 2015-09-24
We extend the recently introduced Bayesian framework `Generative Pulsar Timing Analysis' to incorporate both pulse jitter (high frequency variation in the arrival time of the pulse) and epoch to epoch stochasticity in the shape of the pulse profile. This framework allows for a full timing analysis to be performed on the folded profile data, rather than the site arrival times as is typical in most timing studies. We apply this extended framework both to simulations, and to an 11 yr, 10 cm data set for PSR J1909$-$3744. Using simulations, we show that temporal profile variation can induce timing noise in the residuals that when performing a standard timing analysis is highly covariant with the signal expected from a gravitational wave (GW) background. When working in the profile domain, these variations are de-correlated from the expected GW signal, resulting in significant improvement in the obtained upper limits. Using the PSR J1909$-$3744 data set from the Parkes Pulsar Timing Array project, we find significant evidence for systematic high-frequency profile variation resulting from non-Gaussian noise in the oldest observing system, but no evidence for either detectable pulse jitter, or low-frequency profile shape variation. Using our profile domain framework we therefore obtain upper limits on a red noise process with a spectral index of $\gamma = 13/3$ of $1\times10^{-15}$, consistent with previously published limits.
[84]  oai:arXiv.org:1509.06438  [pdf] - 1374107
Detecting nanohertz gravitational waves with pulsar timing arrays
Comments: 11 pages, 5 figures, review article for the special issue of SCIENCE CHINA Physics, Mechanics & Astronomy to mark the Centenary of General Relativity
Submitted: 2015-09-21
Complementary to ground-based laser interferometers, pulsar timing array experiments are being carried out to search for nanohertz gravitational waves. Using the world's most powerful radio telescopes, three major international collaborations have collected $\sim$10-year high precision timing data for tens of millisecond pulsars. In this paper we give an overview on pulsar timing experiments, gravitational wave detection in the nanohertz regime, and recent results obtained by various timing array projects.
[85]  oai:arXiv.org:1507.06982  [pdf] - 1264177
Limits on Planet Formation Around Young Pulsars and Implications for Supernova Fallback Disks
Comments: 5 pages, 4 figures, accepted to ApJL
Submitted: 2015-07-23
We have searched a sample of 151 young, energetic pulsars for periodic variation in pulse time-of-arrival arising from the influence of planetary companions. We are sensitive to objects with masses two orders of magnitude lower than those detectable with optical transit timing, but we find no compelling evidence for pulsar planets. For the older pulsars most likely to host planets, we can rule out Mercury analogues in one third of our sample and planets with masses $>0.4 M_{\oplus}$ and periods $P_b<1$ yr in all but 5% of such systems. If pulsar planets form primarily from supernova fallback disks, these limits imply that such disks do not form, are confined to $<0.1$ AU radii, are disrupted, or form planets more slowly ($>2$ Myr) than their protoplanetary counterparts.
[86]  oai:arXiv.org:1506.07948  [pdf] - 1259085
Pulsar Observations of Extreme Scattering Events
Comments: Accepted by ApJ
Submitted: 2015-06-25
Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the electron density of order 10s of cm^{-3}. Similar structures have been observed in measurements of pulsar intensity scintillation and group delay. Here we report observations of two ESEs showing increases in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). These allow us to make more complete models of the ESE, including an estimate of the "outer-scale" of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully turbulent on an AU scale. They provide some support for the idea that the structures are extended along the line of sight, such as would be the case for a scattering shell. The dispersion measurements also show a variety of AU scale structures which would not be called ESEs, yet involve electron density variations typical of ESEs and likely have the same origin.
[87]  oai:arXiv.org:1503.08491  [pdf] - 1344168
Frequency-Dependent Dispersion Measures and Implications for Pulsar Timing
Comments: 19 pages, 13 figures; submitted to the Astrophysical Journal
Submitted: 2015-03-29
We analyze the frequency dependence of the dispersion measure (DM), the column density of free electrons to a pulsar, caused by multipath scattering from small scale electron-density fluctuations. The DM is slightly different along each propagation path and the transverse spread of paths varies greatly with frequency, yielding time-of-arrival (TOA) perturbations that scale differently than the inverse square of the frequency, the expected dependence for a cold, unmagnetized plasma. We quantify DM and TOA perturbations analytically for thin phase screens and extended media and verify the results with simulations of thin screens. The rms difference between DMs across an octave band near 1.5~GHz $\sim 4\times10^{-5}\,{\rm pc\ cm^{-3}}$ for pulsars at $\sim 1$~kpc distance. TOA errors from chromatic DMs are of order a few to hundreds of nanoseconds for pulsars with DM $\lesssim 30$~pc~cm$^{-3}$ observed across an octave band but increase rapidly to microseconds or larger for larger DMs and wider frequency ranges. Frequency-dependent DMs introduce correlated noise into timing residuals whose power spectrum is `low pass' in form. The correlation time is of order the geometric mean of the refraction times for the highest and lowest radio frequencies used and thus ranges from days to years, depending on the pulsar. We discuss the implications for methodologies that use large frequency separations or wide bandwidth receivers for timing measurements. Chromatic DMs are partially mitigable by using an additional chromatic term in arrival time models. Without mitigation, our results provide an additional term in the noise model for pulsar timing; they also indicate that in combination with measurement errors from radiometer noise, an arbitrary increase in total frequency range (or bandwidth) will yield diminishing benefits and may be detrimental to overall timing precision.
[88]  oai:arXiv.org:1503.04490  [pdf] - 1232033
Two radio-emission mechanisms in PSR J0901-4624
Comments: Accepted, ApJ Letters. 3 figures, 5 pages
Submitted: 2015-03-15
We have detected sporadic, bright, short-duration radio pulses from PSR J0901$-$4624. These pulses are emitted simultaneously with persistent, periodic emission that dominates the flux density when averaging over many periods of the pulsar. The bright pulses have energies that are consistent with a power-law distribution. The integrated profile of PSR J0901$-$4624 is highly polarized and shows four distinct components. The bright pulses appear to originate near the magnetic pole of the pulsar and have polarization properties unlike that of the underlying emission at the same pulse phase. We conclude that the bright pulses represent a secondary giant-micropulse emission process, possibly from a different region in the pulsar magnetosphere.
[89]  oai:arXiv.org:1503.01841  [pdf] - 1231959
A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars
Comments: 41 pages, 28 figures, accepted for publication in MNRAS
Submitted: 2015-03-05
We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.
[90]  oai:arXiv.org:1412.1599  [pdf] - 1223477
A fast radio burst in the direction of the Carina dwarf spheroidal galaxy
Comments: 5 pages, 3 figures, 1 table, accepted to ApJ Letters. Updated to match published version
Submitted: 2014-12-04, last modified: 2015-01-15
We report the real-time discovery of a fast radio burst (FRB 131104) with the Parkes radio telescope in a targeted observation of the Carina dwarf spheroidal galaxy. The dispersion measure of the burst is 779 cm$^{-3}$ pc, exceeding predictions for the maximum line-of-sight Galactic contribution by a factor of 11. The temporal structure of the burst is characterized by an exponential scattering tail with a timescale of 2.0$^{+0.8}_{-0.5}$ ms at 1582 MHz that scales as frequency to the power $-$4.4$^{+1.6}_{-1.8}$ (all uncertainties represent 95% confidence intervals). We bound the intrinsic pulse width to be $<0.64$ ms due to dispersion smearing across a single spectrometer channel. Searches in 78 hours of follow-up observations with the Parkes telescope reveal no additional sporadic emission and no evidence for associated periodic radio emission. We hypothesize that the burst is associated with the Carina dwarf galaxy. Follow-up observations at other wavelengths are necessary to test this hypothesis.
[91]  oai:arXiv.org:1501.00126  [pdf] - 916316
Understanding pulsar magnetospheres with the SKA
Comments: 21 pages, 3 figures, to be published in: "Advancing Astrophysics with the Square Kilometre Array", Proceedings of Science, PoS(AASKA14)038
Submitted: 2014-12-31
The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.
[92]  oai:arXiv.org:1412.5853  [pdf] - 1223710
The unusual glitch recoveries of the high magnetic field pulsar J1119$-$6127
Comments: 14 pages, 4 figures; Accepted for publication in MNRAS
Submitted: 2014-12-18
Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative changes. Moreover, it was followed by an uncommon glitch relaxation that resulted in a smaller spin-down rate relative to the prediction of the pre-glitch timing model. Here, we present 4 years of new radio timing observations and analyse the total of 16 years of timing data for this source. The new data uncover an ongoing evolution of the spin-down rate, thereby allowing us to exclude permanent changes in the external or internal torque as a stand-alone cause of the peculiar features of the glitch recovery. Furthermore, no additional variations of the radio pulse profile are detected, strengthening the association of the previously observed transient emission features with the glitching activity. A self-consistent measurement of the braking index yields a value $n\simeq2.7$, indicating a trajectory in the $P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in $n$ at the epoch of the 2007 glitch.
[93]  oai:arXiv.org:1408.1694  [pdf] - 910017
A 24-Hour Global Campaign To Assess Precision Timing of the Millisecond Pulsar J1713+0747
Comments: Published in the Astrophysical Journal
Submitted: 2014-08-07, last modified: 2014-12-16
The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.
[94]  oai:arXiv.org:1406.5297  [pdf] - 1215100
Prospects for gravitational-wave detection and supermassive black hole astrophysics with pulsar timing arrays
Comments: Updated to match version accepted to MNRAS. 15 pages, 5 figures
Submitted: 2014-06-20, last modified: 2014-12-15
Large-area sky surveys show that massive galaxies undergo at least one major merger in a Hubble time. Ongoing pulsar timing array (PTA) experiments are aimed at measuring the gravitational wave (GW) emission from binary supermassive black holes (SMBHs) at the centres of galaxy merger remnants. In this paper, using the latest observational estimates for a range of galaxy properties and scaling relations, we predict the amplitude of the GW background generated by the binary SMBH population. We also predict the numbers of individual binary SMBH GW sources. We predict the characteristic strain amplitude of the GW background to lie in the range $5.1\times10^{-16}<A_{\rm yr}<2.4\times10^{-15}$ at a frequency of $(1\,{\rm yr})^{-1}$, with 95% confidence. Higher values within this range, which correspond to the more commonly preferred choice of galaxy merger timescale, will fall within the expected sensitivity ranges of existing PTA projects in the next few years. In contrast, we find that a PTA consisting of at least 100 pulsars observed with next-generation radio telescopes will be required to detect continuous-wave GWs from binary SMBHs. We further suggest that GW memory bursts from coalescing SMBH pairs are not viable sources for PTAs. Both the GW background and individual GW source counts are dominated by binaries formed in mergers between early-type galaxies of masses $\gtrsim5\times10^{10}M_{\odot}$ at redshifts $\lesssim 1.5$. Uncertainties in the galaxy merger timescale and the SMBH mass - galaxy bulge mass relation dominate the uncertainty in our predictions.
[95]  oai:arXiv.org:1410.3323  [pdf] - 1222536
Searching for gravitational wave memory bursts with the Parkes Pulsar Timing Array
Comments: Accepted by MNRAS
Submitted: 2014-10-13
Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.
[96]  oai:arXiv.org:1408.5129  [pdf] - 869631
An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set
Comments: A few typos corrected, to match the published version, 13 pages, 11 figures
Submitted: 2014-08-21, last modified: 2014-09-18
We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude $h_0$ for a range of GW frequencies. For example, at $10^{-8}$ Hz our analysis has excluded with $95\%$ confidence the presence of signals with $h_0\geqslant 1.7\times 10^{-14}$. Our new limits are about a factor of four more stringent than those of Yardley et al. (2010) based on an earlier PPTA data set and a factor of two better than those reported in the recent Arzoumanian et al. (2014) paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of $10^{9} M_{\odot}$ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of $4 \times 10^{-3} {\rm{Mpc}}^{-3} {\rm{Gyr}}^{-1}$ at $95\%$ confidence on the coalescence rate of nearby ($z \lesssim 0.1$) supermassive binary black holes in circular orbits with chirp masses of $10^{10}M_{\odot}$.
[97]  oai:arXiv.org:1408.6933  [pdf] - 1216570
The Three Discrete Nulling Timescales of PSR J1717-4054
Comments: 10 pages, 8 beautiful figures, accepted for publication in MNRAS
Submitted: 2014-08-29
PSR J1717-4054 is one of the small class of pulsars which null on intermediate (~hour) timescales. Such pulsars represent an important link between classical nullers, whose emission vanishes for a few rotations, and the intermittent pulsars which null for months and years. Using the Parkes radio telescope and the Australia Telescope Compact Array, we have studied the emission from PSR J1717-4054 over intervals from single pulses to years. We have identified and characterised nulling at three discrete timescales: the pulsar emits during "active states" separated by nulls lasting many thousands of rotations, while active states themselves are interrupted by nulls with a bimodal distribution of durations; one to two rotations, or tens of rotations. We detect no subpulse structure or flux density variations during active states, and we place stringent limits (<0.1% of the mean active flux density) on pulsed emission during inactive states. Finally, our high-quality data have also allowed us to measure for the first time many important properties of PSR J1717-4054, including its position, spindown rate, spectrum, polarization characteristics, and pulse broadening in the interstellar medium.
[98]  oai:arXiv.org:1407.6271  [pdf] - 881928
Six Faint Gamma-ray Pulsars Seen with the Fermi Large Area Telescope -- Towards a Sample Blending into the Background
Comments: 12 pages, 5 figures. Accepted to A&A
Submitted: 2014-07-23
Context: GeV gamma-ray pulsations from over 140 pulsars have been characterized using the Fermi Large Area Telescope, enabling improved understanding of the emission regions within the neutron star magnetospheres, and the contributions of pulsars to high energy electrons and diffuse gamma rays in the Milky Way. The first gamma-ray pulsars to be detected were the most intense and/or those with narrow pulses. Aims: As the Fermi mission progresses, progressively fainter objects can be studied. In addition to more distant pulsars (thus probing a larger volume of the Galaxy), or ones in high background regions (thus improving the sampling uniformity across the Galactic plane), we detect pulsars with broader pulses or lower luminosity. Adding pulsars to our catalog with inclination angles that are rare in the observed sample, and/or with lower spindown power, will reduce the bias in the currently known gamma-ray pulsar population. Methods: We use rotation ephemerides derived from radio observations to phase-fold gamma rays recorded by the Fermi Large Area Telescope, to then determine the pulse profile properties. Spectral analysis provides the luminosities and, when the signal-to-noise ratio allows, the cutoff energies. We constrain the pulsar distances by different means in order to minimize the luminosity uncertainties. Results: We present six new gamma-ray pulsars with an eclectic mix of properties. Three are young, and three are recycled. They include the farthest, the lowest power, two of the highest duty-cycle pulsars seen, and only the fourth young gamma-ray pulsar with a radio interpulse. We discuss the biases existing in the current gamma-ray pulsar catalog, and steps to be taken to mitigate the bias.
[99]  oai:arXiv.org:1407.0435  [pdf] - 846068
The Role of FAST in Pulsar Timing Arrays
Comments: To be published in a special issue on FAST in RAA
Submitted: 2014-07-01
The Five-hundred-meter Aperture Spherical Telescope (FAST) will become one of the world-leading telescopes for pulsar timing array (PTA) research. The primary goals for PTAs are to detect (and subsequently study) ultra-low-frequency gravitational waves, to develop a pulsar-based time standard and to improve solar system planetary ephemerides. FAST will have the sensitivity to observe known pulsars with significantly improved signal-to-noise ratios and will discover a large number of currently unknown pulsars. We describe how FAST will contribute to PTA research and show that jitter- and timing-noise will be the limiting noise processes for FAST data sets. Jitter noise will limit the timing precision achievable over data spans of a few years while timing noise will limit the precision achievable over many years.
[100]  oai:arXiv.org:1406.4716  [pdf] - 1215051
Limitations in timing precision due to single-pulse shape variability in millisecond pulsars
Comments: 21 pages, 23 figures, accepted by MNRAS
Submitted: 2014-06-18
High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.
[101]  oai:arXiv.org:1406.1363  [pdf] - 832489
The fast C(3P) + CH3OH reaction as an efficient loss process for gas-phase interstellar methanol
Comments: Accepted to RSC Advances
Submitted: 2014-06-05
Rate constants for the C(3P) + CH3OH reaction have been measured in a continuous supersonic flow reactor over the range 50 K to 296 K. C(3P) was created by the in-situ pulsed laser photolysis of CBr4, a multiphoton process which also produced some C(1D), allowing us to investigate simultaneously the low temperature kinetics of the C(1D) + CH3OH reaction. C(1D) atoms were followed by an indirect chemiluminescent tracer method in the presence of excess CH3OH. C(3P) atoms were detected by the same chemiluminescence technique and also by direct vacuum ultra-violet laser induced fluorescence (VUV LIF). Secondary measurements of product H(2S) atom formation have been undertaken allowing absolute H atom yields to be obtained by comparison with a suitable reference reaction. In parallel, statistical calculations have been performed based on ab-initio calculations of the complexes, adducts and transition states (TSs) relevant to the title reaction. By comparison with the experimental H atom yields, the preferred reaction pathways could be determined, placing important constraints on the statistical calculations. The experimental and theoretical work are in excellent agreement, predicting a negative temperature dependence of the rate constant increasing from 2.2 x 10-11 cm3 molecule-1 s-1 at 296 K to 20.0 x 10-11 cm3 molecule-1 s-1 at 50 K. CH3 and HCO are found to be the major products under our experimental conditions. As this reaction is not considered in current astrochemical networks, its influence on interstellar methanol abundances is tested using a model of dense interstellar clouds.
[102]  oai:arXiv.org:1404.5183  [pdf] - 1209067
Binary super-massive black hole environments diminish the gravitational-wave signal in the pulsar timing band
Comments: 15 pages, 7 figures, 1 table, accepted by MNRAS
Submitted: 2014-04-21
We assess the effects of super-massive black hole (SMBH) environments on the gravitational-wave (GW) signal from binary SMBHs. To date, searches with pulsar timing arrays for GWs from binary SMBHs, in the frequency band $\sim1-100$\,nHz, include the assumptions that all binaries are circular and evolve only through GW emission. However, dynamical studies have shown that the only way that binary SMBH orbits can decay to separations where GW emission dominates the evolution is through interactions with their environments. We augment an existing galaxy and SMBH formation and evolution model with calculations of binary SMBH evolution in stellar environments, accounting for non-zero binary eccentricities. We find that coupling between binaries and their environments causes the expected GW spectral energy distribution to be reduced with respect to the standard assumption of circular, GW-driven binaries, for frequencies up to $\sim20$\,nHz. Larger eccentricities at binary formation further reduce the signal in this regime. We also find that GW bursts from individual eccentric binary SMBHs are unlikely to be detectable with current pulsar timing arrays. The uncertainties in these predictions are large, owing to observational uncertainty in SMBH-galaxy scaling relations and the galaxy stellar mass function, uncertainty in the nature of binary-environment coupling, and uncertainty in the numbers of the most massive binary SMBHs. We conclude, however, that low-frequency GWs from binary SMBHs may be more difficult to detect with pulsar timing arrays than currently thought.
[103]  oai:arXiv.org:1401.3003  [pdf] - 800178
The High Time Resolution Universe Pulsar Survey X: Discovery of four millisecond pulsars and updated timing solutions of a further 12
Comments: 21 pages, 9 figures, 2 tables, MNRAS accepted
Submitted: 2014-01-13
We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.
[104]  oai:arXiv.org:1311.3541  [pdf] - 1180727
Evidence of an asteroid encountering a pulsar
Comments: Accepted by The Astrophysical Journal Letters, 12th November 2013
Submitted: 2013-11-14
Debris disks and asteroid belts are expected to form around young pulsars due to fallback material from their original supernova explosions. Disk material may migrate inwards and interact with a pulsar's magnetosphere, causing changes in torque and emission. Long term monitoring of PSR J0738-4042 reveals both effects. The pulse shape changes multiple times between 1988 and 2012. The torque, inferred via the derivative of the rotational period, changes abruptly from September 2005. This change is accompanied by an emergent radio component that drifts with respect to the rest of the pulse. No known intrinsic pulsar processes can explain these timing and radio emission signatures. The data lead us to postulate that we are witnessing an encounter with an asteroid or in-falling debris from a disk.
[105]  oai:arXiv.org:1311.0588  [pdf] - 1180443
The kinematics and orbital dynamics of the PSR B1259$-$63/LS 2883 system from 23 years of pulsar timing
Comments: Accepted to MNRAS. 11 pages, 5 figures
Submitted: 2013-11-04
We present an analysis of $23~$yr of pulse arrival times for PSR B1259$-$63. The pulsar is in a binary orbit about its approximately $20~M_{\odot}$ companion LS$~$2883. Our best-fitting timing solution has none of the pulse-number ambiguities that have plagued previous attempts to model the binary orbit. We measure significant first and second time derivatives of the projected semimajor axis and longitude of the periastron of the orbit. These variations are found to be consistent with the precession of the orbital plane due to classical spin-orbit coupling. The derived moment of inertia and spin of the companion are consistent with the companion rotating at near-breakup velocity. The system configuration is also consistent with the geometry derived from both the polarisation of the radio emission and the eclipse of the pulsar by the equatorial disc of the companion. We find strong evidence for orbital period decay that can be attributed to mass loss from the companion star. We also measure a significant proper motion that locates the birth of the system in the Centaurus OB1 association. By combining proper motion of the pulsar with radial velocity measurements of the companion, we measure the three-dimensional velocity of the system. This velocity is used to constrain the masses of the stars prior to the supernova explosion and the kick the pulsar received at or immediately after the explosion.
[106]  oai:arXiv.org:1310.4569  [pdf] - 1180011
Gravitational-wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution
Comments: Published in Science 18 October 2013
Submitted: 2013-10-16
The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be ${\Omega}_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.
[107]  oai:arXiv.org:1310.2120  [pdf] - 753956
TempoNest: A Bayesian approach to pulsar timing analysis
Comments: 21 pages, 12 Figures
Submitted: 2013-10-08
A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.
[108]  oai:arXiv.org:1310.1852  [pdf] - 745032
Broadband Pulsations from PSR B1821-24: Implications for Emission Models and the Pulsar Population of M28
Comments: 5 figures, 2 tables, accepted for publication in ApJ
Submitted: 2013-10-07
We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density
[109]  oai:arXiv.org:1305.3036  [pdf] - 1166611
Radio properties of the magnetar near Sagittarius A* from observations with the Australia Telescope Compact Array
Comments: 4 pages, 4 figures, published in MNRAS Letters
Submitted: 2013-05-14, last modified: 2013-08-14
We have carried out observations of the newly-discovered magnetar in the direction of Sagittarius A* using the Australia Telescope Compact Array in four frequency bands from 4.5 to 20 GHz. Radio pulsations are clearly detected at all frequencies. We measure the pulsar's dispersion measure to be 1650 +/- 50 pc/cm^3, the highest of any of the known pulsars. Once Faraday rotation has been taken into account, the pulse profile is more than 80% linearly polarized at all frequencies and has a small degree (5%) of circular polarization. The rotation measure of -67000 +/- 500$ rad/m^2 is the largest (in magnitude) ever measured for a pulsar but still a factor 8 smaller than Sgr A* itself. The combination of the dispersion and rotation measures implies an integrated magnetic field strength of -50uG along the line of sight. The flux density appears to have increased by about a factor of two between observations made 30 days apart. This object therefore joins the small class of radio emitting magnetars.
[110]  oai:arXiv.org:1307.7221  [pdf] - 1173020
Dispersion measure variations in a sample of 168 pulsars
Comments: 9 pages, 6 figures, 2 tables, Accepted for publication in MNRAS
Submitted: 2013-07-27
We analyse dispersion measure (DM) variations in six years of radio observations of more than 160 young pulsars, all gamma ray candidates for the Fermi gamma ray telescope mostly located close to the Galactic plane. DMs were fit across 256 MHz of bandwidth for observations centred at 1.4 GHz and across three frequencies -- 0.7 GHz, 1.4 GHz, and 3.1 GHz -- where multifrequency observations were available. Changes in dispersion measure, dDM/dt, were calculated using a weighted linear fit across all epochs of available DMs. DM variations were detected at a 3\sigma level in 11 pulsars, four of which were above 5\sigma: PSRs J0835-4510, J0908-4913, J1824-1945, and J1833-0827. We find that after 28 years of gradual decline, the DM of PSR J0835-4510 is now increasing. The magnitude of variations in three of the four (PSRs J0835-4510, J0908-4913, and J1833-0827) are above what models would predict for an ISM dominated by Kolmogorov turbulence. We attribute this excess as likely due to the pulsar's local environment - the supernova remnants near PSRs J0835-4510 and J1833-0827 and the pulsar wind nebula around PSR J0908-4913. Upper limits were determined for all pulsars without detectable values of dDM/dt, most limits were found to lie above the levels of variations predicted by ISM theory. We find our results to be consistent with scattering estimates from the NE2001 model along these lines of sight.
[111]  oai:arXiv.org:1307.5375  [pdf] - 1172864
Interplanetary spacecraft navigation using pulsars
Comments: accepted by Advances in Space Research
Submitted: 2013-07-20
We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than approx. 20km and velocity measurements with a precision of approx. 0.1m/s.
[112]  oai:arXiv.org:1304.4644  [pdf] - 1166021
A connection between radio state changing and glitch activity in PSR J0742-2822
Comments: Accepted for publication in MNRAS; 6 pages, 6 figures
Submitted: 2013-04-16
PSR J0742-2822 exhibits two distinct emission states that are identified by discrete changes in the observed pulse profile. These changes have previously been shown to correlate with changes in the derivative of the pulse frequency. In this paper, we use observations with the Parkes radio telescope at a centre frequency of 1369 and 3100 MHz to produce high phase resolution polarisation profiles for the two modes and perform a detailed study of the correlation between observed pulse shape and spin-down rate. We find no correlation for at least 200 days prior to a glitch in the pulsar at MJD 55022, following which the correlation becomes strong. This suggests a link between the emission state switching phenomenon and glitch events. We discuss the possibility that emission state switching is driven by the interaction between the magnetosphere and the interior of the neutron star.
[113]  oai:arXiv.org:1212.4360  [pdf] - 619981
Six millisecond pulsars detected by the Fermi Large Area Telescope and the radio/gamma-ray connection of millisecond pulsars
Comments: 19 pages, 10 figures, 4 tables. Published by MNRAS
Submitted: 2012-12-18, last modified: 2013-01-29
We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.
[114]  oai:arXiv.org:1301.6429  [pdf] - 1159306
An Asteroid Belt Interpretation for the Timing Variations of the Millisecond Pulsar B1937+21
Comments: 15 pages, 5 figures, accepted to the Astrophysical Journal
Submitted: 2013-01-27
Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.
[115]  oai:arXiv.org:1211.5887  [pdf] - 1158031
Measurement and correction of variations in interstellar dispersion in high-precision pulsar timing
Comments: 16 pages, 8 figures. Accepted for publication in MNRAS
Submitted: 2012-11-26
Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales. In this paper we demonstrate that without consideration of wavelength- independent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem. Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603-7202 of ~2x10^{-3} cm^{-3}pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.
[116]  oai:arXiv.org:1211.2761  [pdf] - 1157698
Deep X-ray Observations of the Young High-Magnetic-Field Radio Pulsar J1119-6127 and Supernova Remnant G292.2-0.5
Comments: 11 pages, 4 figures, accepted for publication in ApJ
Submitted: 2012-11-12
High-magnetic-field radio pulsars are important transition objects for understanding the connection between magnetars and conventional radio pulsars. We present a detailed study of the young radio pulsar J1119-6127, which has a characteristic age of 1900yr and a spin-down-inferred magnetic field of 4.1e13G, and its associated supernova remnant G292.2-0.5, using deep XMM-Newton and Chandra X-ray Observatory exposures of over 120ks from each telescope. The pulsar emission shows strong modulation below 2.5keV, with a single-peaked profile and a large pulsed fraction of 0.48+/-0.12. Employing a magnetic, partially ionized hydrogen atmosphere model, we find that the observed pulse profile can be produced by a single hot spot of temperature 0.13keV covering about one third of the stellar surface, and we place an upper limit of 0.08keV for an antipodal hot spot with the same area. The nonuniform surface temperature distribution could be the result of anisotropic heat conduction under a strong magnetic field, and a single-peaked profile seems common among high-B radio pulsars. For the associated remnant G292.2-0.5, its large diameter could be attributed to fast expansion in a low-density wind cavity, likely formed by a Wolf-Rayet progenitor, similar to two other high-B radio pulsars.
[117]  oai:arXiv.org:1211.2035  [pdf] - 588272
Detection of 107 glitches in 36 southern pulsars
Comments: Accepted for publication in MNRAS
Submitted: 2012-11-08
Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 1990 and 2011 have been searched for period glitches. A total of 107 glitches were identified in 36 pulsars, where 61 have previously been reported and 46 are new discoveries. Glitch parameters were measured by fitting the timing residual data. Observed relative glitch sizes \Delta\nu_g/\nu range between 10^-10 and 10^-5, where \nu = 1/P is the pulse frequency. We confirm that the distribution of \Delta\nu_g/\nu is bimodal with peaks at approximately 10^-9 and 10^-6. Glitches are mostly observed in pulsars with characteristic ages between 10^3 and 10^5 years, with large glitches mostly occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially also has a bimodal distribution. Large glitches generally have low Q, typically a few per cent, but large Q values are observed in both large and small glitches. Observed time constants for exponential recoveries ranged between 10 and 300 days with some tendency for longer timescales in older pulsars. Shorter timescale recoveries may exist but were not revealed by our data which typically have observation intervals of 2 - 4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent linear increase in \dot\nu in the inter-glitch interval. Where an exponential recovery is also observed, the effects of this are superimposed on the linear increase in \dot\nu. In some cases, the slope of the linear recovery changes at the time of a glitch. The \ddot\nu values characterising the linear changes in \dot\nu are almost always positive and, after subtracting the magnetospheric component of the braking, are approximately proportional to the ratio of |\dot\nu| and the inter-glitch interval, as predicted by vortex-creep models.
[118]  oai:arXiv.org:1210.7179  [pdf] - 1152488
XMM-Newton observation of the very old pulsar J0108-1431
Comments: 11 pages, 11 figures, accepted by ApJ; Figures 2 and 3 are in addition to the ApJ paper version
Submitted: 2012-10-26
We report on an X-ray observation of the 166 Myr old radio pulsar J0108-1431 with XMM-Newton. The X-ray spectrum can be described by a power-law model with a relatively steep photon index Gamma~3 or by a combination of thermal and non-thermal components, e.g., a power-law component with fixed photon index Gamma~2 plus a blackbody component with a temperature of kT=0.11 keV. The two-component model appears more reasonable considering different estimates for the hydrogen column density. The non-thermal X-ray efficiency in the single power-law model is eta^PL (1-10 keV) = L^PL (1-10 keV) / Edot ~ 0.003, higher than in most other X-ray detected pulsars. In the case of the combined model, the non-thermal and thermal X-ray efficiencies are even higher, eta^PL (1-10 keV) ~ eta^bb ~ 0.006. We detected X-ray pulsations at the radio period of P=0.808s with significance of 7sigma. The pulse shape in the folded X-ray lightcurve (0.15-2 keV) is asymmetric, with statistically significant contributions from up to 5 leading harmonics. Pulse profiles at two different energy ranges differ slightly: the profile is asymmetric at low energies, 0.15-1 keV, while at higher energies, 1-2 keV, it has a nearly sinusodial shape. The radio pulse peak leads the 0.15-2 keV X-ray pulse peak by delta phi = 0.06 +/- 0.03.
[119]  oai:arXiv.org:1210.7021  [pdf] - 1152471
Pulse intensity modulation and the timing stability of millisecond pulsars: A case study of PSR J1713+0747
Comments: 9 pages, 9 figures. Astrophysical Journal, accepted
Submitted: 2012-10-25
Most millisecond pulsars, like essentially all other radio pulsars, show timing errors well in excess of what is expected from additive radiometer noise alone. We show that changes in amplitude, shape and pulse phase for the millisecond pulsar J1713+0747 cause this excess error. These changes appear to be uncorrelated from one pulse period to the next. The resulting time of arrival variations are correlated across a wide frequency range and are observed with different backend processors on different days, confirming that they are intrinsic in origin and not an instrumental effect or caused by strongly frequency dependent interstellar scattering. Centroids of single pulses show an rms phase variation \approx 40 microsec, which dominates the timing error and is the same phase jitter phenomenon long known in slower spinning, canonical pulsars. We show that the amplitude modulations of single pulses are modestly correlated with their arrival time fluctuations. We also demonstrate that single-pulse variations are completely consistent with arrival time variations of pulse profiles obtained by integrating N pulses such that the arrival time error decreases proportional to 1/\sqrt{N}. We investigate methods for correcting times of arrival for these pulse shape changes, including multi-component TOA fitting and principal component analysis. These techniques are not found to improve the timing precision of the observations. We conclude that when pulse shape changes dominate timing errors, the timing precision of PSR J1713+0747 can be only improved by averaging over a larger number of pulses.
[120]  oai:arXiv.org:1210.6130  [pdf] - 1152365
The Parkes Pulsar Timing Array Project
Comments: 36 pages, 13 figures, accepted by Publications of the Astronomical Society of Australia (PASA)
Submitted: 2012-10-23
A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.
[121]  oai:arXiv.org:1210.3854  [pdf] - 1152113
Does a 'stochastic' background of gravitational waves exist in the pulsar timing band?
Comments: 14 pages, 9 figures, 2 tables, accepted by ApJ
Submitted: 2012-10-14
We investigate the effects of gravitational waves (GWs) from a simulated population of binary super-massive black holes (SMBHs) on pulsar timing array datasets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing datasets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The "Hellings & Downs" correlations between simulated datasets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a 5-year dataset. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in 5-year pulsar timing array datasets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.
[122]  oai:arXiv.org:1208.3560  [pdf] - 1150745
Development of a pulsar-based timescale
Comments: Accepted for publication in MNRAS
Submitted: 2012-08-17
Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.
[123]  oai:arXiv.org:1205.3089  [pdf] - 511898
Radio Searches of Fermi LAT Sources and Blind Search Pulsars: The Fermi Pulsar Search Consortium
Comments: 2011 Fermi Symposium proceedings - eConf C110509
Submitted: 2012-05-14
We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.
[124]  oai:arXiv.org:1204.6111  [pdf] - 1118320
Optimal Interpolation and Prediction in Pulsar Timing
Comments: 9 pages, accepted by MNRAS
Submitted: 2012-04-27
For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty.
[125]  oai:arXiv.org:1201.6641  [pdf] - 1116285
Limits on the Stochastic Gravitational Wave Background from the North American Nanohertz Observatory for Gravitational Waves
Comments: To be submitted to ApJ
Submitted: 2012-01-31
We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.
[126]  oai:arXiv.org:1110.1271  [pdf] - 1084680
Pulsed Gamma Rays from the Original Millisecond and Black Widow Pulsars: a case for Caustic Radio Emission?
Comments: Accepted for publication in the Astrophysical Journal
Submitted: 2011-10-06, last modified: 2011-11-07
We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nan\c{c}ay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4\sigma$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.
[127]  oai:arXiv.org:1108.0161  [pdf] - 1083099
Discovery of gamma and X-ray pulsations from the young and energetic PSR J1357-6429 with Fermi and XMM-Newton
Comments: Accepted for publication by A&A; 7 pages, 5 figures
Submitted: 2011-07-31
Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.
[128]  oai:arXiv.org:1106.4047  [pdf] - 1077427
Minimum Requirements for Detecting a Stochastic Gravitational Wave Background Using Pulsars
Comments: Submitted to the Astrophysical Journal
Submitted: 2011-06-20
We assess the detectability of a nanohertz gravitational wave (GW) background with respect to additive red and white noise in the timing of millisecond pulsars. We develop detection criteria based on the cross-correlation function summed over pulsar pairs in a pulsar timing array. The distribution of correlation amplitudes is found to be non-Gaussian and highly skewed, which significantly influences detection and false-alarm probabilities. When only white noise and GWs contribute, our detection results are consistent with those found by others. Red noise, however, drastically alters the results. We discuss methods to meet the challenge of GW detection ("climbing mount significance") by distinguishing between GW-dominated and red or white-noise limited regimes. We characterize detection regimes by evaluating the number of millisecond pulsars that must be monitored in a high-cadence, 5-year timing program for a GW background spectrum $h_c(f) = A f^{-2/3}$ with $A = 10^{-15}$ yr$^{-2/3}$. Unless a sample of 20 super-stable millisecond pulsars can be found --- those with timing residuals from red-noise contributions $\sigma_r \lesssim 20$ ns --- a much larger timing program on $\gtrsim 50 - 100$ MSPs will be needed. For other values of $A$, the constraint is $\sigma_r \lesssim 20 {\rm ns} (A/10^{-15} {\rm yr}^{-2/3})$. Identification of suitable MSPs itself requires an aggressive survey campaign followed by characterization of the level of spin noise in the timing residuals of each object. The search and timing programs will likely require substantial fractions of time on new array telescopes in the southern hemisphere as well as on existing ones.
[129]  oai:arXiv.org:1012.3507  [pdf] - 1042608
Pulsars with the Australian Square Kilometre Array Pathfinder
Comments: To appear in proceedings of "Radio Pulsars: An astrophysical key to unlock the secrets of the Universe"
Submitted: 2010-12-15
The Australian Square Kilometre Array Pathfinder (ASKAP) is a 36-element array with a 30-square-degree field of view being built at the proposed SKA site in Western Australia. We are conducting a Design Study for pulsar observations with ASKAP, planning both timing and search observations. We provide an overview of the ASKAP telescope and an update on pulsar-related progress.
[130]  oai:arXiv.org:1010.4794  [pdf] - 275634
Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
Comments: 16 pages and 6 figures. ApJ, accepted
Submitted: 2010-10-22
We investigate rotational spin noise (referred to as timing noise) in non-accreting pulsars: millisecond pulsars, canonical pulsars, and magnetars. Particular attention is placed on quantifying the strength and non-stationarity of timing noise in millisecond pulsars because the long-term stability of these objects is required to detect nanohertz gravitational radiation. We show that a single scaling law is sufficient to characterize timing noise in millisecond and canonical pulsars while the same scaling law underestimates the levels of timing noise in magnetars. The scaling law, along with a detailed study of the millisecond pulsar B1937+21, leads us to conclude that timing noise is latent in most millisecond pulsars and will be measurable in many objects when better arrival time estimates are obtained over long data spans. The sensitivity of a pulsar timing array to gravitational radiation is strongly affected by any timing noise. We conclude that detection of proposed gravitational wave backgrounds will require the analysis of more objects than previously suggested over data spans that depend on the spectra of both the gravitational wave background and of the timing noise. It is imperative to find additional millisecond pulsars in current and future surveys in order to reduce the effects of timing noise.
[131]  oai:arXiv.org:1010.3785  [pdf] - 243797
A Measurement Model for Precision Pulsar Timing
Comments: 30 pages, 9 figures, to be submitted to the Astrophysical Journal
Submitted: 2010-10-18
This paper describes a comprehensive measurement model for the error budget of pulse arrival times with emphasis on intrinsic pulse jitterand plasma propagation effects (particularly interstellar scattering), which are stochastic in time and have diverse dependences on radio frequency. To reduce their contribution, timing measurements can be made over a range of frequencies that depends on a variety of pulsar and instrumentation-dependent factors that we identify. A salient trend for high signal-to-noise measurements of millisecond pulsars is that time-of-arrival precision is limited either by irreducible interstellar scattering or by pulse-phase jitter caused by variable emission within pulsar magnetospheres. A cap on timing errors implies that pulsars must be confined to low dispersion measures (DMs) and observed at high frequencies. Use of wider bandwidths that increase signal-to-noise ratios will degrade timing precision if nondispersive chromatic effects are not mitigated. The allowable region in the DM-frequency plane depends on how chromatic timing perturbations are addressed. Without mitigation, observations at 1.4~GHz or 5~GHz are restricted to $\DM\lesssim 30$ and $\lesssim 100~\DMu$, respectively. With aggressive mitigation of interstellar scattering and use of large telescopes to provide adequate sensitivity at high frequencies (e.g. Arecibo, FAST, phase 1 of the SKA, and the SKA), pulsars with DMs up to 500~$\DMu$ can be used in precision timing applications. We analyze methods that fit arrival times vs. frequency at a given epoch prior to multi-epoch fitting. While the terms of greatest astrophysical interest are achromatic (e.g. orbital and gravitational wave perturbations), measurements may ultimately be limited by similarly achromatic stochasticity in a pulsar's spin rate.
[132]  oai:arXiv.org:0911.5206  [pdf] - 172800
The international pulsar timing array project: using pulsars as a gravitational wave detector
Comments: accepted by Classical and Quantum Gravity. Review talk for the Amaldi8 conference series
Submitted: 2009-11-27
The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.
[133]  oai:arXiv.org:0909.1058  [pdf] - 28042
The North American Nanohertz Observatory for Gravitational Waves
Comments: This document is the NANOGrav consortium's submission to Astro2010's Program Prioritization Panel on Particle Astrophysics and Gravitation
Submitted: 2009-09-05
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) is a consortium of astronomers whose goal is the creation of a galactic scale gravitational wave observatory sensitive to gravitational waves in the nHz-microHz band. It is just one component of an international collaboration involving similar organizations of European and Australian astronomers who share the same goal. Gravitational waves, a prediction of Einstein's general theory of relativity, are a phenomenon of dynamical space-time generated by the bulk motion of matter, and the dynamics of space-time itself. They are detectable by the small disturbance they cause in the light travel time between some light source and an observer. NANOGrav exploits radio pulsars as both the light (radio) source and the clock against which the light travel time is measured. In an array of radio pulsars gravitational waves manifest themselves as correlated disturbances in the pulse arrival times. The timing precision of today's best measured pulsars is less than 100 ns. With improved instrumentation and signal-to-noise it is widely believed that the next decade could see a pulsar timing network of 100 pulsars each with better than 100 ns timing precision. Such a pulsar timing array (PTA), observed with a regular cadence of days to weeks, would be capable of observing supermassive black hole binaries following galactic mergers, relic radiation from early universe phenomena such as cosmic strings, cosmic superstrings, or inflation, and more generally providing a vantage on the universe whose revolutionary potential has not been seen in the 400 years since Galileo first turned a telescope to the heavens.
[134]  oai:arXiv.org:astro-ph/0605145  [pdf] - 1233787
Rocking the Lighthouse: Circumpulsar Asteroids and Radio Intermittency
Comments: 20 pages, 6 figures, submitted to ApJ
Submitted: 2006-05-04
We propose that neutral, circumpulsar debris entering the light cylinder can account for many time-dependent pulsar phenomena that are otherwise difficult to explain. Neutral material avoids propeller ejection and injects sufficient charges -- after heating, evaporation, and ionization -- to alter current flows and pair-production and thus trigger, detune, or extinguish coherent emission. Relevant phenomena, with time scales from seconds to months, include nulls, rotating radio transients (RRATs), rapid changes in pulse profile (``mode changes''), variable subpulse drift rates, quasi-periodic bursts from B1931+24, and torque variations. Over the 10 Myr lifetime of a canonical pulsar with trillion-gauss surface magnetic field, less than a millionth of an Earth mass of material is needed to modulate the Goldreich-Julian current by 100%. Circumpulsar material originates from metal-rich, supernova fallback gas that aggregates into asteroids. Debris disks can inject sufficient material on time scales of interest, yet be too tenuous to form large planets detectable in pulse timing data. Asteroid migration results from collisions and the radiation-driven Yarkovsky and Poynting-Robertson effects. For B1931+24, an asteroid in a $\sim 40$~day elliptical orbit pollutes the magnetosphere stochastically through collisions with other debris. Injection is less likely for hot, young and highly magnetized pulsars or millisecond pulsars that pre-ionize any debris material well outside their small magnetospheres. Injection effects will therefore be most prominent in long-period, cooler pulsars, consistent with the distribution of relevant objects in perid and period derivative. A pulsar's spin history and its radiation-beam orientation may influence whether it displays nulling, RRATs and other effects.
[135]  oai:arXiv.org:astro-ph/0410548  [pdf] - 68397
Magnetospheric Birefringence Induces Unique Polarization Signatures in Neutron-Star Spectra
Comments: 4 pages, 5 figures, changes to reflect accepted version
Submitted: 2004-10-22, last modified: 2006-02-18
We study the propagation of polarization light through the magnetosphere of neutron stars. At intermediate frequencies (the optical through the infrared), both the birefringence induced by the plasma and by quantumelectrodynamics influence the observation polarization of radiation from the surface of the neutron star. Because these two processes compete in this regime, we find that polarization observations can constrain the properties of the neutron-star magnetosphere, specifically the total charge density. We calculate both the phase-resolved and the phase-averaged polarization signatures induced by magnetospheric birefringence.