Pulsar Timing Arrays and Gravity Tests in the Radiative Regime

In this paper, we focus on testing gravity theories in the radiative regime using pulsar timing array observations. After reviewing current techniques to measure the dispersion and alternative polarization of gravitational waves, we extend the framework to the most general situations, where the combinations of a massive graviton and alternative polarization modes are considered. The atlas of the Hellings-Downs functions is completed by the new calculations for these dispersive alternative polarization modes. We find that each mode and corresponding graviton mass introduce characteristic features in the Hellings-Downs function. Thus, in principal, we can not only detect each polarization mode, measure the corresponding graviton mass, but also discriminate the different scenarios. In this way, we can test gravity theories in the radiative regime in a generalized fashion, and such method is a direct experiment, where one can address the gauge symmetry of the gravity theories in their linearised limits. Although current pulsar timing still lacks enough stable pulsars and sensitivity for such practices, we expect that future telescopes with larger collecting area could make such experiments be feasible.

Model-based asymptotically optimal dispersion measure correction for pulsar timing

In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration.

PEACE: Pulsar Evaluation Algorithm for Candidate Extraction -- A software package for post-analysis processing of pulsar survey candidates

Modern radio pulsar surveys produce a large volume of prospective candidates, the majority of which are polluted by human-created radio frequency interference or other forms of noise. Typically, large numbers of candidates need to be visually inspected in order to determine if they are real pulsars. This process can be labor intensive. In this paper, we introduce an algorithm called PEACE (Pulsar Evaluation Algorithm for Candidate Extraction) which improves the efficiency of identifying pulsar signals. The algorithm ranks the candidates based on a score function. Unlike popular machine-learning based algorithms, no prior training data sets are required. This algorithm has been applied to data from several large-scale radio pulsar surveys. Using the human-based ranking results generated by students in the Arecibo Remote Command enter programme, the statistical performance of PEACE was evaluated. It was found that PEACE ranked 68% of the student-identified pulsars within the top 0.17% of sorted candidates, 95% within the top 0.34%, and 100% within the top 3.7%. This clearly demonstrates that PEACE significantly increases the pulsar identification rate by a factor of about 50 to 1000. To date, PEACE has been directly responsible for the discovery of 47 new pulsars, 5 of which are millisecond pulsars that may be useful for pulsar timing based gravitational-wave detection projects.

Application of the Gaussian mixture model in pulsar astronomy -- pulsar classification and candidates ranking for {it Fermi} 2FGL catalog

Machine learning, algorithms to extract empirical knowledge from data, can be used to classify data, which is one of the most common tasks in observational astronomy. In this paper, we focus on Bayesian data classification algorithms using the Gaussian mixture model and show two applications in pulsar astronomy. After reviewing the Gaussian mixture model and the related Expectation-Maximization algorithm, we present a data classification method using the Neyman-Pearson test. To demonstrate the method, we apply the algorithm to two classification problems. Firstly, it is applied to the well known period-period derivative diagram, where we find that the pulsar distribution can be modeled with six Gaussian clusters, with two clusters for millisecond pulsars (recycled pulsars) and the rest for normal pulsars. From this distribution, we derive an empirical definition for millisecond pulsars as $frac{dot{P}}{10^{-17}} leq3.23(frac{P}{100 textrm{ms}})^{-2.34}$. The two millisecond pulsar clusters may have different evolutionary origins, since the companion stars to these pulsars in the two clusters show different chemical composition. Four clusters are found for normal pulsars. Possible implications for these clusters are also discussed. Our second example is to calculate the likelihood of unidentified textit{Fermi} point sources being pulsars and rank them accordingly. In the ranked point source list, the top 5% sources contain 50% known pulsars, the top 50% contain 99% known pulsars, and no known active galaxy (the other major population) appears in the top 6%. Such a ranked list can be used to help the future follow-up observations for finding pulsars in unidentified textit{Fermi} point sources.

The optimal schedule for pulsar timing array observations

