Search

We present high signal-to-noise, full polarization pulse profiles for 40 bright, 'slowly'-rotating (non-recycled) pulsars using the new Ultra-Wideband Low-frequency (UWL; 704-4032 MHz) receiver on the Parkes radio telescope. We obtain updated and accurate interstellar medium parameters towards these pulsars (dispersion measures and Faraday rotation measures), and reveal Faraday dispersion towards PSR J1721-3532 caused by interstellar scattering. We find general trends in the pulse profiles including decreasing fractional linear polarization and increasing degree of circular polarization with increasing frequency, consistent with previous studies, while also revealing new features and frequency evolution. This demonstrates results that can be obtained using UWL monitoring observations of slow pulsars, which are valuable for improving our understanding of pulsar emission and the intervening interstellar medium. The calibrated data products are publicly available., 20 pages, 5 figures, 2 tables, 25 pages of supporting information with 2 figures. Accepted for publication in MNRAS. doi:10.1093/mnras/stab861

We report on a timing programme of 74 young pulsars that have been observed by the Parkes 64-m radio telescope over the past decade. Using modern Bayesian timing techniques, we have measured the properties of 124 glitches in 52 of these pulsars, of which 74 are new. We demonstrate that the glitch sample is complete to fractional increases in spin-frequency greater than $\Delta\nu^{90\%}_{g}/\nu \approx 8.1 \times 10^{-9}$. We measure values of the braking index, $n$, in 33 pulsars. In most of these pulsars, their rotational evolution is dominated by episodes of spin-down with $n > 10$, punctuated by step changes in the spin-down rate at the time of a large glitch. The step changes are such that, averaged over the glitches, the long-term $n$ is small. We find a near one-to-one relationship between the inter-glitch value of $n$ and the change in spin-down of the previous glitch divided by the inter-glitch time interval. We discuss the results in the context of a range of physical models., Comment: 24 pages, 9 figures and 9 tables. Published in MNRAS. (Fixed typo in the abstract and conclusion)

It is very computationally expensive to search for pulsars using time-domain observations, and the volume of data will be enormous with next-generation telescopes such as the Square Kilometre Array. We use artificial neural networks (ANNs), a machine learning method, for the efficient selection of pulsar candidates from radio continuum surveys; this is much cheaper than using time-domain observations. With observed quantities such as radio fluxes, sky position and compactness as inputs, our ANNs output the ‘score’ that indicates the degree of likeliness that an object is a pulsar. We demonstrate ANNs based on existing survey data by the Tata Institute for Fundamental Research (TIFR) Giant Metrewave Radio Telescope (GMRT) Sky Survey (TGSS) and the National Radio Astronomy Observatory (NRAO) Very Large Array (VLA) Sky Survey (NVSS) and we test their performance. The precision, which is the ratio of the number of pulsars classified correctly as pulsars to the number of any objects classified as pulsars, is about $96 {{\ \rm per\ cent}}$. Finally, we apply the trained ANNs to unidentified radio sources and our fiducial ANN with five inputs (the galactic longitude and latitude, the TGSS and NVSS fluxes and compactness) generates 2436 pulsar candidates from 456 866 unidentified radio sources. We need to confirm whether these candidates are truly pulsars by using time-domain observations. More information, such as polarization, will narrow the number of candidates down further.

