Population statistics study of radio and gamma-ray pulsars in the Galactic plane

We present results of our pulsar population synthesis of normal and millisecond pulsars in the Galactic plane. Over the past several years, a program has been developed to simulate pulsar birth, evolution and emission using Monte Carlo techniques. We have added to the program the capability to simulate millisecond pulsars, which are old, recycled pulsars with extremely short periods. We model the spatial distribution of the simulated pulsars by assuming that they start with a random kick velocity and then evolve through the Galactic potential. We use a polar cap/slot gap model for γ-ray emission from both millisecond and normal pulsars. From our studies of radio pulsars that have clearly identifiable core and cone components, in which we fit the polarization sweep as well as the pulse profiles in order to constrain the viewing geometry, we develop a model describing the ratio of radio core-to-cone peak fluxes. In this model, short period pulsars are more cone-dominated than in our previous studies. We present the preliminary results of our recent study and the implications for observing these pulsars with GLAST and AGILE.      

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

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

Region of enhanced flux of cosmic rays with PeV energies toward the pulsars PSR J1840+5640 and LAT PSR J1836+5925

Analysis of the arrival directions of extensive air showers (EASs) detected on the EAS MSU array and the prototype of the EAS-1000 array has revealed a region of enhanced flux of cosmic rays with PeV energies toward the pulsars PSR J1840+5640 and LAT PSR J1836+5925 at a confidence level up to 4.5σ. The first pulsar was discovered almost 30 years ago and is a well-studied old radio pulsar at a distance of 1.7 kpc from the Solar system. The second pulsar belongs to a new class of pulsars discovered by the Fermi Gamma-Ray Observatory whose pulsations are seen neither in the X-ray nor in the radio bands, but only in the gamma-ray energy range (gamma-ray-only pulsars). In our opinion, the existence of a region of enhanced cosmic-ray flux in the data sets obtained on two different arrays suggests that the pulsars can make a noticeable contribution to the flux of Galactic cosmic rays with PeV energies.     

Fe emission and ionized excess absorption in the luminous quasar 3C 109 with XMM–Newton

The Parkes High-Latitude pulsar survey covers a region of the sky enclosed by Galactic longitudes 220° < l < 260° and Galactic latitudes | b | < 60°. The observations have been performed using the 20-cm multibeam receiver on the Parkes 64-m radio telescope. A total of 6456 pointings of 265 s each have been collected. The system adopted provided a sensitivity limit, for long-period pulsars with 5 per cent duty cycles, of ∼0.5 mJy. Data analysis resulted in the detection of 42 pulsars of which 18 were new discoveries. Four of these belong to the class of the millisecond – or recycled – pulsars; three of these four are in binary systems. The double pulsar system J0737−3039 is among those and has been presented elsewhere. Here, we discuss the other discoveries and provide timing parameters for the objects for which we have a phase-connected solution.     

Deep radio observations of the rich globular clusters Terzan 5 and Liller 1

We have obtained deep 6, 20 and 90 cm observation of two of the richest globular clusters in the galaxy, Terzan 5 and Liller 1. In addition to the well known eclipsing pulsar, PSR1744-24A, our 20 cm image of Terzan 5 displays numerous point sources within 30″, or 4 core radii of the cluster center. The density of these objects rises rapidly toward the core, where an elongated region of emission is found. This feature may be due to the presence of a higher density of point sources than can be resolved by the 1″ resolution of our data. Observations at 6 and 90 cm show that the brightest individual sources, as well as the extended emission, possess the steep spectra expected of pulsars. Furthermore, we find that the strength of the sources in the 20 cm image agrees well with the standard pulsar luminosity function. The total luminosity and number of objects observed in Terzan 5 suggest that this cluster may contain 100 pulsars — probably more than any other Galactic cluster. Although we do not resolve individual sources in Liller 1, comparison of our images with other observations suggests that we are again observing the sum of emission from a number of objects. The 90 cm flux densities of Terzan 5 and Liller 1 are 35 and 9 mJy, respectively, implying that a number of bright pulsars in these clusters have been hidden from pulsed searches, perhaps by a combination of dispersion smearing, scattering and doppler broadening.     

GLAST large area telescope multiwavelength planning

Gamma-ray astrophysics depends in many ways on multiwavelength studies. The Gamma-ray Large Area Space Telescope (GLAST) Large Area Telescope (LAT) Collaboration has started multiwavelength planning well before the scheduled 2007 launch of the observatory. Some of the high-priority multiwavelength needs include: (1) availability of contemporaneous radio and X-ray timing of pulsars; (2) expansion of blazar catalogs, including redshift measurements; (3) improved observations of molecular clouds, especially at high galactic latitudes; (4) simultaneous broad-band blazar monitoring; (5) characterization of gamma-ray transients, including gamma ray bursts; (6) radio, optical, X-ray and TeV counterpart searches for reliable and effective sources identification and characterization. Several of these activities are needed to be in place before launch.      

