Application of ‘CLEAN’ in the power spectral analysis of non-uniformly sampled pulsar timing data

Spectral analysis of the residual pulsearrival times of pulsars is a useful tool in understanding the nature of the underlying processes that may be responsible for the timing noise observed from pulsars. Power spectra of pulsar timing residuals may be described by one or a combination of powerlaws. As these spectra are expected to be very steep, it is important to ensure a high dynamic range in the estimation of the spectrum. This is difficult in practice since one is, in general, dealing with timing measurements made at unevenly placed epochs. In this paper, we present a technique based on, ‘CLEAN’ to obtain high dynamic range spectra from unevenly sampled data. We compare the performance of this technique with other techniques including some that were used earlier for estimation of power spectra of pulsar timing residuals.     

New timing parameters and positions for 16 southern radio pulsars

We report timing observations on 16 pulsars obtained using the Ooty Radio Telescope operating at 327 MHz. Using these observations, spanning over a year, we have obtained values of period derivatives for the first time and refined the estimates of the periods for all pulsars in this sample. We also report improved positions for all cases where the earlier position estimates had an uncertainty of a few arcmin.     

Using XMM-Newton to measure the spectrum of the Vela pulsar and its phase variation

We report on our analysis of two XMM-Newton observations of the Vela pulsar performed in December 2000 (total exposure time: 96.5 ks). We succeeded in resolving the pulsar spectrum from the surrounding bright nebular emission taking advantage both of the accurate calibration of the EPIC point spread function and of the Chandra/HRC surface brightness map of the nebula. This made it possible to assess to pulsar spectral shape disentangling its thermal and non-thermal components. Exploiting the photon harvest, we have also been able to perform a phase-resolved study of the pulsar emission.      

Equation of state of neutron star cores and spin down of isolated pulsars

We study possible impact of a softening of the equation of state by a phase transition, or appearance of hyperons, on the spin evolution of isolated pulsars. Numerical simulations are performed using exact 2-D simulations in general relativity. The equation of state of dense matter at supranuclear densities is poorly known. Therefore, the accent is put on the general correlations between evolution and equation of state, and mathematical strictness. General conjectures referring to the structure of the one-parameter families of stationary configurations are formulated. The interplay of the back bending phenomenon and stability with respect to axisymmetric perturbations is described. Changes of pulsar parameters in a corequake following instability are discussed, for a broad choice of phase transitions predicted by different theories of dense matter. The energy release in a corequake, at a given initial pressure, is shown to be independent of the angular momentum of collapsing configuration. This result holds for various types of phases transition, with and without metastability. We critically review observations of pulsars that could be relevant for the detection of the signatures of the phase transition in neutron star cores. This work was partially supported by the Polish MNiI Grant no. 1P03D-008-27.     

The Parkes Multibeam Pulsar Survey – VI. Discovery and timing of 142 pulsars and a Galactic population analysis

We present the discovery and follow-up observations of 142 pulsars found in the Parkes 20-cm multibeam pulsar survey of the Galactic plane. These new discoveries bring the total number of pulsars found by the survey to 742. In addition to tabulating spin and astrometric parameters, along with pulse width and flux density information, we present orbital characteristics for 13 binary pulsars which form part of the new sample. Combining these results from another recent Parkes multibeam survey at high Galactic latitudes, we have a sample of 1008 normal pulsars which we use to carry out a determination of their Galactic distribution and birth rate. We infer a total Galactic population of  30 000 ± 1100  potentially detectable pulsars (i.e. those beaming towards us) having 1.4-GHz luminosities above 0.1 mJy kpc². Adopting the Tauris & Manchester beaming model, this translates to a total of  155 000 ± 6000  active radio pulsars in the Galaxy above this luminosity limit. Using a pulsar current analysis, we derive the birth rate of this population to be  1.4 ± 0.2  pulsars per century. An important conclusion from our work is that the inferred radial density function of pulsars depends strongly on the assumed distribution of free electrons in the Galaxy. As a result, any analyses using the most recent electron model of Cordes & Lazio predict a dearth of pulsars in the inner Galaxy. We show that this model can also bias the inferred pulsar scaleheight with respect to the Galactic plane. Combining our results with other Parkes multibeam surveys we find that the population is best described by an exponential distribution with a scaleheight of 330 pc. Surveys underway at Parkes and Arecibo are expected to improve the knowledge of the radial distribution outside the solar circle, and to discover several hundred new pulsars in the inner Galaxy.     

