Coherent Timing Solution Research Articles

Green Bank Telescope Discovery of the Redback Binary Millisecond Pulsar PSR J0212+5321

We report the discovery of a 2.11 ms binary millisecond pulsar during a targeted search of the redback optical candidate coincident with the γ-ray source 3FGL J0212.5+5320 using the Robert C. Byrd Green Bank Telescope (GBT) with the Breakthrough Listen backend at L band. Over a seven month period, five pointings were made near inferior conjunction of the pulsar in its 20.9 hr orbit, resulting in two detections, lasting 12 and 42 minutes. The pulsar dispersion measure (DM) of 25.7 pc cm−3 corresponds to a distance of 1.15 kpc in the NE2001 Galactic electron density model, consistent with the Gaia parallax distance of 1.16 ± 0.03 kpc for the companion star. We suspect the pulsar experiences wide orbital eclipses, similar to other redbacks, as well as scintillation and DM delays caused by its interaction with its companion and surroundings. Although the pulsar was only detected over ≈3.7% of the orbit, its measured acceleration is consistent with published binary parameters from optical radial velocity spectroscopy and light-curve modeling of the companion star, and it provides a more precise mass ratio and a projected semimajor axis for the pulsar orbit. We also obtained a refined optical photometric orbit ephemeris, and observed variability of the tidally distorted companion over 7 yr. A hard X-ray light curve from NuSTAR shows expected orbit-modulated emission from the intrabinary shock. The pulsar parameters and photometric ephemeris greatly restrict the parameter space required to search for a coherent timing solution including pulsar spin-down rate, either using Fermi γ-rays or further radio pulse detections.

Long term radio and X-ray evolution of the magnetar Swift J1818.0-1607

ABSTRACT We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift J1818.0–1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of −1.37 × 10−11 Hz s−1. The magnitude of the variation in the frequency derivative reduces with time before converging on the mean value. This corresponds to a characteristic age of ∼ 860 years, 2–4 times more than previously estimated. We were able to identify four states in the spin-frequency derivative that were quantified by the amount of modulation about the mean value and the transition between these states seem to be correlated with the change in the radio emission of the magnetar while no correlation is seen in the average radio profile variability on a shorter time-scale (days). The 0.5–12 keV X-ray flux shows a monotonic decrease that can be attributed to thermal emission from a hot spot on the surface of the neutron star that is reducing in size. Such decrease is consistent with what is seen in other magnetars. The potential correlation between the radio emission mode and the behaviour of the spin-down rate hints to a global change in the magnetopshere of the magnetar akin to the correlation seen in a subset of mode-changing radio pulsars and suggests a physical link between the two sub-populations.

Precise optical timing of PSR J1023+0038, the first millisecond pulsar detected with Aqueye+ in Asiago

Abstract We report the first detection of an optical millisecond pulsar with the fast photon counter Aqueye+ in Asiago. This is an independent confirmation of the detection of millisecond pulsations from PSR J1023+0038 obtained with SiFAP at the Telescopio Nazionale Galileo. We observed the transitional millisecond pulsar PSR J1023+0038 with Aqueye+ mounted at the Copernicus telescope in 2018 January. Highly significant pulsations were detected. The rotational period is in agreement with the value extrapolated from the X-ray ephemeris, while the time of passage at the ascending node is shifted by 11.55 ± 0.08 s from the value predicted using the orbital period from the X-rays. An independent optical timing solution is derived over a baseline of a few days that has an accuracy of ∼0.007 in pulse phase (∼12 μs in time). This level of precision is needed to derive an accurate coherent timing solution for the pulsar and to search for possible phase shifts between the optical and X-ray pulses using future simultaneous X-ray and optical observations.

Discovery of a glitch in the accretion-powered pulsar SXP 1062

We present timing analysis of the accretion powered pulsar SXP 1062, based on the observations of \textit{Swift}, \textit{XMM-Newton} and \textit{Chandra} satellites covering a time span of about 2 years. We obtain a phase coherent timing solution which shows that SXP 1062 has been steadily spinning down with a rate $-\,4.29(7) \times 10^{-14}$ Hz s$^{-1}$ leading to a surface magnetic field estimate of about $1.5 \times 10^{14}$ G. We also resolve the binary orbital motion of the system from X-ray data which confirms an orbital period of 656(2) days. On MJD 56834.5, a sudden change in pulse frequency occurs with $\Delta\nu = 1.28(5) \times 10^{-6} $ Hz, which indicates a glitch event. The fractional size of the glitch is $\Delta \nu / \nu \! \sim \! 1.37(6) \times 10^{-3}$ and SXP 1062 continues to spin-down with a steady rate after the glitch. A short X-ray outburst 25 days prior to the glitch does not alter the spin-down of the source; therefore the glitch should be associated with the internal structure of the neutron star. While glitch events are common for isolated pulsars, the glitch of SXP 1062 is the first confirmation of the observability of this type of events among accretion powered pulsars. Furthermore, the value of the fractional change of pulse frequency ensures that we discover the largest glitch reported up to now.

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

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