TeV neutrinos and gamma-rays from pulsars

Recent studies suggest that pulsars could be strong sources of TeV muon neutrinos provided positive ions are accelerated by pulsar polar caps to PeV energies. In such a situation, muon neutrinos are produced through the Δ-resonance in interactions of pulsar-accelerated ions with its thermal radiation field. High-energy gamma-rays should also be produced simultaneously in pulsar environment as both charged and neutral pions are generated in the interactions of energetic hadrons with the ambient photon fields. Here, we estimate TeV gamma-ray flux at the Earth from a few nearby young pulsars. When compared with the observations, we find that proper consideration of the effect of polar cap geometry in flux calculation is important. Incorporating such an effect, we obtain the (revised) event rates at the Earth due to a few potential nearby pulsars. The results suggest that pulsars are unlikely to be detected by the upcoming neutrino telescopes. We also estimate TeV gamma-ray and neutrino fluxes from pulsar nebulae for the adopted model of particle acceleration.     

Populating the Galaxy with pulsars – I. Stellar and binary evolution

The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad hoc manner. Thus, no complete and direct comparison with observations of the pulsar population within the Galactic disc has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed typical millisecond pulsar formation process. Results also focus on reproducing the observed     diagram for Galactic pulsars and how this precludes short time-scales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass.     

Very high energy γ-rays from galactic sources

The field of Very High Energy (VHE) gamma ray astronomy using the Atmospheric Cerenkov Technique has entered an interesting phase with detection of various galactic and extragalactic sources. Among galactic sources, only the Crab nebula has been established as a standard candle. Most observations on pulsars are in agreement as to the necessity for the GeV spectra to steepen at < 200 GeV. While the Imaging method for increase of sensitivity has been successful with many results, an alternate technique —Wavefront Sampling Technique- is also being used by an increasing number of experiments. The recently commissioned experiment at Pachmarhi (PACT) in India is presented as an example of this technique. Preliminary results from this experiment show detection of VHE γ-ray emission from (a) the Crab nebula at a high significance and (b) Crab and Geminga pulsars at > 1.5 TeV which could be the second component of the Outer Gap pulsar models.     

The distribution of space velocities of pulsars

Using reasonable assumptions, we derive the distribution function of pulsar space velocities from the distribution of pulsar heights above the Galactic plane. We find that for 137 pulsars of the north Galactic hemisphere, the space velocity distribution is discrete, and forms two separate groups. We estimate the velocity ranges and the mean pulsar velocity of each group.Translated from Astrofizika, Vol. 37, No. 2, pp. 245–254, April–June, 1994.     

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.     

Detectivity of Unidentified Egret Sources as Gamma-Ray Pulsars with the Magic Telescope

Most of the unidentified gamma ray sources detected near the Galactic plane by EGRET aboard CGRO are expected to be gamma ray pulsars. We present a study about the detectability and identification of some unidentified EGRET sources with the MAGIC telescope. We list some unidentified gamma ray sources from the third EGRET catalogue to be detected with MAGIC taking into account some important conditions such as the variability parameter of the source, spectral index, inclusion in the GeV catalogue (ApJ 488, 1997, p. 872) and possible associations with known X-ray/radio sources located within the error box of the unidentified gamma ray source. We show the required observation time of these gamma ray pulsar candidates to be detected by MAGIC telescope considering reasonable values of cut-off energy. To be more realistic, we have chosen the zenith angle corresponding to the source culmination in the simulation of the effective area A since the observation time is function of the effective area. In addition to this study, it is very important to consider the extrapolated EGRET flux at MAGIC energies above 30 GeV of the gamma ray pulsar candidates taking the MAGIC sensitivity.     

Population studies of the unidentified EGRET sources

The third EGRET catalog contains a large number of unidentified sources. This subset of objects is expected to include known gamma-ray emitters of Galactic origin such as pulsars and supernova remnants, in addition to an extragalactic population of blazars. However, current data allows the intriguing possibility that some of these objects may represent a new class of yet undiscovered gamma-ray sources. Many theoretically motivated candidate emitters (e.g. clumps of annihilating dark matter particles) have been suggested to account for these detections. We take a new approach to determine to what extent this population is Galactic and to investigate the nature of the possible Galactic component. By assuming that galaxies similar to the Milky Way should host comparable populations of objects, we constrain the allowed Galactic abundance and distribution of various classes of gamma-ray sources using the EGRET data set. We find it is highly improbable that a large number of the unidentified sources are members of a Galactic halo population, but that a distribution of the sources entirely in the disk and bulge is plausible. Finally, we discuss the additional constraints and new insights that GLAST will provide.     

候选射电脉冲星的预测

