Pulsar slow glitches in a solid quark star model

A series of five unusual slow glitches of the radio pulsar B1822–09 (PSR J1825–0935) was observed between 1995 and 2005. This is a phenomenon that is understood in a solid quark star model, and reasonable parameters for slow glitches are given in this paper. We propose that, because of increasing shear stress as the pulsar spins down, a slow glitch may occur, beginning with the collapse of a superficial layer of the quark star. This layer of material turns to viscous fluid at first, the viscosity of which helps to deplete the energy released from both the accumulated elastic energy and the gravitation potential. There is then a slow glitch. Numerical calculations show that the slow glitches that have been observed could be reproduced if the effective coefficient of viscosity is ∼10² cm² s⁻¹ and the initial velocity of the superficial layer is of the order of 10⁻¹⁰ cm s⁻¹ in the coordinate rotating frame of the star.     

Spatially resolved STIS spectra of the gravitationally lensed broad absorption line quasar APM08279+5255: the nature of component C and evidence for microlensing

PSR J1806−2125 is a pulsar discovered in the Parkes multibeam pulsar survey with a rotational period of 0.4 s and a characteristic age of 65 kyr. Between MJDs 51462 and 51894 this pulsar underwent an increase in rotational frequency of  Δ ν / ν ≈16×10^-6  . The magnitude of this glitch is ∼2.5 times greater than any previously observed in any pulsar and 16 times greater than the mean glitch size. This Letter gives the parameters of the glitch and compares its properties with those of previously observed events. The existence of such large and rare glitches offers new hope for attempts to observe thermal X-ray emission from the internal heat released following a glitch, and suggests that pulsars which previously have not been observed to glitch may do so on long time-scales .     

Observations of six glitches in PSR B1737−30

Six glitches have been recently observed in the rotational frequency of the young pulsar PSR B1737−30 (J1740−3015) using the 25-m Nanshan telescope of Urumqi Observatory. With a total of 20 glitches in 20 yr, it is one of the most frequently glitching pulsars of the ∼1750 known pulsars. Glitch amplitudes are very variable with fractional increases in rotation rate ranging from 10⁻⁹ to 10⁻⁶. Interglitch intervals are also very variable, but no relationship is observed between interval and the size of the preceding glitch. There is a persistent increase in     , opposite in sign to that expected from slowdown with a positive braking index, which may result from changes in the effective magnetic dipole moment of the star during the glitch.     

Kelvin–Helmholtz instability and circulation transfer at an isotropic–anisotropic superfluid interface in a neutron star

A recent laboratory experiment suggests that a Kelvin–Helmholtz (KH) instability at the interface between two superfluids – one rotating and anisotropic, the other stationary and isotropic – may trigger sudden spin-up of the stationary superfluid. This result suggests that a KH instability at the crust–core (  ¹ S₀–³ P₂  –superfluid) boundary of a neutron star may provide a trigger mechanism for pulsar glitches. We calculate the dispersion relation of the KH instability involving two different superfluids including the normal fluid components and their effects on stability, particularly entropy transport. We show that an entropy difference between the core and crust superfluids reduces the threshold differential shear velocity and threshold crust–core density ratio. We evaluate the wavelength of maximum growth of the instability for neutron star parameters and find the resultant circulation transfer to be within the range observed in pulsar glitches.     

Glitch monitoring in PSRs B1046−58 and B1737−30

Large glitches were recently observed in the spin rates of two pulsars, B1046−58 and B1737−30. The glitches were characterized by fractional increases in rotation rate of 0.77 and  1.44×10⁻⁶  respectively. PSR B1737−30 is the most frequently glitching pulsar and this is the largest glitch so far observed from it. Most of the jump in the spin-down rate accompanying these glitches decayed away on short time-scales of a few days. For PSR B1737−30, there appears to be a cumulative shift in spin-down rate resulting from its frequent glitches. This probably accounts for its braking index of  −4±2  suggested by the available data, while a value of  2.1±0.2  is obtained for B1046−58.     

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.     

