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.     

Radio-Loud and Radio-Quiet Gamma-Ray Pulsars from the Galaxy and the Gould Belt

We present results of a population synthesis study of radio-loud and radio-quiet γ-ray pulsars from the Galactic plane and the Gould Belt. The simulation includes the Parkes multibeam pulsar survey, realistic beam geometries for radio and γ-ray emission from neutron stars and the new electron density model of Cordes and Lazio. Normalizing to the number of radio pulsars observed by a set of nine radio surveys, the simulation suggests a neutron star birth rate of 1.4 neutron stars per century in the Galactic plane. In addition, the simulation predicts 19 radio-loud and 7 radio-quiet γ-ray pulsars from the plane that EGRET should have observed as point sources. Assuming that during the last 5 Myr the Gould Belt produced 100 neutron stars, only 10 of these would be observed as radio pulsars with three radio-loud and four radio-quiet γ-ray pulsars observed by EGRET. These results are in general agreement with the recent number of about 25 EGRET error boxes that contain Parkes radio pulsars. Since the Gould Belt pulsars are relatively close by, the selection of EGRET radio-quiet γ-ray pulsars strongly favors large impact angles, β, in the viewing geometry where the off-beam emission from curvature radiation provides the γ-ray flux. Therefore, the simulated EGRET radio-quiet γ-ray pulsars, being young and nearby, most closely reflect the current shape of the Gould Belt suggesting that such sources may significantly contribute to the EGRET unidentified γ-ray sources correlated with the Gould Belt.     

Contribution to diffuse gamma-rays in the Galactic Centre region from unresolved millisecond pulsars

Diffuse gamma-rays in the Galactic Centre region have been studied. We propose that there exists a population of millisecond pulsars in the Galactic Centre, which emit GeV gamma-rays through synchrotron-curvature radiation as predicted by outer gap models. These GeV gamma-rays from unresolved millisecond pulsars probably contribute to the diffuse gamma-ray spectrum detected by EGRET which displays a break at a few GeV. We have used a Monte Carlo method to obtain simulated samples of millisecond pulsars in the Galactic Centre region covered by EGRET  (∼ 15)  according to the different period and magnetic field distributions from observed millisecond pulsars in the Galactic field and globular clusters, and superposed their synchrotron-curvature spectra to derive the total GeV flux. Our simulated results suggest that there probably exist about 6000 unresolved millisecond pulsars in the region of angular resolution of EGRET, the emissions of which could contribute significantly to the observed diffuse gamma-rays in the Galactic Centre.     

The contribution of unresolved radio-loud AGN to the extragalactic diffuse gamma-ray background

We present a calculation of the blazar contribution to the extragalactic diffuse γ -ray background (EGRB) in the EGRET energy range. Our model is based on inverse-Compton scattering as the dominant γ -ray production process in the jets of flat spectrum radio quasars (FSRQs) and BL Lac objects, and on the unification scheme of radio-loud AGN. According to this picture, blazars represent the beamed fraction of the Fanaroff–Riley radio galaxies (FR galaxies).
The observed log  N –log  S distribution and redshift distribution of both FSRQs and BL Lacs constrain our model. Depending slightly on the evolutionary behaviour of blazars, we find that unresolved AGN underproduce the intensity of the extragalactic background radiation. With our model only 20–40 per cent of the extragalactic background emission can be explained by unresolved blazars if we integrate to a maximum redshift of Z _max=3. For Z _max=5, blazars could account for 40–80 per cent of the EGRB. Roughly 70–90 per cent of the AGN contribution to the EGRB would result from BL Lacs. While the systematic uncertainties in our estimate for the FSRQ contribution appear small, in the case of BL Lacs our model parameters are not consistent with the results from studies in other wavelength regimes, and therefore may have larger systematic uncertainties. Thus we end up with two possibilities, depending on whether we underpredict or overpredict the BL Lac contribution: either unresolved AGN cannot account for the entire EGRB, or unresolved BL Lacs produce the observed background.
We predict a significant flattening of the γ -ray log  N –log  S function in the next two decades of flux below the EGRET threshold.     