Pulse phase-coherent timing and spectroscopy of CXOU J164710.2−45521 outbursts

We present a long-term phase-coherent timing analysis and pulse-phase resolved spectroscopy for the two outbursts observed from the transient anomalous X-ray pulsar CXOU J164710.2−455216. For the first outburst we used 11 Chandra and XMM–Newton observations between 2006 September and 2009 August, the longest baseline yet for this source. We obtain a coherent timing solution with P = 10.61065583(4) s, Ṗ = 9.72(1) × 10−13 s s−1 and P̈ = –1.05(5) × 10−20 s s−2. Under the standard assumptions this implies a surface dipolar magnetic field of ∼1014 G, confirming this source as a standard B magnetar. We also study the evolution of the pulse profile (shape, intensity and pulsed fraction) as a function of time and energy. Using the phase-coherent timing solution we perform a phase-resolved spectroscopy analysis, following the spectral evolution of pulse-phase features, which hints at the physical processes taking place on the star. The results are discussed from the perspective of magnetothermal evolution models and the untwisting magnetosphere model. Finally, we present similar analysis for the second, less intense, 2011 outburst. For the timing analysis we used Swift data together with 2 XMM–Newton and Chandra pointings. The results inferred for both outbursts are compared and briefly discussed in a more general framework.

Optical phase coherent timing of the Crab nebula pulsar with Iqueye at the ESO New Technology Telescope

The Crab nebula pulsar was observed in 2009 January and December with a novel very fast optical photon counter, Iqueye, mounted at the ESO 3.5 m New Technology Telescope. Thanks to the exquisite quality of the Iqueye data, we computed accurate phase coherent timing solutions for the two observing runs and over the entire year 2009. Our statistical uncertainty on the determination of the phase of the main pulse and the rotational period of the pulsar for short (a few days) time intervals are $\approx 1 \, \mu$s and ~0.5 ps, respectively. Comparison with the Jodrell Bank radio ephemerides shows that the optical pulse leads the radio one by ~240 $\mu$s in January and ~160 $\mu$s in December, in agreement with a number of other measurements performed after 1996. A third-order polynomial fit adequately describes the spin-down for the 2009 January plus December optical observations. The phase noise is consistent with being Gaussian distributed with a dispersion $\sigma$ of $\approx 15 \, \mu$s in most observations, in agreement with theoretical expectations for photon noise-induced phase variability.

A coherent acceleration search of the Parkes multibeam pulsar survey – techniques and the discovery and timing of 16 pulsars

A fully coherent acceleration search algorithm has been applied to the Parkes multibeam pulsar survey of the Galactic plane to search for previously undiscovered relativistic binary pulsars. The search has resulted in the discovery of 16 pulsars including a binary millisecond pulsar and an intermittent pulsar. Despite a number of promising candidates there have been no new discoveries of relativistic binary pulsars. Here, we detail the acceleration search performed in our analysis and present coherent timing solutions for each of the pulsars discovered. In light of the lack of discoveries of relativistic binary pulsars, we also discuss the technique of acceleration searching and its effectiveness in finding these systems.

Rotating Radio Transients and their place among pulsars

AbstractSix years ago, the discovery of Rotating Radio Transients (RRATs) marked what appeared to be a new type of sparsely-emitting pulsar. Since 2006, more than 70 of these objects have been discovered in single-pulse searches of archival and new surveys. With a continual inflow of new information about the RRAT population in the form of new discoveries, multi-frequency follow ups, coherent timing solutions, and pulse rate statistics, a view is beginning to form of the place in the pulsar population RRATs hold. Here we review the properties of neutron stars discovered through single pulse searches. We first seek to clarify the definition of the term RRAT, emphasising that “the RRAT population” encompasses several phenomenologies. A large subset of RRATs appears to represent the tail of an extended distribution of pulsar nulling fractions and activity cycles; these objects present several key open questions remaining in this field.

Rotating Radio Transients: new discoveries, timing solutions and musings

