BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

R. N. Manchester

R. N. Manchester appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.IMastro-ph.SRastro-ph.GAastro-ph.COastro-phgr-qcastro-ph.EPhep-ph

Catalog footprint

What is connected

101works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

101 published item(s)

preprint2025arXiv

Frequency evolution of pulsar emission: Further evidence for fan beam model

Aims. We explore frequency-dependent changes in pulsar radio emission by analyzing their profile widths and emission heights, assessing whether the simple radius-to-frequency mapping (RFM) or the fan beam model can describe the data. Methods. Using wideband (704-4032 MHz) Murriyang (Parkes) observations of over 100 pulsars, we measured profile widths at multiple intensity levels, fit Gaussian components, and used aberration-retardation effects to estimate emission altitudes. We compared trends in width evolution and emission height with a fan beam model. Results. Similar to other recent studies, we find that while many pulsars show profiles narrowing with increasing frequency, a substantial fraction show the reverse. The Gaussian decomposition of the profiles reveals that the peak locations of the components vary little with frequency. However, the component widths do, in general, narrow with increasing frequency. This argues that propagation effects are responsible for the width evolution of the profiles rather than emission height. Overall, the evolution of the emission height with frequency is unclear, and clouded by the assumptions in the model. Spin-down luminosity correlates weakly with profile narrowing but not with emission height. Conclusions. The classic picture where pulsars emit at a single emission height which decreases with increasing observing frequency cannot explain the diversity in behavior observed here. Instead, pulsar beams likely originate from extended regions at multiple altitudes, with fan-beam or patchy structures dominating their frequency evolution. Future models must incorporate realistic plasma physics and multi-altitude emission to capture the range of pulsar behaviors.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Pulsar polarization: a broad-band population view with the Parkes Ultra-Wideband receiver

The radio polarization properties of the pulsar population are only superficially captured by the conventional picture of pulsar radio emission. We study the broadband polarization of 271 young radio pulsars, focusing particularly on circular polarization, using high quality observations made with the Ultra-Wideband Low receiver on Murriyang, the Parkes radio telescope. We seek to encapsulate polarization behaviour on a population scale by defining broad categories for frequency- and phase-dependent polarization evolution, studying the co-occurrences of these categorizations and comparing them with average polarization measurements and spin-down energy ($\dot{E}$). This work shows that deviations of the linear polarization position angle (PA) from the rotating vector model (RVM) are linked to the presence of circular polarization features and to frequency evolution of the polarization. Polarization fraction, circular polarization contribution and profile complexity all evolve with $\dot{E}$ across the population, with the profiles of high-$\dot{E}$ pulsars being simple and highly linearly polarized. The relationship between polarization fraction and circular contribution is also seen to evolve such that highly polarized profiles show less variation in circular contribution with frequency than less strongly polarized profiles. This evolution is seen both across the population and across frequency for individual sources. Understanding pulsar radio polarization requires detailed study of individual sources and collective understanding of population-level trends. For the former, we provide visualizations of their phase- and frequency-resolved polarization parameters. For the latter, we have highlighted the importance of including the impact of circular polarization and of $\dot{E}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evaluating the prevalence of spurious correlations in pulsar timing array datasets

Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

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})^α$, 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 $α= -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 $α= -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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two years of pulsar observations with the Ultra-Wideband Receiver on the Parkes radio telescope

The major programme for observing young, non-recycled pulsars with the Parkes telescope has transitioned from a narrow-band system to an ultra-wideband system capable of observing between 704 and 4032 MHz. We report here on the initial two years of observations with this receiver. Results include dispersion measure (DM) and Faraday rotation measure (RM) variability with time, determined with higher precision than hitherto, flux density measurements and the discovery of several nulling and mode changing pulsars. PSR J1703-4851 is shown to be one of a small subclass of pulsars that has a weak and a strong mode which alternate rapidly in time. PSR J1114-6100 has the fourth highest |RM| of any known pulsar despite its location far from the Galactic Centre. PSR J1825-1446 shows variations in both DM and RM likely due to its motion behind a foreground supernova remnant.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Flux density variability of 286 radio pulsars from a decade of monitoring

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Parkes transient events: I. Database of single pulses, initial results and missing FRBs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Results of 12 Years of Pulsar Timing at Nanshan -- I

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 $Δν_g/ν> 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 ($τ_c > 10^5$~yr), the calculated braking indices, $|n|$, are strongly correlated with $τ_c$ and the numbers of positive and negative values of $n$ are almost equal. For young pulsars ($τ_c < 10^5$~yr), there is no correlation between $|n|$ and $τ_c$ and most have $n>0$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Parkes Pulsar Timing Array Project: Second data release

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Timing of young radio pulsars II. Braking indices and their interpretation

In Paper I of this series, we detected a significant value of the braking index ($n$) for 19 young, high-$\dot{E}$ radio pulsars using $\sim$ 10 years of timing observations from the 64-m Parkes radio telescope. Here we investigate this result in more detail using a Bayesian pulsar timing framework to model timing noise and to perform selection to distinguish between models containing exponential glitch recovery and braking index signatures. We show that consistent values of $n$ are maintained with the addition of substantial archival data, even in the presence of glitches. We provide strong arguments that our measurements are unlikely due to exponential recovery signals from unseen glitches even though glitches play a key role in the evolution of a pulsar's spin frequency. We conclude that, at least over decadal time scales, the value of $n$ can be significantly larger than the canonical 3 and discuss the implications for the evolution of pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A pulsar-based timescale from the International Pulsar Timing Array

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An ultra-wide bandwidth (704 to 4032 MHz) receiver for the Parkes radio telescope

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Periodic mode changing in PSR J1048-5832

By analysing the data acquired from the Parkes 64-m radio telescope at 1369 MHz, we report on the phase-stationary non-drift amplitude modulation observed in PSR J1048-5832. The high-sensitivity observations revealed that the central and trailing components of the pulse profile of this pulsar switch between a strong mode and a weak mode periodically. However, the leading component remains unchanged. Polarization properties of the strong and weak modes are investigated. Considering the similarity to mode changing, we argue that the periodic amplitude modulation in PSR J1048$-$5832 is periodic mode changing. The fluctuation spectral analysis showed that the modulation period is very short (~2.1 s or 17 P1), where P1 is the rotation period of the pulsar. We find that this periodic amplitude modulation is hard to explain by existing models that account for the periodic phenomena in pulsars like subpulse drifting.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Characterising the rotational irregularities of the Vela pulsar from 21 yr of phase-coherent timing

Pulsars show two classes of rotational irregularities that can be used to understand neutron-star interiors and magnetospheres: glitches and timing noise. Here we present an analysis of the Vela pulsar spanning nearly 21 yr of observation and including 8 glitches. We identify the relative pulse number of all of the observations between glitches, with the only pulse-number ambiguities existing over glitch events. We use the phase coherence of the timing solution to simultaneously model the timing noise and glitches in a Bayesian framework, allowing us to select preferred models for both. We find the glitches can be described using only permanent and transient changes in spin frequency, i.e., no step changes in frequency derivative. For all of the glitches, we only need two exponentially decaying changes in spin frequency to model the transient components. In contrast to previous studies, we find that the dominant transient components decay on a common $\approx$ 1300 d time scale, and that a larger fraction ( $\gtrsim 25\%$) of glitch amplitudes are associated with these transient components. We also detect shorter-duration transient components of $\approx$ 25 d, as previously observed, but are limited in sensitivity to events with shorter durations by the cadence of our observations. The timing noise is well described by a steep power-law process that is independent of the glitches and subdominant to the glitch recovery. The braking index is constrained to be $<$ 8 with 95% confidence. This methodology can be used to robustly measure the properties of glitches and timing noise in other pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Investigation of nulling and subpulse drifting properties of PSR J1727$-$2739

Aims. To make a detailed study of the nulling and subpulse drifting in PSR J1727$-$2739 for investigation of its radiation properties. Methods. The observations were carried out on 20 March, 2004 using the Parkes 64-m radio telescope, with a central frequency of 1518 MHz. A total of 5568 single pulses were analysed. Results. This pulsar shows well defined nulls with lengths lasting from 6 to 281 pulses and separated by burst phases ranging from 2 to 133 pulses. We estimate a nulling fraction of around 68\%. No emission in the average pulse profile integrated over all null pulses is detected with significance above $3σ$. Most transitions from nulls to bursts are within a few pulses, whereas the transitions from bursts to nulls exhibit two patterns of decay: decrease gradually or rapidly. In the burst phase, we find that there are two distinct subpulse drift modes with vertical spacing between the drift bands of $9.7 \pm 1.6$ and $5.2 \pm 0.9$ pulse periods, while sometimes there is a third mode with no subpulse drifting. Some mode transitions occur within a single burst, while others are separated by nulls. Different modes have different average pulse profiles. Possible physical mechanisms are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Long-term observations of the pulsars in 47 Tucanae. I. A study of four elusive binary systems