We report the discovery of a highly-polarized, highly-variable, steep-spectrum radio source, ASKAP J173608.2-321635, located $\sim$4\,deg from the Galactic center in the Galactic plane. The source was detected six times between 2020 January and 2020 September as part of the Australian Square Kilometre Array Pathfinder Variables and Slow Transients (ASKAP VAST) survey at 888\,MHz. It exhibited a high degree ($\sim 25$\%) of circular polarization when it was visible. We monitored the source with the MeerKAT telescope from 2020 November to 2021 February on a 2--4 week cadence. The source was not detected with MeerKAT before 2021 February 07 when it appeared and reached a peak flux density of 5.6\,mJy. The source was still highly circularly polarized, but also showed up to 80\% linear polarization, and then faded rapidly with a timescale of one day. The rotation measure of the source varied significantly, from $-11.8\pm0.8$\,rad\,m$^{-2}$ to $-64.0\pm1.5$\,rad\,m$^{-2}$, over three days. No X-ray counterpart was found in follow-up \textit{Swift} or \textit{Chandra} observations about a week after the first MeerKAT detection, with upper limits of $\sim 5.0\times10^{31}$\,erg\,s$^{-1}$ (0.3--8\,keV, assuming a distance $\sim10$ kpc). No counterpart is seen in new or archival near-infrared observations down to $J=20.8$\,mag. We discuss possible identifications for ASKAP J173608.2-321635 including a low-mass star/substellar object with extremely low infrared luminosity, a pulsar with scatter-broadened pulses, a transient magnetar, or a Galactic Center Radio Transient: none of these fully explains the observations, which suggests that ASKAP J173608.2-321635 may represent part of a new class of objects being discovered through radio imaging surveys., Comment: 22 pages, 15 figures, Accepted by ApJ

The Variables and Slow Transients Survey (VAST) on the Australian Square Kilometre Array Pathfinder (ASKAP) is designed to detect highly variable and transient radio sources on timescales from 5 s to$\sim\!5$yr. In this paper, we present the survey description, observation strategy and initial results from the VAST Phase I Pilot Survey. This pilot survey consists of$\sim\!162$h of observations conducted at a central frequency of 888 MHz between 2019 August and 2020 August, with a typical rms sensitivity of$0.24\ \mathrm{mJy\ beam}^{-1}$and angular resolution of$12-20$arcseconds. There are 113 fields, each of which was observed for 12 min integration time, with between 5 and 13 repeats, with cadences between 1 day and 8 months. The total area of the pilot survey footprint is 5 131 square degrees, covering six distinct regions of the sky. An initial search of two of these regions, totalling 1 646 square degrees, revealed 28 highly variable and/or transient sources. Seven of these are known pulsars, including the millisecond pulsar J2039–5617. Another seven are stars, four of which have no previously reported radio detection (SCR J0533–4257, LEHPM 2-783, UCAC3 89–412162 and 2MASS J22414436–6119311). Of the remaining 14 sources, two are active galactic nuclei, six are associated with galaxies and the other six have no multi-wavelength counterparts and are yet to be identified.

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., 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., 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., 69 pages, 14 figures, 2 tables; Updates for author correction

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

A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature ($T \sim 1 - 100$ MeV) phase transition with a duration $(\beta/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions., Comment: 9 pages, 4 figures, 2 tables

We report the discovery of the first new pulsar with the Murchison Widefield Array (MWA), PSR J0036$-$1033, a long-period (0.9 s) nonrecycled pulsar with a dispersion measure (DM) of 23.1 ${\rm pc\,cm^{-3}}$. It was found after processing only a small fraction ($\sim$1%) of data from an ongoing all-sky pulsar survey. Follow-up observations have been made with the MWA, the upgraded Giant Metrewave Radio Telescope (uGMRT), and the Parkes 64 m telescopes, spanning a frequency range from $\sim$150 MHz to 4 GHz. The pulsar is faint, with an estimated flux density ($S$) of $\sim$1 mJy at 400 MHz and a spectrum $S(\nu)\,\propto\,\nu^{-2.0 \pm 0.2}$, where $\nu$ is frequency. The DM-derived distance implies that it is also a low-luminosity source ($\sim$ 0.1 ${\rm mJy\,kpc^2}$ at 1400 MHz). The analysis of archival MWA observations reveals that the pulsar's mean flux density varies by up to a factor of $\sim$5-6 on timescales of several weeks to months. By combining MWA and uGMRT data, the pulsar position was determined to arcsecond precision. We also report on polarization properties detected in the MWA and Parkes bands. The pulsar's nondetection in previous pulsar and continuum imaging surveys, the observed high variability, and its detection in a small fraction of the survey data searched to date, all hint at a larger population of pulsars that await discovery in the southern hemisphere, with the MWA and the future low-frequency Square Kilometre Array., Comment: 12 pages, 6 figures, 2 tables, Accepted for publication in The ApJ Letters

We report three new FRBs discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST), namely FRB 181017.J0036+11, FRB 181118 and FRB 181130, through the Commensal Radio Astronomy FAST Survey (CRAFTS). Together with FRB 181123 that was reported earlier, all four FAST-discovered FRBs share the same characteristics of low fluence ($\leq$0.2 Jy ms) and high dispersion measure (DM, $>1000$ \dmu), consistent with the anti-correlation between DM and fluence of the entire FRB population. FRB 181118 and FRB 181130 exhibit band-limited features. FRB 181130 is prominently scattered ($��_s\simeq8$ ms) at 1.25 GHz. FRB 181017.J0036+11 has full-bandwidth emission with a fluence of 0.042 Jy ms, which is one of the faintest FRB sources detected so far. CRAFTS starts to built a new sample of FRBs that fills the region for more distant and fainter FRBs in the fluence-$\rm DM_E$ diagram, previously out of reach of other surveys. The implied all sky event rate of FRBs is $1.24^{+1.94}_{-0.90} \times 10^5$ sky$^{-1}$ day$^{-1}$ at the $95\%$ confidence interval above 0.0146 Jy ms. We also demonstrate here that the probability density function of CRAFTS FRB detections is sensitive to the assumed intrinsic FRB luminosity function and cosmological evolution, which may be further constrained with more discoveries., 9 Pages, 4 Plots and 1 Table. The Astrophysical Journal Letter Accepted

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

