PSR J1939 Research Articles

Estimating 5-year rotation stability of PSR B1937+21 using NICER observations

Using observations of millisecond pulsars (MSPs) can construct a time scale that is independent of terrestrial time standards. We investigate the 5-year rotation stability of PSR B1937 + 21 (PSR J1939 + 2134) using the Neutron Star Interior Composition Explorer (NICER) observations and compare it with that of the radio data from North American Nanohertz Observatory for Gravitational Waves (NANOGrav). We implement a data screening method to improve the accuracy of X-ray pulse times of arrival (TOAs), and develop a data-analysis framework based on an iterative least square technique for the X-ray timing parameter fitting. Finally, we carefully consider the main factors affecting the pulsar rotation stability. We find that the root mean square error (RMS) of the timing residuals of PSR B1937 + 21 based on NICER observations is 4.7 μs (1σ), and the 5-year rotation stability is 1.7 × 10−14 the accuracy of which is comparable to that of NANOGrav data. In order to ensure less than 10−14 of the main factors influencing the rotation stability on 5-year time scale, the errors of R.A. and Dec shall be no more than 2 mas and 6 mas respectively, while the amplitudes of the white noise and red noise must be less than 5 μs and 1 μs, respectively.

Long-term scintillation studies of EPTA pulsars

Context. Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. From the frequency scale of scintillation, one can estimate the geometric time delays from multipath propagation, a source of (typically) unmodeled, correlated noise in pulsar timing. Interstellar scintillation analysis of well-timed pulsars is useful to quantify the effects of time delays and may lead to improved timing precision, enhancing the probability of detecting gravitational waves. Aims. Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars, employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes in the 21-cm and 11-cm bands. Additionally, we aim to identify future possible lines of study using this long-term scintillation data set. Methods. The autocorrelation function of dynamic spectra has been used to estimate the scintillation bandwidth vd and scintillation timescale τd. Results. We present the long-term time series of vd and τd for 13 pulsars. Sanity checks and comparisons indicate that the scintillation parameters of our work and previously published works are mostly consistent. For two pulsars, PSRs J1857+0943 and J1939+2134, we were able to obtain measurements of the vd at both bands, which allowed us to derive the time series of frequency scaling indices with a mean and a standard deviation of 2.82 ± 1.95 and 3.18 ± 0.60, respectively. We found some interesting features which will be studied in more detail in subsequent papers in this series: (i) in the time series of PSR J1939+2134, where vd and τd sharply decrease associated with a sharp increase in the dispersion measure; (ii) PSR J0613-0200 and PSR J0636+5126 show a strong annual variation in the time series of the τd; and (iii) PSR J1939+2134 shows a weak anticorrelation between the scintillation timescale and the dispersion in Westerbork Synthesis Radio Telescope data.

Study of 72 Pulsars Discovered in the PALFA Survey: Timing Analysis, Glitch Activity, Emission Variability, and a Pulsar in an Eccentric Binary

Abstract We present new discoveries and results from long-term timing of 72 pulsars discovered in the Pulsar Arecibo L-band Feed Array (PALFA) survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages ∼30 kyr) with no apparent supernova remnant associations, three mode-changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age ∼1 Gy), and PSR J1954+2529, which likely belongs to a newly emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a nonrecycled neutron star with a 931 ms spin period in an eccentric (e = 0.114) wide (P b = 82.7 days) orbit with a companion of undetermined nature having a minimum mass of ∼0.6 M ⊙. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSPs) have twice the dispersion measure per unit spin period than the known population of MSP in the plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, rotating radio transients, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought.

Space/ground based pulsar timescale for comprehensive PNT system

