The contribution of pulsars to diffuse gamma-rays in the Large Magellanic Cloud

We study the contribution of young pulsars, with characteristic ages of less than 10⁶ yr, to the diffuse γ-ray emission from the Large Magellanic Cloud (LMC). Based on the outer gap model for γ-ray emission proposed by Zhang & Cheng and pulsar properties in the LMC given by Hartmann, Brown & Schnepf, we simulate the properties of the young pulsars in the LMC. We show that γ-rays produced by the pulsars in the LMC may make an important contribution to the diffuse γ-rays in the LMC, especially in the high-energy range. We calculate the γ-ray energy spectrum of the pulsars in the LMC and show that the γ-ray component contributed by the pulsars to the diffuse γ-rays in the high-energy range (above ∼1 GeV) becomes dominant. We expect that none of the young pulsars should be detectable as an individual point source of γ-ray emission by EGRET. We also expect that pulsar contribution above ∼1 GeV in the SMC is very important.     

The efficiency of gamma-ray emission from pulsars

We present a modified scenario of gamma-ray emission from pulsars within the framework of polar cap models. Our model incorporates the possible acceleration of electron–positron pairs created in magnetospheres, and their subsequent contribution to the gamma-ray luminosity L _γ. It also reproduces the empirical trend in L _γ for seven pulsars detected with Compton Gamma-Ray Observatory ( CGRO ) experiments. At the same time it avoids basic difficulties faced by theoretical models when confronted with observational constraints.   We show that the classical and millisecond pulsars form two distinct branches in the L _γ— L _sd diagram (where L _sd is the spin-down luminosity). In particular, we explain why the millisecond pulsar J0437−4715 has not been detected with any of the CGRO instruments despite its very high position in the ranking list of spin-down fluxes (i.e. L _sd/ D ², where D is a distance). The gamma-ray luminosity predicted for this particular object is about one order of magnitude below the upper limit set by EGRET.     

Explosive reconnection in magnetars

The X-ray activity of anomalous X-ray pulsars and soft γ-ray repeaters may result from the heating of their magnetic corona by direct currents dissipated by magnetic reconnection. We investigate the possibility that X-ray flares and bursts observed from anomalous X-ray pulsars and soft γ-ray repeaters result from magnetospheric reconnection events initiated by development of the tearing mode in magnetically dominated relativistic plasma. We formulate equations of resistive force-free electrodynamics, discuss the relation of the latter to ideal electrodynamics, and give examples of both ideal and resistive equilibria. Resistive force-free current layers are unstable towards the development of small-scale current sheets where resistive effects become important. Thin current sheets are found to be unstable due to the development of the resistive force-free tearing mode. The growth rate of the tearing mode is intermediate between the short Alfvén time-scale  τ_A  and a long resistive time-scale  τ_R: Γ∼ 1/(τ_Rτ_A)^1/2  , similar to the case of non-relativistic non-force-free plasma. We propose that growth of the tearing mode is related to the typical rise time of flares, ∼10 ms. Finally, we discuss how reconnection may explain other magnetar phenomena and ways to test 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.     

The cross-correlation between galaxies and groups: probing the galaxy distribution in and around dark matter haloes

Some massive binaries should contain energetic pulsars which inject relativistic leptons from their inner magnetospheres and/or pulsar wind regions. If the binary system is compact enough, then these leptons can initiate inverse Compton (IC) e^± pair cascades in the anisotropic radiation field of a massive star. γ-rays can be produced in the IC cascade during its development in a pulsar wind region and above a shock in a massive star wind region where the propagation of leptons is determined by the structure of a magnetic field around the massive star. For a binary system with specific parameters, we calculate phase-dependent spectra and fluxes of γ-rays escaping as a function of the inclination angle of the system and for different assumptions on injection conditions of the primary leptons (their initial spectra and location of the shock inside the binary). We conclude that the features of γ-ray emission from such massive binaries containing energetic pulsars should allow us to obtain important information on the acceleration of particles by the pulsars, and on interactions of a compact object with the massive star wind. Predicted γ-ray light curves and spectra in the GeV and TeV energy ranges from such binary systems within our Galaxy and Magellanic Clouds should be observed by future AGILE and GLAST satellites and low-threshold Cherenkov telescopes, such as MAGIC, HESS, VERITAS or CANGAROO III.     