We report the first wideband monitoring observations of PSR J1803−3002A, a relatively bright millisecond pulsar (MSP) in the globular cluster NGC 6522 with a spin period of 7.1 ms and no known binary companion. These observations were performed using the Parkes 64 m radio telescope with the Ultra-Wideband Low receiver system, which covers 704–4032 MHz. We confirm that PSR J1803−3002A is an isolated MSP located near the cluster center and probe the emission properties of the pulsar over the wide observed band. The mean pulse profile consists of three components, with the outer components becoming more prominent at higher frequencies, and a mean spectral index for the pulsed emission of −1.66 ± 0.07 over the observed band. The fractional linear and circular polarization increase with increasing frequency, which is unusual for pulsars. We determine a Faraday rotation measure of −107 ± 6 rad m−2 for the pulsar. PSR J1803−3002A is a distant pulsar in the Galactic plane, but our observations show no evidence of pulse broadening due to interstellar scattering. These results demonstrate the power of ultra-wideband receivers and signal processing systems.

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., 18 pages, 7 figures, 5 tables

PSR~J1119$-$6127, a high-magnetic field pulsar detected from radio to high-energy wavelengths, underwent a magnetar-like outburst beginning on July 27, 2016. In this paper, we study the post-outburst multi-wavelength properties of this pulsar from the radio to GeV bands and discuss its similarity with the outburst of the magnetar XTE~J1810$-$197. In phase-resolved spectral analysis of 0.5--10\,keV X-ray data collected in August 2016, the on-pulse and off-pulse spectra are both characterized by two blackbody components and also require a power-law component similar to the hard X-ray spectra of magnetars. This power-law component is no longer distinguishable in data from December, 2016. We likewise find that there was no substantial shift between the radio and X-ray pulse peaks after the 2016 X-ray outburst. The gamma-ray pulsation after the X-ray outburst is confirmed with data taken after 2016 December and the pulse structure and phase difference between the gamma-ray and radio peaks ($\sim$0.4 cycle) are also consistent with those before the X-ray outburst. These multi-wavelength observations suggest that the re-configuration of the global magnetosphere after 2016 magnetar-like outburst at most continued for about six months. We discuss the evolution of the X-ray emission after the 2016 outburst with the untwisting magnetosphere model., 12 pages,11 figures

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)., 18 pages, 8 figures, 5 tables. Accepted for publication by MNRAS. Correcting error in the caption of Figure 8

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

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

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

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

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 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

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