The micro-glitch in PSR B1821−24: a case for a strange pulsar?

The single glitch observed in PSR B1821−24, a millisecond pulsar in M28, is unusual on two counts. First, the magnitude of this glitch is at least an order of magnitude smaller  (Δν/ν∼ 10⁻¹¹)  than the smallest glitch observed to date. Secondly, all other glitching pulsars have strong magnetic fields with   B ≳ 10¹¹ G  and are young, whereas PSR B1821−24 is an old recycled pulsar with a field strength of  2.25 × 10⁹ G  . We have earlier suggested that some of the recycled pulsars could actually be strange quark stars. In this work, we argue that the crustal properties of such a strange pulsar are just right to give rise to a glitch of this magnitude, explaining the scarcity of larger glitches in millisecond pulsars.     

Pulsar braking indices revisited

Using the standard equation for the slowdown of a neutron star, we derive a formula for the braking index via integration rather than the conventional differentiation. The new formula negates the need to measure the second time derivative of the rotation frequency, ν¨ . We show that the method gives similar braking indices for PSR B1509−58 and the Crab pulsar to those already in the literature. We point out that our method is useful for obtaining the braking indices of moderate-aged pulsars without the need for long, phase-connected timing solutions. We applied the method to 20 pulsars and discuss the implications of the results. We find that virtually all the derived braking indices are dominated by the effects of (unseen) glitches, the recovery from which corrupts the value of ν˙ . However, any real, large, positive braking index has implications for magnetic field decay and offers support to recent models of pulsar evolution.     

Flux predictions of high-energy neutrinos from pulsars

Young, rapidly rotating neutron stars could accelerate ions from their surfaces to energies of ∼1 PeV. If protons reach such energies, they will produce pions (with low probability) through resonant scattering with X-rays from the stellar surface. The pions subsequently decay to produce muon neutrinos. Here, we calculate the energy spectrum of muon neutrinos, and estimate the event rates at Earth. The spectrum consists of a sharp rise at ∼50 TeV, corresponding to the onset of the resonance, above which the flux drops with neutrino energy as  ε⁻²_ν  up to an upper energy cut-off that is determined by either kinematics or the maximum energy to which protons are accelerated. We estimate event rates as high as 10–100 km⁻² yr⁻¹ from some candidates, a flux that would be easily detected by IceCube. Lack of detection would allow constraints on the energetics of the poorly understood pulsar magnetosphere.     

Investigating ultra-long gravitational waves with measurements of pulsar rotational parameters

A method is suggested with which to explore the gravitational wave background (GWB) in the frequency range 10⁻¹²–10⁻⁸ Hz. This method is based on the precise measurements of pulsar rotational parameters: the influence of gravitational waves (GWs) in this frequency range will affect these parameters and therefore some conclusions about the energy density of the GWB can be made using analysis of the derivatives of pulsar rotational frequency. The calculated values of the second derivative from a number of pulsars limit the density of the GWB, Ω_gw, as follows:  Ω_gw < 2 × 10⁻⁶  . Also, the time series of the frequency ν of different pulsars in a pulsar array can be cross-correlated pairwise in the same manner as in anomalous residuals analysis, thus providing the possibility of GWB detection in the ultra-low-frequency range.     

Rotation speed and stellar axis inclination from p modes: how CoRoT would see other suns

In the context of future space-based asteroseismic missions, we have studied the problem of extracting the rotation speed and the rotation-axis inclination of solar-like stars from the expected data. We have focused on slow rotators (at most twice solar rotation speed), first, because they constitute the most difficult case and, secondly, because some of the Convection Rotation and planetary Transits ( CoRoT ) main targets are expected to have slow rotation rates. Our study of the likelihood function has shown a correlation between the estimates of inclination of the rotation axis i and the rotational splitting δν of the star. By using the parameters, i and  δν^⋆=δν sin  i   , we propose and discuss new fitting strategies. Monte Carlo simulations have shown that we can extract a mean splitting and the rotation-axis inclination down to solar rotation rates. However, at the solar rotation rate we are not able to correctly recover the angle i , although we are still able to measure a correct  δν^⋆  with a dispersion less than 40 nHz.     