对已知脉冲星大样本的统计和分析表明,绝大多数脉冲星的线偏振大于１０％,谱指数低于－１．比较了脉冲星和射电类星体及致密陡谱源的频谱和偏振特性的差别,以高线偏振、高谱指数、高银纬及在低频下较低流量等作为射电脉冲星候选者的判据,从新近释放的ＶＬＡ巡天和ＷＳＲＴ巡天结果中挑选出一个较粗的候选脉冲星样本．进一步考虑近期在部分天区中已作过高灵敏度的脉冲星巡天,所有落入这些天区内的候选源也都从本样本中排出．以候选源偏振度大小为序,给出了一个共４７颗候选源的子样本．     

On selection effects in pulsar searches

Selection effects are a major source of error in statistical studies of pulsar data since the observed sample is a biased subset of the full galactic pulsar population. It is important to identify all selection effects and make a reasonable model before attempting to determine pulsar properties. Here we discuss a hitherto neglected selection effect which is a function of the periodP of the pulsar. We find that short-P pulsars are more difficult to detect, particularly if their dispersion measures are high. We also discuss a declination-dependent selection effect in the II Molonglo Survey (II MS), and find some evidence for the existence of both selection effects in the pulsar data from this survey. We discuss the implications of these additional selection effects for the recently proposed ‘injection’ of pulsars whereby pulsars seem to switch on only at longerP. Using the II MS data we calculate that the observability of pulsars withP between 0.0 s and 0.5 s is about 18 per cent less with the new selection effects than hitherto believed; the mean correction is 6 per cent forP between 0.5 s and 1.0 s. We conclude that injection is not qualitatively affected by these corrections.     

New pulsar rotation measures and the Galactic magnetic field

We measured a sample of 150 pulsar rotation measures (RMs) using the 20-cm receiver of the Parkes 64-m radio telescope. 46 of the pulsars in our sample have not had their RM values previously published, whereas 104 pulsar RMs have been revised. We used a novel quadratic fitting algorithm to obtain an accurate RM from the calibrated polarization profiles recorded across 256 MHz of receiver bandwidth. The new data are used in conjunction with previously known dispersion measures and the NE2001 electron-density model to study models of the direction and magnitude of the Galactic magnetic field.     

A study of dynamical effects in the observed second time-derivative of the spin or orbital frequencies of pulsars

The observed values of the time-derivatives of the spin or orbital frequency of pulsars are affected by their dynamical properties. We derive thorough analytical expressions for such dynamical contributions in terms of the Galactic coordinates, the proper motion, the pulsar distance, and the radial velocity. We find that the effects of the dynamical terms in the second-derivative of frequencies or parameters based on such second derivatives, e.g., braking index, are usually negligible. However, unique pulsars for which the effects of the dynamical terms are significant can exist. In particular, dynamical effects can make the magnitude of the observed value of the braking index to be in the order of thousand while the true value of it is close to the theoretically expected value three, especially if the pulsars lie close to the Galactic centre. Dynamics can also affect the value of the second derivative of the orbital frequency of a binary pulsar at the first decimal place. We also emphasize the fact that our expressions provide more accurate results than pre-existing approximate ones that exclude some of the terms. Comparison with a set of pulsars showed that the median value of the difference between the results obtained by our method and a pre-existing method is about 50 percent.     

脉冲星的星系磁场

脉冲星的偏振信息是理解脉冲星辐射区的重要手段，利用澳大利亚的64m射电望远镜进行大量的脉冲星观测，得到了一批脉冲星的偏振轮廓和偏振参数，编辑了几乎所有发表的脉冲星轮廓资料，系统总结了脉冲星圆偏辐射的规律，为理论上解释脉冲星辐射这一重要难题提供观测依据和物理限制，利用脉冲星作为探针，研究了银河系磁场结构和模型，确定了银河系BS磁场模型，发现了银河系上下反对称的环向磁场，并首次对星系尺度的发电机类型进行判别，证认出A0型发电机运行于银河系，发现了银晕中的垂直磁场和M31及银盘中的非常延展的磁场，探测到NGC2997星系中由内到外的旋涡磁场，并提出可能有两种发电机在这个星系的不同区域运行。     

Pulsar simulations for the Fermi Large Area Telescope

Pulsars are among the prime targets for the Large Area Telescope (LAT) aboard the recently launched Fermi observatory. The LAT will study the gamma-ray Universe between 20 MeV and 300 GeV with unprecedented detail. Increasing numbers of gamma-ray pulsars are being firmly identified, yet their emission mechanisms are far from being understood. To better investigate and exploit the LAT capabilities for pulsar science, a set of new detailed pulsar simulation tools have been developed within the LAT collaboration. The structure of the pulsar simulator package (PulsarSpectrum) is presented here. Starting from photon distributions in energy and phase obtained from theoretical calculations or phenomenological considerations, gamma-rays are generated and their arrival times at the spacecraft are determined by taking into account effects such as barycentric effects and timing noise. Pulsars in binary systems also can be simulated given orbital parameters. We present how simulations can be used for generating a realistic set of gamma-rays as observed by the LAT, focusing on some case studies that show the performance of the LAT for pulsar observations.