• ### FRB microstructure revealed by the real-time detection of FRB170827(1803.05697)

May 1, 2018 astro-ph.HE
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.
• ### Extreme Scattering Events Towards Two Young Pulsars(1712.00426)

Dec. 1, 2017 astro-ph.HE
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.
• ### Follow up of GW170817 and its electromagnetic counterpart by Australian-led observing programs(1710.05846)

Oct. 16, 2017 astro-ph.SR, astro-ph.HE
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.
• We present results from the first directed search for nontensorial gravitational waves. While general relativity allows for tensorial (plus and cross) modes only, a generic metric theory may, in principle, predict waves with up to six different polarizations. This analysis is sensitive to continuous signals of scalar, vector or tensor polarizations, and does not rely on any specific theory of gravity. After searching data from the first observation run of the advanced LIGO detectors for signals at twice the rotational frequency of 200 known pulsars, we find no evidence of gravitational waves of any polarization. We report the first upper limits for scalar and vector strains, finding values comparable in magnitude to previously-published limits for tensor strain. Our results may be translated into constraints on specific alternative theories of gravity.
• ### Spectral Flattening at Low Frequencies in Crab Giant Pulses(1709.03651)

Sept. 12, 2017 astro-ph.HE
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.
• 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.
• ### Chandra Phase-Resolved Spectroscopy of the High Magnetic Field Pulsar B1509-58(1704.00460)

April 3, 2017 astro-ph.HE
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.
• ### Robust Estimation of Scattering in Pulsar Timing Analysis(1703.02108)

March 6, 2017 astro-ph.IM, astro-ph.HE
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.
• ### Wide-band Profile Domain Pulsar Timing Analysis(1612.05258)

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.
• ### The NANOGrav Nine-Year Data Set: Excess Noise in Millisecond Pulsar Arrival Times(1610.01731)

Dec. 2, 2016 astro-ph.IM, astro-ph.HE
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.
• ### Radio-interferometric monitoring of FRB 131104: A coincident AGN flare, but no evidence for a cosmic fireball(1611.05580)

Nov. 17, 2016 astro-ph.HE
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.
• ### The magnetic field and turbulence of the cosmic web measured using a brilliant fast radio burst(1611.05758)

Nov. 17, 2016 astro-ph.CO, astro-ph.HE
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.
• ### Modelling and mitigating refractive propagation effects in precision pulsar timing observations(1609.07573)

Sept. 24, 2016 astro-ph.IM, astro-ph.HE
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.
• ### Single-Source Gravitational Wave Limits from the J1713+0747 24-hr Global Campaign(1509.05446)

Aug. 20, 2016 astro-ph.IM, astro-ph.HE
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.
• ### The disturbance of a millisecond pulsar magnetosphere(1608.02163)

Aug. 6, 2016 astro-ph.SR, astro-ph.HE
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.
• ### Characterising the rotational irregularities of the Vela pulsar from 21 yr of phase-coherent timing(1604.02661)

April 10, 2016 astro-ph.SR, astro-ph.HE
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.
• ### From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release(1602.05570)

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.
• ### The International Pulsar Timing Array: First Data Release(1602.03640)

Feb. 11, 2016 astro-ph.IM
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.
• ### A pilot ASKAP survey of radio transient events in the region around the intermittent pulsar PSR J1107-5907(1512.02702)

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.
• ### Emission-rotation correlation in pulsars: new discoveries with optimal techniques(1511.05481)

Nov. 18, 2015 astro-ph.HE
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.
• ### Versatile Directional Searches for Gravitational Waves with Pulsar Timing Arrays(1510.08068)

Oct. 27, 2015 astro-ph.IM
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.
• ### Timing analysis for 20 millisecond pulsars in the Parkes Pulsar Timing Array(1510.04434)

Oct. 15, 2015 astro-ph.IM, astro-ph.HE
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.
• ### Gravitational waves from binary supermassive black holes missing in pulsar observations(1509.07320)

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.
• ### Profile Stochasticity in PSR J1909-3744(1509.07505)

Sept. 24, 2015 astro-ph.IM
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.
• ### Pulsar Observations of Extreme Scattering Events(1506.07948)

June 26, 2015 astro-ph.SR, astro-ph.HE
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.