Vela Pulsar Research Articles

Suzaku Detection of Diffuse Hard X-Ray Emission outside Vela X

Abstract Vela X is a large, 3 $^{\circ}$$\times$ 2 $^{\circ}$ , radio-emitting pulsar wind nebula (PWN) powered by the Vela pulsar in the Vela supernova remnant. Using four Suzaku/XIS observations pointed just outside Vela X, we have found hard X-ray emission extending throughout the fields of view. The hard X-ray spectra are well represented by a power-law. The photon index was measured to be constant at $\Gamma$$\sim$ 2.4, similar to that of the southern outer part of Vela X. The power-law flux decreases with increasing distance from the pulsar. These properties lead us to propose that the hard X-ray emission is associated with the Vela PWN. The larger X-ray extension found in this work strongly suggests that distinct populations of relativistic electrons form the X-ray PWN and Vela X, as was recently inferred from multiwavelength spectral modeling of Vela X.

EXTENDED HARD X-RAY EMISSION FROM THE VELA PULSAR WIND NEBULA

The nebula powered by the Vela pulsar is one of the best examples of an evolved pulsar wind nebula, allowing to access the particle injection history and the interaction with the supernova ejecta. We report on the INTEGRAL discovery of extended emission above 18 keV from the Vela nebula. The northern side has no known counterparts and it appears larger and more significant than the southern one, which is in turn partially coincident with the cocoon, the soft X-ray and TeV filament towards the center of the remnant. We also present the spectrum of the Vela nebula in the 18-400 keV energy range as measured by IBIS/ISGRI and SPI on board the INTEGRAL satellite. The apparent discrepancy between IBIS/ISGRI, SPI, and previous measurements is understood in terms of point spread function, supporting the hypothesis of a nebula more diffuse than previously thought. A break at ~25 keV is found in the spectrum within 6' from the pulsar after including the Suzaku XIS data. Interpreted as a cooling break, this points out that the inner nebula is composed by electrons injected in the last ~2000 years. Broad-band modeling also implies a magnetic field higher than 10 microGauss in this region. Finally, we discuss the nature of the northern emission, which might be due to fresh particles injected after the passage of the reverse shock.

PULSAR WIND NEBULAE WITH THICK TOROIDAL STRUCTURE

We investigate a class of pulsar wind nebulae that show synchrotron emission from a thick toroidal structure. The best studied such object is the small radio and X-ray nebula around the Vela pulsar, which can be interpreted as the result of interaction of a mildly supersonic inward flow with the recent pulsar wind. Such a flow near the center of a supernova remnant can be produced in a transient phase when the reverse shock reaches the center of the remnant. Other nebulae with a thick toroidal structure are G106.6+2.9 and G76.9+1.0. Their structure contrasts with young pulsar nebulae like the Crab Nebula and 3C 38, which show a more chaotic, filamentary structure in the synchrotron emission. In both situations, a torus-jet structure is present where the pulsar wind passes through a termination shock, indicating the flow is initially toroidal. We suggest that the difference is due to the Rayleigh-Taylor instability that operates when the outer boundary of the nebula is accelerating into freely expanding supernova ejecta. The instability gives rise to mixing in the Crab and related objects, but is not present in the nebulae with thick toroidal regions.

DISCOVERY OF AN ENERGETIC PULSAR ASSOCIATED WITH SNR G76.9+1.0

