A VLA Search for the Geminga Pulsar: A Bayesian Limit on a Scintillating Source

We derive an upper limit of 3 mJy (95% confidence) for the flux density at 317 MHz of the Geminga pulsar (J0633+1746). Our results are based on 7 hours of fast-sampled VLA data, which we averaged synchronously with the pulse period using a period model based on CGRO/EGRET gamma-ray data. Our limit accounts for the fact that this pulsar is most likely subject to interstellar scintillations on a timescale much shorter than our observing span. Our Bayesian method is quite general and can be applied to calculate the fluxes of other scintillated sources. We also present a Bayesian technique for calculating the flux in a pulsed signal of unknown width and phase. Comparing our upper limit of 3 mJy with the quoted flux density of Geminga at 102 MHz, we calculate a lower limit to its spectral index of 2.7. We discuss some possible reasons for Geminga's weakness at radio wavelengths, and the likelihood that many of the unidentified EGRET sources are also radio-quiet or radio-weak Geminga-like pulsars.Comment: 27 pages, including figures. Published in Ap

On the Mass Distribution and Birth Masses of Neutron Stars

We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as of slow pulsars, which are all believed to be near their birth masses, has a mean of 1.28 M_solar and a dispersion of 0.24 M_solar. These values are consistent with expectations for neutron star formation in core-collapse supernovae. On the other hand, double neutron stars, which are also believed to be near their birth masses, have a much narrower mass distribution, peaking at 1.33 M_solar but with a dispersion of only 0.05 M_solar. Such a small dispersion cannot easily be understood and perhaps points to a particular and rare formation channel. The mass distribution of neutron stars that have been recycled has a mean of 1.48 M_solar and a dispersion of 0.2 M_solar, consistent with the expectation that they have experienced extended mass accretion episodes. The fact that only a very small fraction of recycled neutron stars in the inferred distribution have masses that exceed ~2 M_solar suggests that only a few of these neutron stars cross the mass threshold to form low mass black holes.Comment: replaced to match the ApJ published version; includes most up to date pulsar masses and a revised fitting formul

The Spin Distribution of Fast Spinning Neutron Stars in Low Mass X-Ray Binaries: Evidence for Two Sub-Populations

We study the current sample of rapidly rotating neutron stars in both accreting and non-accreting binaries in order to determine whether the spin distribution of accreting neutron stars in low-mass X-ray binaries can be reconciled with current accretion torque models. We perform a statistical analysis of the spin distributions and show that there is evidence for two sub-populations among low-mass X-ray binaries, one at relatively low spin frequency, with an average of ~300 Hz and a broad spread, and a peaked population at higher frequency with average spin frequency of ~575 Hz. We show that the two sub-populations are separated by a cut-point at a frequency of ~540 Hz. We also show that the spin frequency of radio millisecond pulsars does not follow a log-normal distribution and shows no evidence for the existence of distinct sub-populations. We discuss the uncertainties of different accretion models and speculate that either the accreting neutron star cut-point marks the onset of gravitational waves as an efficient mechanism to remove angular momentum or some of the neutron stars in the fast sub-population do not evolve into radio millisecond pulsars.Comment: Submitted to Ap

Gravitational waves: search results, data analysis and parameter estimation

The Amaldi 10 Parallel Session C2 on gravitational wave (GW) search results, data analysis and parameter estimation included three lively sessions of lectures by 13 presenters, and 34 posters. The talks and posters covered a huge range of material, including results and analysis techniques for ground-based GW detectors, targeting anticipated signals from different astrophysical sources: compact binary inspiral, merger and ringdown; GW bursts from intermediate mass binary black hole mergers, cosmic string cusps, core-collapse supernovae, and other unmodeled sources; continuous waves from spinning neutron stars; and a stochastic GW background. There was considerable emphasis on Bayesian techniques for estimating the parameters of coalescing compact binary systems from the gravitational waveforms extracted from the data from the advanced detector network. This included methods to distinguish deviations of the signals from what is expected in the context of General Relativity

Practical Methods for Continuous Gravitational Wave Detection using Pulsar Timing Data

Gravitational Waves (GWs) are tiny ripples in the fabric of space-time predicted by Einstein's General Relativity. Pulsar timing arrays (PTAs) are well poised to detect low frequency ($10^{-9}$ -- $10^{-7}$ Hz) GWs in the near future. There has been a significant amount of research into the detection of a stochastic background of GWs from supermassive black hole binaries (SMBHBs). Recent work has shown that single continuous sources standing out above the background may be detectable by PTAs operating at a sensitivity sufficient to detect the stochastic background. The most likely sources of continuous GWs in the pulsar timing frequency band are extremely massive and/or nearby SMBHBs. In this paper we present detection strategies including various forms of matched filtering and power spectral summing. We determine the efficacy and computational cost of such strategies. It is shown that it is computationally infeasible to use an optimal matched filter including the poorly constrained pulsar distances with a grid based method. We show that an Earth-term-matched filter constructed using only the correlated signal terms is both computationally viable and highly sensitive to GW signals. This technique is only a factor of two less sensitive than the computationally unrealizable optimal matched filter and a factor of two more sensitive than a power spectral summing technique. We further show that a pairwise matched filter, taking the pulsar distances into account is comparable to the optimal matched filter for the single template case and comparable to the Earth-term-matched filter for many search templates. Finally, using simulated data optimal quality, we place a theoretical minimum detectable strain amplitude of $h>2\times 10^{-15}$ from continuous GWs at frequencies on the order $\sim1/T_{\rm obs}$.Comment: submitted to Ap

Detection and localization of continuous gravitational waves with pulsar timing arrays: the role of pulsar terms

A pulsar timing array is a Galactic-scale detector of nanohertz gravitational waves (GWs). Its target signals contain two components: the `Earth term' and the `pulsar term' corresponding to GWs incident on the Earth and pulsar respectively. In this work we present a Frequentist method for the detection and localization of continuous waves that takes into account the pulsar term and is significantly faster than existing methods. We investigate the role of pulsar terms by comparing a full-signal search with an Earth-term-only search for non-evolving black hole binaries. By applying the method to synthetic data sets, we find that (i) a full-signal search can slightly improve the detection probability (by about five percent); (ii) sky localization is biased if only Earth terms are searched for and the inclusion of pulsar terms is critical to remove such a bias; (iii) in the case of strong detections (with signal-to-noise ratio $\gtrsim$ 30), it may be possible to improve pulsar distance estimation through GW measurements.Comment: 12 pages, 9 figures, typos corrected. To match the published version. Code implementing this method is available at the PPTA Wiki pag