The comprehensive positioning navigation timing (PNT) system in China is a multi-source information fusion system with BeiDou navigation satellite system (BDS) as a core. The high-precision millisecond pulsar timing can enhance the long-term stability of the BDS time benchmark and maintain a space-time benchmark for future deep-space users. In this paper, a ground-based pulsar time service system is proposed for detecting and improving the time benchmark of BDS. The preliminary designs and functions of the system are outlined. At the same time, the method of establishing space and ground-based pulsar time is studied. The ground radio timing data from the international pulsar timing array (IPTA), the X-ray timing data from the neutron star interior composition explorer (NICER) in space, and the simulation data from the 500-meter spherical radio telescope (five-hundred-meter aperture spherical radio telescope, FAST) for three millisecond pulsars are used to analyze the stability of ground/space-based pulsar time. The research results are as follows. The annual stability of the PSR J0437-4715 ground-based pulsar time based on IPTA data is 3.30 × 10<sup>–14</sup>, and the 10-year stability is 1.23 × 10<sup>–15</sup>, respectively. The existence of pulsar red noise can reduce the time stability of the pulsar. The annual stability of the PSR J1939+2134 ground-based pulsar time is 6.51 × 10<sup>–12</sup>. We find that the accuracy of the pulse time of Arrival(TOA) is an important factor that restricts the stability of space-based pulsar time. Based on NICER space X-ray timing data, the stability of the pulsar time for PSR J1824-2452A is 1.36 × 10<sup>–13</sup> in one year. Finally, the simulation analysis of the FAST’s data without considering the influence of red noise is completed, and we find that the PSR J1939+2134 ground-based pulsar time based on the FAST has an annual stability of 2.55 × 10<sup>–15</sup>, a 10-year stability of 1.39 × 10<sup>–16</sup>, and a 20-year stability of 5.08 × 10<sup>–17</sup>. It demonstrates that the powerful pulsar observation capability of FAST will help to improve the accuracy of ground-based pulsar time and enhance the long-term stability of the comprehensive PNT system time benchmark in China.

The 31 yr Rotation History of the Millisecond Pulsar J1939+2134 (B1937+21)

Abstract The timing properties of the millisecond pulsar PSR J1939+2134—very high rotation frequency, very low time derivative of rotation frequency, no timing glitches, and relatively low timing noise—are responsible for its exceptional timing stability over decades. It has been timed by various groups since its discovery, at diverse radio frequencies, using different hardware and analysis methods. Most of the timing data is now available in the public domain in two segments, which have not been combined so far. This work analyzes the combined data by deriving uniform methods of data selection, derivation of Dispersion Measure (DM), accounting for correlation due to “red” noise, etc. The timing noise of this pulsar is very close to a sinusoid, with a period of approximately 31 yr. The main results of this work are: (1) the clock of PSR J1939+2134 is stable at the level of almost one part in 1015 over about 31 yr; (2) the power-law index of the spectrum of electron density fluctuations in the direction of PSR J1939+2134 is 3.86 ± 0.04; (3) a Moon-sized planetary companion, in an orbit of semimajor axis about 11 astronomical units and eccentricity ≈0.2, can explain the timing noise of PSR J1939+2134; (4) precession under electromagnetic torque with very small values of oblateness and wobble angle can also be the explanation but with reduced confidence; and (5) there is an excess timing noise of about 8 μs amplitude during the epochs of steepest DM gradient, of unknown cause.

Physics of radio emission in gamma-ray pulsars

Propagation of radio emission in a pulsar magnetosphere is reviewed. The effects of polarization transfer, induced scattering and reprocessing to high energies are analysed with a special emphasis on the implications for the gamma-ray pulsars. The possibilities of the pulsar plasma diagnostics based on the observed radio pulse characteristics are also outlined. As an example, the plasma number density profiles obtained from the polarization data for the Vela and the gamma-ray millisecond pulsars J1446-4701, J1939+2134 and J1744-1134 are presented. The number densities derived tend to be the highest/lowest when the radio pulse leads/lags the gamma-ray peak. In the PSR J1939+2134, the plasma density profiles for the main pulse and interpulse appear to fit exactly the same curve, testifying to the origin of both radio components above the same magnetic pole and their propagation through the same plasma flow in opposite directions. The millisecond radio pulse components exhibiting flat position angle curves are suggested to result from the induced scattering of the main pulse by the same particles that generate gamma rays. This is believed to underlie the wide-sense radio/gamma-ray correlation in the millisecond pulsars. The radio quietness of young gamma-ray pulsars is attributed to resonant absorption, whereas the radio loudness to the radio beam escape through the periphery of the open field line tube.

MODELING PHASE-ALIGNED GAMMA-RAY AND RADIO MILLISECOND PULSAR LIGHT CURVES