We describe observations of Rotating RAdio Transients (RRATs) that were discovered in a re-analysis of the Parkes Multi-beam Pulsar Survey (PMPS). The sources have now been monitored for sufficiently long to obtain seven new coherent timing solutions, to make a total of 14 now known. Furthermore we announce the discovery of 7 new transient sources, one of which may be extragalactic in origin (with $z\sim0.1$) and would then be a second example of the so-called `Lorimer burst'. The timing solutions allow us to infer neutron star characteristics such as energy-loss rate, magnetic field strength and evolutionary timescales, as well as facilitating multi-wavelength followup by providing accurate astrometry. All of this enables us to consider the question of whether or not RRATs are in any way special, i.e. a distinct and separate population of neutron stars, as has been previously suggested. We see no reason to consider 'RRAT' as anything other than a detection label, the subject of a selection effect in the parameter space searched. However, single-pulse searches can be utilised to great effect to identify pulsars difficult, or impossible, to find by other means, in particular those with long-periods (half of the PMPS RRATs have periods greater than 4 seconds), high-magnetic field strengths ($B\gtrsim 10^{13}$ G) and pulsars approaching the 'death valley'. The detailed nulling properties of such pulsars are unknown but the mounting evidence suggests a broad range of behaviour in the pulsar population. The group of RRATs fit in to the picture where pulsar magnetospheres switch between stable configurations.

Updated phase coherent timing solution of the isolated neutron star RX J0720.4–3125 using recent XMM-Newton and Chandra observations

Since the last phase coherent timing solution of the nearby radio-quiet isolated neutron star RX J0720.4-3125 six new XMM-Newton and three Chandra observations were carried out. The phase coherent timing solutions from previous authors were performed without restricting to a fixed energy band. However, we recently showed that the phase residuals are energy dependent, and thus phase coherent solutions must be computed referring always to the same energy band. We updated the phase coherent timing solution for RX J0720.4-3125 by including the recent XMM-Newton EPIC-pn, MOS1, MOS2 and Chandra ACIS data in the energy range 400-1000~eV. Altogether these observations cover a time span of almost 10~yrs. A further timing solution was obtained including the ROSAT pointed data. In this case, observations cover a time span of $\approx$16~yrs. To illustrate the timing differences between the soft band (120-400~eV) and the hard band (400-1000~eV) a timing solution for the soft band is also presented and the results are verified using a $\mathrm{Z_{n}^{2}}$ test. In contrast to previous work, we obtain almost identical solutions whether or not we include the ROSAT or Chandra data. Thanks to the restriction to the hard band, the data points from EPIC-pn are in better agreement with those from MOS1, MOS2 and Chandra than in previous works. In general the phase residuals are still large and vary with time. In particular, the latest XMM-Newton and Chandra data show that the phase residuals have attained relatively large and negative values.

A Coherent Timing Solution for the Nearby Isolated Neutron Star RX J1308.6+2127/RBS 1223

We present a phase-connected timing solution for the nearby isolated neutron star RX J1308.6+2127 (RBS 1223). From dedicated Chandra observations as well as archival Chandra and XMM-Newton data spanning a period of 5 years, we demonstrate that the 10.31 s pulsations are slowing down steadily at a rate of = 1.120(3) × 10-13 s s-1. Under the assumption that this is due to magnetic dipole torques, we infer a characteristic age of 1.5 Myr and a magnetic field strength of 3.4 × 1013 G. As with RX J0720.4-3125, the only other radio-quiet thermally emitting isolated neutron star for which a timing solution has been derived, the field strength is roughly consistent with what was inferred earlier from the presence of a strong absorption feature in its X-ray spectrum. Furthermore, for both sources the characteristic age is in excess of the cooling age inferred from standard cooling models. The sources differ, however, in their timing noise: while RX J0720.4-3125 showed considerable timing noise, RX J1308.6+2127 appears relatively stable.

A Coherent Timing Solution for the Nearby Isolated Neutron Star RX J0720.4-3125

We present the results of a dedicated effort to measure the spin-down rate of the nearby isolated neutron star RX J0720.4-3125. Comparing arrival times of the 8.39 s pulsations for data from Chandra, we derive an unambiguous timing solution for RX J0720.4-3125 that is accurate to 5 years. Adding data from XMM and ROSAT, the final solution yields = (6.98 ± 0.02) × 10-14 s s-1; for dipole spin-down, this implies a characteristic age of 2 Myr and a magnetic field strength of 2.4 × 1013 G. The phase residuals are somewhat larger than those for purely regular spin-down, but they do not show conclusive evidence for higher order terms or a glitch. From our timing solution as well as recent X-ray spectroscopy, we concur with recent suggestions that RX J0720.4-3125 is most likely an off-beam radio pulsar with a moderately high magnetic field.

Observations of 20 Millisecond Pulsars in 47 Tucanae at 20 Centimeters

We have used a new observing system on the Parkes radio telescope to carry out a series of pulsar observations of the globular cluster 47 Tucanae at 20 cm wavelength. We detected all 11 previously known pulsars and have discovered nine others, all of which are millisecond pulsars in binary systems. We have searched the data for relatively short orbital period systems, and found one pulsar with an orbital period of 96 minutes, the shortest of any known radio pulsar. The increased rate of detections with the new system resulted in improved estimates of the flux density of the previously known pulsars, determination of the orbital parameters of one of them, and a coherent timing solution for another one. Five of the pulsars now known in 47 Tucanae have orbital periods of a few hours and implied companion masses of only ~0.03 M☉. Two of these are eclipsed at some orbital phases, while three are seen at all phases at 20 cm but not always at lower frequencies. Four and possibly six of the other binary systems have longer orbital periods and companion masses ~0.2 M☉, with at least two of them having relatively large orbital eccentricities. All 20 pulsars have rotation periods in the range 2-8 ms.