Freely precessing neutron stars: model and observations

We present a model of a freely precessing neutron star, which is then compared against pulsar observations. The aim is to draw conclusions regarding the structure of the star, and to test theoretical ideas of crust–core coupling and superfluidity. We argue that, on theoretical grounds, it is likely that the core neutron superfluid does not participate in the free precession of the crust. We apply our model to the handful of proposed observations of free precession that have appeared in the literature. Assuming crust-only precession, we find that all but one of the observations are consistent with there being no pinned crustal superfluid at all; the maximum amount of pinned superfluid consistent with the observations is about 10⁻¹⁰ of the total stellar moment of inertia. However, the observations do not rule out the possibility that the crust and neutron superfluid core precess as a single unit. In this case the maximum amount of pinned superfluid consistent with the observations is about 10⁻⁸ of the total stellar moment of inertia. Both of these values are many orders of magnitude less than the 10⁻² value predicted by many theories of pulsar glitches. We conclude that superfluid pinning, at least as it affects free precession, needs to be reconsidered.     

Rotation and activity in the solar-type stars of NGC 2547

We present high-resolution spectroscopy of a sample of 24 solar-type stars in the young (15–40 Myr), open cluster NGC 2547. We use our spectra to confirm cluster membership in 23 of these stars, to determine projected equatorial velocities and chromospheric activity, and to search for the presence of accretion discs. We find examples of both fast ( v _e sin  i >50 km s⁻¹) and slow ( v _e sin  i <10 km s⁻¹) rotators, but no evidence for active accretion in any of the sample. The distribution of projected rotation velocities is indistinguishable from the slightly older IC 2391 and IC 2602 clusters, implying similar initial angular momentum distributions and circumstellar disc lifetimes. The presence of very slow rotators indicates either that long (10–40 Myr) disc lifetimes or internal differential rotation are needed, or that NGC 2547 (and IC 2391/2602) were born with more slowly rotating stars than are presently seen in even younger clusters and associations. The solar-type stars in NGC 2547 follow a similar rotation–activity relationship to that seen in older clusters. X-ray activity increases until a saturation level is reached for v _e sin  i >15–20 km s⁻¹. We are unable to explain why this saturation level, of log( L _x L _bol)≃−3.3, is a factor of 2 lower than in other clusters, but rule out anomalously slow rotation rates or uncertainties in X-ray flux calculations.     

Timing models for the long orbital period binary pulsar PSR B1259–63

The pulsar PSR B1259–63 is in a highly eccentric 3.4-yr orbit with the Be star SS 2883. Timing observations of this pulsar, made over a 7-yr period using the Parkes 64-m radio telescope, cover two periastron passages, in 1990 August and 1994 January. The timing data cannot be fitted by the normal pulsar and Keplerian binary parameters. A timing solution including a (non-precessing) Keplerian orbit and timing noise (represented as a polynomial of fifth order in time) provides a satisfactory fit to the data. However, because the Be star probably has a significant quadrupole moment, we prefer to interpret the data by a combination of timing noise, dominated by a cubic phase term, and ω^. and x ^. terms. We show that the ω^. and x ^. terms are likely to be a result of a precessing orbit caused by the quadrupole moment of the tilted companion star. We further rule out a number of possible physical effects which could contribute to the timing data of PSR B1259–63 on a measurable level.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

The possible role of r-modes in post-glitch relaxation of the Crab pulsar

The loss of angular momentum through gravitational radiation, driven by the excitation of r-modes, is considered for neutron stars that have rotation frequencies lower than the associated critical frequency. We find that for reasonable values of the initial amplitudes of such pulsation modes of the star, being excited at the event of a glitch in a pulsar, the total post-glitch losses correspond to a negligible fraction of the initial rise of the spin frequency in the case of Vela and older pulsars. However, for the Crab pulsar the same effect would result, within a few months, in a decrease in its spin frequency by an amount larger than its glitch-induced frequency increase. This could provide an explanation for the peculiar behaviour observed in the post-glitch relaxations of the Crab pulsar.     

