Search

We have searched for anomalous events using 2,520 hours of archival observations from Murriyang, CSIRO's Parkes radio telescope. These observations were originally undertaken to search for pulsars. We used a machine-learning algorithm based on ResNet and Uniform Manifold Approximation and Projection (UMAP) in order to identify parts of the data stream that potentially contain anomalous signals. Many of these anomalous events are radio frequency interference, which were subsequently filtered using multibeam information. We detected 202 anomalous events and provide their positions and event times. We discuss the possibility that one of the events comes from radio emission from a white dwarf star.The other events are currently of unknown type., Comment: 9 pages, 7 figures

Fast Radio Bursts (FRBs) are short-duration radio transients that occur at random times in host galaxies distributed all over the sky. Large field of view instruments can play a critical role in the blind search for rare FRBs. We present a concept for an all-sky FRB monitor using a compact all-sky phased array (CASPA), which can efficiently achieve an extremely large field of view of $\sim10^4$ square degrees. Such a system would allow us to conduct a continuous, blind FRB search covering the entire southern sky. Using the measured FRB luminosity function, we investigate the detection rate for this all-sky phased array and compare the result to a number of other proposed large field-of-view instruments. We predict a rate of a few FRB detections per week and determine the dispersion measure and redshift distributions of these detectable FRBs. This instrument is optimal for detecting FRBs in the nearby Universe and for extending the high-end of the FRB luminosity function through finding ultraluminous events. Additionally, this instrument can be used to shadow the new gravitational-wave observing runs, detect high energy events triggered from Galactic magnetars and search for other bright, but currently unknown transient signals., Comment: 10 pages, 5 figures, accepted for publication in PASA

