Articles

GRAVITATIONAL WAVES FROM KNOWN PULSARS: RESULTS FROM THE INITIAL DETECTOR ERA

, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2014 April 2 © 2014. The American Astronomical Society. All rights reserved.
, , Citation J. Aasi et al 2014 ApJ 785 119 DOI 10.1088/0004-637X/785/2/119

0004-637X/785/2/119

ABSTRACT

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Pulsars are spinning, magnetized neutron stars with slowly decreasing rotation rates. In the model of a triaxial ellipsoid star, a deformation (possibly from shear strains in the solid part(s) of the star, or from magnetic stresses) can appear as a time-varying quadrupole moment as the star rotates. The observed loss of rotational energy, known as the spin-down luminosity (given by $\dot{E} = I_{zz}\Omega |\dot{\Omega }| = 4\pi ^2I_{zz}f_{\rm rot} |\dot{f}_{\rm rot}|$, where Izz is the moment of inertia around the principal axis (aligned with the rotation axis), frot is the rotation frequency, and $\dot{f}_{\rm rot}$ is the rotational frequency derivative) provides a huge reservoir of energy. Along with magnetic dipole radiation some fraction of this reservoir is potentially dissipated through gravitational wave emission (see Shklovskii 1969; Ostriker & Gunn 1969; Ferrari & Ruffini 1969; Melosh 1969 for four contemporaneous calculations of gravitational wave emission from soon after pulsars were discovered, or e.g., Owen 2006 for a review of more recent emission mechanism calculations). Known pulsars usually have precisely determined frequency evolutions and sky-positions making them ideal targets for gravitational wave detectors. If a pulsar is monitored regularly through electromagnetic observations it can yield a coherent phase model, which allows gravitational wave data to be coherently integrated over months or years.

Since the initial science data runs of the Laser Interferometric Gravitational-wave Observatory (LIGO), Virgo and GEO600, searches have been performed for continuous quasi-monochromatic gravitational wave emission from many known pulsars (Abbott et al. 2004, 2005, 2007a, 2008, 2010; Abadie et al. 2011). Most recently 116 known pulsars were targeted using data from LIGO's fifth science run (S5; Abbott et al. 2010), and the Vela pulsar (J0835−4510) was targeted using data from Virgo's second science run (VSR2). These searches reported no detections, but provided upper limits on the gravitational wave amplitude from the sources and surpassed the so-called spin-down limit (see Section 1.1) for the Crab and Vela pulsars.

We aim here to search for gravitational wave emission from a large selection of stars including some of those with the largest spin-down luminosities. Due to the sensitivity reduction caused at low frequency by seismic noise at the detectors, it is not worthwhile to search for pulsars with rotational frequencies, frot, smaller than about 10 Hz, which corresponds to gravitational wave mass quadrupole emission at frequencies, fgw = 2 frot, smaller than 20 Hz. The exact value of this gravitational wave low-frequency cut-off is rather arbitrary, our choice of taking 20 Hz is motivated by the presence of several noise lines and bumps in Virgo data at lower frequencies. In general, young pulsars, with large spin-down luminosities are searched for at lower frequencies where the Virgo detector has better sensitivity, whereas the search for millisecond pulsars (MSPs) is conducted at higher frequencies where the LIGO detectors are more sensitive. The selection of pulsars will be discussed more fully in Section 2.

1.1. The Signal

The expected quadrupolar gravitational wave signal from a triaxial neutron star 144 steadily spinning about one of its principal axes of inertia is at twice the rotation frequency, with a strain of

Equation (1)

in the detector, where

Equation (2)

is the dimensionless gravitational wave strain amplitude. h0 is dependent on Izz , the fiducial equatorial ellipticity, defined as ε = (Ixx Iyy )/Izz in terms of principal moments of inertia, the rotational frequency, frot, and the distance to the source d. The signal amplitudes in the two polarizations ("+" and "×") depend on the inclination of the star's rotation axis to the line-of-sight, ι, while the detector antenna pattern responses for the two polarization states, F+(t, ψ) and F×(t, ψ), depend on the gravitational wave polarization angle, ψ, as well as the detector location, orientation and source sky position. The gravitational wave phase evolution, ϕ(t), depends on both the intrinsic rotational frequency and frequency derivatives of the pulsar and on Doppler and propagation effects. These extrinsic effects include relativistic modulations caused by the Earth's orbital and rotational motion, the presence of massive bodies in the solar system close to the line-of-sight to the pulsar, the proper motion of the pulsar, and (in the case of a binary system) pulsar orbital motions. We will assume that ϕ(t) is phase-locked to the electromagnetic pulse phase evolution, but with double the value and with an initial phase offset, ϕ0. Given this phase evolution, the four unknown search parameters are simply h0, cos ι, ϕ0 and ψ. The gravitational wave amplitude is related to the star's l = m = 2 mass quadrupole moment via (see, e.g., Owen 2005)

Equation (3)

where Q22 is the slightly non-standard definition of quadrupole moment used in Ushomirsky et al. (2000) and many subsequent papers. This value can be constrained independently of any assumptions about the star's equation of state and moment of inertia.

If we allocate all the spin-down luminosity, $\dot{E}$, to gravitational wave luminosity, $\dot{E}_{\rm gw}$, where

Equation (4)

then we have the canonical "spin-down limit" on gravitational wave strain 145

Equation (5)

where I38 is the star's moment of inertia in the units of 1038 kg m2, and dkpc is the distance to the pulsar in kiloparsecs. The spin-down limit on the signal amplitude corresponds (via Equation (2)) to an upper limit on the star's fiducial ellipticity 146

Equation (6)

Johnson-McDaniel (2013) shows how to relate this to the physical ellipticity of the star's surface for a given equation of state.

A gravitational wave strain upper limit that is below the spin-down limit is an important milestone, as such a measurement is probing uncharted regions of the parameter space. Likewise it directly constrains the fraction of spin-down power that could be due to the emission of gravitational waves, which gives insight into the overall spin-down energy budget.

1.2. The Science Runs

In this paper we have used calibrated data from the Virgo second (Aasi et al. 2012) and fourth science runs (VSR2 and VSR4) and the LIGO sixth science run (S6). Virgo's third science run (VSR3) was relatively insensitive in comparison with VSR4 and has not been included in this analysis. This was partially because Virgo introduced monolithic mirror suspensions prior to VSR4 which improved sensitivity in the low-frequency range. During S6, the two LIGO 4 km detectors at Hanford, Washington (LHO/H1), and Livingston, Louisiana (LLO/L1), were running in an enhanced configuration (Adhikari et al. 2006) over that from the previous S5 run (Abbott et al. 2009). Table 1 shows dates of the runs, the duty factors and science data lengths for each detector that we analyzed.

Table 1. Science Runs

RunDatesDuty FactorData Length
(%)(days)
VSR22009 Jul 7 (20:55 UTC) – 2010 Jan 8 (22:00 UTC)80.4149
VSR42011 Jun 3 (10:27 UTC) – 2011 Sep 5 (13:26 UTC)81.076
S6 Hanford (H1)2009 Jul 7 (21:00 UTC) – 2010 Oct 21 (00:00 UTC)50.6238
S6 Livingston (L1)2009 Jul 7 (21:00 UTC) – 2010 Oct 21 (00:00 UTC)47.9225

Download table as:  ASCIITypeset image

The Virgo and LIGO data used in these analyses have been calibrated through different reconstruction procedures, but both based ultimately on the measured response to actuation controls on the positions of the mirrors that define the interferometers. For Virgo VSR2, the calibration uncertainty was about 5.5% in amplitude and ∼50 mrad (3°) in phase over most of the frequency range (Accadia et al. 2011). For VSR4, the uncertainty amounted to about 7.5% in amplitude and to (40 + 50f) mrad in phase, where f is the frequency in kilohertz, for frequencies up to 500 Hz (Mours & Rolland 2011). For LIGO, the S6 calibration uncertainties over the relevant frequency range (50–1500 Hz) were up to ∼19% in amplitude and ∼170 mrad (10°) in phase for L1, and up to ∼16% in amplitude and ∼120 mrad (7°) for H1 (Bartos et al. 2011). These phase errors are well within the range (i.e., less than 30°as applied in Abbott et al. 2007a) that would cause significant loss in signal power due to decoherence between the pulsar signal and the assumed phase evolution.

1.3. Methods

We used three semi-independent methods (very similar to those used in the Vela pulsar search in Abadie et al. 2011) to search for signals described in Section 1.1. Here, we briefly outline their operation, but for full descriptions we refer the reader to the references below. Two of the search methods work with time domain data that has been heterodyned to remove the signal's phase evolution and then heavily decimated. This leaves a complex data stream in which any signal would only be modulated by the detector's beam pattern. In the first method, this data stream is used to give Bayesian parameter estimates of the unknown signal parameters 147 (Dupuis & Woan 2005). The second method computes the maximum likelihood $\mathcal {F}$-statistic rather than a Bayesian posterior (or in case where ψ and ι are well constrained, the $\mathcal {G}$-statistic) (Jaranowski & Królak 2010). The third method (Astone et al. 2010) makes use of a Short Fourier Transform Database (SFDB) of each detector's data. After the extraction of a small frequency band around the signal's expected frequency, the Doppler effect, Einstein delay and spin-down are removed in the time domain and the data are down-sampled with a re-sampling technique. Two matched filters on the "+" and "×" signal Fourier components are then computed at the five frequencies at which the signal power is spread due to the signal amplitude and phase modulation; they are used to build a detection statistic and to estimate signal parameters in the case of detection. This 5-vector method has been extended over that used in Abadie et al. (2011) to allow for coherent analysis of data from multiple detectors (Astone et al. 2012). Each of these methods can incorporate prior information on the pulsar's inclination and polarization angle. From here on, we will refer to the first method as the Bayesian method, 148 the second as the $\mathcal {F}$/$\mathcal {G}$-statistic method and the third as the 5n-vector method, where n refers to the number of datasets coherently combined.

All three methods apply some data cleaning. The procedure used to obtain the heterodyned data removes extreme outliers by running two passes of a scheme that identifies points with absolute values greater than five times the standard deviation of the dataset. The $\mathcal {F}$/$\mathcal {G}$-statistic method performs further cleaning of this data through the Grubbs test (see Abadie et al. 2011). In the 5n-vector method, after an initial time-domain cleaning before the construction of the SFDB, a further cleaning step is applied on the final down-sampled time series in which the largest outliers belonging to the non-Gaussian tail of the data amplitude distribution are removed.