For the past couple of decades, the Parkes radio telescope has been regularly observing the millisecond pulsars in 47 Tucanae (47 Tuc). This long-term timing program was designed to address a wide range of scientific issues related to these pulsars and the globular cluster where they are located. In this paper, the first of a series, we address one of these objectives: the characterization of four previously known binary pulsars for which no precise orbital parameters were known, namely 47 Tuc P, V, W and X (pulsars 47 Tuc R and Y are discussed elsewhere). We determined the previously unknown orbital parameters of 47 Tuc V and X and greatly improved those of 47 Tuc P and W. For pulsars W and X we obtained, for the first time, full coherent timing solutions across the whole data span, which allowed a much more detailed characterization of these systems. 47 Tuc W, a well-known tight eclipsing binary pulsar, exhibits a large orbital period variability, as expected for a system of its class. 47 Tuc X turns out to be in a wide, extremely circular, 10.9-day long binary orbit and its position is ~3.8 arcmin away from the cluster center, more than three times the distance of any other pulsar in 47 Tuc. These characteristics make 47 Tuc X a very different object with respect to the other pulsars of the cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Proper motions of 15 pulsars: a comparison between Bayesian and frequentist algorithms

We present proper motions for 15 pulsars which are observed regularly by the Nanshan 25-m radio telescope. Two methods, the frequentist method (Coles et al.2011) and the Bayesian (Lentati et al. 2014) method, are used and the results are compared. We demonstrate that the two methods can be applied to young pulsar data sets that exhibit large amounts of timing noise with steep spectral exponents and give consistent results. The measured positions also agree with very-long-baseline interferometric positions. Proper motions for four pulsars are obtained for the first time, and improved values are obtained for five pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The disturbance of a millisecond pulsar magnetosphere

Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The International Pulsar Timing Array: First Data Release

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $λ^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gravitational waves from binary supermassive black holes missing in pulsar observations

Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pulsar Observations of Extreme Scattering Events

Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the electron density of order 10s of cm^{-3}. Similar structures have been observed in measurements of pulsar intensity scintillation and group delay. Here we report observations of two ESEs showing increases in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). These allow us to make more complete models of the ESE, including an estimate of the "outer-scale" of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully turbulent on an AU scale. They provide some support for the idea that the structures are extended along the line of sight, such as would be the case for a scattering shell. The dispersion measurements also show a variety of AU scale structures which would not be called ESEs, yet involve electron density variations typical of ESEs and likely have the same origin.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pulsars and Gravity

Pulsars are wonderful gravitational probes. Their tiny size and stellar mass give their rotation periods a stablility comparable to that of atomic frequency standards. This is especially true of the rapidly rotating "millisecond pulsars" (MSPs). Many of these rapidly rotating pulsars are in orbit with another star, allowing pulsar timing to probe relativistic perturbations to the orbital motion. Pulsars have provided the most stringent tests of theories of relativistic gravitation, especially in the strong-field regime, and have shown that Einstein's general theory of relativity is an accurate description of the observed motions. Many other gravitational theories are effectively ruled out or at least severely constrained by these results. MSPs can also be used to form a "Pulsar Timing Array" (PTA). PTAs are Galactic-scale interferometers that have the potential to directly detect nanohertz gravitational waves from astrophysical sources. Orbiting super-massive black holes in the cores of distant galaxies are the sources most likely to be detectable. Although no evidence for gravitational waves has yet been found in PTA data sets, the latest limits are seriously constraining current ideas on galaxy and black-hole evolution in the early Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Single-pulse radio observations of the Galactic Center magnetar PSR J1745-2900

In this paper, we report radio observations of the Galactic Center magnetar PSR J1745-2900 at six epochs between June and October, 2014. These observations were carried out using the new Shanghai Tian Ma Radio Telescope at a frequency of 8.6 GHz. Both the flux density and integrated profile of PSR J1745-2900 show dramatic changes from epoch to epoch showing that the pulsar was in its "erratic" phase. On MJD 56836, the flux density of this magnetar was about 8.7 mJy, which was ten times large than that reported at the time of discovery, enabling a single-pulse analysis. The emission is dominated by narrow "spiky" pulses which follow a log-normal distribution in peak flux density. From 1913 pulses, we detected 53 pulses whose peak flux density is ten times greater than that of the integrated profile. They are concentrated in pulse phase at the peaks of the integrated profile. The pulse widths at the 50% level of these bright pulses was between 0.2 to 0.9 deg, much narrower than that of integrated profile (~12 deg). The observed pulse widths may be limited by interstellar scattering. No clear correlation was found between the widths and peak flux density of these pulses and no evidence was found for subpulse drifting. Relatively strong spiky pulses are also detected in the other five epochs of observation, showing the same properties as that detected in MJD 56836. These strong spiky pulses cannot be classified as "giant" pulses but are more closely related to normal pulse emission.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Parkes multibeam pulsar survey: VII. Timing of four millisecond pulsars and the underlying spin period distribution of the Galactic millisecond pulsar population

We present timing observations of four millisecond pulsars discovered in the Parkes 20-cm multibeam pulsar survey of the Galactic plane. PSRs J1552-4937 and J1843-1448 are isolated objects with spin periods of 6.28 and 5.47 ms respectively. PSR J1727-2946 is in a 40-day binary orbit and has a spin period of 27 ms. The 4.43-ms pulsar J1813-2621 is in a circular 8.16-day binary orbit around a low-mass companion star with a minimum companion mass of 0.2 solar masses. Combining these results with detections from five other Parkes multibeam surveys, gives a well-defined sample of 56 pulsars with spin periods below 20 ms. We develop a likelihood analysis to constrain the functional form which best describes the underlying distribution of spin periods for millisecond pulsars. The best results were obtained with a log-normal distribution. A gamma distribution is less favoured, but still compatible with the observations. Uniform, power-law and Gaussian distributions are found to be inconsistent with the data. Galactic millisecond pulsars being found by current surveys appear to be in agreement with a log-normal distribution which allows for the existence of pulsars with periods below 1.5 ms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Timing analysis for 20 millisecond pulsars in the Parkes Pulsar Timing Array

We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with $M_p=1.44\pm0.07$ $M_\odot$ and $M_p=1.47\pm0.03$ $M_\odot$ respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16% level of precision, $D=156.79\pm0.25$ pc, one of the most fractionally precise distance measurements of any star to date.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Versatile Directional Searches for Gravitational Waves with Pulsar Timing Arrays

By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "phased-up" to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set