We report the discovery of PSR J2022+3842, a 24 ms radio and X-ray pulsar in the supernova remnant G76.9+1.0, in observations with the Chandra X-ray telescope, the Robert C. Byrd Green Bank Radio Telescope, and the Rossi X-ray Timing Explorer (RXTE). The pulsar's spin-down rate implies a rotation-powered luminosity Edot = 1.2 x 10^{38} erg/s, a surface dipole magnetic field strength B_s = 1.0 x 10^{12} G, and a characteristic age of 8.9 kyr. PSR J2022+3842 is thus the second-most energetic Galactic pulsar known, after the Crab, as well as the most rapidly-rotating young, radio-bright pulsar known. The radio pulsations are highly dispersed and broadened by interstellar scattering, and we find that a large (delta-f / f ~= 1.9 x 10^{-6}) spin glitch must have occurred between our discovery and confirmation observations. The X-ray pulses are narrow (0.06 cycles FWHM) and visible up to 20 keV, consistent with magnetospheric emission from a rotation-powered pulsar. The Chandra X-ray image identifies the pulsar with a hard, unresolved source at the midpoint of the double-lobed radio morphology of SNR G76.9+1.0 and embedded within faint, compact X-ray nebulosity. The spatial relationship of the X-ray and radio emissions is remarkably similar to extended structure seen around the Vela pulsar. The combined Chandra and RXTE pulsar spectrum is well-fitted by an absorbed power-law model with column density N_H = (1.7\pm0.3) x 10^{22} cm^{-2} and photon index Gamma = 1.0\pm0.2; it implies that the Chandra point-source flux is virtually 100% pulsed. For a distance of 10 kpc, the X-ray luminosity of PSR J2022+3842 is L_X(2-10 keV) = 7.0 x 10^{33} erg s^{-1}. Despite being extraordinarily energetic, PSR J2022+3842 lacks a bright X-ray wind nebula and has an unusually low conversion efficiency of spin-down power to X-ray luminosity, L_X/Edot = 5.9 x 10^{-5}.

THE BRIGHTEST GAMMA-RAY FLARING BLAZAR IN THE SKY: AGILE AND MULTI-WAVELENGTH OBSERVATIONS OF 3C 454.3 DURING 2010 NOVEMBER

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

FERMI GAMMA-RAY SPACE TELESCOPE OBSERVATIONS OF THE GAMMA-RAY OUTBURST FROM 3C454.3 IN NOVEMBER 2010

The flat-spectrum radio quasar 3C 454.3 underwent an extraordinary 5-day gamma-ray outburst in November 2010 where the daily flux measured with the Fermi Large Area Telescope (LAT) at photon energies E>100 MeV reached (66+/-2) x 10^-6 ph cm^-2 s^-1. This is a factor of 3 higher than its previous maximum flux recorded in December 2009 and ~5 times brighter than the Vela pulsar, which is normally the brightest source in the gamma-ray sky. The 3-hr peak flux was (85+/-5) x 10^-6 ph cm^-2 s^-1, corresponding to an apparent isotropic luminosity of 2.1+/-0.2 10^50 erg s^-1, the highest ever recorded for a blazar. In this paper, we investigate the features of this exceptional event in the gamma-ray band of the Fermi-LAT. In contrast to previous flares of the same source observed with the Fermi-LAT, clear spectral changes are observed during the flare.

A study of 315 glitches in the rotation of 102 pulsars

The rotation of more than 700 pulsars has been monitored using the 76-m Lovell Telescope at Jodrell Bank. Here we report on a new search for glitches in the observations, revealing 128 new glitches in the rotation of 63 pulsars. Combining these new data with those already published we present a database containing 315 glitches in 102 pulsars. The database was used to study the glitch activity among the pulsar population, finding that it peaks for pulsars with a characteristic age tau_c ~ 10kyr and decreases for longer values of tau_c, disappearing for objects with tau_c > 20Myr. The glitch activity is also smaller in the very young pulsars (tau_c <~ 1kyr). The cumulative effect of glitches, a collection of instantaneous spin up events, acts to reduce the regular long term spindown rate |nudot| of the star. The percentage of |nudot| reversed by glitch activity was found to vary between 0.5% and 1.6% for pulsars with spindown rates |nudot| between 10^(-14) and 3.2*10^(-11) Hz/s, decreasing to less than 0.01% at both higher and lower spindown rates. These ratios are interpreted in terms of the amount of superfluid involved in the generation of glitches. In this context the activity of the youngest pulsar studied, the Crab pulsar, may be explained by quake-like activity within the crust. Pulsars with low spindown rates seem to exhibit mostly small glitches, matching well the decrease of their crustal superfluid. Through the analysis of glitch sizes it was found that the particular glitching behaviour of PSR J0537-6910 and the Vela pulsar may be shared by most Vela-like pulsars. These objects present most of their glitches with characteristic frequency and frequency derivative jumps, occurring at regular intervals of time. Their behaviour is different from other glitching pulsars of similar characteristic age.

Nuclear constraints on the inner edge of neutron star crusts