The birthrate of magnetars

Magnetars, neutron stars with ultrastrong magnetic fields  ( B ∼ 10¹⁴−10¹⁵G)  , manifest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray pulsars. This study estimates the birthrate of magnetars to be ∼0.22 per century with a Galactic population comprising ∼17 objects. A population synthesis was carried out based on the five anomalous X-ray pulsars detected in the ROSAT All Sky Survey by comparing their number to that of massive OB stars in a well-defined volume. Additionally, the group of seven X-ray dim isolated neutron stars detected in the same survey were found to have a birthrate of ∼2 per century with a Galactic population of ∼22 000 objects.     

Source population synthesis and the Galactic diffuse gamma-ray emission

Population synthesis is used to study the contribution from unresolved sources to the Galactic ridge emission measured by EGRET. Synthesized source counts are compared with the 3rd EGRET catalogue at low and high latitudes. For pulsar-like populations, 5–10% of the emission >100 MeV comes from sources below the EGRET threshold. A steeper luminosity function can increase this to 20% without violating EGRET source statistics. Less luminous populations can produce much higher values without being detected. Since the unresolved source spectrum is different from the interstellar spectrum, it could provide an explanation of the observed MeV and GeV excesses above the predictions, and we give an explicit example of how this could work.     

A Survey of Unidentified Egret Sources at TeV Energies

The Whipple Observatory 10 m γ-ray telescope has been used to survey the error boxes of 24 EGRET unidentified sources in an attempt to find counterparts at energies of 350 GeV and above. In no case is a statistically significant signal found in the EGRET error box which implies that, at least for this sample, the γ-ray spectra of these sources steepen between 100 MeV and 350 GeV.     

The gamma-ray conversion efficiency of rotation-powered pulsars

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

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.     

Almost analytic models of ultramagnetized neutron star envelopes

Recent ROSAT measurements show that the X-ray emission from isolated neutron stars is modulated at the stellar rotation period. To interpret these measurements, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We present nearly analytic models of the thermal structure of the envelopes of ultramagnetized neutron stars. Specifically, we examine the limit in which only the ground Landau level is filled. We use the models to estimate the amplitude of modulation expected from non-uniformities in the surface temperatures of strongly magnetized neutron stars. In addition, we estimate cooling rates for stars with fields B  ∼ 10¹⁵ − 10¹⁶ G, which are relevant to models that invoke 'magnetars' to account for soft γ-ray emission from some repeating sources.     

Inverse Compton e± pair cascade model for the γ-ray production in massive binary LSI +61° 303

We apply an inverse Compton   e ^±  pair cascade model for γ-ray production in the massive binary system LSI +61° 303 assuming that electrons are accelerated already inside the inner part of the jet launched by the compact object. γ-ray spectra, affected by the cascade process, and lower energy spectra, from the synchrotron cooling of the highest energy electrons in the jet, are calculated as a function of the phase of this binary system. γ-ray spectra expected in such a model have different shape than those ones produced by electrons in the jet directly to observer. Moreover, the model predicts clear anticorrelation between γ-ray fluxes in the GeV (1–10 GeV) and TeV (>200 GeV) energy ranges with the peak of the TeV emission at the phase ∼0.5 (the peak half-width ranges between the phases ∼0.4–0.9, for the inclination of the binary system equal to 60°, and ∼0.4–0.1 for 30°). The fine features of TeV γ-ray emission (fluxes and spectral shapes) as a function of the phase of the binary system are consistent with recent observations reported by the MAGIC collaboration. Future simultaneous observations in the GeV energies (by the GLAST and AGILE telescopes) and in the TeV energies (by the MAGIC and VERITAS telescopes) should test other predictions of the considered model supporting or disproving the hypothesis of acceleration of electrons already in the inner part of the microquasar jets.     

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).     