We have incorporated some limits from the previous LIGO S5 results (Abbott et al. 2010) as priors in the Bayesian analysis. However, the S6/VSR2,4 phase models were produced with updated pulsar ephemerides resulting in an unknown phase offset between them and the S5 results. We have, therefore, simply used the S5 marginalized posterior on h0 and cos ι, p(h0, cos ι), as our prior for the new results. In the case of glitching pulsars (see Section 2), we used the same approach and (incoherently) combined the separate coherent analyses produced between glitches. In the case of the $\mathcal {F}$/$\mathcal {G}$-statistic method, the results from different detectors or different inter-glitch periods are combined incoherently by adding the respective statistics. Also, for the 5n-vector method, results from different inter-glitch periods are incoherently combined by summing the corresponding statistics. Our reasons for not coherently combining the data over glitches are twofold. The first is that we do not know how a glitch would effect the relative phase offset between the electromagnetic pulses and the gravitational wave signal. The second reason is a practical consideration based on the timing solutions we have for our pulsars. For three of the four glitching pulsars (all except J0537−6910) in this analysis we have separate timing solutions for each inter-glitch period. These separate timing solutions, as provided by the pulsar timing software TEMPO(2), do not give an epoch defined at a fixed pulse phase (i.e., the epoch is not given as the time of the peak of a particular pulse), so there is some unknown phase offset between the separate solutions. However, if this phase offset were known (e.g., by going back to the original pulsar pulse time of arrival data) the gain in sensitivity would still be minimal: for the Vela pulsar and J1813−1246 the data from VSR4 (which was after the glitches in these pulsars) was much more sensitive, so completely dominates the result; for J1952+3252 the post-glitch data contains the latter part of S6, which was more sensitive and would again dominate the results. For J0537−6910 the epoch for each inter-glitch timing solution is defined (see Appendix), so the electromagnetic phase could be tracked over the glitches, but again our result is dominated by the longest and most sensitive inter-glitch period.

Even without a detection, all three methods can be used to produce upper limits on the gravitational wave amplitude from the pulsars. Here, we will quote 95% confidence upper limits on the amplitude. In the Bayesian method, an upper limit on the h0 posterior (after marginalization over the orientation parameters) is found by calculating the upper bound, from zero, on the integral over this posterior that encloses 95% of the probability. In the $\mathcal {F}$/$\mathcal {G}$-statistic method, a frequentist upper limit is calculated through Monte-Carlo simulations, which find the value of h0 for which 95% of trials exceed the maximum likelihood statistic. 149 The 5n-vector method computes an upper limit on the H0 posterior, given the actual value of the detection statistic, and the marginalization over the other parameters is implicitly done in the Monte Carlo simulation used to compute the likelihood. The amplitude, H0, is linked to the classical h0, given by Equation (2), by the relation $H_0=({h_0}/{2})\sqrt{1+6\cos {}^2\iota +\cos {}^4\iota }$ (see Equation (A5) in Abadie et al. 2011). H0 is the strain amplitude of a linearly polarized signal with polarization angle ψ = 0. In order to convert an upper limit on H0 to an upper limit on h0, we use the previous equation replacing the coefficient on the right hand side with its mean value over the distribution of cos ι used in the upper limit procedure. This is justified by the fact that the posterior distribution of H0 is not dependent on ι. The three methods have been tested with hardware and software simulated signal injections to check that they can recover the expected signal model (see, e.g., Abadie et al. 2011). In the Bayesian analysis these upper limits are really 95% credibility, or degrees-of-belief, values, whereas for the frequentist analysis these are 95% confidence values. These are both asking different questions and in general should not be expected to produce identical results. A brief discussion of this is given in the first search for a pulsar in LIGO data in Abbott et al. (2004), whilst a more technical discussion of the differences between the upper limits can be found in Röver et al. (2011).

2. PULSAR SELECTION

The sensitivity of the Virgo and LIGO detectors allows us to target pulsars with frot > 10 Hz. Currently the Australia Telescope National Facility (ATNF) pulsar catalog (Manchester et al. 2005) contains data for 368 pulsars (out of a total of 2264) consistent with this criterion. 150 The majority of these (∼90%) are recycled MSPs that have been spun-up to high rotation frequencies by accretion from a binary companion which may still be present (see, e.g., Lorimer 2008, for an overview of MSPs and binary pulsars). MSPs spin down slowly (with $\dot{f}_{\rm rot}$ between approximately −10−14 and −10−17 Hz s−1) and have characteristic ages 151 greater than a few times 108 yr, implying a comparatively weak surface polar magnetic field (108 G ≲ Bs ≲ 109 G, via the relation for an orthogonal rotator with radius 10 km and Izz = I38 of $B_s = 3.3\!\times \!10^{19}(|\dot{f}_{\rm rot}|/f_{\rm rot} ^3)^{1/2}$ G) compared to "normal" pulsars. About 10% are young pulsars with $\dot{f}_{\rm rot}$ between approximately −10−10 and −10−12 Hz s−1, characteristic ages of between ∼1000 and a few tens of thousands of years, and therefore with the large implied surface magnetic fields of "normal" pulsars, Bs ∼ 1012 G. They are situated towards the low-frequency end of our sensitivity range.

Young pulsars have large spin-downs and relatively low frequencies, so in general have the highest gravitational wave spin-down limits, see Equation (5). This makes them particularly important targets as the limits can be within reach of current detectors. Equations (3) and (6) show that to produce emission at around the spin-down limit the required mass quadrupole/ellipticity would have to be large, at a level consistent with only the most exotic neutron star equations of state (see the discussion in Section 4). Such strong emission is unlikely, but its detection would have profound implications. Young pulsars also often show rotational anomalies such as glitches and timing noise (see, e.g., Helfand et al. 1980). The underlying causes of such phenomena are still quite uncertain, and gravitational wave data would be a powerful constraint. For the MSPs, the spin-down limits are generally several orders of magnitude below those for the young pulsars. They are located, however, in a more sensitive frequency range.

2.1. Electromagnetic Pulsar Observations

For this analysis, we have obtained ephemerides using radio, X-ray and γ-ray observations. The radio telescope observations have come from a variety of sources: the 12.5-m telescope and Lovell telescope at Jodrell Bank in the UK, the 26-m telescope at Hartebeesthoek in South Africa, the 15-m eXperimental Development Model (XDM) telescope in South Africa, the Giant Metrewave Radio Telescope (GMRT) in India, the Robert C. Byrd Green Bank Radio Telescope in the US, the Parkes radio telescope in Australia, the Nançay Decimetric Radio Telescope in France and the Hobart radio telescope in Australia. High energy X-ray and γ-ray timings have been obtained from the Rossi X-ray Timing Explorer (RXTE) and the Fermi Large Area Telescope (LAT).

In total, for this analysis, we collected timing solutions for 179 pulsars. This selection includes 73 pulsars that have not been previously studied. However, for five of the pulsars targeted in the S3/S4 analysis (Abbott et al. 2007a) and another eleven of the pulsars targeted in the S5 analysis (Abbott et al. 2010), new coherent timing solutions were not available, so these stars 152 have not been included in this search.

2.1.1. High Interest Targets

As discussed in Abbott et al. (2008), 153 due to our ignorance of the correct neutron star equation of state there is a large uncertainty in the moments of inertia for our targets, from 1 to 3 × 1038 kg m2. Therefore, the canonical spin-down limit estimates could be increased by a factor of ∼1.7. Also, there are uncertainties in some pulsar distance measurements of up to a factor of two which could further increase or decrease the spin-down limit. We therefore identified all sources that were within a factor of four of the canonical spin-down limit as worthy of special attention. Seven of the pulsars for which we have obtained timing solutions beat, or approach to within a factor of four, this limit. The electromagnetic observation epochs for each pulsar (which include each inter-glitch epoch for pulsars that glitched during the analysis) are given in Table 2.

Table 2. Electromagnetic Observation Epochs for the High Interest Pulsars

MJD and Date 
J0534+2200 (Crab pulsar) 
54997 (2009 Jun 15) – 55814 (2011 Sep 10) 
J0537−6910 (N157B) 
54897 (2009 Mar 7) – 55041 (2009 Jul 29) 
55045 (2009 Aug 2) – 55182 (2009 Dec 17) 
55185 (2009 Dec 20) – 55263 (2010 Mar 8) 
55275 (2010 Mar 20) – 55445 (2010 Sep 6) 
55458 (2010 Sep 19) – 55503 (2010 Nov 3) 
J0835−4510 (Vela pulsar) 
54983 (2009 Jun 1) – 55286 (2010 Mar 31) 
55713 (2011 Jun 1) – 55827 (2011 Sep 23) 
J1813−1246 
54693 (2008 Aug 15) – 55094 (2009 Sep 20) 
55094 (2009 Sep 20) – 55828 (2011 Sep 24) 
J1833−1034 (G21.5−0.9) 
55041 (2009 Jul 29) – 55572 (2011 Jan 11) 
J1913+1011 
54867 (2009 Feb 5) – 55899 (2011 Dec 4) 
J1952+3252 (CTB 80) 
54589 (2008 May 3) – 55325 (2010 May 9) 
55331 (2010 May 15) – 55802 (2011 Aug 29) 

Download table as:  ASCIITypeset image

Further details of these observations are given below:

  • J0534+2200 (the Crab pulsar). We have used the Jodrell Bank Monthly Ephemeris (Lyne et al. 1993) to track the phase of the Crab pulsar over the period of our runs. This ephemeris has timing solutions using the DE200 solar system ephemeris and the Barycentric Dynamical Time (TDB) time coordinate system. During S6/VSR2,4 the pulsar did not show signs of any timing glitches.
  • J0537−6910 (N157B). Long-term X-ray timing has been performed with the RXTE (Middleditch et al. 2006). Recent data covering S6 shows four glitches over the span of our science runs and the ephemerides for each inter-glitch epoch are given in the Appendix. The timing solutions used the DE200 solar system ephemeris (see Marshall et al. 1998) and the TDB time coordinate system. Several more glitches have been observed since the end of our science runs, but we do not report on them here.
  • J0835−4510 (the Vela pulsar). Radio observations over the period of VSR2 were taken with the Hobart radio telescope in Tasmania and the Hartebeesthoek 26-m radio telescope in South Africa (Abadie et al. 2011). Radio timing over the VSR4 run was performed with the XDM telescope and the 26-m telescope at Hartebeesthoek. The timing solutions have used the DE405 solar system ephemeris and the Barycentric Coordinate Time (TCB) time coordinate system. Vela was observed to glitch on 2010 July 31 (Buchner 2010), between VSR2 and VSR4, but it has not glitched since then.
  • J1813−1246. This pulsar was discovered in a search of gamma-ray data from the Fermi LAT (Abdo et al. 2009), and using the unbinned maximum likelihood methods of Ray et al. (2011) timing measurements were made covering all our runs. It was observed to glitch once during this time on 2009 September 20. Pre-and-post glitch timing solutions have been produced using the DE405 solar system ephemeris and the TDB time coordinate system.
  • J1833−1034 (G21.5−0.9). The period from the start of S6/VSR2 until 2011 January is covered by observations made with the GMRT (Roy et al. 2012). During this period, one glitch was observed, with a best fit epoch of 2009 November 6 (MJD 55142±2). To remove its effect, an ephemeris fit was performed on timing data excluding 80 days after the glitch. The timing solution uses the DE405 solar system ephemeris and the TDB time coordinate system.
  • J1913+1011. This pulsar was observed at Jodrell Bank and showed no timing anomalies over the science runs. The timing solution uses the DE405 solar system ephemeris and the TDB time coordinate system.
  • J1952+3252 (CTB 80). This pulsar was observed over the whole of our science runs at Nançay and Jodrell Bank. It glitched on 2010 May 11 (MJD 55327), between the end of S6/VSR2 and the start of VSR4. Phase incoherent pre- and post-glitch timing solutions have been produced using the DE405 solar system ephemeris and the TCB time ephemeris. The solution include fits to the timing noise using the tempo2 Fitwaves method described in Hobbs et al. (2006).

