Statistical studies of pulsar glitches

Shemar & Lyne have previously presented observations and an analysis of 32 glitches and their subsequent relaxations observed in a total of 15 pulsars. These data are brought together in this paper with those published by other authors. We show quantitatively how glitch activity decreases linearly with decreasing rate of slow-down. As indicated previously from studies of the Vela pulsar, the analysis suggests that 1.7 per cent of the moment of inertia of a typical neutron star is normally contained in pinned superfluid which releases its excess angular momentum at the time of a glitch. There is a broad range of glitch amplitude and there is a strong indication that pulsars with large magnetic fields suffer many small glitches while others show a smaller number of large glitches. Transient effects following glitches are very marked in young pulsars and decrease linearly with decreasing rate of slow-down, suggesting that the amount of loosely pinned superfluid decreases with age. We suggest that the low braking index of the Vela and Crab pulsars cannot be caused by a decreasing moment of inertia and should be attributed to step increases in the effective magnetic moment of the neutron star at the glitches.     

Self-similar pressure oscillations in neutron star envelopes as probes of neutron star structure

We study eigenmodes of acoustic oscillations of high multipolarity l ∼ 100–1000 and high frequency (∼100 kHz), localized in neutron star envelopes. We show that the oscillation problem is self-similar. Once the oscillation spectrum is calculated for a given equation of state (EOS) in the envelope and given stellar mass M and radius R , it can be rescaled to a star with any M and R (but the same EOS in the envelope). For l ≳ 300, the modes can be subdivided into the outer and inner ones. The outer modes are mainly localized in the outer envelope. The inner modes are mostly localized near the neutron drip point, being associated with the softening of the EOS after the neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and accreted matter and show that the spectra of the inner modes are essentially different. A detection and identification of high-frequency pressure modes would allow one to infer M and R and determine also the EOS in the envelope (accreted or ground state) providing a new, potentially powerful method to explore the main parameters and internal structure of neutron stars.     

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 .     

The evolutionary status of W Ursae Majoris-type systems

Well-determined physical parameters of 130 W Ursae Majoris (W UMa) systems were collected from the literature. Based on these data, the evolutionary status and dynamical evolution of W UMa systems are investigated. It is found that there is no evolutionary difference between W- and A-type systems in the   M – J   diagram, which is consistent with the results derived from the analysis of observed spectral type and of   M – R   and   M – L   diagrams of W UMa systems.   M – R   and   M – L   diagrams of W- and A-type systems indicate that a large amount of energy should be transferred from the more massive to the less massive component, so that they are not in thermal equilibrium and undergo thermal relaxation oscillation. Moreover, the distribution of angular momentum, together with the distribution of the mass ratio, suggests that the mass ratio of the observed W UMa systems decreases with decreasing total mass. This could be the result of the dynamical evolution of W UMa systems, which suffer angular momentum loss and mass loss as a result of the magnetic stellar wind. Consequently, the tidal instability forces these systems towards lower q values and finally to rapidly rotating single stars.     

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.     

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.     

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.     

Stellar velocity profiles and line strengths out to four effective radii in the early-type galaxies NGC 3379 and 821

We use the integral-field spectrograph SAURON to measure the stellar line-of-sight velocity distribution and absorption line strengths out to four effective radii ( R _e) in the early-type galaxies NGC 3379 and 821. With our newly developed observing technique, we can now probe these faint regions in galaxies that were previously not accessible with traditional long-slit spectroscopy. We make optimal use of the large field-of-view and high throughput of the spectrograph: by adding the signal of all ∼1400 lenslets into one spectrum, we obtain sufficient signal-to-noise in a few hours of observing time to reliably measure the absorption line kinematics and line strengths out to large radius.
We find that the line strength gradients previously observed within 1 R _e remain constant out to at least 4 R _e, which puts constraints on the merger histories of these galaxies. The stellar halo populations are old and metal poor. By constructing orbit-based Schwarzschild dynamical models, we find that dark matter is necessary to explain the observed kinematics in NGC 3379 and 821, with 30–50 per cent of the total matter being dark within 4 R _e. The radial anisotropy in our best-fitting halo models is less than in our models without halo, due to differences in orbital structure. The halo also has an effect on the  Mg  b – V _esc  relation: its slope is steeper when a dark matter halo is added to the model.     

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.     