The contribution to the Galactic diffuse gamma-ray spectrum from unresolved rotation-powered pulsars

We consider the contribution to the Galactic diffuse γ-ray emission from unresolved γ-ray pulsars. Based on the thick outer gap model, Monte Carlo methods are used to simulate the properties (period, distance, magnetic field, etc.) of the Galactic population of rotation-powered pulsars the gamma-ray flux of which is lower than the threshold sensitivity of the EGRET detector on the Compton Gamma-Ray Observatory . Furthermore, the contribution to the Galactic diffuse γ-ray spectrum from the unresolved γ-ray pulsars is calculated. Our results indicate that the unresolved γ-ray pulsars contribute ∼5 to ∼10 per cent to the measured Galactic diffuse γ-ray emission if the birth rate of neutron stars in the Galaxy is 1 to 2 per century, and that these pulsars contribute significantly to the observed Galactic diffuse γ-ray emission above 1 GeV. Comparing the model spectrum with the observed spectrum, we show that the unresolved γ-ray pulsars contribute very little to the diffuse emission at lower energies but can account for ∼50 per cent of the observed spectrum above 1 GeV if the product of the birth rate of neutron stars and the γ-ray beaming fraction is about unity. Such a large pulsar contribution can explain the difference (∼60 per cent) between the intensity of the Galactic diffuse emission as measured by EGRET above 1 GeV and model predictions based on cosmic ray–matter interaction only.     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the gamma-ray burst (GRB) prompt emission site from the explosion centre R , by determining the location of the electron's self-absorption frequency in the GRB prompt optical-to-X/γ-ray spectral energy distribution, assuming that the optical and the γ-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely the bulk Lorentz factor of the emitting source Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 300  . For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R -constraints for this non-detection sample and those for the detection sample.     

Gamma-rays and cosmic rays from a pulsar in Cygnus OB2

We argue that γ-ray sources observed in the direction of the Cygnus OB2 association in the GeV and TeV energy range are due to a pulsar that was created by a supernova a few tens of thousands of years ago. The GeV emission is produced by a middle-aged pulsar, a factor of 2 older than the Vela pulsar. The TeV emission is produced by high-energy hadrons and/or leptons accelerated in pulsar wind nebulae. We suggest, moreover, that the excess of cosmic rays at ∼10¹⁸ eV observed from the direction of the Cygnus region can also be related to the appearance of this very energetic pulsar in the Cyg OB2 association. Some of the relativistic hadrons, captured in strong magnetic fields of a high-density region of Cyg OB2, produce neutrons and γ-rays in collisions with matter. These neutrons can arrive from Cyg OB2, creating an excess of cosmic rays.     

One-dimensional electric field structure of an outer gap accelerator – II. γ-ray production resulting from inverse Compton scattering

We study the structure of a stationary and axisymmetric charge-deficient region (or a potential gap) in the outer magnetosphere of a spinning neutron star. A large electric field along the magnetic field lines is created in this potential gap and accelerates migratory electrons (e⁻) and/or positrons (e⁺) to ultrarelativistic energies. Assuming that the gap is immersed in a dense soft photon field, these relativistic e^± radiate γ -ray photons via inverse Compton (IC) scattering. These γ -rays, in turn, produce yet more radiating particles by colliding with ambient soft photons, leading to a pair-production cascade in the gap. The replenished charges partially screen the longitudinal electric field, which is self-consistently solved together with the distribution of e^± and γ -ray photons. It is demonstrated that the voltage drop in the gap is not more than 10¹⁰ V when the background X-ray radiation is as luminous as 10³⁷ erg s⁻¹. However, this value increases with decreasing X-ray luminosity and attains 10¹² V when the X-ray radiation is 10³⁶ erg s⁻¹. In addition, we find useful expressions of the spatial distribution of the particle fluxes and longitudinal electric field, together with the relationship between the voltage drop and the current density. Amazingly, these expressions are valid not only when IC scattering dominates but also when curvature radiation dominates.     