Spectral variation in the X-ray pulsar GX 1+4 during a low-flux episode

The X-ray pulsar GX 1+4 was observed with the RXTE satellite for a total of 51 ks between 1996 July 19 and 21. During this period the flux decreased smoothly from an initial mean level of ≈6×10³⁶ erg s⁻¹ to a minimum of ≈4×10³⁵ erg s⁻¹ (2–60 keV, assuming a source distance of 10 kpc) before partially recovering towards the initial level at the end of the observation.
BATSE pulse timing measurements indicate that a torque reversal took place approximately 10 d after this observation. Both the mean pulse profile and the photon spectrum varied significantly. The observed variation in the source may provide important clues as to the mechanism of torque reversals.
The single best-fitting spectral model was based on a component originating from thermal photons with kT ₀≈1 keV Comptonized by a plasma of temperature kT ≈7 keV. Both the flux modulation with phase during the brightest interval and the evolution of the mean spectra over the course of the observation are consistent with variations in this model component; with, in addition, a doubling of the column density n _H contributing to the mean spectral change.
A strong flare of duration ≲50 s was observed during the interval of minimum flux, with the peak flux ≈20 times the mean level. Although beaming effects are likely to mask the true variation in M ˙ thought to give rise to the flare, the timing of a modest increase in flux prior to the flare is consistent with dual episodes of accretion resulting from successive orbits of a locally dense patch of matter in the accretion disc.     

Red supergiants in the LMC — II. Spectrophotometry and model atmospheres

We have analysed polarization data for a large number of isolated pulsars to investigate the evolution of pulsar radio beams. Assuming that a circular beam is directed along the axis of a dipolar magnetic field, we demonstrate that the distribution of magnetic inclination angles for the parent population of all pulsars is not flat but highly concentrated towards small inclination angles and that, consequently, the average beaming fraction is only ∼ 10 per cent. Furthermore, we find that there is a tendency for the beam axis to align with the rotational axis on a time-scale of ∼ 10⁷ yr. This has interesting consequences for statistical studies of the pulsar population. Finally, the luminosity of pulsars is shown to be independent of the impact parameter, which indicates that pulsar beams are sharp-edged and have a relatively flat integrated intensity distribution.     

Scintillation parameters for 49 pulsars

We report on multi-epoch, multifrequency observations of 64 pulsars with high spectral and time resolution. Scintillation parameters were obtained for 49 pulsars, including 13 millisecond pulsars. Scintillation speeds were derived for all 49, which doubles the number of pulsars with speeds measured in this way. There is excellent agreement between the scintillation speed and proper motion for the millisecond pulsars in our sample using the simple assumption of a mid-placed scattering screen. This indicates that the scaleheight of scattering electrons is similar to that of the dispersing electrons. In addition, we present observations of the Vela pulsar at 14 and 23 GHz, and show that the scintillation bandwidth scales as ν^3.93 over a factor of 100 in observing frequency. We show that for PSR J0742−2822, and perhaps PSR J0837−4135, the Gum nebula is responsible for the high level of turbulence along their lines of sight, contrary to previous indications. There is a significant correlation between the scintillation speeds and the product of the pulsar's period and period derivative for the 'normal' pulsars. However, we believe this to be caused by selection effects both in pulsar detection experiments and in the choice of pulsars used in scintillation studies.     

The gamma-ray conversion efficiency of rotation-powered pulsars

We use a new self-consistent model to derive the conversion efficiency from rotation power to γ-ray power for pulsars (ε_th). Our result indicates that ε_th∝τ^6/7 P ², where τ and P are the characteristic age and period of the pulsar, which shows that although the efficiency increases with the characteristic age of the pulsar, it also depends on the pulsar period. We test our model results with the survey of high-energy γ-rays from pulsars by EGRET. Our model not only successfully explains the efficiency of the confirmed γ-ray pulsars but also explains why the γ-ray efficiency of millisecond pulsars is so low.