Radio astronomy file formats are now required to store wide frequency bandwidths and multiple simultaneous receiver beams and must be able to account for versatile observing modes and numerous calibration strategies. The need to capture and archive high-time and high frequency-resolution data, along with the comprehensive metadata that fully describe the data, implies that a new data format and new processing software are required. This requirement is suited to a well-defined, hierarchically-structured and flexible file format. In this paper we present the Spectral-Domain Hierarchical Data Format (`SDHDF') -- a new file format for radio astronomy data, in particular for single dish or beam-formed data streams. Since 2018, SDHDF has been the primary format for data products from the spectral-line and continuum observing modes at Murriyang, the CSIRO Parkes 64-m radio telescope, and we demonstrate that this data format can also be used to store observations of pulsars and fast radio bursts., Comment: Supplementary material (SDHDF definition): https://ars.els-cdn.com/content/image/1-s2.0-S2213133724000192-mmc1.pdf

The Parkes 20 cm Multibeam pulsar surveys have discovered nearly half of the known pulsars and revealed many distant pulsars with high dispersion measures. Using a sample of 1,301 pulsars from these surveys, we have explored the spatial distribution and birth rate of normal pulsars. The pulsar distances used to calculate the pulsar surface density are estimated from the YMW16 electron-density model. When estimating the impact of the Galactic background radiation on our survey, we projected pulsars in the Galaxy onto the Galactic plane, assuming that the flux density distribution of pulsars is uniform in all directions, and utilized the most up-to-date background temperature map. We also used an up-to-date version of the ATNF Pulsar Catalogue to model the distribution of pulsar flux densities at 1400 MHz. We derive an improved radial distribution for the pulsar surface density projected on to the Galactic plane, which has a maximum value at $\sim$4 kpc from the Galactic Centre. We also derive the local surface density and birthrate of pulsars, obtaining 47 $\pm$ 5 $\mathrm{kpc^{-2}}$ and $\sim$ 4.7 $\pm$ 0.5 $\mathrm{kpc^{-2}\ Myr^{-1}}$, respectively. For the total number of potentially detectable pulsars in the Galaxy, we obtain (1.1 $\pm$ 0.2) $\times$ $10^{4}$ and (1.1 $\pm$ 0.2) $\times$ $10^{5}$ before and after applying the TM98 beaming correction model. The radial distribution function is used to estimate the proportion of pulsars in each spiral arm and the Galactic centre.

Understanding the jitter noise resulting from single-pulse phase and shape variations is important for the detection of gravitational waves using pulsar timing array. We presented measurements of jitter noise and single-pulse variability of 12 millisecond pulsars that are part of the International Pulsar Timing Array sample using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We found that the levels of jitter noise can vary dramatically among pulsars. A moderate correlation with a correlation coefficient of 0.57 between jitter noise and pulse width is detected. To mitigate jitter noise, we performed matrix template matching using all four Stokes parameters. Our results revealed a reduction in jitter noise ranging from 6.7\% to 39.6\%. By performing longitude-resolved fluctuation spectrum analysis, we identified periodic intensity modulations in 10 pulsars. In PSR J0030+0451, we detected single-pulses with energies more than 10 times the average pulse energy, suggesting the presence of giant pulses. We also observed a periodic mode-changing phenomenon in PSR J0030+0451. We examined the achievable timing precision by selecting a sub-set of pulses with a specific range of peak intensity, but no significant improvement in timing precision is achievable., Comment: 14 pages, 9 figures, Accepted for publication in ApJ

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA., Comment: 21 pages, 9 figures, submitted to ApJ

Black widows are millisecond pulsars ablating their companions. The material blown from the companion blocks the radio emission, resulting in radio eclipses. The properties of the eclipse medium are poorly understood. Here, we present direct evidence of the existence of magnetic fields in the eclipse medium of the black widow PSR J2051$-$0827 using observations made with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We detect a regular decrease in rotation measure (RM) in the egress of eclipse, changing from $60\,\rm rad\,m^{-2}$ to $-28.7\,\rm rad\,m^{-2}$. The RM gradually changes back to normal when the line-of-sight moves away from the eclipse. The estimated line-of-sight magnetic field strength in the eclipse medium is $\sim 0.1$ G. The RM reversal could be caused by a change of the magnetic field strength along the line of sight due to binary orbital motion. The RM reversal phenomenon has also been observed in some repeating fast radio bursts (FRBs), and the study of spider pulsars may provide additional information about the origin of FRBs., Comment: 7 pages, 3 figures, accept for publication in ApJ

Rotating Radio Transients (RRATs) are neutron stars that emit sporadic radio bursts. We detected 1955 single pulses from RRAT J1913+1330 using the 19-beam receiver of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). These pulses were detected in 19 distinct clusters, with 49.4% of them occurring with a waiting time of one rotation period. The energy distribution of these individual pulses exhibited a wide range, spanning three orders of magnitude, reminiscent of repeating fast radio bursts (FRBs). Furthermore, we observed abrupt variations in pulse profile, width, peak flux, and fluence between adjacent sequential pulses. These findings suggest that this RRAT could be interpreted as a pulsar with extreme pulse-to-pulse modulation. The presence of sequential pulse trains during active phases, along with significant pulse variations in profile, fluence, flux, and width, should be intrinsic to a subset of RRATs. Our results indicate that J1913+1330 represents a peculiar source that shares certain properties with populations of nulling pulsars, giant pulses, and FRBs from different perspectives. The dramatic pulse-to-pulse variation observed in J1913+1330 could be attributed to unstable pair creation above the polar cap region and the variation of the site where streaming pairs emit coherently. Exploring a larger sample of RRATs exhibiting similar properties to J1913+1330 has the potential to significantly advance our understanding of pulsars, RRATs, and FRBs., Comment: 12 pages, 10 figures, 1 table, accepted for publication in ApJ

The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.

Flux densities are basic observation parameters to describe pulsars. In the most updated pulsar catalog, 24% of the listed radio pulsars have no flux density measurement at any frequency. Here, we report the first flux density measurements, spectral indices, pulse profiles, and correlations of the spectral index with pulsar parameters for 19 pulsars employing the Ultra-Wideband Low (UWL) receiver system installed on the Parkes radio telescope. The results for spectral indices of 17 pulsars are in the range between -0.6 and -3.10. The polarization profiles of thirteen pulsars are shown. There is a moderate correlation between the spectral index and spin frequency. For most pulsars detected, the S/N ratio of pulse profile is not high, so DM, Faraday rotation measure (RM), and polarization can not be determined precisely. Twenty-nine pulsars were not detected in our observations. We discuss the possible explanations for why these pulsars were not detected., Comment: 19 pages, 6 figures, 5 tables, accepted for publication in RAA

We report the phase-connected timing ephemeris, polarization pulse profiles, Faraday rotation measurements, and Rotating-Vector-Model (RVM) fitting results of twelve millisecond pulsars (MSPs) discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal radio Astronomy FAST survey (CRAFTS). The timing campaigns were carried out with FAST and Arecibo over three years. Eleven of the twelve pulsars are in neutron star - white dwarf binary systems, with orbital periods between 2.4 and 100 d. Ten of them have spin periods, companion masses, and orbital eccentricities that are consistent with the theoretical expectations for MSP - Helium white dwarf (He WD) systems. The last binary pulsar (PSR J1912$-$0952) has a significantly smaller spin frequency and a smaller companion mass, the latter could be caused by a low orbital inclination for the system. Its orbital period of 29 days is well within the range of orbital periods where some MSP - He WD systems have shown anomalous eccentricities, however, the eccentricity of PSR J1912$-$0952 is typical of what one finds for the remaining MSP - He WD systems., Comment: 11 pages, 5 figures, MNRAS accepted

We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $\tau\propto\Delta{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering due to multi-path propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of $\sim100$ m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz., Comment: 8 pages, 4 figures, MNRAS accepted

The environment of actively repeating fast radio bursts (FRBs) has been shown to be complex and varying. The recently localized FRB 20190520B is extremely active, has the largest confirmed host dispersion measure, and is only the second FRB source associated with a compact, persistent radio source (PRS). The main tracer of the magneto-ionic environments is the rotation measure (RM), a path-integral of the line-of-sight component of magnetic field strength (B) and electron density, which does not allow a direct probe of the B-field configuration. Here we report direct evidence for a B-field reversal based on the observed sign change and extreme variation of FRB 20190520B's RM, which changed from $\sim10000$ rad m$^{-2}$ to $\sim-16000$ rad m$^{-2}$ between June 2021 and January 2022. Such extreme RM reversal has never been observed before in any FRB nor in any astronomical object. The implied short-term change of the B-field configuration in or around the FRB could be due to the vicinity of massive black holes, or a magnetized companion star in binary systems, or a young supernova remnant along the line of sight., Comment: It was merged into arXiv:2202.11112. Science380,599-603(2023). DOI:10.1126/science.abo6526

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^{\alpha}$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $\alpha = -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $\alpha = -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability., Comment: 17 pages, 12 figures, accepted in MNRAS

Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating bodies and also reveal new relativistic effects that have been expected but are now observed for the first time. These include effects of light propagation in strong gravitational fields which are currently not testable by any other method. We observe retardation and aberrational light-bending that allow determination of the pulsar's spin direction. In total, we have detected seven post-Keplerian (PK) parameters, more than for any other binary pulsar. For some of these effects, the measurement precision is so high that for the first time we have to take higher-order contributions into account. These include contributions of A's effective mass loss (due to spin-down) to the observed orbital period decay, a relativistic deformation of the orbit, and effects of the equation of state of super-dense matter on the observed PK parameters via relativistic spin-orbit coupling. We discuss the implications of our findings, including those for the moment of inertia of neutron stars. We present the currently most precise test of general relativity's (GR's) quadrupolar description of gravitational waves, validating GR's prediction at a level of $1.3 \times 10^{-4}$ (95% conf.). We demonstrate the utility of the Double Pulsar for tests of alternative theories by focusing on two specific examples and discuss some implications for studies of the interstellar medium and models for the formation of the Double Pulsar. Finally, we provide context to other types of related experiments and prospects for the future., Comment: 56 pages, 23 Figures. Published by Physical Review X. Uploaded ArXiv version is authors' (pre-proof) version with abbreviated abstract. Resubmission for updated acknowledgement information. For final published version see https://link.aps.org/doi/10.1103/PhysRevX.11.041050

We report on an {unusually} bright observation of PSR J1720$-$0533 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulsar is in a black widow system that {was discovered by the Commensal Radio Astronomy FAST Survey (CRAFTS). By coincidence, a bright scintillation maximum was simultaneous with the eclipse in our observation which allowed for precise measurements of flux density variations, as well as dispersion measure (DM) and polarization.} We found that there are quasi-periodic pulse emission variations with a modulation period of $\sim$ {22\,s} during the ingress of the eclipse, which could be caused by plasma lensing. {No such periodic modulation was found during the egress of the eclipse. } {The linear polarization of the pulsar disappears before the eclipse, even before there is a visually obvious change in DM. We also found that the pulse scattering maybe play an important role in the eclipse of PSR J1720$-$0533.}, Comment: 7 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

We report the follow-up of 10 pulsars discovered by the Five-hundred-meter Aperture Spherical radio-Telescope (FAST) during its commissioning. The pulsars were discovered at a frequency of 500-MHz using the ultra-wide-band (UWB) receiver in drift-scan mode, as part of the Commensal Radio Astronomy FAST Survey (CRAFTS). We carried out the timing campaign with the 100-m Effelsberg radio-telescope at L-band around 1.36 GHz. Along with 11 FAST pulsars previously reported, FAST seems to be uncovering a population of older pulsars, bordering and/or even across the pulsar death-lines. We report here two sources with notable characteristics. PSR J1951$+$4724 is a young and energetic pulsar with nearly 100% of linearly polarized flux density and visible up to an observing frequency of 8 GHz. PSR J2338+4818, a mildly recycled pulsar in a 95.2-d orbit with a Carbon-Oxygen white dwarf (WD) companion of $\gtrsim 1\rm{M}_{\odot}$, based on estimates from the mass function. This system is the widest WD binary with the most massive companion known to-date. Conspicuous discrepancy was found between estimations based on NE2001 and YMW16 electron density models, which can be attributed to under-representation of pulsars in the sky region between Galactic longitudes $70^o

We have searched for weak fast radio burst (FRB) events using a database containing 568,736,756 transient events detected using the Parkes radio telescope between 1997 and 2001. In order to classify these pulses, and to identify likely FRB candidates, we used a machine learning algorithm based on ResNet. We identified 81 new candidate FRBs and provide details of their positions, event times, and dispersion measures. These events were detected in only one beam of the Parkes multibeam receiver. We used a relatively low S/N cutoff threshold when selecting these bursts and some have dispersion measures only slightly exceeding the expected Galactic contribution. We therefore present these candidate FRBs as a guide for follow-up observations in the search for repeating FRBs., Comment: 8 pages, 5 figures, 2 tables, accepted by MNRAS

A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background., Comment: 10 pages, 5 figures

The event rate, energy distribution, and time-domain behaviour of repeating fast radio bursts (FRBs) contains essential information regarding their physical nature and central engine, which are as yet unknown. As the first precisely-localized source, FRB 121102 has been extensively observed and shows non-Poisson clustering of bursts over time and a power-law energy distribution. However, the extent of the energy distribution towards the fainter end was not known. Here we report the detection of 1652 independent bursts with a peak burst rate of 122~hr^{-1}, in 59.5 hours spanning 47 days. A peak in the isotropic equivalent energy distribution is found to be ~4.8 x 10^{37} erg at 1.25~GHz, below which the detection of bursts is suppressed. The burst energy distribution is bimodal, and well characterized by a combination of a log-normal function and a generalized Cauchy function. The large number of bursts in hour-long spans allow sensitive periodicity searches between 1 ms and 1000 s. The non-detection of any periodicity or quasi-periodicity poses challenges for models involving a single rotating compact object. The high burst rate also implies that FRBs must be generated with a high radiative efficiency, disfavoring emission mechanisms with large energy requirements or contrived triggering conditions., Comment: 69 pages, 14 figures, 2 tables; Updates for author correction

The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 \pm 0.2 M_\odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_\odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments., Comment: 44 pages, 32 figures, 7 tables. Accepted for publication in MNRAS; v2 includes corrected typos in parameter tables and reference added

The Parkes telescope has been monitoring 286 radio pulsars approximately monthly since 2007 at an observing frequency of 1.4 GHz. The wide dispersion measure (DM) range of the pulsar sample and the uniformity of the observing procedure make the data-set extremely valuable for studies of flux density variability and the interstellar medium. Here, we present flux density measurements and modulation indices of these pulsars over this period. We derive the structure function from the light curves and discuss the contributions to it from measurement noise, intrinsic variability and interstellar scintillation. Despite a large scatter, we show that the modulation index is inversely correlated with DM, and can be generally described by a power-law with an index of $\sim-0.7$ covering DMs from $\sim10$ to 1000 cm$^{-3}$ pc. We present refractive timescales and/or lower limits for a group of 42 pulsars. These often have values significantly different from theoretical expectations, indicating the complex nature of the interstellar medium along individual lines of sight. In particular, local structures and non-Kolmogorov density fluctuations are likely playing important roles in the observed flux density variation of many of these pulsars., Comment: 27 pages, 46 figures, accepted for publication in MNRAS

Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no convincing detections we place upper limits on the GW amplitude for different event durations. From these bounds we place limits on the cosmic string tension of G mu ~ 10^{-5}, and highlight that this bound is independent from those obtained using other techniques. We discuss the physical implications of our results and the prospect of probing cosmic strings in the era of Square Kilometre Array (SKA)., Comment: Accepted for publication by MNRAS

Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterise noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify "exponential dip" events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713$+$0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods., Comment: 18 pages, 7 figures, 5 tables

We have observed PSR~B1534+12 (J1537+1155), a pulsar with a neutron star companion, using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We found that this pulsar shows two distinct emission states: a weak state with a wide pulse profile and a burst state with a narrow pulse profile. The weak state is always present. We cannot, with our current data, determine whether the pulse energy of the weak state follows a normal or a log-normal distribution. The burst state energy distribution follows a power-law. The amplitude of the single pulse emission in the burst state varies significantly; the peak flux intensity of the brightest pulse is 334 times stronger than that of the average pulse. We also examined the timing precision achievable using only bright pulses, which showed no demonstrable improvement because of pulse jitter and therefore quantified the jitter noise level for this pulsar., Comment: 6 pages, 6 figures, accepted for publication in The Astrophysical Journal Letters

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar., Comment: 16 pages, 6 figures, accepted for publication in A&A

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS)., Comment: 18 pages, 8 figures, 5 tables. Accepted for publication by MNRAS. Correcting error in the caption of Figure 8

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., Comment: 7 pages, 7 figures, accepted for publication in ApJ

We have used the Nanshan 25-m Radio Telescope at Xinjiang Astronomical Observatory to obtain timing observations of 87 pulsars from 2002 July to 2014 March. Using the "Cholesky" timing analysis method we have determined positions and proper motions for 48 pulsars, 24 of which are improved positions compared to previously published values. We also present the first published proper motions for nine pulsars and improved proper motions for 21 pulsars using pulsar timing and position comparison method. The pulsar rotation parameters are derived and are more accurate than previously published values for 36 pulsars. Glitches are detected in three pulsars: PSRs J1722$-$3632, J1852$-$0635 and J1957+2831. For the first two, the glitches are large, with $\Delta\nu_g/\nu > 10^{-6}$, and they are the first detected glitches in these pulsars. PSR J1722$-$3632 is the second oldest pulsar with large glitch. For the middle-age pulsars ($\tau_c > 10^5$~yr), the calculated braking indices, $|n|$, are strongly correlated with $\tau_c$ and the numbers of positive and negative values of $n$ are almost equal. For young pulsars ($\tau_c < 10^5$~yr), there is no correlation between $|n|$ and $\tau_c$ and most have $n>0$., Comment: 21 pages, 8 figures, accepted for publication by ApJ

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference., Comment: 15 pages, 4 Figures, 6 Tables, 2 Appendix, accepted for publication in ApJS

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., Comment: 26 pages, 6 figures. Accepted for publication in PASA. The data release can be downloaded from https://doi.org/10.25919/5db90a8bdeb59

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., Comment: submitted to PASA

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., Comment: accepted by MNRAS

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., Comment: Submitted to MNRAS and in review, 23 pages, 5 figures

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., Comment: 8 pages, accepted for publication in MNRAS

We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pulsar timing alone. These include: the orbital inclination $i$; the longitude of the ascending node $\Omega$; and the pulsar system transverse velocity. We use the annual variations to resolve the previous ambiguity in the sense of the inclination angle. Using the correct sense, and a prior probability distribution given by a constraint from pulsar timing ($i=73\pm3^\circ$), we find $\Omega=24.8\pm1.8^\circ$ and we estimate the pulsar distance to be $D=10^{+4}_{-3}$ kpc. This then gives us an estimate of this pulsar's proper motion of $\mu_{\alpha}\cos{\delta}=2.9\pm1.0$ mas yr$^{-1}$ in right ascension and $\mu_{\delta}=1.8\pm0.6$ mas yr$^{-1}$ in declination. Finally, we obtain measurements of the spatial structure of the interstellar electron density fluctuations, including: the spatial scale and anisotropy of the diffraction pattern; the distribution of scattering material along the line of sight; and spatial variation in the strength of turbulence from epoch to epoch. We find that the scattering is dominated by a thin screen at a distance of $(0.724\pm0.008)D$, with an anisotropy axial ratio $A_{\rm r} = 2.14\pm0.11$., Comment: 17 pages, 8 figures, 2 tables. Accepted for publication in MNRAS

PSR B1828-11 is a young pulsar once thought to be undergoing free precession and recently found instead to be switching magnetospheric states in tandem with spin-down changes. Here we show the two extreme states of the mode-changing found for this pulsar and comment briefly on its interpretation., Comment: 4 pages, in proceedings of "Pulsar Astrophysics -- The Next 50 Years," eds. Weltevrede et al

The young pulsar PSR B1828-11 has long been known to show correlated shape and spin-down changes with timescales of roughly 500 and 250 days, perhaps associated with large-scale magnetospheric switching. Here we present multi-hour observations with the Parkes and Green Bank Telescopes at multiple phases across the roughly 500-day cycle and show that the pulsar undergoes mode-changing between two stable, extreme profile states. The fraction of time spent in each profile state naturally accounts for the observed overall "shape parameter" (defined to be 0 for wide profiles and 1 for narrow ones); this and the variable rate of the mode transitions are directly related to the spin-down changes. We observe that the mode transition rate could plausibly function as an additional parameter governing the chaotic behaviour in this object which was proposed earlier by Seymour and Lorimer. Free precession is not needed to account for the variations., Comment: 13 pages, accepted to MNRAS

We report the discovery of a new fast radio burst (FRB), FRB~010312, in archival data from a 1.4\,GHz survey of the Magellanic Clouds using the multibeam receiver on the Parkes 64\,m-diameter radio telescope. These data sets include the Lorimer burst (FRB~010724), which it pre-dates and which we also re-detect. The new burst has a much higher dispersion measure of 1187\,cm$^{-3}$pc. The burst is one of the broadest found to date, the second earliest FRB known, and the ninth FRB discovered with a dispersion measure larger than 1000\,cm$^{-3}\,$pc. Our discovery indicates that there are likely to be more burst events still to be found in the existing Parkes data archive., Comment: 4 pages, 1 figure, 1 table. Accepted for publication in MNRAS Letters

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., Comment: Accepted for publication by the Monthly Notices of the Royal Astronomical Society

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., Comment: 9 papes, 3 figures, ApJ accepted

We have observed the remnant of supernova SN~1987A (SNR~1987A), located in the Large Magellanic Cloud (LMC), to search for periodic and/or transient radio emission with the Parkes 64\,m-diameter radio telescope. We found no evidence of a radio pulsar in our periodicity search and derived 8$\sigma$ upper bounds on the flux density of any such source of $31\,\mu$Jy at 1.4~GHz and $21\,\mu$Jy at 3~GHz. Four candidate transient events were detected with greater than $7\sigma$ significance, with dispersion measures (DMs) in the range 150 to 840\,cm$^{-3}\,$pc. For two of them, we found a second pulse at slightly lower significance. However, we cannot at present conclude that any of these are associated with a pulsar in SNR~1987A. As a check on the system, we also observed PSR~B0540$-$69, a young pulsar which also lies in the LMC. We found eight giant pulses at the DM of this pulsar. We discuss the implications of these results for models of the supernova remnant, neutron star formation and pulsar evolution., Comment: 7 pages, 3 figures, 2 tables. Accepted for publication in MNRAS

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation., Comment: Published in ApJ (2018 April 5); 13 pages, 4 figures

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., Comment: journal version

In Dai et al. (2016) we developed a formalism for computing variance images from standard, interferometric radio images containing time and frequency information. Variance imaging with future radio continuum surveys allows us to identify radio pulsars and serves as a complement to conventional pulsar searches which are most sensitive to strictly periodic signals. Here, we carry out simulations to predict the number of pulsars that we can uncover with variance imaging on future continuum surveys. We show that the Australian SKA Pathfinder (ASKAP) Evolutionary Map of the Universe (EMU) survey can find $\sim30$ normal pulsars and $\sim40$ millisecond pulsars (MSPs) over and above the number known today, and similarly an all-sky continuum survey with SKA-MID can discover $\sim140$ normal pulsars and $\sim110$ MSPs with this technique. Variance imaging with EMU and SKA-MID will detect pulsars with large duty cycles and is therefore a potential tool for finding MSPs and pulsars in relativistic binary systems. Compared with current pulsar surveys at high Galactic latitudes in the southern hemisphere, variance imaging with EMU and SKA-MID will be more sensitive, and will enable detection of pulsars with dispersion measures between $\sim10$ and 100 cm$^{-3}$ pc., Comment: 8 pages, 5 figures, accepted for publication in MNRAS

During February 2016, CSIRO Astronomy and Space Science and the Max-Planck-Institute for Radio Astronomy installed, commissioned and carried out science observations with a phased array feed (PAF) receiver system on the 64m diameter Parkes radio telescope. Here we demonstrate that the PAF can be used for pulsar observations and we highlight some unique capabilities. We demonstrate that the pulse profiles obtained using the PAF can be calibrated and that multiple pulsars can be simultaneously observed. Significantly, we find that an intrinsic polarisation leakage of -31dB can be achieved with a PAF beam offset from the centre of the field of view. We discuss the possibilities for using a PAF for future pulsar observations and for searching for fast radio bursts with the Parkes and Effelsberg telescopes., Comment: 10 pages, 8 figures, 2 tables. It has been accepted for publication in PASA

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., Comment: accepted by MNRAS,22 pages, 5 figures

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., Comment: 12 pages, 6 figures, Accepted to MNRAS

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., Comment: 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

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., Comment: 22 pages, submitted to MNRAS