For several of these pulsars potential constraints on their orientations (the inclination ι and polarization angle ψ 154 ) are available from observations of their pulsar wind nebulae (Ng & Romani 2004, 2008). These are listed in Table 3 where the uncertainties used are estimated from the systematic and statistical values given in Ng & Romani (2004, 2008), and the mean angle value is used if multiple fits are given (e.g., fits to the inner and outer tori of the Crab pulsar wind nebula). We briefly discussed how these constraints are used in the analyses in Section 1.3.

Table 3. Implied Orientations of Pulsars from their Pulsar Wind Nebulae Observations (Ng & Romani 2004, 2008)

Pulsarιψ
J0534+2200 (Crab pulsar)62fdg2 ± 1fdg935fdg2 ± 1fdg5
J0537−691092fdg8 ± 0fdg941fdg0 ± 2fdg2
J0835−4510 (Vela pulsar)63fdg6 ± 0fdg640fdg6 ± 0fdg1
J1833−103485fdg4 ± 0fdg345° ± 1° 
J1952+3252 a  ⋅⋅⋅−11fdg5 ± 8fdg6

Note. a The polarization angle is not taken from a fit to the pulsar wind nebula, but instead is the average of the angle calculated from proper motion measurements and Hα observations of a bow shock (Ng & Romani 2004).

Download table as:  ASCIITypeset image

For J0534+2200 and J0537−6910, the Bayesian method also makes use of results from the LIGO S5 run (Abbott et al. 2010) as a prior on the h0 and cos ι parameters. During S5, both of these pulsars glitched, and the data for each inter-glitch period was analyzed independently. Results were also produced assuming that the data could be analyzed coherently over the glitches. To avoid the assumptions about coherence over the glitches, we have used the independent inter-glitch results that gave the lowest h0 as the prior for the current analysis (see Table 3 of Abbott et al. 2010).

3. RESULTS

None of the searches yielded evidence of a gravitational wave signal, and upper limits have been placed on signal strengths. These limits are subject to the uncertainties in the amplitude calibration, as discussed in Section 1.2. For the joint results, which combine data from multiple detectors, the sensitivity is often dominated by the most sensitive instrument. Therefore, we expect the amplitude uncertainty due to calibration uncertainties to also be dominated by the most sensitive instrument. So, following the calibration error given in Section 1.2, below ∼50 Hz we have an amplitude uncertainty of ∼6%, and above that we have uncertainty of ∼20%. The phase uncertainties are small enough to have a negligible contribution to the possible amplitude uncertainty.

3.1. Data Selection

As discussed in Section 2.1, for a few pulsars the electromagnetic observations did not always span the S6/VSR2,4 runs completely, and some pulsars glitched during the runs. As a result, we deal with these instances separately. In most cases, we can use all the data coherently, but, in other cases, sections of data must be combined incoherently. The relative sensitivities of the detectors at the pulsar frequencies also dictate whether we have used Virgo-only, LIGO-only or Virgo and LIGO data (see Figure 1). For J0537−6910, only the LIGO data has been used because of its better sensitivity at the corresponding frequency, and results from each inter-glitch period have been combined incoherently. For J0835−4510 (the Vela pulsar), a glitch occurred just prior to VSR4 and we had no phase-connected timing solution between VSR2 (Abadie et al. 2011) and VSR4 epochs. The VSR2 results and VSR4 data have therefore been incoherently combined. For J1813−1246, the results from the pre- and post-glitch periods using all data from S6 and VSR2,4 have been combined incoherently. For J1833−1034, only VSR2 data up to the time of the observed glitch has been used.

Figure 1.

Figure 1.  $h_0^{95\%}$ upper limits (given by ⋆) for 195 pulsars from the LIGO and Virgo S3/S4, S5, S6, VSR2, and VSR4 runs. The curves give estimated relative strain sensitivities of these runs and potential future science runs. The sensitivities are based on the harmonic mean of the observation time (T) weighted one-sided power spectral densities Sn from all detectors operating during the given run, and are given by $10.8\sqrt{S_n/T}$, where the scale factor of 10.8 is given in Dupuis & Woan (2005). The AdV/aLIGO curve assumes a joint analysis of two equally sensitive advanced LIGO detectors and the advanced Virgo detector operating at their full design sensitivities with one year of coherent integration (the sensitivity curves are those given in Aasi et al. 2013c). The ▾ give the spin-down limits for all (non-Globular Cluster) pulsars, based on values taken from the ATNF catalog and assuming the canonical moment of inertia. The ⋆ show the observational upper limits from Tables 56, and 7, with the seven high interest pulsars represented by the larger, lighter colored stars. Results for pulsars using the previous S3/S4 and S5 data are given by the small lighter colored stars.

Standard image High-resolution image

The parameters and results (from the three different analyses discussed in Section 1.3) for the seven pulsars highlighted in Section 2.1 are given in Tables 4 and 5, respectively. Table 5 gives the 95% upper limit on the gravitational wave amplitude, $h_0^{95\%}$, the equivalent limits (via Equation (3)) on the stars fiducial ellipticity, ε, and mass quadrupole moment, Q22, the ratio of the limit to the spin-down limit, $h_0^{95\%}/h_0^{\rm sd}$, and the limit on the gravitational wave luminosity compared to the total spin-down luminosity. This final value is given in the form of the percentage of the spin-down luminosity required to produce a gravitational wave at the amplitude limit (it can be seen from Equation (4) that this is just the square of the ratio $h_0^{95\%}/h_0^{\rm sd}$). For those pulsars with constrained orientations (see Table 3) the results with and without the constraints are also given in Table 5. Despite the very tight constraints given in Table 3 it should be noted that these results would only show minor differences if the angle (ι and ψ) errors were expanded to several times the given values. This is because the recovered posterior probability distributions on these parameters are slowly and smoothly varying over their parameter ranges. A brief discussion of the differences between the upper limits from the different methods is given in Section 1.3.

Table 4. The Properties of the Pulsars of High Interest

Pulsarαδ frot fgw $\dot{f}_{\rm rot}$ d $\dot{E}$ a hsd a
(Hz)(Hz)(Hz s−1)(kpc)(W)
J0534+2200 (Crab)05h34m31fs9722°00'52farcs0729.7259.44−3.7 × 10−10 2.0 b 4.6 × 1031 1.4 × 10−24
J0537−6910 (N157B)05h37m47fs36−69°10'20farcs4061.97123.94−2.0 × 10−10 50.0 c 4.9 × 1031 3.0 × 10−26
J0835−4510 (Vela)08h35m20fs61−45°10'34farcs8811.1922.39−1.6 × 10−11 0.29 d 6.9 × 1029 3.3 × 10−24
J1813−124618h13m23fs74−12°46'00farcs8620.8041.60−7.6 × 10−12 1.9 e 6.2 × 1029 2.6 × 10−25
J1833−1034 (G21.5−0.9)18h33m33fs61−10°34'16farcs6116.1632.33−5.3 × 10−11 4.8 f 3.4 × 1030 3.0 × 10−25
J1913+101119h13m20fs3410°11'23farcs1127.8555.70−2.6 × 10−12 4.5 g 2.8 × 1029 2.3 × 10−25
J1952+3252 (CTB 80)19h52m58fs1132°52'41farcs2425.3050.59−3.7 × 10−12 3.0 g 3.7 × 1029 1.0 × 10−25

Notes. aThe spin-down luminosity, $\dot{E}$, and spin-down gravitational wave amplitude limit, hsd, both assume a canonical moment of inertia of Izz = 1038 kg m2. bSee Appendix of Kaplan et al. (2008). cPietrzyński et al. (2013). dDodson et al. (2003). eThis distance is the average of the two estimates from Wang (2011), which allow a distance between ∼0.9–3.5 kpc. fTian & Leahy (2008). gThe distance is taken from the ATNF pulsar catalog (Manchester et al. 2005).

Download table as:  ASCIITypeset image

Table 5. Upper Limits for the High Interest Pulsars a

Analysis $h_0^{95\%}$ ε Q22 (kg m2) $h_0^{95\%}/h_0^{\rm sd}$ $\dot{E}_{\rm gw}/\dot{E}$ %
J0534+2200 (Crab)
Bayesian1.6 (1.4) × 10−25 8.6 (7.5) × 10−5 6.6 (5.8) × 1033 0.11 (0.10)1.2 (1.0)
$\mathcal {F}$/$\mathcal {G}$-statistic2.3 (1.8) × 10−25 12.3 (9.6) × 10−5 11.6 (7.4) × 1033 0.16 (0.13)2.6 (1.7)
5n-vector1.8 (1.6) × 10−25 9.7 (8.6) × 10−5 7.4 (6.6) × 1033 0.12 (0.11)1.4 (1.2)
J0537−6910
Bayesian3.8 (4.4) × 10−26 1.2 (1.4) × 10−4 0.9 (1.0) × 1034 1.4 (1.7)200 (290)
$\mathcal {F}$/$\mathcal {G}$-statistic1.1 (1.0) × 10−25 3.4 (3.1) × 10−4 2.6 (2.4) × 1034 4.1 (3.9)1700 (1500)
5n-vector4.5 (6.7) × 10−26 1.4 (2.1) × 10−4 1.1 (1.6) × 1034 1.6 (2.4)260 (580)
J0835−4510 (Vela)
Bayesian1.1 (1.0) × 10−24 6.0 (5.5) × 10−4 4.7 (4.2) × 1034 0.33 (0.30)11 (9.0)
$\mathcal {F}$/$\mathcal {G}$-statistic4.2 (9.0) × 10−25 2.3 (4.9) × 10−4 1.8 (3.8) × 1034 0.13 (0.27)1.7 (7.3)
5n-vector1.1 (1.1) × 10−24 6.0 (6.0) × 10−4 4.7 (4.7) × 1034 0.33 (0.33)11 (11)
J1813−1246
Bayesian3.4 × 10−25 3.5 × 10−4 2.7 × 1034 1.3170
$\mathcal {F}$/$\mathcal {G}$-statistic7.1 × 10−25 7.4 × 10−4 5.7 × 1034 2.7730
5n-vector4.8 × 10−25 4.9 × 10−4 3.8 × 1034 1.8320
J1833−1034
Bayesian1.3 (1.4) × 10−24 5.7 (6.1) × 10−3 4.4 (4.7) × 1035 4.3 (4.6)1800 (2100)
$\mathcal {F}$/$\mathcal {G}$-statistic1.2 (1.2) × 10−24 5.2 (5.2) × 10−3 4.0 (4.0) × 1035 3.9 (3.9)1500 (1500)
5n-vector1.4 (2.0) × 10−24 6.1 (8.7) × 10−3 4.7 (6.7) × 1035 4.6 (6.6)2100 (4400)
J1913+1011
Bayesian1.6 × 10−25 2.2 × 10−4 1.7 × 1034 2.9840
$\mathcal {F}$/$\mathcal {G}$-statistic2.9 × 10−25 4.1 × 10−4 3.1 × 1034 5.32800
5n-vector2.5 × 10−25 3.4 × 10−4 2.7 × 1034 4.52000
J1952+3252
Bayesian2.7 (2.5) × 10−25 3.0 (2.8) × 10−4 2.3 (2.1) × 1034 2.6 (2.5)680 (630)
$\mathcal {F}$/$\mathcal {G}$-statistic6.0 × 10−25 6.7 × 10−4 5.1 × 1034 5.83400
5n-vector3.1 (3.2) × 10−25 3.4 (3.5) × 10−4 2.6 (2.7) × 1034 3.0 (3.1)900 (960)