Searching for hidden Wolf–Rayet stars in the Galactic plane – 15 new Wolf–Rayet stars

We report the discovery of 15 previously unknown Wolf–Rayet (WR) stars found as part of an infrared (IR) broad-band study of candidate WR stars in the Galaxy. We have derived an empirically based selection algorithm which has selected ∼5000 WR candidate stars located within the Galactic plane drawn from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (mid-IR) and Two-Micron All-Sky Survey (near-IR) catalogues. Spectroscopic follow-up of 184 of these reveals 11 nitrogen-rich (WN) and four carbon-rich (WC) WR stars. Early WC subtypes are absent from our sample and none shows evidence for circumstellar dust emission. Of the candidates which are not WR stars, ∼120 displayed hydrogen emission-line features in their spectra. Spectral features suggest that the majority of these are in fact B supergiants/hypergiants, ∼40 of these are identified Be/B[e] candidates.
Here, we present the optical spectra for six of the newly detected WR stars, and the near-IR spectra for the remaining nine of our sample. With a WR yield rate of ∼7 per cent and a massive star detection rate of ∼65 per cent, initial results suggest that this method is one of the most successful means for locating evolved, massive stars in the Galaxy.     

TeV gamma-rays from accreting magnetars in massive binaries

This paper focuses on neutron stars (NS) of the magnetar type inside massive binary systems. We determine the conditions under which the matter from the stellar wind can penetrate the inner magnetosphere of the magnetar. At a certain distance from the NS surface, the magnetic pressure can balance the gravitational pressure of the accreting matter, creating a very turbulent, magnetized transition region. It is suggested that this region provides good conditions for the acceleration of electrons to relativistic energies. These electrons lose energy due to the synchrotron process and inverse Compton (IC) scattering of the radiation from the nearby massive stellar companion, producing high-energy radiation from X-rays up to ∼TeV γ-rays. The primary γ-rays can be further absorbed in the stellar radiation field, developing an IC  e^±  pair cascade. We calculate the synchrotron X-ray emission from primary electrons and secondary  e^±  pairs and the IC γ-ray emission from the cascade process. It is shown that quasi-simultaneous observations of the TeV γ-ray binary system LSI +61 303 in the X-ray and TeV γ-ray energy ranges can be explained with such an accreting magnetar model.     

Unidentified EGRET sources and the extragalactic gamma-ray background

The large majority of EGRET point sources remain to this day without an identified low-energy counterpart. Whatever the nature of the EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse gamma-ray background: if most unidentified objects are extragalactic, faint unresolved sources of the same class contribute to the background, as a distinct extragalactic population; on the other hand, if most unidentified sources are Galactic, their counterparts in external galaxies will contribute to the unresolved emission from these systems. Understanding this component of the gamma-ray background, along with other guaranteed contributions from known sources, is essential in any attempt to use gamma-ray observations to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the extragalactic gamma-ray background is likely to be important, and we find that it is. Additionally, we comment on how the anticipated GLAST measurement of the diffuse gamma-ray background will change, depending on the nature of the majority of these sources.     

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.     

Studies of accretion flows around rotating black holes – I. Particle dynamics in a pseudo-Kerr potential