New limits on the orbital parameters of 1E 1048.1–5937 and 1E 2259+586 from RXTE observations

We report on two Rossi X-ray Timing Explorer ( RXTE ) observations of the anomalous X-ray pulsars 1E 1048.1–5937 and 1E 2259+586. Both sources have continued their almost constant spin-down during 1995/96. We carried out a search for orbital Doppler shifts, in their observed spin frequencies, deriving stringent limits on the projected semi-axis. Unless these systems have unlikely small inclinations, main-sequence companions can be excluded. If 1E 1048.1–5937 and 1E 2259+586 are indeed binary systems, their companion stars must be either white dwarfs, or helium-burning stars with M ≲ 0.8 M⊙, possibly underfilling their Roche lobe.     

A search for giant pulses in Vela-like pulsars

We have carried out a survey for 'giant pulses' in six young, Vela-like pulsars. In no cases did we find single pulses with flux densities more than 10 times the mean flux density. However, in PSR  B1706–44  we have detected giant micro-pulses very similar to those seen in the Vela pulsar. In PSR  B1706–44  these giant micro-pulses appear on the trailing edge of the profile and have an intrinsic width of ∼1 ms. The cumulative probability distribution of their intensities is best described by a power law. If the power law continues to higher intensities, then  3.7×10⁶  rotations are required to obtain a pulse with 20× the mean pulse flux. This number is similar to the giant pulse rate in PSR B1937+21 and PSR  B1821–24  but significantly higher than that for the Crab.     

Narrow-band CCD photometry of giant H ii regions

We have obtained accurate CCD narrow-band Hβ and Hα photometry of giant H  ii regions (GEHRs) in M33, NGC 6822 and M101. Comparison with previous determinations of emission-line fluxes shows large discrepancies; their probable origins are discussed. Combining our new photometric data with global velocity dispersion ( σ ) derived from emission linewidths, we review the     relation. A re-analysis of the properties of the GEHRs included in our sample shows that age spread and the superposition of components in multiple regions introduce a considerable spread in the regression. Combining the information available in the literature regarding ages of the associated clusters, evolutionary footprints on the interstellar medium, and kinematical properties of the knots that build up the multiple GEHRs, we find that a subsample – which we refer to as young and single GEHRs – do follow a tight relation in the     plane.     

Kinematical studies of the low-mass X-ray binary GR Mus (XB 1254–690)

We present simultaneous high-resolution optical spectroscopy and X-ray data of the X-ray binary system GR Mus (XB 1254–690), obtained over a full range of orbital phases. The X-ray observations are used to re-establish the orbital ephemeris for this source. The optical data include the first spectroscopic detection of the donor star in this system through the use of the Doppler Tomography technique on the Bowen fluorescence blend (∼4630–4650 Å). In combination with an estimate for the orbital parameters of the compact object using the wings of the He  ii λ4686 emission line, dynamical mass constraints of  1.20 ≤ M _X /M_⊙≤ 2.64  for the neutron star and  0.45 ≤ M ₂/M_⊙≤ 0.85  for the companion are derived.     

M4–18: the planetary nebula and its WC10 central star

We present a detailed analysis of the planetary nebula M4–18 (G146.7+07.6) and its WC10-type Wolf–Rayet (WR) central star, based on high‐quality optical spectroscopy (WHT/UES, INT/IDS, WIYN/DensPak) and imaging ( HST /WFPC2). From a non-LTE model atmosphere analysis of the stellar spectrum, we derive T _eff=31 kK,     v _∞=160 km s⁻¹ and abundance number ratios of H/He<0.5, C/He=0.60 and O/He=0.10. These parameters are remarkably similar to those of He 2–113 ([WC10]). Assuming an identical stellar mass to that determined by De Marco et al. for He 2–113, we obtain a distance of 6.8 kpc to M4–18 [ E ( B−V )=0.55 mag from nebular and stellar techniques]. This implies that the planetary nebula of M4–18 has a dynamical age of ∼3100 yr, in contrast to ≥270 yr for He 2–113. This is supported by the much higher electron density of the latter. These observations may be reconciled with evolutionary predictions only if [WC]-type stars exhibit a range in stellar masses.
Photoionization modelling of M4–18 is carried out using our stellar WR flux distribution, together with blackbody and Kurucz energy distributions obtained from Zanstra analyses. We conclude that the ionizing energy distribution from the WR model provides the best consistency with the observed nebular properties, although discrepancies remain.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.     