Notes. Limits with constrained orientations are given in parentheses. aDetector calibration errors mean that for pulsars with fgw below and above 50 Hz (see Table 4) there are ∼6% and ∼20% uncertainties respectively on these limits.

Download table as:  ASCIITypeset image

One of the new targets, J1824−2452I (which is an interesting pulsar that is seen to switch between being accretion and rotation powered; Papitto et al. 2013), had a coherent timing solution that covered 2006, so S5 data from the LIGO detectors has been reanalyzed for this result. For all other pulsars, we have used only the VSR2 and VSR4 data if fgw < 40 Hz, and have coherently combined VSR2, VSR4 and S6 data from H1 and L1 for pulsars with fgw > 40 Hz. All the available science mode data (i.e., when the detectors were operating in a stable state) has been used, with details given in Table 1.

For the 19 pulsars with fgw < 40 Hz the results can be found in Table 6. Because of their low frequencies, none of these pulsars had been targeted before.

Table 6. Limits on the Gravitational Wave Amplitude for Known Pulsars with fgw < 40 Hz Using VSR2,4 Data

Pulsar frot fgw $\dot{f}_{\rm rot}$ d $h_0^{\rm sd}$ $h_0^{95\%}$ ε Q22 $h_0^{95\%}/h_0^{\rm sd}$
(Hz)(Hz)(Hz s−1)(kpc)(kg m2)
J0106+485512.0324.05−6.2 × 10−14 7.37.9 × 10−27 7.4 × 10−25 8.9 × 10−3 6.9 × 1035 94
J0609+213017.9535.91−7.6 × 10−17 1.89.1 × 10−28 4.4 × 10−25 5.9 × 10−4 4.6 × 1034 490
J1528−314616.4432.88−6.7 × 10−17 1.01.6 × 10−27 6.0 × 10−25 5.2 × 10−4 4.0 × 1034 360
J1718−382513.3926.78−2.4 × 10−12 4.28.0 × 10−26 7.6 × 10−25 4.2 × 10−3 3.3 × 1035 9.5
J1747−295810.1220.24−6.3 × 10−12 2.52.6 × 10−25 1.8 × 10−24 1.0 × 10−2 7.9 × 1035 7.0
J1748−2446J12.4524.89−2.0 × 10−16 5.55.8 × 10−28 a 1.3 × 10−24 1.1 × 10−2 8.3 × 1035 2200
J1753−191415.8831.77−4.9 × 10−16 2.81.6 × 10−27 5.5 × 10−25 1.4 × 10−3 1.1 × 1035 340
J1753−224010.5121.02−6.9 × 10−17 3.56.0 × 10−28 2.2 × 10−24 1.6 × 10−2 1.3 × 1036 3700
J1809−191712.0824.17−3.7 × 10−12 3.71.2 × 10−25 1.2 × 10−24 7.2 × 10−3 5.6 × 1035 9.9
J1828−110113.8827.76−2.9 × 10−12 7.35.0 × 10−26 1.0 × 10−24 9.2 × 10−3 7.1 × 1035 20
J1831−095214.8729.73−1.8 × 10−12 4.36.5 × 10−26 5.7 × 10−25 2.6 × 10−3 2.0 × 1035 8.7
J1833−082711.7223.45−1.3 × 10−12 4.55.9 × 10−26 b 1.6 × 10−24 1.2 × 10−2 9.6 × 1035 27
J1856+024512.3624.72−9.5 × 10−12 10.36.9 × 10−26 10.0 × 10−25 1.6 × 10−2 1.2 × 1036 15
J1904+041214.0728.13−2.8 × 10−17 4.02.8 × 10−28 9.2 × 10−25 4.4 × 10−3 3.4 × 1035 3300
J1915+160616.9433.88−2.5 × 10−15 7.11.4 × 10−27 b 3.9 × 10−25 2.3 × 10−3 1.7 × 1035 280
J1928+174614.5529.10−2.8 × 10−12 8.14.3 × 10−26 6.8 × 10−25 6.2 × 10−3 4.8 × 1035 16
J1954+283610.7921.57−2.5 × 10−12 1.62.4 × 10−25 1.2 × 10−24 4.0 × 10−3 3.1 × 1035 5.1
J2043+274010.4020.80−1.3 × 10−13 1.18.1 × 10−26 1.7 × 10−24 4.1 × 10−3 3.2 × 1035 21
J2235+150616.7333.46−2.9 × 10−17 1.19.2 × 10−28 b 4.2 × 10−25 4.0 × 10−4 3.1 × 1034 450

Notes. Detector calibration errors mean that for these pulsars there are ∼6% uncertainties on these limits. aThe pulsar's spin-down is calculated using a characteristic spin-down age of 109 yr. bThe pulsar's spin-down is corrected for proper motion effects.

Download table as:  ASCIITypeset image

Results for pulsars with fgw > 40 Hz using S6 and VSR2,4 are shown in Table 7. Distances to pulsars in Terzan 5 (with designations J1748−2446) are assumed to be 5.5 kpc (Ortolani et al. 2007) rather than the value of 8.7 kpc given in the ATNF catalog, and distances to the pulsars in M28 (with designations J1824−2452) are assumed to be 5.5 kpc (Harris 1996) rather than the distance of 4.9 kpc given in Abbott et al. (2010). Unless otherwise specified in the table for all other pulsars we use the distance values given by the DIST value in the ATNF catalog (Manchester et al. 2005), which generally are dispersion measure calculations from the electron density distribution model of Taylor & Cordes (1993). For the 16 pulsars where new timing solutions were not available during the most recent runs (see Section 2.1), we include the results from the LIGO S3/S4 (Abbott et al. 2007a) and S5 analysis (Abbott et al. 2010).

Table 7. Limits on the Gravitational Wave Amplitude for Known Pulsars with fgw > 40 Hz