We report detection of the binary companion to the millisecond pulsar J2302+4442 based on the deep observations performed with the Gran Telescopio Canarias. The observations revealed an optical source with $r'$ = 23.33$\pm$0.02 and $i'$=23.08$\pm$0.02, whose position coincides with the pulsar radio position. By comparing the source colour and magnitudes with the white dwarf cooling predictions, we found that it likely represents a He or CO-core white dwarf and estimated its mass of 0.52$^{+0.25}_{-0.19}$ M$_{\odot}$ and effective temperature of 6300$^{+1000}_{-800}$ K. Combining our results with the radio timing measurements, we set constraints on the binary system inclination angle $i=73^{+6}_{-5}$ degrees and the pulsar mass $m_p=3.1^{+2.7}_{-2.0}$ M$_{\odot}$., Accepted for publication in MNRAS

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

With the largest dish Five-hundred-meter Aperture Spherical radio Telescope (FAST), both the mean and single pulses of PSR B2016$+$28, especially including the single-pulse structure, are investigated in detail in this study. The mean pulse profiles at different frequencies can be well fitted in a conal model, and the peak separation of intensity-dependent pulse profiles increases with intensity. The integrated pulses are obviously frequency dependent (pulse width decreases by $\sim\,20\%$ as frequency increases from 300 MHz to 750 MHz), but the structure of single pulses changes slightly (the corresponding correlation scale decreases by only $\sim\,1\%$). This disparity between mean and single pulses provides independent evidence for the existence of the RS-type vacuum inner gap, indicating a strong bond between particles on the pulsar surface. Diffused drifting sub-pulses are analyzed. The results show that the modulation period along pulse series ($P_3$) is positively correlated to the separation between two adjacent sub-pulses ($P_2$). This correlation may hint a rough surface on the pulsar, eventually resulting in the irregular drift of sparks. All the observational results may have significant implications in the dynamics of pulsar magnetosphere and are discussed extensively in this paper., Sci. China-Phys. Mech. Astron. 62, 959505 (2019)

We report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age between 28-10 kyr. The radio spectral index for this SNR PWN pulsar system is flat (-0.29 $\pm$ 0.01) consistent with other similar objects. We infer that the putative pulsar has a radio spectral index of -1.8, which is typical for Galactic pulsars. We searched for dispersion measures (DMs) up to 1000 cm/pc^3 but found no convincing candidates with a S/N greater than 8. We produce a polarisation map for this PWN at 5500 MHz and find a mean fractional polarisation of P $\sim 23$ percent. The X-ray power-law spectrum (Gamma $\sim 2$) is indicative of non-thermal synchrotron emission as is expected from PWN-pulsar system. Finally, we detect DEM S5 in Infrared (IR) bands. Our IR photometric measurements strongly indicate the presence of shocked gas which is expected for SNRs. However, it is unusual to detect such IR emission in a SNR with a supersonic bow-shock PWN. We also find a low-velocity HI cloud of $\sim 107$ km/s which is possibly interacting with DEM S5. SNR DEM S5 is the first confirmed detection of a pulsar-powered bow shock nebula found outside the Galaxy., 20 pages, 15 figures, accepted for publication in MNRAS

We present the earliest X-ray observations of the 2018 outburst of XTE J1810-197, the first outburst since its 2003 discovery as the prototypical transient and radio-emitting anomalous X-ray pulsar (AXP). The Monitor of All-sky X-ray Image (MAXI) detected XTE J1810-197 immediately after a November 20-26 visibility gap, contemporaneous with its reactivation as a radio pulsar, first observed on December 8. On December 13 the Nuclear Spectroscopic Telescope Array (NUSTAR) detected X-ray emission up to at least 30 keV, with a spectrum well-characterized by a blackbody plus power-law model with temperature kT = 0.74+/-0.02 keV and photon index Gamma = 4.4+/-0.2 or by a two-blackbody model with kT = 0.59+/-0.04 keV and kT = 1.0+/-0.1 keV, both including an additional power-law component to account for emission above 10 keV, with Gamma_h = -0.2+/-1.5 and Gamma_h = 1.5+/-0.5, respectively. The latter index is consistent with hard X-ray flux reported for the non-transient magnetars. In the 2-10 keV bandpass, the absorbed flux is 2E-10 erg/s/cm^2, a factor of 2 greater than the maximum flux extrapolated for the 2003 outburst. The peak of the sinusoidal X-ray pulse lags the radio pulse by approx. 0.13 cycles, consistent with their phase relationship during the 2003 outburst. This suggests a stable geometry in which radio emission originates on magnetic field lines containing currents that heat a spot on the neutron star surface. However, a measured energy-dependent phase shift of the pulsed X-rays suggests that all X-ray emitting regions are not precisely co-aligned., 7 pages, 4 Figues, 1 Table, Latex, emulateapj style; Added Table/Figure for the 3-30 keV spectral results. Astrophysical Journal Letters, in Press

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

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models., Comment: 13pages with 12 figures