Since the discovery of the first eight gamma-ray millisecond pulsars (MSPs) by the Fermi Large Area Telescope, this population has been steadily expanding. Four of the more recent detections, PSR J00340534, PSR J1939+2134 (B1937+21; the first MSP ever discovered), PSR J1959+2048 (B1957+20; the first discovery of a black widow system), and PSR J2214+3000, exhibit a phenomenon not present in the original discoveries: nearly phase-aligned radio and gamma-ray light curves (LCs). To account for the phase alignment, we explore models where both the radio and gamma-ray emission originate either in the outer magnetosphere near the light cylinder or near the polar caps. Using a Markov Chain Monte Carlo technique to search for best-fit model parameters, we obtain reasonable LC fits for the first three of these MSPs in the context of altitude-limited outer gap (alOG) and two-pole caustic (alTPC) geometries (for both gamma-ray and radio emission). These models differ from the standard outer gap (OG)/two-pole caustic (TPC) models in two respects: the radio emission originates in caustics at relatively high altitudes compared to the usual conal radio beams, and we allow both the minimum and maximum altitudes of the gamma-ray and radio emission regions to vary within a limited range (excluding the minimum gamma-ray altitude of the alTPC model, which is kept constant at the stellar radius, and that of the alOG model, which is set to the position-dependent null charge surface altitude). Alternatively, phase-aligned solutions also exist for emission originating near the stellar surface in a slot gap scenario (low-altitude slot gap (laSG) models). We find that the alTPC models provide slightly better LC fits than the alOG models, and both of these give better fits than the laSG models (for the limited range of parameters considered in the case of the laSG models). Thus, our fits imply that the phase-aligned LCs are likely of caustic origin, produced in the outer magnetosphere, and that the radio emission for these pulsars may come from close to the light cylinder. In addition, we were able to constrain the minimum and maximum emission altitudes with typical uncertainties of 30% of the light cylinder radius. Our results therefore describe a third gamma-ray MSP subclass, in addition to the two previously found by Venter et al.: those with LCs fit by standard OG/TPC models and those with LCs fit by pair-starved polar cap models.

Stokes tomography of radio pulsar magnetospheres - II. Millisecond pulsars

The radio polarization characteristics of millisecond pulsars (MSPs) differ significantly from those of non-recycled pulsars. In particular, the position angle (PA) swings of many MSPs deviate from the S-shape predicted by the rotating vector model, even after relativistic aberration is accounted for, indicating that they have non-dipolar magnetic geometries, likely due to a history of accretion. Stokes tomography uses phase portraits of the Stokes parameters as a diagnostic tool to infer a pulsar's magnetic geometry and orientation. This paper applies Stokes tomography to MSPs, generalizing the technique to handle interpulse emission. We present an atlas of look-up tables for the Stokes phase portraits and PA swings of MSPs with current-modified dipole fields, filled core and hollow cone beams, and two empirical linear polarization models. We compare our look-up tables to data from 15 MSPs and find that the Stokes phase portraits for a current-modified dipole approximately match several MSPs whose PA swings are flat or irregular and cannot be reconciled with the standard axisymmetric rotating vector model. PSR J1939+2134 and PSR J0437$-$4715 are modelled in detail. The data from PSR J1939+2134 at 0.61\,GHz can be fitted well with a current-modified dipole at $(\alpha, i) = (22 \pm 2^\circ, 80 \pm 1^\circ)$ and emission altitude 0.4 $r_\text{LC}$. The fit is less accurate for PSR J1939+2134 at 1.414\,GHz, and for PSR J0437$-$4715 at 1.44\,GHz, indicating that these objects may have a more complicated magnetic field geometry, such as a localized surface anomaly or a polar magnetic mountain.

A search for periodic gravitational waves from three recycled pulsars using the AURIGA detector

<i>Aims. <i/>We report on a search for continuous gravitational wave emission from three recycled radio pulsars, performed by using the data of the resonant detector AURIGA. Given the spin rate of the selected targets – the isolated pulsar PSR J1939+2134 and the binary pulsars PSR J0024-7204J and PSR J0218+4232 – the expected frequency of the emitted gravitational waves falls in the high sensitivity range of the detector.<i>Methods. <i/>The main topic is the method, meaning that the statistical analysis is performed by implementing a slightly modified version of the Feldman and Cousins Unified Approach.<i>Results. <i/>By using ephemerides provided by suitable radio observations of the targets, we were able to demodulate the Doppler shifts within a coherence time of 1 day, and then incoherently sum 10 daily spectra collected from December 8th to December 17th, 2006. We have found upper limits on the gravitational wave amplitudes in the order of a few units of 10<sup>-23<sup/> at 90% Confidence Level (C.L.), which translate to limits in the ellipticity of the targeted pulsars of at 90% C.L.<i>Conclusions. <i/>The same framework can then be applied to data coming from most sensitive experiments as VIRGO or LIGO; moreover, an application to recently discovered transients in X-ray pulsars is discussed.

Upper limits on the strength of periodic gravitational waves from PSR J1939+2134

The first science run of the LIGO and GEO gravitational wave detectors presented the opportunity to test methods of searching for gravitational waves from known pulsars. Here we present new direct upper limits on the strength of waves from the pulsar PSR J1939+2134 using two independent analysis methods, one in the frequency domain using frequentist statistics and one in the time domain using Bayesian inference. Both methods show that the strain amplitude at Earth from this pulsar is less than a few times 10−22.