Pulsar frot fgw $\dot{f}_{\rm rot}$ d $h_0^{\rm sd}$ Prior $h_0^{95\%}$ S6/VSR2,4 $h_0^{95\%}$ ε Q22 $h_0^{95\%}/h_0^{\rm sd}$
(Hz)(Hz)(Hz s−1)(kpc)(kg m2)
J0023+0923327.85655.69−1.3 × 10−15 0.91.7 × 10−27  ⋅⋅⋅6.9 × 10−26 1.5 × 10−7 1.1 × 1031 41
J0024−7204C173.71347.421.6 × 10−15 4.06.1 × 10−28 a 5.8 × 10−25 1.7 × 10−25 5.5 × 10−6 4.2 × 1032 290
J0024−7204D186.65373.301.6 × 10−16 4.01.8 × 10−28 a 4.5 × 10−26 2.9 × 10−26 7.8 × 10−7 6.1 × 1031 160
J0024−7204E282.78565.56−7.8 × 10−15 4.01.1 × 10−27 a 9.9 × 10−26 9.7 × 10−26 1.1 × 10−6 8.9 × 1031 92
J0024−7204F381.16762.32−9.3 × 10−15 4.09.9 × 10−28 a 8.8 × 10−26 6.7 × 10−26 4.4 × 10−7 3.4 × 1031 67
J0024−7204G247.50495.002.6 × 10−15 4.06.6 × 10−28 a 9.9 × 10−26 8.2 × 10−26 1.3 × 10−6 9.8 × 1031 120
J0024−7204H311.49622.993.0 × 10−16 4.02.0 × 10−28 a 6.5 × 10−26 5.2 × 10−26 5.1 × 10−7 3.9 × 1031 260
J0024−7204I286.94573.893.9 × 10−15 4.07.4 × 10−28 a 5.2 × 10−26 5.8 × 10−26 6.7 × 10−7 5.2 × 1031 79
J0024−7204J476.05952.092.4 × 10−15 4.04.5 × 10−28 a 1.0 × 10−25 6.6 × 10−26 2.8 × 10−7 2.1 × 1031 150
J0024−7204L230.09460.18−3.6 × 10−15 4.08.0 × 10−28 b 5.8 × 10−26 4.2 × 10−26 7.4 × 10−7 5.7 × 1031 52
J0024−7204M271.99543.97−4.3 × 10−15 4.08.0 × 10−28 b 6.2 × 10−26 7.0 × 10−26 9.0 × 10−7 7.0 × 1031 88
J0024−7204N327.44654.892.5 × 10−15 4.05.6 × 10−28 a 8.4 × 10−26 5.1 × 10−26 4.5 × 10−7 3.5 × 1031 91
J0024−7204O378.31756.62−4.2 × 10−15 4.06.7 × 10−28 a 9.3 × 10−26 5.8 × 10−26 3.8 × 10−7 3.0 × 1031 86
J0024−7204Q247.94495.89−3.9 × 10−15 4.08.0 × 10−28 b 5.8 × 10−26 5.1 × 10−26 7.9 × 10−7 6.1 × 1031 64
J0024−7204R287.32574.64−4.6 × 10−15 4.08.0 × 10−28 b 5.6 × 10−26 3.6 × 10−26 4.1 × 10−7 3.2 × 1031 45
J0024−7204S353.31706.61−5.6 × 10−15 4.08.0 × 10−28 b 6.9 × 10−26 6.1 × 10−26 4.6 × 10−7 3.6 × 1031 76
J0024−7204T131.78263.56−2.1 × 10−15 4.08.0 × 10−28 b 3.3 × 10−26 3.8 × 10−26 2.0 × 10−6 1.6 × 1032 47
J0024−7204U230.26460.53−4.9 × 10−15 4.09.3 × 10−28 a 5.7 × 10−26 4.1 × 10−26 7.2 × 10−7 5.6 × 1031 43
J0024−7204Y455.24910.47−7.2 × 10−15 4.08.0 × 10−28 b 9.4 × 10−26 7.2 × 10−26 3.3 × 10−7 2.6 × 1031 90
J0030+0451205.53411.06−4.3 × 10−16 0.34.1 × 10−27 a  ⋅⋅⋅7.2 × 10−26 1.1 × 10−7 8.7 × 1030 17
J0034−0534532.711065.43−1.5 × 10−16 1.04.4 × 10−28 a  ⋅⋅⋅1.8 × 10−25 1.5 × 10−7 1.1 × 1031 410
J0218+4232430.46860.92−1.4 × 10−14 5.87.9 × 10−28 a 1.5 × 10−25 8.7 × 10−26 6.5 × 10−7 5.0 × 1031 110
J0307+7443316.85633.70−1.7 × 10−15 0.6 c 3.3 × 10−27  ⋅⋅⋅8.6 × 10−26 1.1 × 10−7 8.9 × 1030 26
J0340+41303.09606.18−6.5 × 10−16 2.74.4 × 10−28  ⋅⋅⋅5.6 × 10−26 3.8 × 10−7 3.0 × 1031 130
J0407+160738.9177.82−1.2 × 10−16 4.13.5 × 10−28 6.2 × 10−26 5.1 × 10−26 3.2 × 10−5 2.5 × 1033 140
J0437−4715173.69347.38−4.1 × 10−16 0.27.8 × 10−27 a 5.7 × 10−25 1.2 × 10−25 1.6 × 10−7 1.2 × 1031 16
J0605+3757 d 366.58733.15−6.5 × 10−16 0.7 c 1.6 × 10−27  ⋅⋅⋅9.3 × 10−26 1.1 × 10−7 8.6 × 1030 59
J0610−2100258.98517.96−2.3 × 10−16 5.61.3 × 10−28 a  ⋅⋅⋅7.9 × 10−26 1.6 × 10−6 1.2 × 1032 590
J0613−0200326.60653.20−9.4 × 10−16 0.91.5 × 10−27 a 1.1 × 10−25 5.6 × 10−26 1.1 × 10−7 8.6 × 1030 37
J0614−3329317.59635.19−1.7 × 10−15 3.06.4 × 10−28  ⋅⋅⋅8.5 × 10−26 5.9 × 10−7 4.5 × 1031 130
J0621+100234.6669.31−5.5 × 10−17 1.95.4 × 10−28 a 1.5 × 10−25 9.6 × 10−26 3.6 × 10−5 2.7 × 1033 180
J0711−6830182.12364.23−2.9 × 10−16 1.09.8 × 10−28 a 5.0 × 10−26 3.4 × 10−26 2.5 × 10−7 2.0 × 1031 35
J0737−3039A44.0588.11−3.4 × 10−15 1.16.5 × 10−27 a 7.8 × 10−26 6.0 × 10−26 8.0 × 10−6 6.2 × 1032 9.2
J0751+1807287.46574.92−6.3 × 10−16 0.43.0 × 10−27 a 1.6 × 10−25 1.1 × 10−25 1.2 × 10−7 9.6 × 1030 36
J0900−314490.01180.02−4.0 × 10−16 0.82.1 × 10−27  ⋅⋅⋅1.8 × 10−25 4.3 × 10−6 3.3 × 1032 88
J1012+5307190.27380.54−4.1 × 10−16 0.71.7 × 10−27 a 6.9 × 10−26 4.6 × 10−26 2.1 × 10−7 1.6 × 1031 27
J1017−7156427.62855.24−3.1 × 10−16 8.18.5 × 10−29  ⋅⋅⋅1.0 × 10−25 1.1 × 10−6 8.3 × 1031 1200
J1022+100160.78121.56−1.6 × 10−16 0.52.5 × 10−27 4.5 × 10−26 4.8 × 10−26 1.6 × 10−6 1.2 × 1032 19
J1024−0719193.72387.431.3 × 10−16 0.51.3 × 10−27 a 5.0 × 10−26 4.6 × 10−26 1.4 × 10−7 1.1 × 1031 35
J1038+003234.6669.32−7.8 × 10−17 2.45.1 × 10−28  ⋅⋅⋅1.2 × 10−25 5.5 × 10−5 4.3 × 1033 230
J1045−4509133.79267.59−3.1 × 10−16 0.25.3 × 10−27 a 4.3 × 10−26 3.0 × 10−26 9.3 × 10−8 7.2 × 1030 5.7
J1231−1411271.45542.911.6 × 10−15 0.54.3 × 10−27 a  ⋅⋅⋅7.5 × 10−26 1.1 × 10−7 8.4 × 1030 17
J1300+1240160.81321.62−7.9 × 10−16 0.63.0 × 10−27 a  ⋅⋅⋅4.9 × 10−26 2.7 × 10−7 2.1 × 1031 16
J1301+0833542.381084.76−3.1 × 10−15 0.92.1 × 10−27  ⋅⋅⋅1.1 × 10−25 7.9 × 10−8 6.1 × 1030 51
J1435−6100 e 106.98213.95−2.8 × 10−16 3.34.0 × 10−28 2.7 × 10−25  ⋅⋅⋅1.8 × 10−5 1.4 × 1033 670
J1453+1902172.64345.29−3.2 × 10−16 0.91.2 × 10−27 a  ⋅⋅⋅1.4 × 10−25 1.0 × 10−6 8.0 × 1031 120
J1455−3330125.20250.40−2.5 × 10−16 0.71.5 × 10−27 a 5.1 × 10−26 3.6 × 10−26 4.1 × 10−7 3.1 × 1031 24
J1518+0204A180.06360.13−2.9 × 10−15 8.04.0 × 10−28 b  ⋅⋅⋅7.8 × 10−26 4.6 × 10−6 3.5 × 1032 190
J1518+490424.4348.86−1.3 × 10−17 0.78.5 × 10−28 a  ⋅⋅⋅2.9 × 10−25 8.1 × 10−5 6.2 × 1033 340
J1537+115526.3852.76−1.6 × 10−15 1.06.3 × 10−27 a  ⋅⋅⋅2.4 × 10−25 8.4 × 10−5 6.5 × 1033 39
J1600−3053277.94555.88−6.5 × 10−16 2.45.1 × 10−28 a 5.6 × 10−26 6.7 × 10−26 4.9 × 10−7 3.8 × 1031 130
J1603−720267.38134.75−5.4 × 10−17 1.64.4 × 10−28 a 2.3 × 10−26 2.3 × 10−26 1.9 × 10−6 1.5 × 1032 51
J1614−2230317.38634.763.9 × 10−16 1.85.0 × 10−28 a  ⋅⋅⋅6.4 × 10−26 2.7 × 10−7 2.1 × 1031 130
J1623−263190.29180.57−5.1 × 10−15 2.22.7 × 10−27 a 5.7 × 10−26 5.1 × 10−26 3.3 × 10−6 2.5 × 1032 19
J1629−6902 e 166.65333.30−2.8 × 10−16 1.47.7 × 10−28 3.2 × 10−25  ⋅⋅⋅3.8 × 10−6 2.9 × 1032 420
J1630+3734 d 301.38602.75−8.9 × 10−16 0.9 c 1.5 × 10−27  ⋅⋅⋅9.0 × 10−26 2.2 × 10−7 1.7 × 1031 60
J1640+2224316.12632.25−1.6 × 10−16 1.24.9 × 10−28 a 6.7 × 10−26 5.1 × 10−26 1.4 × 10−7 1.1 × 1031 110
J1641+3627A96.36192.72−1.5 × 10−15 7.54.3 × 10−28 b  ⋅⋅⋅5.0 × 10−26 9.5 × 10−6 7.3 × 1032 120
J1643−1224216.37432.75−8.5 × 10−16 0.43.8 × 10−27 a 4.4 × 10−26 3.6 × 10−26 7.6 × 10−8 5.9 × 1030 9.4
J1701−3006A190.78381.57−3.0 × 10−15 6.94.7 × 10−28 b 5.8 × 10−26 3.6 × 10−26 1.6 × 10−6 1.3 × 1032 78
J1701−3006B f 278.25556.502.7 × 10−14 6.94.7 × 10−28 b 7.6 × 10−26  ⋅⋅⋅1.6 × 10−6 1.2 × 1032 160
J1701−3006C f 131.36262.721.1 × 10−15 6.94.7 × 10−28 b 3.5 × 10−26  ⋅⋅⋅3.3 × 10−6 2.6 × 1032 76
J1709+2313215.93431.85−6.9 × 10−17 1.82.5 × 10−28 a  ⋅⋅⋅9.3 × 10−26 8.6 × 10−7 6.7 × 1031 370
J1713+0747218.81437.62−3.9 × 10−16 1.11.0 × 10−27 a 4.5 × 10−26 3.5 × 10−26 1.8 × 10−7 1.4 × 1031 34
J1719−1438172.71345.41−2.3 × 10−16 1.65.8 × 10−28  ⋅⋅⋅1.6 × 10−25 2.0 × 10−6 1.6 × 1032 270
J1721−2457285.99571.98−2.4 × 10−16 1.64.7 × 10−28 a  ⋅⋅⋅5.5 × 10−26 2.5 × 10−7 1.9 × 1031 120
J1730−2304123.11246.22−3.1 × 10−16 0.52.5 × 10−27 5.8 × 10−26 4.5 × 10−26 3.6 × 10−7 2.8 × 1031 18
J1731−1847426.52853.04−4.6 × 10−15 4.06.6 × 10−28  ⋅⋅⋅1.3 × 10−25 6.6 × 10−7 5.1 × 1031 190
J1732−5049188.23376.47−5.0 × 10−16 1.87.3 × 10−28 5.3 × 10−26 4.6 × 10−26 5.6 × 10−7 4.3 × 1031 63
J1738+0333170.94341.87−6.7 × 10−16 2.08.1 × 10−28 a  ⋅⋅⋅1.1 × 10−25 1.8 × 10−6 1.4 × 1032 140
J1741+1351266.87533.74−2.2 × 10−15 1.41.6 × 10−27  ⋅⋅⋅1.1 × 10−25 5.1 × 10−7 3.9 × 1031 67
J1744−1134245.43490.85−4.3 × 10−16 0.42.5 × 10−27 a 1.1 × 10−25 6.3 × 10−26 1.0 × 10−7 8.0 × 1030 25
J1745−095251.61103.22−7.6 × 10−17 2.44.1 × 10−28 a  ⋅⋅⋅6.0 × 10−26 1.3 × 10−5 9.8 × 1032 150
J1748−2446A86.48172.96−1.4 × 10−15 5.55.8 × 10−28 b 3.9 × 10−26 3.2 × 10−26 5.6 × 10−6 4.3 × 1032 55
J1748−2446C118.54237.08−1.9 × 10−15 5.55.8 × 10−28 b 5.0 × 10−26 3.8 × 10−26 3.5 × 10−6 2.7 × 1032 65
J1748−2446D212.14424.27−3.4 × 10−15 5.55.8 × 10−28 b 6.8 × 10−26 5.3 × 10−26 1.5 × 10−6 1.2 × 1032 91
J1748−2446E455.00910.00−7.2 × 10−15 5.55.8 × 10−28 b 9.0 × 10−26 7.3 × 10−26 4.6 × 10−7 3.6 × 1031 130
J1748−2446F180.50361.00−2.9 × 10−15 5.55.8 × 10−28 b 8.3 × 10−26 7.4 × 10−26 3.0 × 10−6 2.3 × 1032 130
J1748−2446G46.1492.29−7.3 × 10−16 5.55.8 × 10−28 b 5.9 × 10−26 4.5 × 10−26 2.8 × 10−5 2.1 × 1033 78
J1748−2446H203.01406.02−3.2 × 10−15 5.55.8 × 10−28 b 7.8 × 10−26 5.8 × 10−26 1.8 × 10−6 1.4 × 1032 99
J1748−2446I104.49208.98−1.7 × 10−15 5.55.8 × 10−28 b 3.6 × 10−26 3.7 × 10−26 4.5 × 10−6 3.4 × 1032 64
J1748−2446K336.74673.48−5.3 × 10−15 5.55.8 × 10−28 b 6.8 × 10−26 5.9 × 10−26 6.8 × 10−7 5.3 × 1031 100
J1748−2446L445.49890.99−7.1 × 10−15 5.55.8 × 10−28 b 1.4 × 10−25 1.1 × 10−25 7.5 × 10−7 5.8 × 1031 200
J1748−2446M280.15560.29−4.4 × 10−15 5.55.8 × 10−28 b 1.0 × 10−25 9.5 × 10−26 1.6 × 10−6 1.2 × 1032 160
J1748−2446N115.38230.76−1.8 × 10−15 5.55.8 × 10−28 b 5.8 × 10−26 6.4 × 10−26 6.3 × 10−6 4.8 × 1032 110
J1748−2446O596.441192.87−9.4 × 10−15 5.55.8 × 10−28 b 2.6 × 10−25 2.6 × 10−25 9.6 × 10−7 7.4 × 1031 450
J1748−2446P f 578.501157.00−8.7 × 10−14 5.55.8 × 10−28 b 1.6 × 10−25  ⋅⋅⋅6.1 × 10−7 4.8 × 1031 267
J1748−2446Q355.64711.29−5.6 × 10−15 5.55.8 × 10−28 b 8.8 × 10−26 9.4 × 10−26 9.7 × 10−7 7.5 × 1031 160
J1748−2446R198.86397.73−3.2 × 10−15 5.55.8 × 10−28 b 8.1 × 10−26 5.1 × 10−26 1.7 × 10−6 1.3 × 1032 87
J1748−2446S163.49326.98−2.6 × 10−15 5.55.8 × 10−28 b 4.5 × 10−26 5.3 × 10−26 2.6 × 10−6 2.0 × 1032 91
J1748−2446T141.15282.29−2.2 × 10−15 5.55.8 × 10−28 b 5.1 × 10−26 3.0 × 10−26 2.0 × 10−6 1.5 × 1032 52
J1748−2446U304.03608.06−4.8 × 10−15 5.55.8 × 10−28 b  ⋅⋅⋅1.4 × 10−25 1.9 × 10−6 1.5 × 1032 230
J1748−2446V482.51965.02−7.6 × 10−15 5.55.8 × 10−28 b 1.3 × 10−25 1.3 × 10−25 7.2 × 10−7 5.6 × 1031 220
J1748−2446W237.80475.60−3.8 × 10−15 5.55.8 × 10−28 b 9.4 × 10−26 9.5 × 10−26 2.2 × 10−6 1.7 × 1032 160
J1748−2446X333.42666.83−5.3 × 10−15 5.55.8 × 10−28 b 8.3 × 10−26 4.8 × 10−26 5.7 × 10−7 4.4 × 1031 83
J1748−2446Y488.24976.49−7.7 × 10−15 5.55.8 × 10−28 b 2.1 × 10−25 2.1 × 10−25 1.2 × 10−6 9.0 × 1031 370
J1748−2446Z406.08812.15−6.4 × 10−15 5.55.8 × 10−28 b 8.5 × 10−26 8.6 × 10−26 6.8 × 10−7 5.2 × 1031 150
J1748−2446aa172.77345.54−2.7 × 10−15 5.55.8 × 10−28 b 2.3 × 10−25 1.5 × 10−25 6.6 × 10−6 5.1 × 1032 260
J1748−2446ab195.32390.65−3.1 × 10−15 5.55.8 × 10−28 b 4.7 × 10−26 4.0 × 10−26 1.4 × 10−6 1.1 × 1032 69
J1748−2446ac196.58393.17−3.1 × 10−15 5.55.8 × 10−28 b 7.2 × 10−26 5.6 × 10−26 1.9 × 10−6 1.5 × 1032 97
J1748−2446ad f 716.361432.701.7 × 10−14 5.55.8 × 10−28 b 1.8 × 10−25  ⋅⋅⋅4.5 × 10−7 3.5 × 1031 300
J1748−2446ae273.33546.66−4.3 × 10−15 5.55.8 × 10−28 b 6.6 × 10−26 6.8 × 10−26 1.2 × 10−6 9.2 × 1031 120
J1748−2446af302.63605.26−4.8 × 10−15 5.55.8 × 10−28 b 1.1 × 10−25 6.0 × 10−26 8.5 × 10−7 6.5 × 1031 100
J1748−2446ag224.82449.64−3.6 × 10−15 5.55.8 × 10−28 b 9.4 × 10−26 5.1 × 10−26 1.3 × 10−6 1.0 × 1032 87
J1748−2446ah201.40402.81−3.2 × 10−15 5.55.8 × 10−28 b 5.5 × 10−26 4.0 × 10−26 1.3 × 10−6 10.0 × 1031 69
J1748−2446ai47.1194.21−7.5 × 10−16 5.55.8 × 10−28 b  ⋅⋅⋅1.6 × 10−25 9.2 × 10−5 7.1 × 1033 270
J1751−2857255.44510.87−7.3 × 10−16 1.49.5 × 10−28  ⋅⋅⋅6.8 × 10−26 3.6 × 10−7 2.8 × 1031 72
J1756−225135.1470.27−1.3 × 10−15 2.91.7 × 10−27 9.7 × 10−26 7.4 × 10−26 4.1 × 10−5 3.2 × 1033 45
J1757−5322 e 112.74225.48−3.3 × 10−16 1.41.0 × 10−27 3.7 × 10−25  ⋅⋅⋅9.4 × 10−6 7.3 × 1032 360
J1801−1417275.85551.713.1 × 10−16 1.84.8 × 10−28 a 6.2 × 10−26 7.3 × 10−26 4.1 × 10−7 3.1 × 1031 150
J1801−3210134.16268.331.8 × 10−17 5.05.8 × 10−29  ⋅⋅⋅3.5 × 10−26 2.3 × 10−6 1.8 × 1032 610
J1802−212479.07158.13−4.4 × 10−16 3.35.7 × 10−28 a  ⋅⋅⋅4.7 × 10−26 6.0 × 10−6 4.6 × 1032 83
J1803−30140.82281.63−2.2 × 10−15 7.84.1 × 10−28 b 5.5 × 10−26 4.1 × 10−26 3.8 × 10−6 3.0 × 1032 100
J1804−073543.2986.58−6.9 × 10−16 8.43.8 × 10−28 b 8.4 × 10−26 8.8 × 10−26 9.4 × 10−5 7.2 × 1033 230
J1804−2717107.03214.06−4.7 × 10−16 1.21.4 × 10−27 2.4 × 10−26 2.2 × 10−26 5.3 × 10−7 4.1 × 1031 15
J1807−2459A326.86653.71−5.2 × 10−15 2.71.2 × 10−27 b 1.5 × 10−25 9.6 × 10−26 5.8 × 10−7 4.5 × 1031 81
J1810+1744601.411202.82−1.6 × 10−15 2.55.3 × 10−28  ⋅⋅⋅1.8 × 10−25 2.9 × 10−7 2.3 × 1031 340
J1810−200530.4760.93−5.0 × 10−17 4.02.6 × 10−28 a 2.2 × 10−25 1.6 × 10−25 1.6 × 10−4 1.3 × 1034 630
J1811−2405375.86751.71−1.9 × 10−15 1.71.1 × 10−27  ⋅⋅⋅8.5 × 10−26 2.4 × 10−7 1.9 × 1031 80
J1823−3021A183.82367.65−2.9 × 10−15 7.94.1 × 10−28 b 4.0 × 10−26 3.5 × 10−26 1.9 × 10−6 1.5 × 1032 86
J1824−2452A327.41654.81−1.7 × 10−13 5.53.4 × 10−27 a 7.9 × 10−26 5.5 × 10−26 6.7 × 10−7 5.2 × 1031 16
J1824−2452B f 152.75305.505.6 × 10−15 5.55.8 × 10−28 b 4.2 × 10−26  ⋅⋅⋅2.3 × 10−6 1.8 × 1032 72
J1824−2452C f 240.48480.96−9.8 × 10−15 5.55.8 × 10−28 b 6.6 × 10−26  ⋅⋅⋅1.5 × 10−6 1.2 × 1032 110
J1824−2452E f 184.53369.063.7 × 10−15 5.55.8 × 10−28 b 7.5 × 10−26  ⋅⋅⋅2.9 × 10−6 2.2 × 1032 130
J1824−2452F f 407.97815.94−1.6 × 10−15 5.55.8 × 10−28 b 9.8 × 10−26  ⋅⋅⋅7.6 × 10−7 5.9 × 1031 170
J1824−2452G f 169.23338.46−5.2 × 10−15 5.55.8 × 10−28 b 7.3 × 10−26  ⋅⋅⋅3.3 × 10−6 2.6 × 1032 130
J1824−2452H f 216.01432.02−3.6 × 10−15 5.55.8 × 10−28 b 8.2 × 10−26  ⋅⋅⋅2.3 × 10−6 1.8 × 1032 140
J1824−2452I g 254.33508.67−5.4 × 10−15 5.55.8 × 10−28 b 2.2 × 10−25  ⋅⋅⋅4.4 × 10−6 3.4 × 1032 370
J1824−2452J f 247.54495.084.7 × 10−15 5.55.8 × 10−28 b 1.1 × 10−25  ⋅⋅⋅2.3 × 10−6 1.7 × 1032 180
J1841+013033.5967.18−9.2 × 10−15 3.24.2 × 10−27 1.6 × 10−25 1.3 × 10−25 8.7 × 10−5 6.7 × 1033 31
J1843−1113541.811083.62−2.8 × 10−15 2.09.3 × 10−28 1.6 × 10−25 1.1 × 10−25 1.8 × 10−7 1.4 × 1031 120
J1853+1303244.39488.78−5.1 × 10−16 1.67.3 × 10−28 a  ⋅⋅⋅8.5 × 10−26 5.4 × 10−7 4.2 × 1031 120
J1857+0943186.49372.99−6.1 × 10−16 0.91.6 × 10−27 a 7.3 × 10−26 5.7 × 10−26 3.5 × 10−7 2.7 × 1031 35
J1903+0327465.14930.27−3.8 × 10−15 6.53.6 × 10−28 a  ⋅⋅⋅1.6 × 10−25 1.1 × 10−6 8.8 × 1031 450
J1905+0400264.24528.48−2.8 × 10−16 1.36.3 × 10−28 a 7.4 × 10−26 5.0 × 10−26 2.3 × 10−7 1.8 × 1031 80
J1909−3744339.32678.63−1.9 × 10−16 1.34.7 × 10−28 a 8.2 × 10−26 5.9 × 10−26 1.5 × 10−7 1.2 × 1031 120
J1910+1256200.66401.32−3.4 × 10−16 1.95.4 × 10−28 a  ⋅⋅⋅7.7 × 10−26 8.8 × 10−7 6.8 × 1031 140
J1910−5959A306.17612.33−1.9 × 10−16 4.51.4 × 10−28 a 7.7 × 10−26 5.5 × 10−26 6.2 × 10−7 4.8 × 1031 390
J1910−5959B119.65239.30−1.9 × 10−15 4.57.1 × 10−28 b 3.8 × 10−26 2.5 × 10−26 1.9 × 10−6 1.4 × 1032 35
J1910−5959C189.49378.981.1 × 10−18 4.51.4 × 10−29 a 4.4 × 10−26 3.2 × 10−26 9.5 × 10−7 7.4 × 1031 2300
J1910−5959D110.68221.35−1.8 × 10−15 4.57.1 × 10−28 b 3.1 × 10−26 2.1 × 10−26 1.8 × 10−6 1.4 × 1032 29
J1910−5959E218.73437.47−3.5 × 10−15 4.57.1 × 10−28 b 4.8 × 10−26 3.6 × 10−26 8.0 × 10−7 6.1 × 1031 50
J1911+0101A e 276.36552.715.0 × 10−16 7.44.3 × 10−28 b 6.0 × 10−25  ⋅⋅⋅1.4 × 10−5 1.1 × 1033 1400
J1911+0101B e 185.72371.456.9 × 10−17 7.44.3 × 10−28 b 7.4 × 10−25  ⋅⋅⋅3.8 × 10−5 2.9 × 1033 1700
J1911+1347216.17432.34−7.9 × 10−16 1.69.6 × 10−28 7.0 × 10−26 4.8 × 10−26 3.9 × 10−7 3.0 × 1031 50
J1911−1114275.81551.61−4.7 × 10−16 1.66.7 × 10−28 a 5.7 × 10−26 6.3 × 10−26 3.1 × 10−7 2.4 × 1031 94
J1918−0642130.79261.58−4.0 × 10−16 1.41.0 × 10−27 a  ⋅⋅⋅4.0 × 10−26 7.7 × 10−7 6.0 × 1031 39
J1939+2134641.931283.86−4.3 × 10−14 5.01.3 × 10−27 a 1.8 × 10−25 1.3 × 10−25 3.6 × 10−7 2.8 × 1031 96
J1944+0907192.86385.71−3.6 × 10−16 1.38.6 × 10−28 a  ⋅⋅⋅5.5 × 10−26 4.4 × 10−7 3.4 × 1031 64
J1955+2908163.05326.10−7.5 × 10−16 5.43.2 × 10−28 a 7.0 × 10−26 5.4 × 10−26 2.6 × 10−6 2.0 × 1032 170
J1959+2048622.121244.24−4.4 × 10−15 1.51.4 × 10−27 a  ⋅⋅⋅1.5 × 10−25 1.4 × 10−7 1.1 × 1031 110
J2007+272240.8281.64−1.6 × 10−15 6.87.4 × 10−28  ⋅⋅⋅7.1 × 10−26 6.9 × 10−5 5.3 × 1033 96
J2010−1323191.45382.90−1.8 × 10−16 1.36.0 × 10−28  ⋅⋅⋅6.3 × 10−26 5.2 × 10−7 4.0 × 1031 100
J2017+0603345.28690.56−9.6 × 10−16 1.31.0 × 10−27  ⋅⋅⋅1.3 × 10−25 3.4 × 10−7 2.6 × 1031 130
J2019+2425254.16508.32−1.7 × 10−16 0.97.2 × 10−28 a 9.2 × 10−26 5.6 × 10−26 1.9 × 10−7 1.4 × 1031 79
J2033+17168.10336.19−2.3 × 10−16 1.46.8 × 10−28 a 7.5 × 10−26 8.0 × 10−26 9.2 × 10−7 7.1 × 1031 120
J2043+1711420.19840.38−7.3 × 10−16 1.19.4 × 10−28 a  ⋅⋅⋅7.3 × 10−26 1.1 × 10−7 8.5 × 1030 78
J2051−0827221.80443.59−6.1 × 10−16 1.31.0 × 10−27 a 7.5 × 10−26 5.3 × 10−26 3.3 × 10−7 2.5 × 1031 51
J2124−3358202.79405.59−4.4 × 10−16 0.34.0 × 10−27 a 4.9 × 10−26 3.9 × 10−26 6.7 × 10−8 5.2 × 1030 9.9
J2129−5721268.36536.72−1.5 × 10−15 0.44.7 × 10−27 a 6.2 × 10−26 5.2 × 10−26 6.8 × 10−8 5.3 × 1030 11
J2140−2310A90.75181.50−1.4 × 10−15 9.23.5 × 10−28 b  ⋅⋅⋅5.9 × 10−26 1.6 × 10−5 1.2 × 1033 170
J2145−075062.30124.59−1.0 × 10−16 0.61.8 × 10−27 a 3.8 × 10−26 2.9 × 10−26 1.0 × 10−6 7.9 × 1031 16
J2214+3000320.59641.18−1.5 × 10−15 1.31.3 × 10−27  ⋅⋅⋅7.2 × 10−26 2.2 × 10−7 1.7 × 1031 54
J2215+5135383.20766.40−4.1 × 10−15 3.38.0 × 10−28  ⋅⋅⋅1.6 × 10−25 8.7 × 10−7 6.7 × 1031 200
J2229+2643335.82671.631.7 × 10−16 1.44.1 × 10−28 a 9.9 × 10−26 6.5 × 10−26 2.0 × 10−7 1.5 × 1031 160
J2241−5236457.31914.62−1.4 × 10−15 0.72.1 × 10−27  ⋅⋅⋅8.9 × 10−26 6.9 × 10−8 5.3 × 1030 42
J2302+4442192.59385.18−5.1 × 10−16 0.81.8 × 10−27  ⋅⋅⋅4.5 × 10−26 2.2 × 10−7 1.7 × 1031 26
J2317+1439290.25580.51−1.3 × 10−16 1.92.8 × 10−28 a 8.8 × 10−26 5.6 × 10−26 3.0 × 10−7 2.3 × 1031 200
J2322+2057207.97415.94−1.8 × 10−16 0.89.6 × 10−28 a 1.1 × 10−25 5.4 × 10−26 2.3 × 10−7 1.8 × 1031 57