Hot accretion discs with thermal Comptonization and advection in luminous black hole sources

We solve for the structure of a hot accretion disc with unsaturated thermal Comptonization of soft photons and with advection, generalizing the classical model of Shapiro et al. The upper limit on the accretion rate due to advection constrains the luminosity to ≲ 0.15 y^3/5 α^7/5 of the Eddington limit, where y and α are the Compton and viscosity parameters, respectively. The characteristic electron temperature and Thomson optical depth of the inner flow at accretion rates within an order of magnitude of that upper limit are ∼ 10⁹ K and ∼ 1, respectively. The resulting spectra are then in close agreement with the X-ray and soft γ-ray spectra from black hole binaries in the hard state and Seyferts. At low accretion rates, bremsstrahlung becomes the dominant radiative process.     

Gamma-ray emission from individual classical novae

Classical novae are important producers of radioactive nuclei, such as ⁷Be, ¹³N, ¹⁸F, ²²Na and ²⁶Al. The disintegration of these nuclei produces positrons (except for ⁷Be) that through annihilation with electrons produce photons of energies 511 keV and below. Furthermore, ⁷Be and ²²Na decay producing photons with energies of 478 and 1275 keV, respectively, well in the γ-ray domain. Therefore, novae are potential sources of γ-ray emission. We have developed two codes in order to analyse carefully the γ-ray emission of individual classical novae: a hydrodynamical one, which follows both the accretion and the explosion stages, and a Monte Carlo one, able to treat both production and transfer of γ-ray photons. Both codes have been coupled in order to simulate realistic explosions. The properties of γ-ray spectra and γ-ray light curves (for the continuum and for the lines at 511, 478 and 1275 keV) have been analysed, with a special emphasis on the difference between carbon–oxygen and oxygen–neon novae. Predictions of detectability of individual novae by the future SPI spectrometer on board the INTEGRAL satellite are made. Concerning ²⁶Al, its decay produces photons of 1809 keV but this occurs on a time-scale much longer than the typical time interval between nova outbursts in the Galaxy, making it undetectable in individual novae. The accumulated emission of ²⁶Al from many Galactic novae has not been modelled in this paper.     

Emergence of magnetic field due to spin-polarized baryon matter in neutron stars

A model of the ferromagnetic origin of magnetic fields of neutron stars is considered. In this model, the magnetic phase transition occurs inside the core of neutron stars soon after formation. However, owing to the high electrical conductivity the core magnetic field is initially fully screened. We study how this magnetic field emerges for an outside observer. After some time, the induced field that screens the ferromagnetic field decays enough to uncover a detectable fraction of the ferromagnetic field. We calculate the time-scale of decay of the screening field and study how it depends on the size of the ferromagnetic core. We find that the same fractional decay of the screening field occurs earlier for larger cores. We conjecture that weak fields of millisecond pulsars, B ∼10⁸–10⁹ G, could be identified with ferromagnetic fields of unshielded fraction ε ∼10⁻⁴–10⁻³ resulting from the decay of screening fields by a factor 1− ε in ∼10⁸ yr since their birth.     

GeV emission from neutron-rich internal shocks of some long γ-ray bursts