We show that the widely used parabolic approximation to the Equation of State (EOS) of asymmetric nuclear matter leads systematically to significantly higher core-crust transition densities and pressures. Using an EOS for neutron-rich nuclear matter constrained by the isospin diffusion data from heavy-ion reactions in the same sub-saturation density range as the neutron star crust, the density and pressure at the inner edge separating the liquid core from the solid crust of neutron stars are determined to be 0.040 fm − 3 ≤ ρ t ≤ 0.065 fm − 3 and 0.01 MeV/fm 3 ≤ P t ≤ 0.26 MeV/fm 3 , respectively. Implications of these constraints on the Vela pulsar are discussed.

Evolving into Magnetars from Normal Pulsars with a Low Braking Index

Anomalous x-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are believed to be candidates for magnetars, and they are powered by the decay of ultra-strong magnetic fields of > 1014 G. From the modified spin-down relation of pulsars Ṗ ∝ P2−n, we find that the Vela pulsar would evolve into the classes of magnetars under some assumptions that pulsars lose their rotational energy only by magnetic dipole radiation and the braking index is a constant. Our rough calculation indicates that only pulsars with n ∼ 1.3 – 1.6 can evolve into magnetars. Pulsars like Vela with a low braking index may be the progenitors of AXPs and SGRs. Regarding the mechanism evolved into magnetars, we suggest that pulsars' surface magnetic field component may be increased by frequent glitches.

FERMILARGE AREA TELESCOPE OBSERVATIONS OF THE VELA PULSAR

The Vela pulsar is the brightest persistent source in the GeV sky and thus is the traditional first target for new gamma-ray observatories. We report here on initial Fermi Large Area Telescope observations during verification phase pointed exposure and early sky survey scanning. We have used the Vela signal to verify Fermi timing and angular resolution. The high quality pulse profile, with some 32,400 pulsed photons at E>0.03 GeV, shows new features, including pulse structure as fine as 0.3ms and a distinct third peak, which shifts in phase with energy. We examine the high energy behavior of the pulsed emission; initial spectra suggest a phase-averaged power law index of Gamma=1.51{+0.05/-0.04} with an exponential cut-off at E_c=2.9+/-0.1 GeV. Spectral fits with generalized cut-offs of the form e^{-(E/E_c)^b} require b<1, which is inconsistent with magnetic pair attenuation, and thus favor outer magnetosphere emission models. Finally, we report on upper limits to any unpulsed component, as might be associated with a surrounding synchrotron wind nebula (PWN).

Towards explanation of the X-ray-radio correlation in the Vela pulsar

Recent observations of the Vela pulsar have revealed a peculiar connection of its emission in the soft X-ray and radio ranges. We suggest the model of the radio pulse formation in the Vela pulsar, develop the theory of the radio photon reprocessing to high energies, and on this basis interpret the observed X-ray-radio connection. The processes of spontaneous and induced scattering of radio waves off the spiralling particles and their observational consequences are examined. The particles are assumed to acquire relativistic gyration energies due to resonant absorption of the radio emission in the outer magnetosphere of a pulsar. The spectral and angular distributions of the spontaneously scattered power are analysed and compared with the characteristics of the particle synchrotron emission. The consequences of intensity transfer from the radio beam to the background in the course of induced scattering are studied as well. It is demonstrated that the induced scattering can account for the basic features of the Vela's radio profile and its pulse-to-pulse fluctuations. In particular, it can explain a greater role of the leading component and its earlier arrival in stronger pulses. The studies of the radio photon reprocessing to high energies in application to the Vela pulsar show that the scattered and synchrotron spectra peak at ∼0.8 and 0.2 keV, respectively, with the corresponding luminosities of 10 29 and 10 31 erg s -1 . The reprocessing inside the light cylinder results in the profile component, which stretches from the radio pulse position to ∼30° earlier in pulse phase. The synchrotron emission beyond the light cylinder presumably contributes to the component that lags the radio pulse by ∼90°. Within the framework of our model, the observed X-ray-radio connection is explained in terms of the interplay between the processes of induced and spontaneous scattering of the radio pulse.

Pulse Profiles, Spectra, and Polarization Characteristics of Nonthermal Emissions from the Crab‐like Pulsars