Notes. Detector calibration errors mean that for pulsars with fgw below and above 50 Hz there are ∼6% and ∼20% uncertainties on these limits respectively. aThe pulsar's spin-down is corrected for proper motion effects. bThe pulsar's spin-down is calculated using a characteristic spin-down age of 109 yr. cThe pulsar's distance is calculated using the NE2001 model of Cordes & Lazio (2002). dRecently discovered pulsar with timing solution from after the end of S6 or VSR4. As this is an MSP the timing solution should extrapolate back well over the search period and it is unlikely that any glitches occurred. eResults from S3/S4 (Abbott et al. 2007a). fResults from S5 (Abbott et al. 2010). gNew result from S5 using data covering 2006.

Download table as:  ASCIITypeset images: 1 2 3

The gravitational wave amplitude upper limits as a function of frequency are plotted in Figures 1 and 2 (showing a version just containing the seven high interest pulsars), which also show bands giving the expected sensitivity of the analysis. The upper limits in histogram form for all pulsars can be seen in Figure 3. The histograms show that the distribution of h0 upper limits is peaked just below 10−25, corresponding to equivalent peaks on ε and Q22 of ∼10−6 and ∼10−32 kg m2. The spin-down limit ratios shows that we are within a factor of 100 for just over half of the pulsars. It is interesting to see that due to the shape of the detector sensitivity curves the lower frequency young pulsars (analyzed only with Virgo data) have the highest amplitude limits, but as several have high spin-down luminosities they have an approximately uniform spread in spin-down limit ratios.