In the neutron-rich internal shocks model for γ-ray bursts (GRBs), the Lorentz factors (LFs) of ion shells are variable, and so are the LFs of accompanying neutron shells. For slow neutron shells with a typical LF of approximate tens, the typical β-decay radius is  ∼10¹⁴–10¹⁵ cm  . As GRBs last long enough  [ T ₉₀ > 14(1 + z ) s]  , one earlier but slower ejected neutron shell will be swept successively by later ejected ion shells in the range  ∼10¹³–10¹⁵ cm  , where slow neutrons have decayed significantly. Part of the thermal energy released in the interaction will be given to the electrons. These accelerated electrons will mainly be cooled by the prompt soft γ-rays and give rise to GeV emission. This kind of GeV emission is particularly important for some very long GRBs and is detectable for the upcoming satellite Gamma-Ray Large Area Space Telescope (GLAST).     

One-dimensional electric field structure of an outer gap accelerator – I. γ-ray production resulting from curvature radiation

We study the structure of a stationary and axisymmetric charge-deficient region (or potential gap) in the outer magnetosphere of a spinning neutron star. Assuming the existence of global current flow patterns in the magnetosphere, the charge depletion causes a large electric field along the magnetic field lines. This longitudinal electric field accelerates migratory electrons and/or positrons to ultrarelativistic energies. These relativistic electrons/positrons radiate γ -ray photons by curvature radiation. These γ -rays, in turn, produce yet more radiating particles by colliding with ambient X-ray photons, leading to a pair production cascade in the gap. The replenished charges partially screen the longitudinal electric field, which is self-consistently solved together with the distribution of e^± and γ -ray photons. We find the voltage drop in the gap as a function of the soft photon luminosity. It is demonstrated that the voltage drop is less than 3×10¹³ V when the background X-ray radiation is as luminous as Vela . However, this value increases with decreasing X-ray luminosity and attains 3×10¹⁵ V when the X-ray luminosity is as low as L _X=10³¹ erg s⁻¹.     

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.     

TeV γ-rays from old supernova remnants

We study the emission from an old supernova remnant (SNR) with an age of around 10⁵ yr and that from a giant molecular cloud (GMC) encountered by the SNR. When the SNR age is around 10⁵ yr, proton acceleration is efficient enough to emit TeV γ-rays both at the shock of the SNR and that in the GMC. The maximum energy of primarily accelerated electrons is so small that TeV γ-rays and X-rays are dominated by hadronic processes,  π⁰  -decay and synchrotron radiation from secondary electrons, respectively. However, if the SNR is older than several 10⁵ yr, there are few high-energy particles emitting TeV γ-rays because of the energy-loss effect and/or the wave-damping effect occurring at low-velocity isothermal shocks. For old SNRs or SNR–GMC interacting systems capable of generating TeV γ-ray emitting particles, we calculated the ratio of TeV γ-ray (1–10 TeV) to X-ray (2–10 keV) energy flux and found that it can be more than  ∼10²  . Such a source showing large flux ratio may be a possible origin of recently discovered unidentified TeV sources.     

High-energy γ-rays from globular clusters

It is expected that specific globular clusters (GCs) can contain up to a hundred of millisecond pulsars. These pulsars can accelerate leptons at the shock waves originated in collisions of the pulsar winds and/or inside the pulsar magnetospheres. Energetic leptons diffuse gradually through the GC Comptonizing stellar and microwave background radiation. We calculate the GeV–TeV γ-ray spectra for different models of injection of leptons and parameters of the GCs assuming reasonable, of the order of 1 per cent, efficiency of energy conversion from the pulsar winds into the relativistic leptons. It is concluded that leptons accelerated in the GC cores should produce well localized γ-ray sources which are concentric with these GCs. The results are shown for four specific GCs (47 Tuc, Ter 5, M13 and M15), in which significant population of millisecond pulsars have been already discovered. We argue that the best candidates, which might be potentially detected by the present Cherenkov telescopes and the planned satellite telescopes (AGILE, GLAST), are 47 Tuc on the Southern hemisphere, and M13 on the Northern hemisphere. We conclude that detection (or non-detection) of GeV–TeV γ-ray emission from GCs by these instruments put important constraints on the models of acceleration of leptons by millisecond pulsars.     