We discuss nonthermal emission mechanisms of the Crab-like pulsars with both a two-dimensional electrodynamic study and a three-dimensional model. We investigate the emission process in the outer gap accelerator. In the two-dimensional electrodynamic study, we solve the Poisson equation of the accelerating electric field in the outer gap and the equation of motion of the primary particles with the synchrotron and the curvature radiation processes and the pair-creation process. We show a solved gap structure that produces a γ-ray spectrum consistent with EGRET observations. Based on the two-dimensional model, we construct a three-dimensional emission model to calculate the synchrotron and the inverse Compton processes of the secondary pairs produced outside the outer gap. We calculate the pulse profiles, the phase-resolved spectra, and the polarization characteristics in optical through γ-ray bands for comparison with the observation of the Crab pulsar and PSR B0540-69. For the Crab pulsar, we find that the outer gap geometry extending from near the stellar surface to near the light cylinder produces a complex morphology change of the pulse profiles as a function of the photon energy. This predicted morphology change is quite similar to that of the observations. The calculated phase-resolved spectra are consistent with the data from the optical to the γ-ray bands. We demonstrate that the 10%-20% of the polarization degree in the optical emissions from the Crab pulsar and the Vela pulsar is explained by the synchrotron emissions from the particle gyration motion. For PSR B0540-69, the observed pulse profile with a single broad pulse is reproduced for an emission region thicker and an inclination angle between the rotational axis and the magnetic axis smaller than the Crab pulsar.

Angular velocity glitches of the Vela pulsar in a model with “centrifugal buoyancy” of vortices

The forces acting on the solid crust of a differentially rotating neutron star are examined when a nonuniform excess of chemical potential exists. The resultant of the external forces, a stress force, is expressed in terms of a centrifugal buoyancy force and the deformation of the star’s crust under the action of this force is calculated. It is shown that there is a region within the star where the resulting stresses lead to fracture of the crust when the difference in the angular velocities of the superfluid and normal components reaches a critical value. The “centrifugal buoyancy” mechanism for generating a glitch is used to estimate the parameters of glitches in the Vela pulsar.

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.

Gamma-Ray Emission from Pulsar Outer Magnetospheres

We investigate a stationary pair production cascade in the outer magnetosphere of an isolated, spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate gamma-rays via curvature and inverse-Compton processes. Some of such gamma-rays collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved together with the energy distribution of particles and gamma-rays at each point along the field lines. By solving the set of Maxwell and Boltzmann equations, we demonstrate that an external injection of charged particles at nearly Goldreich-Julian rate does not quench the gap but shifts its position and that the particle energy distribution cannot be described by a power-law. The injected particles are accelerated in the gap and escape from it with large Lorentz factors. We show that such escaping particles migrating outside of the gap contribute significantly to the gamma-ray luminosity for young pulsars and that the soft gamma-ray spectrum between 100 MeV and 3 GeV observed for the Vela pulsar can be explained by this component. We also discuss that the luminosity of the gamma-rays emitted by the escaping particles is naturally proportional to the square root of the spin-down luminosity.

Interaction of pulsar winds with interstellar medium: numerical simulation

A time-dependent numerical simulation of relativistic wind interaction with interstellar medium was performed. The winds are ejected from the magnetosphere of rotation-powered pulsars. The particle flux in the winds is isotropic while the energy flux is strongly anisotropic. Pure hydrodynamical winds and magnetized winds with magnetization in the range 3 × 10 −3 ‐10 −1 were modelled. The numerical solutions reproduce the most spectacular features observed in the central part of the plerions: toroidal structure and jet-like features. A comparison of the solutions with the observations confirms that magnetization of the wind from a pulsar in the Crab nebula is of the order of 3 × 10 −3 . The plasma flow in the nebula is variable. The variability has different time-scales. Non-periodic 2‐3 yr scale variations are superimposed on the quasi-periodical variations with period ∼40 yr. An increase of the wind’s magnetization results in a decrease of the size of the synchrotron nebula. The size of the Vela nebula can be reproduced provided that the wind from Vela is more highly magnetized then the wind from Crab. Nevertheless, the magnetization of the wind from the Vela pulsar is most likely below 0.1. Ke yw ords: MHD ‐ shock waves ‐ pulsars: general ‐ ISM: individual: Crab nebula ‐ jets and outflows ‐ supernova remnants.

Tachyonic synchrotron radiation from γ-ray pulsars