We study the formation of the absorption features, called the cyclotron–annihilation lines, in the γ-spectra of the neutron stars (pulsars), owing to the fundamental quantum-electrodynamic effect of the one–photon pair creation in magnetized vacuum. As a result, we substantiate a new method for the determination of the neutron star magnetic fields B based on measuring the interval between the main annihilation and the first cyclotron–annihilation absorption lines. It is found that these lines may be easily resolved, and, consequently, the method is surely applicable if the following conditions are satisfied. (i) A γ-source has to be compact enough and located near a star, but not close to its magnetic poles. For instance, it may be a disc in the plane of a star magnetic equator with latitudinal angular width less than     and radial extent up to 25 per cent of the star radius. (ii) The source is to produce detectable γ-radiation at large angles ≳60° to the local magnetic field. Being situated in a closed field line region and having a broad radiation pattern, such a source is not what is usually considered in the context of the polar cap and outer gap models of the pulsar γ-emission dealing with open field lines only. (iii) Magnetic field strength must lie in a certain narrow interval with the centre at  ∼(3–4) × 10¹²  G. Its width depends on the star orientation and disc radial extend and in the most favourable case is about 20–30 per cent of its lower boundary. Finally, the influence of the star rotation on this method employment is considered and new possibilities arising from forthcoming polarization observations are briefly discussed.     

Absorption of nuclear γ-rays on the starlight radiation in FR I sources: the case of Centaurus A

Several BL Lac objects are confirmed sources of variable and strongly Doppler-boosted TeV emission produced in the nuclear portions of their relativistic jets. It is more than probable that also many of the Fanaroff–Riley type I (FR I) radio galaxies, believed to be the parent population of BL Lacs, are TeV sources, for which Doppler-hidden nuclear γ-ray radiation may be only too weak to be directly observed. Here we show, however, that about 1 per cent of the total time-averaged TeV radiation produced by the active nuclei of low-power FR I radio sources is inevitably absorbed and re-processed by photon–photon annihilation on the starlight photon field, and the following emission of the created and quickly isotropized electron–positron pairs. In the case of the radio galaxy Centaurus A, we found that the discussed mechanism can give a distinctive observable feature in the form of an isotropic γ-ray halo. It results from the electron–positron pairs injected to the interstellar medium of the inner parts of the elliptical host by the absorption process, and upscattering starlight radiation via the inverse-Compton process mainly to the GeV–TeV photon energy range. Such a galactic γ-ray halo is expected to possess a characteristic spectrum peaking at ∼0.1 TeV photon energies, and the photon flux strong enough to be detected by modern Cherenkov Telescopes and, in the future, by GLAST. These findings should apply as well to the other nearby FR I sources.     

Be phenomenon in open clusters: results from a survey of emission-line stars in young open clusters

Emission-line stars in young open clusters are identified to study their properties, as a function of age, spectral type and evolutionary state. 207 open star clusters were observed using the slitless spectroscopy method and 157 emission stars were identified in 42 clusters. We have found 54 new emission-line stars in 24 open clusters, out of which 19 clusters are found to house emission stars for the first time. About 20 per cent clusters harbour emission stars. The fraction of clusters housing emission stars is maximum in both the 0–10 and 20–30 Myr age bin (∼40 per cent each). Most of the emission stars in our survey belong to Classical Be class (∼92 per cent) while a few are Herbig Be stars (∼6 per cent) and Herbig Ae stars (∼2 per cent). The youngest clusters to have Classical Be stars are IC 1590, NGC 637 and 1624 (all 4 Myr old) while NGC 6756 (125–150 Myr) is the oldest cluster to have Classical Be stars. The Classical Be stars are located all along the main sequence (MS) in the optical colour–magnitude diagrams (CMDs) of clusters of all ages, which indicates that the Be phenomenon is unlikely due to core contraction near the turn-off. The distribution of Classical Be stars as a function of spectral type shows peaks at B1–B2 and B6–B7 spectral types. The Be star fraction [N(Be)/N(B+Be)] is found to be less than 10 per cent for most of the clusters and NGC 2345 is found to have the largest fraction (∼26 per cent). Our results indicate there could be two mechanisms responsible for the Classical Be phenomenon. Some are born Classical Be stars (fast rotators), as indicated by their presence in clusters younger than 10 Myr. Some stars evolve to Classical Be stars, within the MS lifetime, as indicated by the enhancement in the fraction of clusters with Classical Be stars in the 20–30 Myr age bin.