The Parkes Southern Pulsar Survey — II. Final results and population analysis

A survey of the entire southern sky for millisecond and low-luminosity pulsars using the ATNF Parkes radio telescope has now been completed. The survey detected 298 pulsars, of which 101 were previously unknown. The new pulsars include 17 millisecond pulsars. This is the largest sample of both normal and millisecond pulsars detected in any survey. Combining our sample with other recent surveys in the Northern Hemisphere, we present a statistical study of the populations of both normal and millisecond pulsars. We find that the improved statistics allow us to estimate the number and birth-rate of both types of pulsar down to a 400-MHz luminosity limit of 1 mJy kpc². The local surface densities of potentially observable normal pulsars and millisecond pulsars are both about 30 kpc⁻², corresponding to ∼ 30000 potentially observable pulsars of each type in the Galaxy. Once beaming effects are taken into consideration we estimate that the active population of normal pulsars is ∼ 160000. Although there is evidence for flattening of the luminosity function of normal pulsars, this is not evident for millisecond pulsars which probably have a substantial population with luminosities below 1 mJy kpc². After correcting for beaming effects, we estimate that a normal pulsar is born with a luminosity greater than 1 mJy kpc² between once every 60 and 330 yr in the Galaxy. The birth-rate of millisecond pulsars is at least 3 × 10⁻⁶ yr⁻¹ above the same luminosity limit. Modelling the observed transverse speeds of millisecond pulsars using a dynamical simulation, we find their mean birth velocity to be 130 ± 30 km s⁻¹, significantly lower than that of the normal pulsars.     

Precision timing measurements of PSR J1012+5307

We present results and applications of high-precision timing measurements of the binary millisecond pulsar J1012+5307. Combining our radio timing measurements with results based on optical observations, we derive complete 3D velocity information for this system. Correcting for Doppler effects, we derive the intrinsic spin parameters of this pulsar and a characteristic age of 8.6±1.9 Gyr . Our upper limit for the orbital eccentricity of only 8×10⁻⁷ (68 per cent confidence level) is the smallest ever measured for a binary system. We demonstrate that this makes the pulsar an ideal laboratory in which to test certain aspects of alternative theories of gravitation. Our precision measurements suggest deviations from a simple pulsar spin-down timing model, which are consistent with timing noise and the extrapolation of the known behaviour of slowly rotating pulsars.     

A deep search for pulsar wind nebulae using pulsar gating

Using the Australia Telescope Compact Array (ATCA) we have imaged the fields around five promising pulsar candidates to search for radio pulsar wind nebulae (PWNe). We have used the ATCA in its pulsar-gating mode; this enables an image to be formed containing only off-pulse visibilities, thereby dramatically improving the sensitivity to any underlying PWN. Data from the Molonglo Observatory Synthesis Telescope were also used to provide sensitivity on larger spatial scales. This survey found a faint new PWN around PSR B0906−49; here we report on non-detections of PWNe towards PSRs B1046−58, B1055−52, B1610−50 and J1105−6107. Our radio observations of the field around PSR B1055−52 argue against previous claims of an extended X-ray and radio PWN associated with the pulsar. If these pulsars power unseen, compact radio PWNe, upper limits on the radio flux indicate that a fraction of less than 10⁻⁶ of their spin-down energy is used to power this emission. Alternatively, PSRs B1046−58 and B1610−50 may have relativistic winds similar to other young pulsars and the unseen PWN may be resolved and fainter than our surface brightness sensitivity threshold. We can then determine upper limits on the local interstellar medium (ISM) density of 2.2×10⁻³ and 1×10⁻² cm⁻³, respectively. Furthermore, we derive the spatial velocities of these pulsars to be ∼450 km s⁻¹ and thus rule out the association of PSR B1610−50 with supernova remnant (SNR) G332.4+00.1 (Kes 32). Strong limits on the ratio of unpulsed to pulsed emission are also determined for three pulsars.