Superluminal radiation emitted by electrons orbiting in strong magnetic fields is investigated. We show that electrons gyrating in the surface fields of rotation-powered neutron stars can radiate superluminal quanta (tachyons) via the synchrotron mechanism. The tachyonic luminosity of γ-ray pulsars is inferred from COMPTEL and EGRET observations, and so is the magnetospheric electron population generating this radiation. In the surface fields, electromagnetic synchrotron radiation in the γ-ray band is suppressed by a quantum cutoff, but not so tachyonic γ-radiation. This provides an exceptional opportunity to search for tachyon radiation, unspoiled by electromagnetic emission. Estimates of the superluminal power radiated and the tachyonic count rates are obtained for each of the seven established γ-ray pulsars, the Crab and Vela pulsars, as well as PSR B1706−44, Geminga, PSR B1055−52, B1951+32, and B1509−58. Detection mechanisms such as tachyonic ionization and Compton scattering are analyzed with regard to superluminal γ-rays.

The Very Young Radio Pulsar J1357-6429

We report the discovery of a radio pulsar with a characteristic age of 7300 years, making it one of the 10 apparently youngest Galactic pulsars known. PSR J1357-6429, with a spin period of P = 166 ms and spin-down luminosity of 3.1e36 ergs/s, was detected during the Parkes multibeam survey of the Galactic plane. We have measured a large rotational glitch in this pulsar, with Delta P/P = -2.4e-6, similar in magnitude to those experienced occasionally by the Vela pulsar. At a nominal distance of only ~ 2.5 kpc, based on the measured free electron column density of 127 pc/cc and the electron distribution model of Cordes & Lazio, this may be, after the Crab, the nearest very young pulsar known. The pulsar is located near the radio supernova remnant candidate G309.8-2.6.

Improving the angular resolution of EGRET and new limits on supersymmetric dark matter near the galactic center

Using the EGRET data and an improved point source analysis, including an energy-dependent point spread function and an unbinned maximum likelihood technique, we have been able to place considerably lower limits on the γ-ray flux from the galactic center region. We also test this method on known sources, the Crab and Vela pulsars. In both cases, we find that our method improves the angular precision of EGRET data over the 3EG catalog. This new limit on γ-rays from the galactic center can be used to test models of annihilating supersymmetric dark matter and galactic halo profiles. We find that the present EGRET data can limit many supersymmetric models if the density of the galactic dark matter halo is cuspy or spiked toward the galactic center. We also discuss the ability of GLAST to test these models.

Observations of PSR J2021+3651 and its X‐Ray Pulsar Wind Nebula G75.2+0.1

We present results from X-ray and radio observations of the recently discovered young Vela-like pulsar PSR J2021+3651, which is coincident with the EGRET gamma-ray source GeV 2020+3658. A 19.0-ks Chandra ACIS-S observation has revealed a ~20'' x 10'' pulsar wind nebula that is reminiscent of the equatorial tori seen around some young pulsars, along with thermal emission from an embedded point source (kT = 0.15 +/- 0.02 keV). We name the nebula G75.2+0.1. Its spectrum is well fit by an absorbed power-law model with photon index 1.7 +/- 0.3, hydrogen column density nH = 7.8 +/- 1.7 x 10^21 cm^-2, and an unabsorbed 0.3-10.0 keV flux of 1.9 +/- 0.3 x 10^-12 erg cm^-2 s^-1. We have spatially fit G75.2+0.1 with a model that assumes a toroidal morphology, and from this we infer that the torus is highly inclined 83 deg +/- 1 deg to the line of sight. A 20.8-ks Chandra observation in continuous-clocking mode reveals a possible pulse detection, with a pulsed fraction of ~37% and an H-test probability of occuring by chance of 1.2 x 10^-4. Timing observations with the Arecibo radio telescope spanning two years show that PSR J2021+3651 glitched sometime between MJDs 52616 and 52645 with parameters delta(v)/v = (2.587 +/- 0.002) x 10^-6 and delta(dot(v))/v = (6.2 +/- 0.3) x 10^-3, similar to those of the largest glitches observed in the Vela pulsar. PSR J2021+3651 is heavily scattered (T_sc = 17.7 ms +/- 0.9 ms at 1 GHz) and exhibits a significant amount of timing noise.