We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude $h_0$ for a range of GW frequencies. For example, at $10^{-8}$ Hz our analysis has excluded with $95\%$ confidence the presence of signals with $h_0\geqslant 1.7\times 10^{-14}$. Our new limits are about a factor of four more stringent than those of Yardley et al. (2010) based on an earlier PPTA data set and a factor of two better than those reported in the recent Arzoumanian et al. (2014) paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of $10^{9} M_{\odot}$ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of $4 \times 10^{-3} {\rm{Mpc}}^{-3} {\rm{Gyr}}^{-1}$ at $95\%$ confidence on the coalescence rate of nearby ($z \lesssim 0.1$) supermassive binary black holes in circular orbits with chirp masses of $10^{10}M_{\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Binary super-massive black hole environments diminish the gravitational-wave signal in the pulsar timing band

We assess the effects of super-massive black hole (SMBH) environments on the gravitational-wave (GW) signal from binary SMBHs. To date, searches with pulsar timing arrays for GWs from binary SMBHs, in the frequency band $\sim1-100$\,nHz, include the assumptions that all binaries are circular and evolve only through GW emission. However, dynamical studies have shown that the only way that binary SMBH orbits can decay to separations where GW emission dominates the evolution is through interactions with their environments. We augment an existing galaxy and SMBH formation and evolution model with calculations of binary SMBH evolution in stellar environments, accounting for non-zero binary eccentricities. We find that coupling between binaries and their environments causes the expected GW spectral energy distribution to be reduced with respect to the standard assumption of circular, GW-driven binaries, for frequencies up to $\sim20$\,nHz. Larger eccentricities at binary formation further reduce the signal in this regime. We also find that GW bursts from individual eccentric binary SMBHs are unlikely to be detectable with current pulsar timing arrays. The uncertainties in these predictions are large, owing to observational uncertainty in SMBH-galaxy scaling relations and the galaxy stellar mass function, uncertainty in the nature of binary-environment coupling, and uncertainty in the numbers of the most massive binary SMBHs. We conclude, however, that low-frequency GWs from binary SMBHs may be more difficult to detect with pulsar timing arrays than currently thought.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Limitations in timing precision due to single-pulse shape variability in millisecond pulsars

High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Searching for gravitational wave memory bursts with the Parkes Pulsar Timing Array

Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Role of FAST in Pulsar Timing Arrays

The Five-hundred-meter Aperture Spherical Telescope (FAST) will become one of the world-leading telescopes for pulsar timing array (PTA) research. The primary goals for PTAs are to detect (and subsequently study) ultra-low-frequency gravitational waves, to develop a pulsar-based time standard and to improve solar system planetary ephemerides. FAST will have the sensitivity to observe known pulsars with significantly improved signal-to-noise ratios and will discover a large number of currently unknown pulsars. We describe how FAST will contribute to PTA research and show that jitter- and timing-noise will be the limiting noise processes for FAST data sets. Jitter noise will limit the timing precision achievable over data spans of a few years while timing noise will limit the precision achievable over many years.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of the Radio Remnant of Supernova 1987A: Morphological Changes from Day 7000

We present radio imaging observations of supernova remnant 1987A at 9 GHz, taken with the Australia Telescope Compact Array over 21 years from 1992 to 2013. By employing a Fourier modeling technique to fit the visibility data, we show that the remnant structure has evolved significantly since day 7000 (mid-2006): the emission latitude has gradually decreased, such that the overall geometry has become more similar to a ring structure. Around the same time, we find a decreasing trend in the east-west asymmetry of the surface emissivity. These results could reflect the increasing interaction of the forward shock with material around the circumstellar ring, and the relative weakening of the interaction with the lower-density material at higher latitudes. The morphological evolution caused an apparent break in the remnant expansion measured with a torus model, from a velocity of 4600+150-200 km/s between day 4000 and 7000 to 2400+100-200 km/s after day 7000. However, we emphasize that there is no conclusive evidence for a physical slowing of the shock at any given latitude in the expanding remnant, and that a change of radio morphology alone appears to dominate the evolution. This is supported by our ring-only fits which show a constant expansion of 3890+/-50 km/s without deceleration between days 4000 and 9000. We suggest that once the emission latitude no longer decreases, the expansion velocity obtained from the torus model should return to the same value as that measured with the ring model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational-wave Limits from Pulsar Timing Constrain Supermassive Black Hole Evolution

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be $Ω_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-resolution radio observations of SNR 1987A at high frequencies

We present new imaging observations of the remnant of Supernova (SN) 1987A at 44 GHz, performed in 2011 with the Australia Telescope Compact Array (ATCA). The $0\farcs35\times0\farcs23$ resolution of the diffraction-limited image is the highest achieved to date in high-dynamic range. We also present a new ATCA image at 18 GHz derived from 2011 observations, which is super-resolved to $0\farcs25$. The flux density is 40$\pm$2 mJy at 44 GHz and 81$\pm$6 mJy at 18 GHz. At both frequencies, the remnant exhibits a ring-like emission with two prominent lobes, and an east-west brightness asymmetry that peaks on the eastern lobe. A central feature of fainter emission appears at 44 GHz. A comparison with previous ATCA observations at 18 and 36 GHz highlights higher expansion velocities of the remnant eastern side. The 18-44 GHz spectral index is $α=-0.80$ ($S_ν\proptoν^α$). The spectral index map suggests slightly steeper values at the brightest sites on the eastern lobe, whereas flatter values are associated with the inner regions. The remnant morphology at 44 GHz generally matches the structure seen with contemporaneous X-ray and H$α$ observations. Unlike the H$α$ emission, both the radio and X-ray emission peaks on the eastern lobe. The regions of flatter spectral index align and partially overlap with the optically-visible ejecta. Simple free-free absorption models suggest that emission from a pulsar wind nebula or a compact source inside the remnant may now be detectable at high frequencies, or at low frequencies if there are holes in the ionised component of the ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interplanetary spacecraft navigation using pulsars

We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than approx. 20km and velocity measurements with a precision of approx. 0.1m/s.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

PSR J1723-2837: An Eclipsing Binary Radio Millisecond Pulsar

We present a study of PSR J1723-2837, an eclipsing, 1.86 ms millisecond binary radio pulsar discovered in the Parkes Multibeam survey. Radio timing indicates that the pulsar has a circular orbit with a 15 hr orbital period, a low-mass companion, and a measurable orbital period derivative. The eclipse fraction of ~15% during the pulsar's orbit is twice the Roche lobe size inferred for the companion. The timing behavior is significantly affected by unmodeled systematics of astrophysical origin, and higher-order orbital period derivatives are needed in the timing solution to account for these variations. We have identified the pulsar's (non-degenerate) companion using archival ultraviolet, optical, and infrared survey data and new optical photometry. Doppler shifts from optical spectroscopy confirm the star's association with the pulsar and indicate a pulsar-to-companion mass ratio of 3.3 +/- 0.5, corresponding to a companion mass range of 0.4 to 0.7 Msun and an orbital inclination angle range of between 30 and 41 degrees, assuming a pulsar mass range of 1.4-2.0 Msun. Spectroscopy indicates a spectral type of G for the companion and an inferred Roche-lobe-filling distance that is consistent with the distance estimated from radio dispersion. The features of PSR J1723-2837 indicate that it is likely a "redback" system. Unlike the five other Galactic redbacks discovered to date, PSR J1723-2837 has not been detected as a gamma-ray source with Fermi. This may be due to an intrinsic spin-down luminosity that is much smaller than the measured value if the unmeasured contribution from proper motion is large.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Six millisecond pulsars detected by the Fermi Large Area Telescope and the radio/gamma-ray connection of millisecond pulsars

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The double pulsar: evidence for neutron star formation without an iron core-collapse supernova

The double pulsar system PSR J0737-3039A/B is a double neutron star binary, with a 2.4-hour orbital period, which has allowed measurement of relativistic orbital perturbations to high precision. The low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, point to a different evolutionary scenario compared to most other known double neutron star systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 90 +/- 11 +/- 5 deg. Furthermore, we find a mean 95% upper limit on the misalignment between its spin and orbital angular momentum axes of 3.2 deg, assuming that the observed emission comes from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The International Pulsar Timing Array

The International Pulsar Timing Array (IPTA) is an organisation whose raison d'etre is to facilitate collaboration between the three main existing PTAs (the EPTA in Europe, NANOGrav in North America and the PPTA in Australia) in order to realise the benefits of combined PTA data sets in reaching the goals of PTA projects. Currently, shared data sets for 39 pulsars are available for IPTA-based projects. Operation of the IPTA is administered by a Steering Committee consisting of six members, two from each PTA, plus the immediate past Chair in a non-voting capacity. A Constitution and several Agreements define the framework for the collaboration. Web pages provide information both to members of participating PTAs and to the general public. With support from an NSF PIRE grant, the IPTA facilitates the organisation of annual Student Workshops and Science Meetings. These are very valuable both in training new students and in communicating current results from IPTA-based research.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The kinematics and orbital dynamics of the PSR B1259$-$63/LS 2883 system from 23 years of pulsar timing

We present an analysis of $23~$yr of pulse arrival times for PSR B1259$-$63. The pulsar is in a binary orbit about its approximately $20~M_{\odot}$ companion LS$~$2883. Our best-fitting timing solution has none of the pulse-number ambiguities that have plagued previous attempts to model the binary orbit. We measure significant first and second time derivatives of the projected semimajor axis and longitude of the periastron of the orbit. These variations are found to be consistent with the precession of the orbital plane due to classical spin-orbit coupling. The derived moment of inertia and spin of the companion are consistent with the companion rotating at near-breakup velocity. The system configuration is also consistent with the geometry derived from both the polarisation of the radio emission and the eclipse of the pulsar by the equatorial disc of the companion. We find strong evidence for orbital period decay that can be attributed to mass loss from the companion star. We also measure a significant proper motion that locates the birth of the system in the Centaurus OB1 association. By combining proper motion of the pulsar with radial velocity measurements of the companion, we measure the three-dimensional velocity of the system. This velocity is used to constrain the masses of the stars prior to the supernova explosion and the kick the pulsar received at or immediately after the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The radio remnant of Supernova 1987A at high frequencies and high resolution

As the remnant of Supernova (SN) 1987A has been getting brighter over time, new observations at high frequencies have allowed imaging of the radio emission at unprecedented detail. We present a new radio image at 44 GHz of the supernova remnant (SNR), derived from observations performed with the Australia Telescope Compact Array (ATCA) in 2011. The diffraction-limited image has a resolution of $349\times225$ mas, which is the highest achieved to date in high-dynamic range images of the SNR. We also present a new image at 18 GHz, also derived from ATCA observations performed in 2011, which is super-resolved to $0./!/!^{\prime\prime}25$. The new 44 and 18 GHz images yield the first high-resolution spectral index map of the remnant. The comparison of the 44 GHz image with contemporaneous X-ray and H$α$ observations allows further investigations of the nature of the remnant asymmetry and sheds more light into the progenitor hypotheses and SN explosion. In light of simple free-free absorption models, we discuss the likelihood of detecting at 44 GHz the possible emission originating from a pulsar wind nebula (PWN) or a compact source in the centre of the remnant.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Shapiro delay detection in the binary system hosting the millisecond pulsar PSR J1910-5959A

PSR J1910-5959A is a binary pulsar with a helium white dwarf companion located about 6 arcmin from the center of the globular cluster NGC6752. Based on 12 years of observations at the Parkes radio telescope, the relativistic Shapiro delay has been detected in this system. We obtain a companion mass Mc = 0.180+/-0.018Msun (1sigma) implying that the pulsar mass lies in the range 1.1Msun <= Mp <= 1.5Msun. We compare our results with previous optical determinations of the companion mass, and examine prospects for using this new measurement for calibrating the mass-radius relation for helium white dwarfs and for investigating their evolution in a pulsar binary system. Finally we examine the set of binary systems hosting a millisecond pulsar and a low mass helium white dwarf for which the mass of both stars has been measured. We confirm that the correlation between the companion mass and the orbital period predicted by Tauris & Savonije reproduces the observed values but find that the predicted Mp - Pb correlation over-estimates the neutron star mass by about 0.5Msun in the orbital period range covered by the observations. Moreover, a few systems do not obey the observed Mp - Pb correlation. We discuss these results in the framework of the mechanisms that inhibit the accretion of matter by a neutron star during its evolution in a low-mass X-ray binary.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chandra observations of the old pulsar PSR B1451-68

We present 35 ks Chandra ACIS observations of the 42 Myr old radio pulsar PSR B1451-68. A point source is detected 0.32" +/- 0.73" from the expected radio pulsar position. It has ~200 counts in the 0.3-8 keV energy range. We identify this point source as the X-ray counterpart of the radio pulsar. PSR B1451-68 is located close to a 2MASS point source, for which we derive 7% as the upper limit on the flux contribution to the measured pulsar X-ray flux. The pulsar spectrum can be described by either a power-law model with photon index Gamma=2.4 (+0.4/-0.3) and a unrealistically high absorbing column density N(H)= (2.5 (+1.2/-1.3)) * 10^(21) cm^-2, or by a combination of a kT=0.35 (+0.12/-0.07) keV blackbody and a Gamma = 1.4 +/- 0.5 power-law component for N(H)[DM]= 2.6 * 10^(20) cm^-2, estimated from the pulsar dispersion measure. At the parallactic, Lutz-Kelker bias corrected distance of 480 pc, the non-thermal X-ray luminosities in the 0.3-8 keV energy band are either Lx(nonth)= (11.3 +/- 1.7) * 10^(29) erg/s or Lx(nonth)= (5.9 (+4.9/-5.0)) * 10^(29) erg/s, respectively. This corresponds to non-thermal X-ray efficiencies of either eta(nonth)= Lx(nonth) / (dE/dt) ~ 0.005 or 0.003, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Changes in Polarization Position Angle across the Eclipse in the Double Pulsar System

We investigate the changes in polarization position angle in radiation from pulsar A around the eclipse in the Double Pulsar system PSR J0737-3039A/B at the 20 cm and 50 cm wavelengths using the Parkes 64-m telescope. The changes are ~2σ during and shortly after the eclipse at 20 cm but less significant at 50 cm. We show that the changes in position angle during the eclipse can be modelled by differential synchrotron absorption in the eclipse regions. Position angle changes after the eclipse are interpreted as Faraday rotation in the magnetotail of pulsar B. Implied charge densities are consistent with the Goldreich-Julian density, suggesting that the particle energies in the magnetotail are mildly relativistic.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining the optical emission from the double pulsar system J0737-3039

We present the first optical observations of the unique system J0737-3039 (composed of two pulsars, hereafter PSR-A and PSR-B). Ultra-deep optical observations, performed with the High Resolution Camera of the Advanced Camera for Surveys on board the Hubble Space Telescope could not detect any optical emission from the system down to m_F435W=27.0 and m_F606W=28.3. The estimated optical flux limits are used to constrain the three-component (two thermal and one non-thermal) model recently proposed to reproduce the XMM-Newton X-ray spectrum. They suggest the presence of a break at low energies in the non-thermal power law component of PSR-A and are compatible with the expected black-body emission from the PSR-B surface. The corresponding efficiency of the optical emission from PSR-A's magnetosphere would be comparable to that of other Myr-old pulsars, thus suggesting that this parameter may not dramatically evolve over a time-scale of a few Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Detection of 107 glitches in 36 southern pulsars

Timing observations from the Parkes 64-m radio telescope for 165 pulsars between 1990 and 2011 have been searched for period glitches. A total of 107 glitches were identified in 36 pulsars, where 61 have previously been reported and 46 are new discoveries. Glitch parameters were measured by fitting the timing residual data. Observed relative glitch sizes Δν_g/νrange between 10^-10 and 10^-5, where ν= 1/P is the pulse frequency. We confirm that the distribution of Δν_g/νis bimodal with peaks at approximately 10^-9 and 10^-6. Glitches are mostly observed in pulsars with characteristic ages between 10^3 and 10^5 years, with large glitches mostly occurring in the younger pulsars. Exponential post-glitch recoveries were observed for 27 large glitches in 18 pulsars. The fraction Q of the glitch that recovers exponentially also has a bimodal distribution. Large glitches generally have low Q, typically a few per cent, but large Q values are observed in both large and small glitches. Observed time constants for exponential recoveries ranged between 10 and 300 days with some tendency for longer timescales in older pulsars. Shorter timescale recoveries may exist but were not revealed by our data which typically have observation intervals of 2 - 4 weeks. For most of the 36 pulsars with observed glitches, there is a persistent linear increase in \dotνin the inter-glitch interval. Where an exponential recovery is also observed, the effects of this are superimposed on the linear increase in \dotν. In some cases, the slope of the linear recovery changes at the time of a glitch. The \ddotνvalues characterising the linear changes in \dotνare almost always positive and, after subtracting the magnetospheric component of the braking, are approximately proportional to the ratio of |\dotν| and the inter-glitch interval, as predicted by vortex-creep models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Development of a pulsar-based timescale

Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Does a 'stochastic' background of gravitational waves exist in the pulsar timing band?

We investigate the effects of gravitational waves (GWs) from a simulated population of binary super-massive black holes (SMBHs) on pulsar timing array datasets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing datasets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The "Hellings & Downs" correlations between simulated datasets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a 5-year dataset. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in 5-year pulsar timing array datasets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measurement and correction of variations in interstellar dispersion in high-precision pulsar timing

Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales. In this paper we demonstrate that without consideration of wavelength- independent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem. Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603-7202 of ~2x10^{-3} cm^{-3}pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measurement of the electron density and magnetic field of the solar wind using millisecond pulsars

The magnetic field of the solar wind near the Sun is very difficult to measure directly. Measurements of Faraday rotation of linearly polarized radio sources occulted by the solar wind provide a unique opportunity to estimate this magnetic field, and the technique has been widely used in the past. However Faraday rotation is a path integral of the product of electron density and the projection of the magnetic field on the path. The electron density near the Sun can be measured by several methods, but it is quite variable. Here we show that it is possible to measure the path integrated electron density and the Faraday rotation simultaneously at 6-10 $R_\odot$ using millisecond pulsars as the linearly polarized radio source. By analyzing the Faraday rotation measurements with and without the simultaneous electron density observations we show that these observations significantly improve the accuracy of the magnetic field estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Interpolation and Prediction in Pulsar Timing

For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

PSR B0826-34: Sometimes a rotating radio transient

We report on the detection of sporadic, strong single pulses co-existing with a periodic weak emission in the duration of weak mode of PSR B0826-34. The intensities and durations of these pulses are comparable with that of the sub-pulses in the strong mode, and these pulses are distributed within the phase ranges of the main-pulse and interpulse of the strong-mode average profile. These results suggest that there are most possibly sporadic, very short timescale turn-on of strong-mode emission during the weak-mode state of the pulsar. The emission features of the bursts of strong pulses of PSR B0826-34 during its weak-mode state are similar to those of the rotating radio transients (RRATs). PSR B0826-34 is the second pulsar known which oscillates between pulsar-like and RRAT-like modes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pulsars at Parkes

The first pulsar observations were made at Parkes on March 8, 1968, just 13 days after the publication of the discovery paper by Hewish and Bell. Since then, Parkes has become the world's most successful pulsar search machine, discovering nearly two thirds of the known pulsars, among them many highly significant objects. It has also led the world in pulsar polarisation and timing studies. In this talk I will review the highlights of pulsar work at Parkes from those 1968 observations to about 2006 when the Parkes Multibeam Pulsar Survey was essentially completed and the Parkes Pulsar Timing Array project was established.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiation properties of extreme nulling pulsar J1502-5653

We report on radiation properties of extreme nulling pulsar J1502-5653, by analyzing the data acquired from the Parkes 64-m telescope at 1374 MHz. The radio emission from this pulsar exhibits sequences of several tens to several hundreds consecutive burst pulses, separated by null pulses, and the appearance of the emission seems quasi-periodic. The null fraction from the data is estimated to be 93.6%. No emission is detected in the integrated profile of all null pulses. Systematic modulations of pulse intensity and phase are found at the beginning of burst-pulse sequences just after null. The intensity usually rises to a maximum for the first few pulses, then declines exponentially afterwards, and becomes stable after few tens of pulse periods. The peak phase appears at later longitudes for the first pulse, then drifts to earlier longitudes rapidly, and then systematic drifting gradually vanishes while the intensity becomes stable. In this pulsar, the intensity variation and phase modulation of pulses are correlated in a short duration after the emission starts following a null. Observed properties of the pulsar are compared with other nulling pulsars published previously, and the possible explanation for phase modulation is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Parkes Pulsar Timing Array Project

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Timing the main-sequence-star binary pulsar J1740-3052

PSR J1740-3052 is a young pulsar in orbit around a companion that is most likely a B-type main-sequence star. Since its discovery more than a decade ago, data have been taken at several frequencies with instruments at the Green Bank, Parkes, Lovell, and Westerbork telescopes. We measure scattering timescales in the pulse profiles and dispersion measure changes as a function of binary orbital phase and present evidence that both of these vary as would be expected due to a wind from the companion star. Using pulse arrival times that have been corrected for the observed periodic dispersion measure changes, we find a timing solution spanning 1997 November to 2011 March. This includes measurements of the advance of periastron and the change in the projected semimajor axis of the orbit and sets constraints on the orbital geometry. From these constraints, we estimate that the pulsar received a kick of at least ~50 km/s at birth. A quasi-periodic signal is present in the timing residuals with a period of 2.2 times the binary orbital period. The origin of this signal is unclear.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

X-ray Observations of High-B Radio Pulsars

The study of high-magnetic-field pulsars is important for examining the relationships between radio pulsars, magnetars, and X-ray-isolated neutron stars (XINSs). Here we report on X-ray observations of three such high-magnetic-field radio pulsars. We first present the results of a deep XMM-Newton observation of PSR J1734-3333, taken to follow up on its initial detection in 2009. The pulsar's spectrum is well fit by a blackbody with a temperature of 300 +/- 60 eV, with bolometric luminosity L_bb = 2.0(+2.2 -0.7)e+32 erg/s = 0.0036E_dot for a distance of 6.1 kpc. We detect no X-ray pulsations from the source, setting a 1 sigma upper limit on the pulsed fraction of 60% in the 0.5-3 keV band. We compare PSR J1734-3333 to other rotation-powered pulsars of similar age and find that it is significantly hotter, supporting the hypothesis that the magnetic field affects the observed thermal properties of pulsars. We also report on XMM-Newton and Chandra observations of PSRs B1845-19 and J1001-5939. We do not detect either pulsar, setting 3 sigma upper limits on their blackbody temperatures of 48 and 56 eV, respectively. Despite the similarities in rotational properties, these sources are significantly cooler than all but one of the XINSs, which we attribute to the two groups having been born with different magnetic fields and hence evolving differently.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

XMM-Newton observation of the very old pulsar J0108-1431

We report on an X-ray observation of the 166 Myr old radio pulsar J0108-1431 with XMM-Newton. The X-ray spectrum can be described by a power-law model with a relatively steep photon index Gamma~3 or by a combination of thermal and non-thermal components, e.g., a power-law component with fixed photon index Gamma~2 plus a blackbody component with a temperature of kT=0.11 keV. The two-component model appears more reasonable considering different estimates for the hydrogen column density. The non-thermal X-ray efficiency in the single power-law model is eta^PL (1-10 keV) = L^PL (1-10 keV) / Edot ~ 0.003, higher than in most other X-ray detected pulsars. In the case of the combined model, the non-thermal and thermal X-ray efficiencies are even higher, eta^PL (1-10 keV) ~ eta^bb ~ 0.006. We detected X-ray pulsations at the radio period of P=0.808s with significance of 7sigma. The pulse shape in the folded X-ray lightcurve (0.15-2 keV) is asymmetric, with statistically significant contributions from up to 5 leading harmonics. Pulse profiles at two different energy ranges differ slightly: the profile is asymmetric at low energies, 0.15-1 keV, while at higher energies, 1-2 keV, it has a nearly sinusodial shape. The radio pulse peak leads the 0.15-2 keV X-ray pulse peak by delta phi = 0.06 +/- 0.03.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chandra Observations of the High-Magnetic-Field Radio Pulsar J1718-3718

High-magnetic-field pulsars represent an important class of objects for studying the relationship between magnetars and radio pulsars. Here we report on four Chandra observations of the high-magnetic-field pulsar J1718-3718 ($B=7.4\times10^{13}$\,G) taken in 2009 as well as on a re-analysis of 2002 Chandra observations of the region. We also report an improved radio position for this pulsar based on ATCA observations. We detect X-ray pulsations at the pulsar's period in the 2009 data, with pulsed fraction 52%$\pm$13% in the 0.8--2.0 keV band. We find that the X-ray pulse is aligned with the radio pulse. The data from 2002 and 2009 show consistent spectra and fluxes; a merged overall spectrum is well fit by a blackbody of temperature $186^{+19}_{-18}$\,eV, slightly higher than predicted by standard cooling models, however, the best-fit neutron star atmosphere model is consistent with standard cooling. We find the bolometric luminosity $L^{\infty}_{\rm bb}=4^{+5}_{-2} \times10^{32}$\,erg\,s$^{-1}\sim0.3\dot{E}$, for a distance of 4.5 kpc. We compile measurements of the temperatures of all X-ray detected high-$B$ pulsars as well as those of low-$B$ radio pulsars and find evidence for the former being on average hotter than the latter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraining the coalescence rate of supermassive black-hole binaries using pulsar timing

Pulsar timing observations are used to place constraints on the rate of coalescence of supermassive black-hole (SMBH) binaries as a function of mass and redshift. In contrast to the indirect constraints obtained from other techniques, pulsar timing observations provide a direct constraint on the black-hole merger rate. This is possible since pulsar timing is sensitive to the gravitational waves (GWs) emitted by these sources in the final stages of their evolution. We find that upper bounds calculated from the recently published Parkes Pulsar Timing Array data are just above theoretical predictions for redshifts below 10. In the future, with improved timing precision and longer data spans, we show that a non-detection of GWs will rule out some of the available parameter space in a particular class of SMBH binary merger models. We also show that if we can time a set of pulsars to 10ns timing accuracy, for example, using the proposed Square Kilometre Array, it should be possible to detect one or more individual SMBH binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Discovery of gamma and X-ray pulsations from the young and energetic PSR J1357-6429 with Fermi and XMM-Newton

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-Precision Timing of 5 Millisecond Pulsars: Space Velocities, Binary Evolution and Equivalence Principles

We present high-precision timing of five millisecond pulsars (MSPs) carried out for more than seven years; four pulsars are in binary systems and one is isolated. We are able to measure the pulsars' proper motions and derive an estimate for their space velocities. The measured two-dimensional velocities are in the range 70-210 km/s, consistent with those measured for other MSPs. We also use all the available proper motion information for isolated and binary MSPs to update the known velocity distribution for these populations. As found by earlier works, we find that the velocity distribution of binary and isolated MSPs are indistinguishable with the current data. Four of the pulsars in our observing program are highly recycled with low-mass white dwarf companions and we are able to derive accurate binary parameters for these systems. For three of these binary systems we are able to place initial constraints on the pulsar masses with best-fit values in the range 1.0-1.6 M_sun. The implications of the results presented here to our understanding of binary pulsar evolution are discussed. The updated parameters for the binary systems studied here, together with recently discovered similar systems, allowed us to update previous limits on the the violation of the strong equivalence principle through the parameter |Delta| to 4.6x10^-3 (95% confidence) and the violation of Lorentz-invariance/momentum-conservation through the parameter |hat{alpha}_3| to 5.5x10^-20 (95% confidence).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observations of Energetic High Magnetic Field Pulsars with the Fermi Large Area Telescope

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On detection of the stochastic gravitational-wave background using the Parkes pulsar timing array

We search for the signature of an isotropic stochastic gravitational-wave background in pulsar timing observations using a frequency-domain correlation technique. These observations, which span roughly 12 yr, were obtained with the 64-m Parkes radio telescope augmented by public domain observations from the Arecibo Observatory. A wide range of signal processing issues unique to pulsar timing and not previously presented in the literature are discussed. These include the effects of quadratic removal, irregular sampling, and variable errors which exacerbate the spectral leakage inherent in estimating the steep red spectrum of the gravitational-wave background. These observations are found to be consistent with the null hypothesis, that no gravitational-wave background is present, with 76 percent confidence. We show that the detection statistic is dominated by the contributions of only a few pulsars because of the inhomogeneity of this data set. The issues of detecting the signature of a gravitational-wave background with future observations are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polarization observations of 20 millisecond pulsars

Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pulsar timing analysis in the presence of correlated noise

Pulsar timing observations are usually analysed with least-square-fitting procedures under the assumption that the timing residuals are uncorrelated (statistically "white"). Pulsar observers are well aware that this assumption often breaks down and causes severe errors in estimating the parameters of the timing model and their uncertainties. Ad hoc methods for minimizing these errors have been developed, but we show that they are far from optimal. Compensation for temporal correlation can be done optimally if the covariance matrix of the residuals is known using a linear transformation that whitens both the residuals and the timing model. We adopt a transformation based on the Cholesky decomposition of the covariance matrix, but the transformation is not unique. We show how to estimate the covariance matrix with sufficient accuracy to optimize the pulsar timing analysis. We also show how to apply this procedure to estimate the spectrum of any time series with a steep red power-law spectrum, including those with irregular sampling and variable error bars, which are otherwise very difficult to analyse.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pulsar Timing Arrays and their Applications

Millisecond pulsars are intrinsically very stable clocks and precise measurement of their observed pulse periods can be used to study a wide variety of astrophysical phenomena. In particular, observations of a large sample of millisecond pulsars at regular intervals, constituting a Pulsar Timing Array (PTA), can be used as a detector of low-frequency gravitational waves and to establish a standard of time independent of terrestrial atomic timescales. Three major timing array projects have been established: The European Pulsar Timing Array (EPTA), the North American pulsar timing array (NANOGrav) and the Parkes Pulsar Timing Array (PPTA). Results from the PPTA project are described in some detail and future prospects for PTA projects are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radio Properties of Rotating Radio Transients I: searches for periodicities and randomness in pulse arrival times

We have analysed the long- and short-term time dependence of the pulse arrival times and the pulse detection rates for eight Rotating Radio Transient (RRAT) sources from the Parkes Multi-beam Pulsar Survey (PMPS). We find significant periodicities in the individual pulse arrival times from six RRATs. These periodicities range from 30 minutes to 2100 days and from one to 16 independent (i.e. non-harmonically related) periodicities are detected for each RRAT. In addition, we find that pulse emission is a random (i.e. Poisson) process on short (hour-long) time scales but that most of the objects exhibit longer term (months-years) non-random behaviour. We find that PSRs J1819-1458 and J1317-5759 emit more doublets (two consecutive pulses) and triplets (three consecutive pulses) than is expected in random pulse distributions. No evidence for such an excess is found for the other RRATs. There are several different models for RRAT emission depending on both extrinsic and intrinsic factors which are consistent with these properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotation measure variations for 20 millisecond pulsars

We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Parkes Observatory Pulsar Data Archive

The Parkes pulsar data archive currently provides access to 144044 data files obtained from observations carried out at the Parkes observatory since the year 1991. Around 10^5 files are from surveys of the sky, the remainder are observations of 775 individual pulsars and their corresponding calibration signals. Survey observations are included from the Parkes 70cm and the Swinburne Intermediate Latitude surveys. Individual pulsar observations are included from young pulsar timing projects, the Parkes Pulsar Timing Array and from the PULSE@Parkes outreach program. The data files and access methods are compatible with Virtual Observatory protocols. This paper describes the data currently stored in the archive and presents ways in which these data can be searched and downloaded.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

29 Glitches Detected at Urumqi Observatory

Glitches detected in pulsar timing observations at the Nanshan radio telescope of Urumqi Observatory between 2002 July and 2008 December are presented. In total, 29 glitches were detected in 19 young pulsars, with this being the first detection of a glitch in 12 of these pulsars. Fractional glitch amplitudes range from a few parts in 10^{-11} to 3.9 x 10^{-6}. Three "slow" glitches are identified in PSRs J0631+1036, B1822-09 and B1907+10. Post-glitch recoveries differ greatly from pulsar to pulsar and for different glitches in the same pulsar. Most large glitches show some evidence for exponential post-glitch recovery on timescales of 100 -- 1000 days, but in some cases, e.g., PSR B1758-23, there is little or no recovery. Beside exponential recoveries, permanent increases in slowdown rate are seen for the two large glitches in PSRs B1800-21 and B1823-13. These and several other pulsars also show a linear increase in nudot following the partial exponential recovery, which is similar to the Vela pulsar post-glitch behaviour. Analysis of the whole sample of known glitches show that fractional glitch amplitudes are correlated with characteristic age with a peak at about 10^5 years, but there is a spread of two or three orders of magnitude at all ages. Glitch activity is positively correlated with spin-down rate, again with a wide spread of values.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detection of Gravitational Waves using Pulsar Timing

Pulsars are very stable clocks in space which have many applications to problems in physics and astrophysics. Observations of double-neutron-star binary systems have given the first observational evidence for the existence of gravitational waves (GWs) and shown that Einstein's general theory of relativity is an accurate description of gravitational interactions in the regime of strong gravity. Observations of a large sample of pulsars spread across the celestial sphere forming a "Pulsar Timing Array" (PTA), can in principle enable a positive detection of the GW background in the Galaxy. The Parkes Pulsar Timing Array (PPTA) is making precise timing measurements of 20 millisecond pulsars at three radio frequencies and is approaching the level of timing precision and data spans which are needed for GW detection. These observations will also allow us to establish a "Pulsar Timescale" and to detect or limit errors in the Solar System ephemerides used in pulsar timing analyses. Combination of PPTA data with that of other groups to form an International Pulsar Timing Array (IPTA) will enhance the sensitivity to GWs and facilitate reaching other PTA goals. The principal source of GWs at the nanoHertz frequencies to which PTAs are sensitive is believed to be super-massive binary black holes in the cores of distant galaxies. Current results do not signficantly limit models for formation of such black-hole binary systems, but in a few years we expect that PTAs will either detect GWs or seriously constrain current ideas about black-hole formation and galaxy mergers. Future instruments such as the Square Kilometre Array (SKA) should not only detect GWs from astrophysical sources but also enable detailed studies of the sources and the gravitational theories used to account for the GW emission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discovery of a highly energetic pulsar associated with IGR J14003-6326 in a young uncataloged Galactic supernova remnant G310.6-1.6

We report the discovery of 31.18 ms pulsations from the INTEGRAL source IGR J14003-6326 using the Rossi X-ray Timing Explorer (RXTE). This pulsar is most likely associated with the bright Chandra X-ray point source lying at the center of G310.6-1.6, a previously unrecognised Galactic composite supernova remnant with a bright central non-thermal radio and X-ray nebula, taken to be the pulsar wind nebula (PWN). PSR J1400-6325 is amongst the most energetic rotation-powered pulsars in the Galaxy, with a spin-down luminosity of Edot = 5.1E+37 erg.s-1. In the rotating dipole model, the surface dipole magnetic field strength is B_s = 1.1E+12 G and the characteristic age tau_c = P/2Pdot = 12.7 kyr. The high spin-down power is consistent with the hard spectral indices of the pulsar and the nebula of 1.22 +/- 0.15 and 1.83 +/- 0.08, respectively, and a 2-10 keV flux ratio F_PWN/F_PSR ~ 8. Follow-up Parkes observations resulted in the detection of radio emission at 10 and 20 cm from PSR J1400-6325 at a dispersion measure of ~ 560 cm-3 pc, which implies a relatively large distance of 10 +/- 3 kpc. However, the resulting location off the Galactic Plane of ~ 280 pc would be much larger than the typical thickness of the molecular disk, and we argue that G310.6-1.6 lies at a distance of ~ 7 kpc. There is no gamma-ray counterpart to the nebula or pulsar in the Fermi data published so far. A multi-wavelength study of this new composite supernova remnant, from radio to very-high energy gamma-rays, suggests a young (< 1000 yr) system, formed by a sub-energetic (~ 1E+50 ergs), low ejecta mass (M_ej ~ 3 Msun) SN explosion that occurred in a low-density environment (n_0 ~ 0.01 cm-3).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measuring the mass of solar system planets using pulsar timing

High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Observations and Modelling of Relativistic Spin Precession in PSR J1141-6545

Observations of the binary pulsar PSR J1141-6545 using the Parkes radio telescope over 9.3 years show clear time-variations in pulse width, shape and polarization. We interpret these variations in terms of relativistic precession of the pulsar spin axis about the total angular momentum vector of the system. Over the nine years, the pulse width at the 50% level has changed by more than a factor of three. Large variations have also been observed in the 1400-MHz mean flux density. The pulse polarization has been monitored since 2004 April using digital filterbank systems and also shows large and systematic variations in both linear and circular polarization. Position angle variations, both across the pulse profile and over the data span, are complex, with major differences between the central and outer parts of the pulse profile. Modelling of the observed position angle variations by relativistic precession of the pulsar spin axis shows that the spin-orbit misalignment angle is about 110 deg and that the precessional phase has passed through 180 deg during the course of our observations. At the start of our observations, the line-of-sight impact parameter was about 4 deg in magnitude and it reached a minimum very close to 0 deg around early 2007, consistent with the observed pulse width variations. We have therefore mapped approximately one half of the emission beam, showing that it is very asymmetric with respect to the magnetic axis. The derived precessional parameters imply that the pre-supernova star had a mass of about 2 Msun and that the supernova recoil kick velocity was relatively small. With the reversal in the rate of change of the impact parameter, we predict that over the next decade we will see a reversed "replay" of the variations observed in the past decade.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Observations of radio pulses from CU Virginis

The magnetic chemically peculiar star CU Virginis is a unique astrophysical laboratory for stellar magnetospheres and coherent emission processes. It is the only known main sequence star to emit a radio pulse every rotation period. Here we report on new observations of the CU Virginis pulse profile in the 13 and 20\,cm radio bands. The profile is known to be characterised by two peaks of 100$\%$ circularly polarised emission that are thought to arise in an electron-cyclotron maser mechanism. We find that the trailing peak is stable at both 13 and 20\,cm, whereas the leading peak is intermittent at 13\,cm. Our measured pulse arrival times confirm the discrepancy previously reported between the putative stellar rotation rates measured with optical data and with radio observations. We suggest that this period discrepancy might be caused by an unknown companion or by instabilities in the emission region. Regular long-term pulse timing and simultaneous multi-wavelength observations are essential to clarify the behaviour of this emerging class of transient radio source.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pulsars with the Australian Square Kilometre Array Pathfinder

The Australian Square Kilometre Array Pathfinder (ASKAP) is a 36-element array with a 30-square-degree field of view being built at the proposed SKA site in Western Australia. We are conducting a Design Study for pulsar observations with ASKAP, planning both timing and search observations. We provide an overview of the ASKAP telescope and an update on pulsar-related progress.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Searches for gravitational waves from known pulsars with S5 LIGO data

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Sensitivity of the Parkes Pulsar Timing Array to Individual Sources of Gravitational Waves

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Wide radio beams from gamma-ray pulsars

We investigate the radio and gamma-ray beaming properties of normal and millisecond pulsars by selecting two samples from the known populations. The first, Sample G, contains pulsars which are detectable in blind searches of gamma-ray data from the Fermi Large Area Telescope. The second, Sample R, contains pulsars detectable in blind radio searches which have spin-down luminosities Edot > 10^{34} erg/s. We analyse the fraction of the gamma-ray-selected Sample G which have detectable radio pulses and the fraction of the radio-selected Sample R which have detectable gamma-ray pulses. Twenty of our 35 Sample G pulsars have already observed radio pulses. This rules out low-altitude polar-cap beaming models if, as is currently believed, gamma-ray beams are generated in the outer magnetosphere and are very wide. We further find that, for the highest-Edot pulsars, the radio and gamma-ray beams have comparable beaming factors, i.e., the beams cover similar regions of the sky as the star rotates. For lower-Edot gamma-ray emitting pulsars, the radio beams have about half of the gamma-ray sky coverage. These results suggest that, for high-Edot young and millisecond pulsars, the radio emission originates in wide beams from regions high in the pulsar magnetosphere, probably close to the null-charge surface and to the gamma-ray emitting regions. Furthermore, it suggests that for these high-Edot pulsars, as in the gamma-ray case, features in the radio profile represent caustics in the emission beam pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High resolution 36 GHz imaging of the Supernova Remnant of SN1987A

The aftermath of supernova (SN) 1987A continues to provide spectacular insights into the interaction between a SN blastwave and its circumstellar en- vironment. We here present 36 GHz observations from the Australia Telescope Compact Array of the radio remnant of SN 1987A. These new images, taken in 2008 Apr and 2008 Oct, substantially extend the frequency range of an ongo- ing monitoring and imaging program conducted between 1.4 and 20 GHz. Our 36.2 GHz images have a diffraction-limited angular resolution of 0.3-0.4 arcseconds, which covers the gap between high resolution, low dynamic range VLBI images of the remnant and low resolution, high dynamic range images at frequencies between 1 and 20 GHz. The radio morphology of the remnant at 36 GHz is an elliptical ring with enhanced emission on the eastern and western sides, similar to that seen previously at lower frequencies. Model fits to the data in the Fourier domain show that the emitting region is consistent with a thick inclined torus of mean radius 0.85 arcsec, and a 2008 Oct flux density of 27 +/- 6 mJy at 36.2 GHz. The spectral index for the remnant at this epoch, determined between 1.4 GHz and 36.2 GHz, is -0.83. There is tentative evidence for an unresolved central source with flatter spectral index.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High Resolution X-ray Imaging of Supernova Remnant 1987A

We report observations of the remnant of Supernova 1987A with the High Resolution Camera (HRC) onboard the Chandra X-ray Observatory. A direct image from the HRC resolves the annular structure of the X-ray remnant, confirming the morphology previously inferred by deconvolution of lower resolution data from the Advanced CCD Imaging Spectrometer. Detailed spatial modeling shows that the a thin ring plus a thin shell gives statistically the best description of the overall remnant structure, and suggests an outer radius 0.96" +/- 0.05" +/- 0.03" for the X-ray-emitting region, with the two uncertainties corresponding to the statistical and systematic errors, respectively. This is very similar to the radius determined by a similar modeling technique for the radio shell at a comparable epoch, in contrast to previous claims that the remnant is 10-15% smaller at X-rays than in the radio band. The HRC observations put a flux limit of 0.010 cts/s (99% confidence level, 0.08-10 keV range) on any compact source at the remnant center. Assuming the same foreground neutral hydrogen column density as towards the remnant, this allows us to rule out an unobscured neutron star with surface temperature T^\infty > 2.5MK observed at infinity, a bright pulsar wind nebula or a magnetar.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Multi-frequency Radio Measurements of SN 1987A over 22 Years

We present extensive observations of the radio emission from the remnant of SN 1987A made with the Australia Telescope Compact Array (ATCA), since the first detection of the remnant in 1990. The radio emission has evolved in time providing unique information on the interaction of the supernova shock with the circumstellar medium. We particularly focus on the monitoring observations at 1.4, 2.4, 4.8 and 8.6 GHz, which have been made at intervals of 4-6 weeks. The flux density data show that the remnant brightness is now increasing exponentially, while the radio spectrum is flattening. The current spectral index value of -0.68 represents an 18+/-3% increase over the last 8 years. The exponential trend in the flux is also found in the ATCA imaging observations at 9 GHz, which have been made since 1992, approximately twice a year, as well as in the 843 MHz data set from the Molonglo Observatory Synthesis Telescope from 1987 to March 2007. Comparisons with data at different wavelengths (X-ray, Hα) are made. The rich data set that has been assembled in the last 22 years forms a basis for a better understanding of the evolution of the supernova remnant.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

PSRCHIVE and PSRFITS: Definition of the Stokes Parameters and Instrumental Basis Conventions

This paper defines the mathematical convention adopted to describe an electromagnetic wave and its polarisation state, as implemented in the PSRCHIVE software and represented in the PSRFITS definition. Contrast is made between the convention that has been widely accepted by pulsar astronomers and the IAU/IEEE definitions of the Stokes parameters. The former is adopted as the PSR/IEEE convention, and a set of useful parameters are presented for describing the differences between the PSR/IEEE standard and the conventions (either implicit or explicit) that form part of the design of observatory instrumentation. To aid in the empirical determination of instrumental convention parameters, well-calibrated average polarisation profiles of PSR J0304+1932 and PSR J0742-2822 are presented at radio wavelengths of approximately 10, 20, and 40 cm.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Pulsar Timing with the Parkes Radio Telescope for the Fermi Mission

We report here on two years of timing of 168 pulsars using the Parkes radio telescope. The vast majority of these pulsars have spin-down luminosities in excess of 10^34 erg/s and are prime target candidates to be detected in gamma-rays by the Fermi Gamma-Ray Space Telescope. We provide the ephemerides for the ten pulsars being timed at Parkes which have been detected by Fermi in its first year of operation. These ephemerides, in conjunction with the publicly available photon list, can be used to generate gamma-ray profiles from the Fermi archive. We will make the ephemerides of any pulsars of interest available to the community upon request. In addition to the timing ephemerides, we present the parameters for 14 glitches which have occurred in 13 pulsars, seven of which have no previously known glitch history. The Parkes timing programme, in conjunction with Fermi observations, is expected to continue for at least the next four years.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status Update of the Parkes Pulsar Timing Array

The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a successful detection of gravitational waves within a decade, or alternatively to rule out all current predictions for gravitational wave backgrounds formed by supermassive black-hole mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Timing Observations of Rotating Radio Transients

We present radio timing measurements of six rotating radio transient (RRAT) sources discovered in the Parkes Multibeam Pulsar Survey. These provide four new phase-connected timing solutions and two updated ones, making a total of seven of the original 11 reported RRATs now with high-precision rotational and astrometric parameters. Three of these seven RRATs have magnetic fields greater than 10^13 G, with spin-down properties similar to those of the magnetars and X-ray detected isolated neutron stars. Another two of these RRATs have long periods and large characteristic ages, and lie near the `death-line' for radio pulsar emission. The remaining two RRATs with timing solutions have properties typical of the bulk of the pulsar population. The new solutions offer insights into what might be responsible for the unusual emission properties. We demonstrate that the RRATs have significantly longer periods and higher magnetic fields than normal radio pulsars, and find no correlation with other spin-down parameters. These solutions also provide precise positions, which will facilitate follow-up studies at high energies, crucial for relating these sources with other neutron star populations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Timing stability of millisecond pulsars and prospects for gravitational-wave detection

Analysis of high-precision timing observations of an array of approx. 20 millisecond pulsars (a so-called "timing array") may ultimately result in the detection of a stochastic gravitational-wave background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a five-year timescale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within five years to a decade are good.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Gravitational wave detection using pulsars: status of the Parkes Pulsar Timing Array project

The first direct detection of gravitational waves may be made through observations of pulsars. The principal aim of pulsar timing array projects being carried out worldwide is to detect ultra-low frequency gravitational waves (f ~ 10^-9 to 10^-8 Hz). Such waves are expected to be caused by coalescing supermassive binary black holes in the cores of merged galaxies. It is also possible that a detectable signal could have been produced in the inflationary era or by cosmic strings. In this paper we review the current status of the Parkes Pulsar Timing Array project (the only such project in the Southern hemisphere) and compare the pulsar timing technique with other forms of gravitational-wave detection such as ground- and space-based interferometer systems.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

High-Resolution Timing Observations of Spin-Powered Pulsars with the AGILE Gamma-Ray Telescope

AGILE is a small gamma-ray astronomy satellite mission of the Italian Space Agency dedicated to high-energy astrophysics launched in 2007 April. Its 1 microsecond absolute time tagging capability coupled with a good sensitivity in the 30 MeV-30 GeV range, with simultaneous X-ray monitoring in the 18-60 keV band, makes it perfectly suited for the study of gamma-ray pulsars following up on the CGRO/EGRET heritage. In this paper we present the first AGILE timing results on the known gamma-ray pulsars Vela, Crab, Geminga and B 1706-44. The data were collected from 2007 July to 2008 April, exploiting the mission Science Verification Phase, the Instrument Timing Calibration and the early Observing Pointing Program. Thanks to its large field of view, AGILE collected a large number of gamma-ray photons from these pulsars (about 10,000 pulsed counts for Vela) in only few months of observations. The coupling of AGILE timing capabilities, simultaneous radio/X-ray monitoring and new tools aimed at precise photon phasing, exploiting also timing noise correction, unveiled new interesting features at sub-millisecond level in the pulsars' high-energy light-curves.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Precision timing of PSR J0437-4715: an accurate pulsar distance, a high pulsar mass and a limit on the variation of Newton's gravitational constant

Analysis of ten years of high-precision timing data on the millisecond pulsar PSR J0437-4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, Pbdot, determined at the 1.5% level of precision (Dk = 157.0 +/- 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 Solar System ephemerides. The precise measurement of Pbdot allows a limit on the variation of Newton's gravitational constant, |Gdot/G| < 23 x 10^{-12} 1/yr. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a(Sun)/c| < 1.5 x 10^{-18} 1/s at 95% confidence, and derive a pulsar mass, m(psr) = 1.76 +/- 0.20 M, one of the highest estimates so far obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Pulsar rotation measures and the large-scale structure of Galactic magnetic field

The large-scale magnetic field of our Galaxy can be probed in three dimensions using Faraday rotation of pulsar signals. We report on the determination of 223 rotation measures from polarization observations of relatively distant southern pulsars made using the Parkes radio telescope. Combined with previously published observations these data give clear evidence for large-scale counterclockwise fields (viewed from the north Galactic pole) in the spiral arms interior to the Sun and weaker evidence for a counterclockwise field in the Perseus arm. However, in interarm regions, including the Solar neighbourhood, we present evidence that suggests that large-scale fields are clockwise. We propose that the large-scale Galactic magnetic field has a bisymmetric structure with reversals on the boundaries of the spiral arms. Streaming motions associated with spiral density waves can directly generate such a structure from an initial inwardly directed radial field. Large-scale fields increase toward the Galactic Center, with a mean value of about 2~$μ$G in the Solar neighbourhood and 4~$μ$G at a Galactocentric radius of 3 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

PSR J1833-1034: Discovery of the Central Young Pulsar in the Supernova Remnant G21.5-0.9

We have discovered the pulsar associated with the supernova remnant G21.5-0.9. PSR J1833-1034, with spin period P = 61.8 ms and dispersion measure 169 pc/cc, is very faint, with pulse-averaged flux density of approximately 70 microJy at a frequency of 1.4 GHz, and was first detected in a deep search with the Parkes telescope. Subsequent observations with Parkes and the Green Bank Telescope have confirmed this detection, and yield a period derivative dP/dt = 2.02e-13. These spin parameters imply a characteristic age tau_c = 4.8 kyr and a spin-down luminosity dE/dt = 3.3e37 ergs/s, the latter value exceeded only by the Crab pulsar among the rotation-powered pulsars known in our Galaxy. The pulsar has an unusually steep radio spectrum in the 0.8-2.0 GHz range, with power law index ~3.0, and a narrow single-peaked pulse profile with full-width at half maximum of 0.04P. We have analyzed 350 ks of archival Chandra X-ray Observatory High Resolution Camera (HRC) data, and find a point-like source of luminosity approximately 3e-5dE/dt, offset from the center of an elliptical region of size ~7"x5" and luminosity approximately 1e-3dE/dt within which likely lies the pulsar wind termination shock. We have searched for X-ray pulsations in a 30 ks HRC observation without success, deriving a pulsed fraction upper limit for a sinusoidal pulse shape of about 70% of the pulsar flux. We revisit the distance to G21.5-0.9 based on HI and CO observations, arguing that it is 4.7+-0.4 kpc. We use existing X-ray and radio observations of the pulsar wind nebula, along with the measured properties of its engine and a recent detection of the supernova remnant shell, to argue that G21.5-0.9 and PSR J1833-1034 are much younger than tau_c, and likely their true age is ~<1000 yr. In that case, the initial spin period of the pulsar was >~55 ms.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

A Double-Pulsar System - A Rare Laboratory for Relativistic Gravity and Plasma Physics

The clock-like properties of pulsars moving in the gravitational fields of their unseen neutron-star companions have allowed unique tests of general relativity and provided evidence for gravitational radiation. We report here the detection of the 2.8-sec pulsar J0737-3039B as the companion to the 23-ms pulsar J0737-3039A in a highly-relativistic double-neutron-star system, allowing unprecedented tests of fundamental gravitational physics. We observe a short eclipse of J0737-3039A by J0737-3039B and orbital modulation of the flux density and pulse shape of J0737-3039B, probably due to the influence of J0737-3039A's energy flux upon its magnetosphere. These effects will allow us to probe magneto-ionic properties of a pulsar magnetosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Parkes multibeam pulsar survey: IV. Discovery of 180 pulsars and parameters for 281 previously known pulsars

The Parkes multibeam pulsar survey has led to the discovery of more than 700 pulsars. In this paper, we provide timing solutions, flux densities and pulse profiles for 180 of these new discoveries. Two pulsars, PSRs J1736-2843 and J1847-0130 have rotational periods P > 6s and are therefore among the slowest rotating radio pulsars known. Conversely, with P = 1.8ms, PSR J1843-1113 has the third shortest period of pulsars currently known. This pulsar and PSR J1905+0400 (P = 3.8ms) are both solitary. We also provide orbital parameters for a new binary system, PSR J1420-5625, which has P = 34ms, an orbital period of 40 days and a minimum companion mass of 0.4 solar masses. The 10 degree-wide strip along the Galactic plane that was surveyed is known to contain 264 radio pulsars that were discovered prior to the multibeam pulsar survey. We have redetected almost all of these pulsars and provide new dispersion measure values and flux densities at 20cm for the redetected pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The Parkes Multibeam Pulsar Survey

The Parkes multibeam pulsar survey is a high-frequency, fast-sampled survey of the Galactic Plane, expected to discover at least 500 new pulsars. To date, over 200 pulsars have been found, including several young pulsars and at least one with a very high magnetic field. Seven of the new stars are in binary systems; this number includes one probable double-neutron-star system, and one pulsar with an extremely massive companion.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

RXTE Absolute Timing Results for the Pulsars B1821-24 and B1509-58

Observations with the Rossi X-ray Timing Explorer and the Jodrell Bank, Parkes, and Green Bank telescopes have enabled us to determine the time delay between radio and X-ray pulses in the two isolated pulsars B1821-24 and B1509-58. For the former we find that the narrow X-ray and radio pulse components are close to being coincident in time, with the radio peak leading by 0.02 period (60 +/- 20 microsec), while the wide X-ray pulse component lags the last of the two wider radio components by about 0.08 period. For the latter pulsar we find, using the standard value for the dispersion measure, that the X-ray pulse lags the radio by about 0.27 period, with no evidence for any energy-dependence in the range 2-100 keV. However, uncertainties in the history of the dispersion measure for this pulsar make a comparison with previous results difficult. It is clear that there are no perceptable variations in either the lag or the dispersion measure at time scales of a year or less.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Timing Measurements and Their Implications for Four Binary Millisecond Pulsars

We present timing observations of four millisecond pulsars, using data obtained over three years at the ATNF Parkes and NRAL Jodrell Bank radio telescopes. Astrometric, spin, and binary parameters are updated, and substantially improved for three pulsars, PSRs J0613-0200, J1045-4509 and J1643-1224. We have measured the time variation of the projected semi-major axis of the PSR J0437-4715 orbit due to its proper motion, and use it to constrain the inclination of the orbit and the mass of the companion. Some evidence is found for changes in the dispersion measures of PSRs J1045-4509 and J1643-1224. Limits are placed on the existence of planetary mass companions, ruling out companions with masses and orbits similar to the terrestrial planets of the solar system for eight pulsars.

0 citations0 reviewsNo code linkedOpen access