The broad-band noise characteristics of selected anomalous X-ray pulsars and soft gamma-ray repeaters

We present the broad-band noise structure of selected anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) in the 2–60 keV energy band. We have analysed Rossi X-Ray Timing Explorer Proportional Counter Array archival light curves for four AXPs and one SGR. We detect that the persistent emission of these sources shows band-limited noise at low frequencies in the range 0.005–0.05 Hz varying from 2.5 to 70 per cent integrated rms in times of prolonged quiescence and following outbursts. We discovered band-limited red noise in 1E 2259+586 only for ∼2 yr after its major 2002 outburst. The system shows no broad-band noise otherwise. Although this rise in noise in 1E 2259+586 occurred following an outburst which included a rotational glitch, the other glitching AXPs showed no obvious change in broad-band noise, thus it does not seem that this noise is correlated with glitches. The only source that showed significant variation in broad-band noise was 1E   1048.1−5937  , where the noise gradually rose for 1.95 yr at a rate of ∼3.6 per cent per year. For this source the increases in broad-band noise was not correlated with the large increases in persistent and pulsed flux, or its two short SGR-like bursts. This rise in noise did commence after a long burst, however, given the sparsity of this event, and the possibility that similar bursts went unnoticed the trigger for the rise is noise in 1E   1048.1−5937  is not as clear as for 1E 2259+586. The other three sources indicate a persistent band-limited noise at low levels in comparison.     

The symmetry energy and incompressibility constrained by the observations of glitching pulsars

We investigate the masses of glitching pulsars in order to constrain their equation of state(EOS). The observations of glitches(sudden jumps in rotational frequency) may provide information on the interior physics of neutron stars. With the assumption that glitches are triggered by superfluid neutrons, the masses of glitching neutron stars can be estimated using observations of maximum glitches.Together with the observations of thermal emission from glitching pulsars Vela and J1709–4429, the slope of symmetry energy and incompressibility of nuclear matter at saturation density can be constrained.The slope of symmetry energy L should be larger than 67 MeV while the lower limit of incompressibility for symmetric nuclear matter K_0 is 215 MeV. We also obtain a relationship between L and K_0:6.173 MeV + 0.283 K_0≤ L ≤ 7.729 MeV + 0.291 K_0. The restricted EOSs are consistent with the observations of 2-solar-mass neutron stars and gravitational waves from a binary neutron star inspiral.     

Timing characteristics of the hard X-ray emission from bright X-ray pulsars based on INTEGRAL data

We review the results of a timing analysis of the observations for ten bright X-ray pulsars (with fluxes >100 mCrab in the 20–100 keV energy band) that fell within the INTEGRAL field of view from 2003 to 2007. The dependence of the pulse profile on the energy and intrinsic source luminosity has been investigated; particular attention has been paid to searching for changes in the pulse profile near the cyclotron frequency. The dependence of the pulsed fraction for X-ray pulsars on their luminosity and energy band has been studied in detail for the first time.     

On the aberration–retardation effects in pulsars

The magnetospheric locations of pulsar radio emission region are not well known. The actual form of the so-called radius-to-frequency mapping should be reflected in the aberration–retardation (A/R) effects that shift and/or delay the photons depending on the emission height in the magnetosphere. Recent studies suggest that in a handful of pulsars the A/R effect can be discerned with respect to the peak of the central core emission region. To verify these effects in an ensemble of pulsars, we launched a project analysing multifrequency total intensity pulsar profiles obtained from the new observations from the Giant Meterwave Radio Telescope (GMRT), Arecibo Observatory (AO) and archival European Pulsar Network (EPN) data. For all these profiles, we measure the shift of the outer cone components with respect to the core component, which is necessary for establishing the A/R effect. Within our sample of 23 pulsars, seven show the A/R effects, 12 of them (doubtful cases) show a tendency towards this effect, while the remaining four are obvious counterexamples. The counterexamples and doubtful cases may arise from uncertainties in the determination of the location of the meridional plane and/or the core emission component. Hence, it appears that the A/R effects are likely to operate in most pulsars from our sample. We conclude that in cases where those effects are present the core emission has to originate below the conal emission region.     

Slow glitches in the pulsar B1822-09