A search for radio pulsars around low-mass white dwarfs

Low-mass white dwarfs can be produced either in low-mass X-ray binaries by stable mass transfer to a neutron star, or in a common envelope phase with a heavier white dwarf companion. We have searched eight low-mass white dwarf candidates recently identified in the Sloan Digital Sky Survey for radio pulsations from pulsar companions, using the Green Bank Telescope at 340 MHz. We have found no pulsations down to flux densities of 0.6–0.8 mJy kpc⁻² and conclude that a given low-mass helium-core white dwarf has a probability of  <0.18 ± 0.05  of being in a binary with a radio pulsar.     

Constraints on the mass and abundance of black holes in the Galactic halo: the high-mass limit

We establish constraints on the mass and abundance of black holes in the Galactic halo by determining their impact on globular clusters, which are conventionally considered to be little evolved. Using detailed Monte Carlo simulations and simple evolutionary models, we argue that black holes with masses M _bh≳(1–3)×10⁶ M_⊙ can comprise no more than a fraction f _bh≈0.17 of the total halo density at Galactocentric radius R ≈8 kpc. This bound arises from requiring stability of the cluster mass function. A more restrictive bound may be derived if we demand that the probability of destruction of any given, low-mass M _c≈(2.5–7.5)×10⁴ M_⊙] globular cluster not exceed 50 per cent; this bound is f _bh≲0.025–0.05 at R ≈8 kpc. This constraint improves those based on disc heating and dynamical friction arguments as well as current lensing results. At smaller radius the constraint on f _bh strengthens, while at larger radius an increased fraction of black holes is allowed.     

The evolution of the cluster X-ray scaling relations in the Wide Angle ROSAT Pointed Survey sample at 0.6 < z < 1.0

The X-ray properties of a sample of 11 high-redshift  (0.6 < z < 1.0)  clusters observed with Chandra and/or XMM–Newton are used to investigate the evolution of the cluster scaling relations. The observed evolution in the normalization of the   L – T , M – T , M _g– T   and M – L relations is consistent with simple self-similar predictions, in which the properties of clusters reflect the properties of the Universe at their redshift of observation. Under the assumption that the model of self-similar evolution is correct and that the local systems formed via a single spherical collapse, the high-redshift L – T relation is consistent with the high- z clusters having virialized at a significantly higher redshift than the local systems. The data are also consistent with the more realistic scenario of clusters forming via the continuous accretion of material.
The slope of the L – T relation at high redshift  ( B = 3.32 ± 0.37)  is consistent with the local relation, and significantly steeper than the self-similar prediction of   B = 2  . This suggests that the same non-gravitational processes are responsible for steepening the local and high- z relations, possibly occurring universally at   z ≳ 1  or in the early stages of the cluster formation, prior to their observation.
The properties of the intracluster medium at high redshift are found to be similar to those in the local Universe. The mean surface-brightness profile slope for the sample is  β= 0.66 ± 0.05  , the mean gas mass fractions within   R _{2500( z )}  and   R _{200( z )}  are  0.069 ± 0.012  and  0.11 ± 0.02  , respectively, and the mean metallicity of the sample is  0.28 ± 0.11 Z_⊙  .     

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.     

The Be/X-ray transient V0332+53: evidence for a tilt between the orbit and the equatorial plane?

We present optical and infrared observations of BQ Cam, the optical counterpart to the Be/X-ray transient system V0332+53. BQ Cam is shown to be an O8–9Ve star, which places V0332+53 at a distance of ∼7 kpc. H α spectroscopy and infrared photometry are used to discuss the evolution of the circumstellar envelope. Owing to the low inclination of the system, parameters are strongly constrained. We find strong evidence for a tilt of the orbital plane with respect to the circumstellar disc (presumably on the equatorial plane). Even though the periastron distance is only ≈10 R _*, during the present quiescent state the circumstellar disc does not extend to the distance of periastron passage. Under these conditions, X-ray emission is effectively prevented by centrifugal inhibition of accretion. The circumstellar disc is shown to be optically thick at optical and infrared wavelengths, which, together with its small size, is taken as an indication of tidal truncation.