The globular cluster Omega Centauri is the most massive and luminous cluster in the Galaxy. The $\gamma$-ray source FL8Y J1326.7$-$4729 is coincident with the core of the cluster, leading to speculation that hitherto unknown radio pulsars or annihilating dark matter may be present in the cluster core. Here we report on the discovery of five millisecond pulsars in Omega Centauri following observations with the Parkes radio telescope. Four of these pulsars are isolated with spin periods of 4.1, 4.2, 4.6 and 6.8 ms. The fifth has a spin period of 4.8 ms and is in an eclipsing binary system with an orbital period of 2.1 hours. Deep radio continuum images of the cluster centre with the Australian Telescope Compact Array reveal a small population of compact radio sources making it likely that other pulsars await discovery. We consider it highly likely that the millisecond pulsars are the source of the $\gamma$-ray emission. The long-term timing of these pulsars opens up opportunities to explore the dynamics and interstellar medium of the cluster., Comment: 8 pages, 2 figures, accepted for publication in the Astrophysical Journal Letters

The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state., Comment: 7 pages, 6 figures, accepted by ApJ

This work is devoted to investigate the flow behavior and microstructure evolution of AlSi10Mg by additive manufacturing in hot deformation. Tensile test is carried out under different temperatures, i.e. 200 ∘ C, 250 ∘ C, 300 ∘ C, 350 ∘ C and 400 ∘ C, with various strain rates of 0.004/s, 0.002/s and 0.0004/s. Theoretically, the modified Arrhenius-type model of additive manufacturing materials at high temperature is established. Experimentally, scanning electron microscope and optical microscope are used to analyze the mechanism of hot deformation. It is found that at 200 ∘ C both the dendritic eutectic Si and acicular eutectic Si precipitate, forming cellular structure and strengthening phases. And then, the acicular eutectic Si is partially dissolved into Al matrix, but there is no obvious growth of dendritic Si at 300 ∘ C. With the temperature increasing, the cellular structure and melt pool boundary gradually disappear. For 400 ∘ C temperature, like homogenization process, the precipitation of saturated solid solution forms the dispersed phase, and the stress-strain curve shows a weak hardening. Meanwhile, the statistical result shows that that the constitutive model agrees well with the experimental results at high temperature. This study may provide guidance for the improvement of additive manufacturing material properties by post-treatment.

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., Comment: 29 pages. 16 figures. Accepted for publication in MNRAS

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., 25 pages, 6 tables, 5 figures. Accepted for publication in MNRAS

Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line-of-sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 years for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of the Earth, pulsar, and scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our $>$5000 observations and we show that they originate from scattering by thin screens located at distances $D_1 = 89.8 \pm 0.4$ pc and $D_2 = 124 \pm 3$ pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle $i = 137.1 \pm 0.3^\circ$, and longitude of ascending node, $\Omega=206.3\pm0.4^\circ$. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modelling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars, and provides parameters essential for testing theories of gravity and constraining neutron star masses., Comment: 19 pages, 10 figures. Accepted for publication in ApJ