Figure 2.

Figure 2. Zoomed version of Figure 1 focusing on the seven high interest pulsars. The outlier at ∼32 Hz is J1833−1034 for which only VSR2 data was used.

Standard image High-resolution image
Figure 3.

Figure 3. Upper limits in histogram form for all pulsars for h0, ε, Q22 and the spin-down limit ratio. The gray shaded area represents results from the S6/VSR2,4 analysis combining all detectors, the blue shaded area represents results from the VSR2,4-only analyses. These also contain the seven high interest pulsars for which the Bayesian method values have been plotted based on no assumptions about the pulsar orientations. Previous upper limits from the S5 analysis are given by the unfilled histogram.

Standard image High-resolution image

4. DISCUSSION

We have seen no credible evidence for gravitational wave emission from any known pulsar, but have been able to place upper limits on the gravitational wave amplitude from an unprecedented number of pulsars. In this work we have produced entirely new results for 73 pulsars and updated the results of previously searches for 106 pulsars, with results from a further 16 from previous analyses being reproduced here. A total of 195 pulsars have now been targeted over the lifetimes of the first generation of interferometric gravitational wave detectors.

4.1. Quadrupole Estimates

As discussed in the introduction, we have targeted the gravitational-wave signature of the time-varying l = m = 2 quadrupole moment. There is great uncertainty, however, as to whether neutron stars can form and sustain sufficient elastic deformations to give an observable quadrupole, and this, in turn, makes it difficult to model a realistic source population. The recent work by Johnson-McDaniel & Owen (2013) (also see, e.g., Owen 2005; Pitkin 2011) on the maximum sustainable quadrupole for a variety of neutron star equations of state indicates that relatively large quadrupoles can indeed be sustained. Johnson-McDaniel & Owen (2013) find that solid quark stars could sustain quadrupoles of up to 1037 kg m2 (or fiducial ellipticities of order 0.1), hybrid stars could sustain quadrupoles of up to 1035 kg m2 (or fiducial ellipticities of order 1 × 10−3), while for normal neutron stars the stiffest equations of state allow quadrupoles of ∼1 × 1033 kg m2 (or fiducial ellipticities of ∼1 × 10−5). It is worth noting that these are maximum allowable quadrupoles, and it is still unknown whether they are realized in nature for reasons described in Abbott et al. (2007b).

A mass quadrupole may also be generated by distortional pressure from the star's magnetic field (see, e.g., Bonazzola & Gourgoulhon 1996; Cutler 2002; Ciolfi et al. 2010). The external dipole field of a pulsar is usually estimated from its rotational spin-down, assuming this is due to magnetic dipole radiation (equivalent to the gravitational wave spin-down limit that we define). As discussed in Section 2, this gives external surface dipole field strengths of ∼109 G for MSPs and ∼1012 G for normal pulsars. Internal fields of this magnitude are too small to induce mass quadrupoles that would be currently observable, but the field strengths of some magnetars are at a suitable level (though rotating too slowly to be detectable sources for ground-based gravitational wave detectors). Unfortunately, internal field strengths and configurations are not well understood, and the mechanisms for burying fields beneath the surface are uncertain. Studies of one young pulsar with a braking index of n ≈ 1 (Espinoza et al. 2011) may point towards an evolving and increasing external magnetic field, with an internal field leaking out over time, but recently other mechanisms have been proposed to explain the evolution of the field that do not rely on an increasing magnetic field (e.g., Ho & Andersson 2012; Çalişkan et al. 2013). Mastrano & Melatos (2012) discuss the prospects of constraining field strength and configuration for recycled MSPs using gravitational wave data. This is also discussed in Pitkin (2011), who shows limits that could be obtained on fully poloidal or toroidal field configurations. Further estimates of the quadrupoles that can be generated by internal magnetic fields for a given equation of state are given in Haskell et al. (2008, 2009), Akgun & Wasserman (2007).

4.2. High Interest Pulsars

For the seven high interest pulsars the results are all close to (or beat) the spin-down limits. In particular, our upper limits are significantly below the spin-down limit for the Crab and Vela pulsars, further improving over past results. The mass quadrupole limits are generally within 1034–1035 kg m2, with the Crab pulsar upper limit slightly lower at ∼7 × 1033 kg m2. Therefore, for these stars to emit gravitational waves at current sensitivities the emission would most likely have to come from a quark star or one with a hybrid core, whilst the Crab pulsar is about an order of magnitude above the maximum quadrupoles expected for purely crustal emission. However, for advanced detectors the sensitivity for Crab pulsar would be consistent with most optimistic predictions for normal neutron stars. For J0537−6910, which has a quadrupole limit close to the Crab pulsar, future prospects may not be so good for reaching the most optimistic prediction for normal neutron stars. This is due to the requirement for phase coherent timing, which for these analyses relied on the no-longer-operational RXTE.

For the Crab and Vela pulsars, our results now limit the gravitational wave emission to contribute ≲ 1% and ≲ 10% of their respective spin-down luminosities, with an improvement of about a factor of 4 for Vela with respect to previous results. These limits are compatible with the observed braking indices of the pulsars, which are n = 2.51 and n ≈ 1.4 respectively (see, e.g., Palomba 2000).