In order to maximize the sensitivity of pulsar timing arrays to a stochastic gravitational wave background, we present computational techniques to optimize observing schedules. The techniques are applicable to both single and multi-telescope experiments. The observing schedule is optimized for each telescope by adjusting the observing time allocated to each pulsar while keeping the total amount of observing time constant. The optimized schedule depends on the timing noise characteristics of each individual pulsar as well as the performance of instrumentation. Several examples are given to illustrate the effects of different types of noise. A method to select the most suitable pulsars to be included in a pulsar timing array project is also presented.

Testing Gravity Theories In The Radiative Regime Using Pulsar Timing Arrays

General relativity has predicted the existence of gravitational waves (GW), which are waves of the distortions of space-time with two degrees of polarization and the propagation speed of light. Alternative theories predict more polarizations, up to a maximum of six, and possible deviation of propagation speed from the light speed. The present paper reviews recent proposals to test the gravity theories in the radiation regime by observing GWs using pulsar timing arrays.

Gravitational wave astronomy of single sources with a pulsar timing array

Abbreviated: We investigate the potential of detecting the gravitational wave from individual binary black hole systems using pulsar timing arrays (PTAs) and calculate the accuracy for determining the GW properties. This is done in a consistent analysis, which at the same time accounts for the measurement of the pulsar distances via the timing parallax. We find that, at low redshift, a PTA is able to detect the nano-Hertz GW from super massive black hole binary systems with masses of $sim10^8 - 10^{10},M_{sun}$ less than $sim10^5$,years before the final merger, and those with less than $sim10^3 - 10^4$ years before merger may allow us to detect the evolution of binaries. We derive an analytical expression to describe the accuracy of a pulsar distance measurement via timing parallax. We consider five years of bi-weekly observations at a precision of 15,ns for close-by ($sim 0.5 - 1$,kpc) pulsars. Timing twenty pulsars would allow us to detect a GW source with an amplitude larger than $5times 10^{-17}$. We calculate the corresponding GW and binary orbital parameters and their measurement precision. The accuracy of measuring the binary orbital inclination angle, the sky position, and the GW frequency are calculated as functions of the GW amplitude. We note that the pulsar term, which is commonly regarded as noise, is essential for obtaining an accurate measurement for the GW source location. We also show that utilizing the information encoded in the GW signal passing the Earth also increases the accuracy of pulsar distance measurements. If the gravitational wave is strong enough, one can achieve sub-parsec distance measurements for nearby pulsars with distance less than $sim 0.5 - 1$,kpc.

Detecting massive gravitons using pulsar timing arrays

Massive gravitons are features of some alternatives to general relativity. This has motivated experiments and observations that, so far, have been consistent with the zero mass graviton of general relativity, but further tests will be valuable. A basis for new tests may be the high sensitivity gravitational wave experiments that are now being performed, and the higher sensitivity experiments that are being planned. In these experiments it should be feasible to detect low levels of dispersion due to nonzero graviton mass. One of the most promising techniques for such a detection may be the pulsar timing program that is sensitive to nano-Hertz gravitational waves. Here we present some details of such a detection scheme. The pulsar timing response to a gravitational wave background with the massive graviton is calculated, and the algorithm to detect the massive graviton is presented. We conclude that, with 90% probability, massles gravitons can be distinguished from gravitons heavier than $3times 10^{-22}$ eV (Compton wave length $lambda_{rm g}=4.1 times 10^{12}$ km), if biweekly observation of 60 pulsars are performed for 5 years with pulsar RMS timing accuracy of 100 ns. If 60 pulsars are observed for 10 years with the same accuracy, the detectable graviton mass is reduced to $5times 10^{-23}$ eV ($lambda_{rm g}=2.5 times 10^{13}$ km); for 5-year observations of 100 or 300 pulsars, the sensitivity is respectively $2.5times 10^{-22}$ ($lambda_{rm g}=5.0times 10^{12}$ km) and $10^{-22}$ eV ($lambda_{rm g}=1.2times 10^{13}$ km). Finally, a 10-year observation of 300 pulsars with 100 ns timing accuracy would probe graviton masses down to $3times 10^{-23}$ eV ($lambda_{rm g}=4.1times 10^{13}$ km).

Lower bounds of altitudes for pulsar $gamma$-ray radiation

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $gamma-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $gamma$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $gamma$-B optical depth for $gamma$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.