This work is devoted to study the mechanical behaviors of polylactic acid (PLA) materials from additive manufacturing, and an orthotropic model is established to predict the mechanical properties under arbitrary printing orientation. Firstly, the morphology of PLA material is analyzed by using scanning electron microscopy, from which the orthotropic behavior of PLA material is obtained. Three printing planes are adopted, and on each printing plane different printing angles may be selected. The mechanical parameters, including Young's modulus, yielding stress, and Poisson's ratio, for material under different printing directions are determined via quasi-static experiments. Secondly, the orthotropic constitutive model of PLA materials under different printing angles is thus obtained, and the prediction method of orthotropic mechanical properties is built based on the coordinate transformation matrix, where the orthotropic coordinate transformation matrix is acquired by attitude angles (i.e., Euler angle, the rotation angle of the local coordinate system relative to the global coordinate system). Finally, the theoretical prediction method was verified, and high-quality printing methods were recommended. In addition, the obtained results of the model show that: for PLA material, the orthotropic hypothesis model is superior to the transverse isotropic hypothesis one. This present method is not only suitable for predicting the constitutive model of printed specimens in any direction but also for other materials of fused deposition modeling.

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

We present the spin frequency and profile evolution of the radio pulsar J1119$-$6127 following magnetar-like X-ray bursts from the system in 2016 July. Using data from the Parkes radio telescope, we observe a smooth and fast spin-down process subsequent to the X-ray bursts resulting in a net change in the pulsar rotational frequency of $\Delta\nu\approx-4\times10^{-4}$\,Hz. During the transition, a net spin-down rate increase of $\Delta\dot\nu\approx-1\times10^{-10}$\,Hz\,s$^{-1}$ is observed, followed by a return of $\dot{\nu}$ to its original value. In addition, the radio pulsations disappeared after the X-ray bursts and reappeared about two weeks later with the flux density at 1.4\,GHz increased by a factor of five. The flux density then decreased and undershot the normal flux density followed by a slow recovery back to normal. The pulsar's integrated profile underwent dramatic and short-term changes in total intensity, polarization and position angle. Despite the complex evolution, we observe correlations between the spin-down rate, pulse profile shape and radio flux density. Strong single pulses have been detected after the X-ray bursts with their energy distributions evolving with time. The peculiar but smooth spin frequency evolution of PSR~J1119$-$6127 accompanied by systematic pulse profile and flux density changes are most likely to be a result of either reconfiguration of the surface magnetic fields or particle winds triggered by the X-ray bursts. The recovery of spin-down rate and pulse profile to normal provides us the best case to study the connection between high magnetic-field pulsars and magnetars., Comment: Accepted for publication in MNRAS on 2018 July 24

A new benzocoumarin derivative, dendrocoumarin (1), along with itolide A were isolated from the stems of Dendrobium nobile. The structure of the new compound 1 elucidated on the basis of NMR and mass spectroscopic data. Compounds 1 and 2 showed broad spectrum antibacterial activity against five terrestrial pathogenic bacteria.

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities., Comment: Fixed a typo in Fig. 5; Published in Phys. Rev. D as Editor's Suggestion article, 17 pages, 7 figures, 4 tables

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

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

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

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

A pulsar timing array (PTA) refers to a program of regular, high-precision timing observations of a widely distributed array of millisecond pulsars. Here we review the status of the three primary PTA projects and the joint International Pulsar Timing Array project. We discuss current results related to ultra-low-frequency gravitational wave searches and highlight opportunities for the near future., Comment: 14 pages, 3 figures; National Science Review

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., 10 pages, 8 figures, 2 tables. It has been accepted for publication in PASA

We identified a highly-polarized, steep-spectrum radio source in a deep image with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope at 888 MHz. After considering and rejecting a stellar origin for this source, we discovered a new millisecond pulsar (MSP) using observations from the Parkes radio telescope. This pulsar has period 2.77 ms and dispersion measure 228.27 pc/cm**3. Although this pulsar does not yet appear to be particularly remarkable, the short spin period, wide profile and high dispersion measure do make it relatively hard to discover through traditional blind periodicity searches. Over the course of several weeks we see changes in the barycentric period of this pulsar that are consistent with orbital motion in a binary system, but the properties of any binary need to be confirmed by further observations. While even a deep ASKAP survey may not identify large numbers of new MSPs compared to the existing population, it would be competitive with existing all-sky surveys and could discover interesting new MSPs at high Galactic latitude without the need for computationally-expensive all-sky periodicity searches., ApJ, in press