Given various assumptions about the magnetic field discussed above, our results constrain the internal field of the Crab pulsar to be less than ∼1016 G (e.g., Cutler 2002). For the other high interest pulsars, the limits on the magnetic field would be even higher than this, so we have not included them here.

Johnson-McDaniel (2013) relates the limits on the l = m = 2 quadrupole moment from the gravitational wave emission to the physical surface deformation of a star for a variety of equations of state, which can be compared to the oblateness due to rotation (note that there is no particular reason to expect a relation between these quantities). His results showed that previous gravitational wave limits for the Crab pulsar constrained the surface deformation from the l = m = 2 quadrupole to be well below the rotational deformation for all equations of state and neutron star masses. Our new results slightly improve these limits, with the physical surface deformation limited to less than ∼30 cm, maximized over masses and equations of state. For the Vela pulsar, our new results limit non-axisymmetric quadrupole deformations to be ≲ 100 cm, which is smaller than the expected rotational oblateness for equations of state with large radii.

For PSR J0537−6910, the quality of S6 data at the corresponding frequency was relatively poor, and the upper limits are no better than those produced during S5 (Abbott et al. 2010). If this pulsar were, however, at the upper end of the moment of inertia range (∼3 × 1038 kg m2) the spin-down limit would be increased by a factor of ∼1.7, and we would now fall below it. 155

4.3. Other Highlights

Several other pulsar upper limits are within a factor of 10 of their spin-down limits. For the MSPs, three upper limits are within a factor of ten of the spin-down limit: J1045−4509, a factor of 6; J1643−1224, a factor of 10; and J2124−3358, also a factor of 10. The upper limit that corresponds to the smallest ellipticity/mass quadrupole is from J2124−3358 with ε = 6.7 × 10−8 and Q22 = 5.2 × 1030 kg m2. Although this value is currently above the spin-down limit, it is well within allowable maximum deformations for all neutron star equations of state (see, e.g., Pitkin 2011). The gravitational wave spin-down limits for these pulsars require quadrupoles that are well within reasonable theoretical ranges, so they will make intriguing targets for the advanced generation of detectors.

For the young pulsars only targeted with Virgo VSR2 and VSR4 data a further five are within a factor of 10 of the spin-down limit (see Table 6). All of these would be required to have an exotic equation of state to be observed at around their spin-down limits in future detectors.

4.4. Future Prospects

The search results described in this paper assume that the pulsar gravitational-wave phase evolution is very well known and tied very closely to the observed electromagnetic phase. However, precession (e.g., Zimmermann & Szedenits 1979; Jones & Andersson 2002) or other models (Jones 2010) could give emission at both the rotation frequency and twice the rotation frequency. Additionally, as discussed in Abbott et al. (2008), emission may be offset from the electromagnetic phase model. We therefore will be applying methods to search for gravitational waves from known pulsars at multiple harmonics and with narrow bandwidths around the observed electromagnetic values in archival and future datasets.

We look forward to the era of the advanced LIGO (aLIGO) (Harry & LIGO Scientific Collaboration 2010) and Virgo (AdV) (Acernese et al. 2009; Accadia et al. 2012) detectors (see Aasi et al. 2013c, for estimates of the aLIGO and AdV observation schedule and sensitivity evolution), as well as the KAGRA detector (Somiya 2012). Ongoing radio pulsar surveys are discovering new objects that will be targeted with future detectors. Currently, the High Time Resolution Universe survey with the Parkes and Effelsberg telescopes (Keith et al. 2010) has discovered 29 new MSPs (Keith 2013; Ng et al. 2013) and could discover up to ∼75 once complete. The high sensitivity Arecibo PALFA survey is discovering new pulsars (Lazarus 2013) and making use of distributed computing through Einstein@home (Allen et al. 2013). The Green Bank Drift-scan survey and the Green Bank North Celestial Cap survey are also discovering new and interesting sources (Lynch et al. 2013). Many interesting high energy pulsars, undetectable in the radio frequency band, are also being detected by the Fermi LAT (Saz Parkinson et al. 2013). Fermi is also providing targets to facilitate radio searches which are finding many new MSPs. In addition, new analyses of archive data, such as using Einstein@home to search through Parkes Multi-beam Pulsar Survey data, are still yielding new results (Knispel et al. 2013). In the near future, there are exciting prospects from the Low Frequency Array, which could detect the majority of radio pulsars within ∼2 kpc, giving of order 1000 new pulsars (van Leeuwen & Stappers 2010; Stappers et al. 2011), and perform deep searches for pulsars in globular clusters.

Finally, we should emphasize that known pulsar searches are not the only searches looking for gravitational waves from rotating, galactic neutron stars. There have been, or are under way, several directed searches looking for sources of unknown frequency and spin-down in particular objects, e.g., globular clusters, supernova remnants (e.g., Abadie et al. 2010; Chung et al. 2011), the Galactic center (Aasi et al. 2013a), and low-mass X-ray binaries. There are also several semi-coherent, all-sky, wide-frequency band searches (e.g., Abadie et al. 2012; Aasi et al. 2013b). Very similar pipelines will be used during the advanced detector era, yielding signal candidates, performing follow-ups and, in case of detection, source parameter estimation.

The authors gratefully acknowledge the support of the United States National Science Foundation for the construction and operation of the LIGO Laboratory, the Science and Technology Facilities Council of the United Kingdom, the Max-Planck-Society, and the State of Niedersachsen/Germany for support of the construction and operation of the GEO600 detector, and the Italian Istituto Nazionale di Fisica Nucleare and the French Centre National de la Recherche Scientifique for the construction and operation of the Virgo detector. The authors also gratefully acknowledge the support of the research by these agencies and by the Australian Research Council, the International Science Linkages program of the Commonwealth of Australia, the Council of Scientific and Industrial Research of India, the Istituto Nazionale di Fisica Nucleare of Italy, the Spanish Ministerio de Economía y Competitividad, the Conselleria d'Economia Hisenda i Innovació of the Govern de les Illes Balears, the Foundation for Fundamental Research on Matter supported by the Netherlands Organisation for Scientific Research, the Polish Ministry of Science and Higher Education, the FOCUS Programme of Foundation for Polish Science, the Royal Society, the Scottish Funding Council, the Scottish Universities Physics Alliance, the National Aeronautics and Space Administration, OTKA of Hungary, the Lyon Institute of Origins (LIO), the National Research Foundation of Korea, Industry Canada and the Province of Ontario through the Ministry of Economic Development and Innovation, the National Science and Engineering Research Council Canada, the Carnegie Trust, the Leverhulme Trust, the David and Lucile Packard Foundation, the Research Corporation, and the Alfred P. Sloan Foundation. The Nançay Radio Observatory is operated by the Paris Observatory, associated with the French Centre National de la Recherche Scientifique. LIGO Document No. LIGO-P1200104.

APPENDIX: EPHEMERIS FOR J0537−6910

Over the span of S6, VSR2, and VSR4 RXTE made observations of J0537−6910. It was observed to glitch four times during this period and phase connected ephemerides were produced for each inter-glitch segment. These ephemerides, given in Table 8, use a DE200 sky position of α = 05h37m47fs36 and δ = −69°10'20farcs4 (Wang et al. 2001).

Table 8.  RXTE Ephemerides for J0537−6910 During the Period of S6, VSR2, and VSR4

EM Observation SpanEpoch t0 frot $\dot{f}_{\rm rot}$ $\ddot{f}_{\rm rot}$
(MJD and Date)(MJDTDB)(Hz)(10−10 Hz s−1)(10−20 Hz s−2)
54897 (2009 Mar 7) – 55041 (2009 Jul 29)54966.426661602261.9765822442−1.99483000.960
55045 (2009 Aug 2) – 55182 (2009 Dec 17)55140.037772713361.9736036203−1.99399242.090
55185 (2009 Dec 20) – 55263 (2010 Mar 8)55221.056291231861.9722204755−1.99519732.090
55275 (2010 Mar 20) – 55445 (2010 Sep 6)55313.648883824461.9706584446−1.99533300.673
55458 (2010 Sep 19) – 55503 (2010 Nov 3)55475.685920861561.9678758925−1.99525740.673

Download table as:  ASCIITypeset image

Footnotes

  • 144 

    We use "triaxial neutron star" as shorthand for a star with some asymmetry with respect to its rotation axis and therefore a triaxial moment of inertia ellipsoid.

  • 145 

    As noted in Johnson-McDaniel (2013), the versions of this equation given inline in the first paragraph of Abbott et al. (2008), as Equation (1) in Abbott et al. (2010) and as Equation (14) in Abadie et al. (2011) are incorrect and should have I38 substituted for $I_{38}^{1/2}$.

  • 146 

    Again, as noted in Johnson-McDaniel (2013), the versions of this equation given inline in Section 3 of Abbott et al. (2008) and as Equation (7) in Abbott et al. (2010) are incorrect and should have I38 substituted for $I_{38}^{-1}$.

  • 147 

    For this analysis the parameter posterior distributions were recreated using a Markov chain Monte Carlo (Abbott et al. 2010). For each pulsar five independent chains were produced with 50,000 burn-in samples and 200,000 posterior samples in each. The chains were thinned using the autocorrelation length to give uncorrelated samples, and to test for convergence, the chains were then examined by eye, and a Gelman-Rubins test was performed (see, e.g., Brooks & Roberts 1998).

  • 148 

    For this analysis, the results were produced with version 6.16 of the LSC Algorithm Library Suite (LALSuite) https://www.lsc-group.phys.uwm.edu/daswg/projects/lalsuite.html.

  • 149 

    The $\mathcal {F}$/$\mathcal {G}$-statistic is most suitable for signal detection, whilst the upper limit derived from it here is mainly given for completeness. In the future a more sophisticated method, such as that of Feldman & Cousins (1998), may be used to produce frequentist confidence intervals for this analysis.

  • 150 
  • 151 

    Characteristic age is given by $\tau = -(1/(n-1))(f_{\rm rot}/\dot{f}_{\rm rot})$, which, for a magnetic dipole braking index of n = 3, gives $\tau = -f_{\rm rot}/(2\dot{f}_{\rm rot})$, and for purely gravitational wave (quadrupole) spin-down would be n = 5, giving $\tau = -f_{\rm rot}/(4\dot{f}_{\rm rot})$ (a "gravitar," Palomba 2005; Knispel & Allen 2008).

  • 152 

    The five additional pulsars targeted in S3/S4 were J1435−6100, J1629−6902, J1757−5322, J1911+0101A and J1911+0101B and the eleven additional pulsars targeted in S5 were J1701−3006B, J1701−3006C, J1748−2446P, J1748−2446ad, J1824−2452B, J1824−2452C, J1824−2452E, J1824−2452F, J1824−2452G, J1824−2452H, J1824−2452J.

  • 153 

    Note that Johnson-McDaniel (2013) computes even larger potential moments of inertia at ∼5 × 1038 kg m2 for some solid quark stars.

  • 154 

    In Ng & Romani (2008) the inclination is denoted by ζ and the position angle Ψ is equivalent to our polarization angle. Our searches are insensitive to rotations of 90° in the polarization angle, so our quoted values are rotated into the range −45° < ψ < 45°.

  • 155 

    The distance to the Large Magellanic Cloud is known to ∼2% (Pietrzyński et al. 2013), so does not significantly contribute to the uncertainty on the spin-down limit.

Please wait… references are loading.
10.1088/0004-637X/785/2/119