The pulsar B1822-09 (J1825-0935) experienced a series of five unusual, slow glitches over the 1995–2004 interval. The results of further study of this unusual glitch phenomenon are presented. It is also reported the detection a new glitch of typical signature that occurred in the pulsar period in 2006 January.      

Impact of neutron star oscillations on the accelerating electric field in the polar cap of pulsar

Pulsar “standard model”, that considers a pulsar as a rotating magnetized conducting sphere surrounded by plasma, is generalized to the case of oscillating star. We developed an algorithm for calculation of the Goldreich-Julian charge density for this case. We consider distortion of the accelerating zone in the polar cap of pulsar by neutron star oscillations. It is shown that for oscillation modes with high harmonic numbers (l,m) changes in the Goldreich-Julian charge density caused by pulsations of neutron star could lead to significant altering of an accelerating electric field in the polar cap of pulsar. In the moderately optimistic scenario, that assumes excitation of the neutron star oscillations by glitches, it could be possible to detect altering of the pulsar radioemission due to modulation of the accelerating field. This work was partially supported by RFBR grant 04-02-16720, and by the grants N.Sh.-5218.2006.2 and RNP-2.1.1.5940.     

Simulations of glitches in isolated pulsars

Many radio pulsars exhibit glitches wherein the star's spin rate increases fractionally by ∼10⁻¹⁰–10⁻⁶. Glitches are ascribed to variable coupling between the neutron star crust and its superfluid interior. With the aim of distinguishing among different theoretical explanations for the glitch phenomenon, we study the response of a neutron star to two types of perturbations to the vortex array that exists in the superfluid interior: (1) thermal motion of vortices pinned to inner crust nuclei, initiated by sudden heating of the crust, (e.g., a starquake), and (2) mechanical motion of vortices (e.g., from crust cracking by superfluid stresses). Both mechanisms produce acceptable fits to glitch observations in four pulsars, with the exception of the 1989 glitch in the Crab pulsar, which is best fitted by the thermal excitation model. The two models make different predictions for the generation of internal heat and subsequent enhancement of surface emission. The mechanical glitch model predicts a negligible temperature increase. For a pure and highly conductive crust, the thermal glitch model predicts a surface temperature increase of as much as ∼2 per cent, occurring several weeks after the glitch. If the thermal conductivity of the crust is lowered by a high concentration of impurities, however, the surface temperature increases by ∼10 per cent about a decade after a thermal glitch. A thermal glitch in an impure crust is consistent with the surface emission limits following the 2000 January glitch in the Vela pulsar. Future surface emission measurements coordinated with radio observations will constrain glitch mechanisms and the conductivity of the crust.     

The variation of the angle between the magnetic and rotation axes in a collapsing star

We derive the formula for the variation of the magnetic axis in a collapsing oblique rotator. Our results of calculation show that the case of the magnetic axis coincident with the rotation axis is unstable, and the angle between the two axes will grow on any small perturbation.     

The population of pulsars with interpulses and the implications for beam evolution

The observed fraction of pulsars with interpulses, their period distribution and the observed pulse width versus pulse period correlation are shown to be inconsistent with a model in which the angle α between the magnetic axis and the rotation axis is random. This conclusion appears to be unavoidable, even when non-circular beams are considered. Allowing the magnetic axis to align from a random distribution at birth with a time-scale of  ∼7 × 10⁷ yr  can, however, explain those observations well. The time-scale derived is consistent with that obtained via independent methods. The probability that a pulsar beam intersects the line of sight is a function of the angle α and therefore beam evolution has important consequences for evolutionary models and for estimations of the total number of neutron stars. The validity of the standard formula for the spin-down rate, which is independent of α, appears to be questionable.     

Vortex–interface interactions and generation of glitches in pulsars

We show that the crust–core interface in neutron stars acts as a potential barrier to the peripheral neutron vortices approaching the interface in the model in which these are coupled to the proton vortex clusters. This elementary barrier arises because of the interaction of vortex magnetic flux with the Meissner currents set up by the crustal magnetic field at the interface. The dominant part of the force is derived from the cluster–interface interaction. As a result of the stopping of the continuous neutron vortex current through the interface, angular momentum is stored in the superfluid layers in the vicinity of the crust–core interface during the interglitch period. Discontinuous annihilation of proton vortices at the boundary restores the neutron vortex current and spins up the observable crust on short time-scales, leading to a glitch in the spin characteristics of a pulsar.