ABSTRACT Observations by the Fermi Large Area Telescope of γ-ray millisecond pulsar (MSP) light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger ...pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ∼10 GeV. Fermi has also uncovered many new MSPs, impacting Galactic stellar population models. We investigate the contribution of Galactic MSPs to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day MSPs. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow (BW) and redback (RB) systems. Since MSPs are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of MSPs are only modest contributors around a few tens of GeV to the lepton fluxes measured by the Alpha Magnetic Spectrometer, PAMELA, and Fermi, after which this component cuts off. The contribution by BWs and RBs may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.
In more than four years of observation the Large Area Telescope on board the Fermi satellite has identified pulsed γ-ray emission from more than 80 young or middle-aged pulsars, in most cases ...providing light curves with high statistics. Fitting the observed profiles with geometrical models can provide estimates of the magnetic obliquity α and of the line of sight angle ζ, yielding estimates of the radiation beaming factor and radiated luminosity. Using different γ-ray emission geometries (Polar Cap, Slot Gap, Outer Gap, One Pole Caustic) and core plus cone geometries for the radio emission, we fit γ-ray light curves for 76 young or middle-aged pulsars and we jointly fit their γ-ray plus radio light curves when possible. We find that a joint radio plus γ-ray fit strategy is important to obtain (α,ζ) estimates that can explain simultaneously detectable radio and γ-ray emission: when the radio emission is available, the inclusion of the radio light curve in the fit leads to important changes in the (α,ζ) solutions. The most pronounced changes are observed for Outer Gap and One Pole Caustic models for which the γ-ray only fit leads to underestimated α or ζ when the solution is found to the left or to the right of the main α-ζ plane diagonal respectively. The intermediate-to-high altitude magnetosphere models, Slot Gap, Outer Gap, and One pole Caustic, are favoured in explaining the observations. We find no apparent evolution of α on a time scale of 106 years. For all emission geometries our derived γ-ray beaming factors are generally less than one and do not significantly evolve with the spin-down power. A more pronounced beaming factor vs. spin-down power correlation is observed for Slot Gap model and radio-quiet pulsars and for the Outer Gap model and radio-loud pulsars. The beaming factor distributions exhibit a large dispersion that is less pronounced for the Slot Gap case and that decreases from radio-quiet to radio-loud solutions. For all models, the correlation between γ-ray luminosity and spin-down power is consistent with a square root dependence. The γ-ray luminosities obtained by using the beaming factors estimated in the framework of each model do not exceed the spin-down power. This suggests that assuming a beaming factor of one for all objects, as done in other studies, likely overestimates the real values. The data show a relation between the pulsar spectral characteristics and the width of the accelerator gap. The relation obtained in the case of the Slot Gap model is consistent with the theoretical prediction.
ABSTRACT Nearly one-third of the sources listed in the Third Fermi Large Area Telescope (LAT) catalog (3FGL) remain unassociated. It is possible that predicted and even unanticipated gamma-ray source ...classes are present in these data waiting to be discovered. Taking advantage of the excellent spectral capabilities achieved by the Fermi LAT, we use machine-learning classifiers (Random Forest and XGBoost) to pinpoint potentially novel source classes in the unassociated 3FGL sample outside the Galactic plane. Here we report a total of 34 high-confidence Galactic candidates at . The currently favored standard astrophysical interpretations for these objects are pulsars or low-luminosity globular clusters hosting millisecond pulsars (MSPs). Yet these objects could also be interpreted as dark matter annihilation taking place in ultra-faint dwarf galaxies or dark matter subhalos. Unfortunately, Fermi LAT spectra are not sufficient to break degeneracies between the different scenarios. Careful visual inspection of archival optical images reveals no obvious evidence for low-luminosity globular clusters or ultra-faint dwarf galaxies inside the 95% error ellipses. If these are pulsars, this would bring the total number of MSPs at to 106, down to an energy flux 4.0 × 10−12 erg cm−2 s−1 between 100 MeV and 100 GeV. We find this number to be in excellent agreement with predictions from a new population synthesis of MSPs that predicts 100-126 high-latitude 3FGL MSPs depending on the choice of high-energy emission model. If, however, these are dark matter substructures, we can place upper limits on the number of Galactic subhalos surviving today and on dark matter annihilation cross sections. These limits are beginning to approach the canonical thermal relic cross section for dark matter particle masses below ∼100 GeV in the bottom quark ( ) annihilation channel.
Thanks to the huge amount of γ-ray pulsar photons collected by the Fermi Large Area Telescope since its launch in June 2008, it is now possible to constrain γ-ray geometrical models by comparing ...simulated and observed light-curve morphological characteristics. We assumed vacuum-retarded dipole (VRD) pulsar magnetic field and tested simulated and observed morphological light-curve characteristics in the framework of two pole emission geometries, Polar Cap (PC) and Slot Gap (SG), and one pole emission geometries, traditional Outer Gap (OG) and One Pole Caustic (OPC). Radio core plus cone emission was assumed for the pulsars of the simulated sample. We compared simulated and observed recurrence of class shapes and peak multiplicity, peak separation, radio-lag distributions, and trends of peak separation and radio lag as a function of observable and non-observable pulsar parameters. We studied how pulsar morphological characteristics change in multi-dimensional observable and non-observable pulsar parameter space. The PC model gives the poorest description of the LAT pulsar light-curve morphology. The OPC best explains both the observed γ-ray peak multiplicity and shape classes. The OPC and SG models describe the observed γ-ray peak-separation distribution for low- and high-peak separations, respectively. This suggests that the OPC geometry best explains the single-peak structure but does not manage to describe the widely separated peaks predicted in the framework of the SG model as the emission from the two magnetic hemispheres. The OPC radio-lag distribution shows higher agreement with observations suggesting that assuming polar radio emission, the γ-ray emission regions are likely to be located in the outer magnetosphere. Alternatively, the radio emission altitude could be higher that we assumed. We compared simulated non-observable parameters with the same parameters estimated for LAT pulsars in the framework of the same models. The larger agreement between simulated and LAT estimations in the framework of the OPC suggests that the OPC model best predicts the observed variety of profile shapes. The larger agreement obtained between observations and the OPC model predictions jointly with the need to explain the abundant 0.5 separated peaks with two-pole emission geometries, calls for thin OPC gaps to explain the single-peak geometry but highlights the need of two-pole caustic emission geometry to explain widely separated peaks.
With the large sample of young γ-ray pulsars discovered by the Fermi Large Area Telescope (LAT), population synthesis has become a powerful tool for comparing their collective properties with model ...predictions. We synthesised a pulsar population based on a radio emission model and four γ-ray gap models (Polar Cap, Slot Gap, Outer Gap, and One Pole Caustic). Applying γ-ray and radio visibility criteria, we normalise the simulation to the number of detected radio pulsars by a select group of ten radio surveys. The luminosity and the wide beams from the outer gaps can easily account for the number of Fermi detections in 2 years of observations. The wide slot-gap beam requires an increase by a factor of ~10 of the predicted luminosity to produce a reasonable number of γ-ray pulsars. Such large increases in the luminosity may be accommodated by implementing offset polar caps. The narrow polar-cap beams contribute at most only a handful of LAT pulsars. Using standard distributions in birth location and pulsar spin-down power (Ė), we skew the initial magnetic field and period distributions in a an attempt to account for the high ĖFermi pulsars. While we compromise the agreement between simulated and detected distributions of radio pulsars, the simulations fail to reproduce the LAT findings: all models under-predict the number of LAT pulsars with high Ė, and they cannot explain the high probability of detecting both the radio and γ-ray beams at high Ė. The beaming factor remains close to 1.0 over 4 decades in Ė evolution for the slot gap whereas it significantly decreases with increasing age for the outer gaps. The evolution of the enhanced slot-gap luminosity with Ė is compatible with the large dispersion of γ-ray luminosity seen in the LAT data. The stronger evolution predicted for the outer gap, which is linked to the polar cap heating by the return current, is apparently not supported by the LAT data. The LAT sample of γ-ray pulsars therefore provides a fresh perspective on the early evolution of the luminosity and beam width of the γ-ray emission from young pulsars, calling for thin and more luminous gaps.
Cyclotron decay and absorption rates have been well studied in the literature, focusing primarily on spectral, angular, and polarization dependence. Astrophysical applications usually do not require ...retention of information on the electron spin state, and these are normally averaged in obtaining the requisite rates. In magnetic fields, higher order quantum processes such as Compton scattering become resonant at the cyclotron frequency and its harmonics, with the resonances being formally divergent. Such divergences are usually eliminated by accounting for the finite lifetimes of excited Landau states. This practice requires the use of spin-dependent cyclotron rates in order to obtain accurate determinations of process rates very near cyclotronic resonances, the phase-space domain most relevant for certain applications to pulsar models. This paper develops previous results in the literature to obtain compact analytic expressions for cyclotron decay rates/widths in terms of a series of Legendre functions of the second kind; these expressions can be used expediently in astrophysical models. The rates are derived using two popular eigenstate formalisms, namely, that due to Sokolov & Ternov and that due to Johnson & Lippmann. These constitute two sets of eigenfunctions of the Dirac equation that diagonalize different operators and accordingly yield different spin-dependent cyclotron rates. This paper illustrates the attractive Lorentz transformation characteristics of the Sokolov & Ternov formulation, which is another reason why it is preferable when electron spin information must be explicitly retained.
We present results of a population synthesis of millisecond pulsars from the Galactic disk. Excluding globular clusters, we model the spatial distribution of millisecond pulsars by assuming their ...birth in the Galactic disk with a random kick velocity and evolve them to the present within the Galactic potential. We assume that normal and millisecond pulsars are standard candles described with a common radio luminosity model that invokes a new relationship between radio core and cone emission suggested by recent studies. In modeling the radio emission beams, we explore the relativistic effects of time delay, aberration, and sweep-back of the open field lines. While these effects are essential to understanding pulse profiles, the phase-averaged flux is adequately described without a relativistic model. We use a polar cap acceleration model for the gamma -ray emission. We present the preliminary results of our recent study and the implications for observing millisecond pulsars with GLAST and AGILE.
Non-thermal quiescent X-ray emission extending between 10 keV and around 150 keV has been seen in about 10 magnetars by RXTE, INTEGRAL, Suzaku, NuSTAR and Fermi-GBM. For inner magnetospheric models ...of such hard X-ray signals, inverse Compton scattering is anticipated to be the most efficient process for generating the continuum radiation, because the scattering cross section is resonant at the cyclotron frequency. We present hard X-ray upscattering spectra for uncooled monoenergetic relativistic electrons injected in inner regions of pulsar magnetospheres. These model spectra are integrated over bundles of closed field lines and obtained for different observing perspectives. The spectral turnover energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the turnovers inferred in magnetar hard X-ray tails. Electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulse phases. Our spectral computations use a new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.
We present results of a pulsar population synthesis study that incorporates a number of recent developments and some significant improvements over our previous study. We have included the results of ...the Parkes multibeam pulsar survey in our select group of nine radio surveys, doubling our sample of radio pulsars. More realistic geometries for the radio and gamma-ray beams are included in our Monte Carlo computer code, which simulates the characteristics of the Galactic population of radio and gamma-ray pulsars. We adopted with some modifications the radio-beam geometry of Arzoumanian, Chernoff, and Cordes. For the gamma-ray beam, we have assumed the slot gap geometry described in the work of Muslimov and Harding. To account for the shape of the distribution of radio pulsars in the image-P diagram, we continue to find that decay of the magnetic field on a timescale of 2.8 Myr is needed. With all nine surveys, our model predicts that EGRET should have seen seven radio-quiet (below the sensitivity of these radio surveys) and 19 radio-loud gamma-ray pulsars. AGILE (nominal sensitivity map) is expected to detect 13 radio-quiet and 37 radio-loud gamma-ray pulsars, while GLAST, with greater sensitivity, is expected to detect 276 radio-quiet and 344 radio-loud gamma-ray pulsars. When the Parkes multibeam pulsar survey is excluded, the ratio of radio-loud to radio-quiet gamma-ray pulsars decreases, especially for GLAST. The decrease for EGRET is 45%, implying that some fraction of EGRET unidentified sources are radio-loud gamma-ray pulsars. In the radio geometry adopted, short-period pulsars are core dominated. Unlike the EGRET gamma-ray pulsars, our model predicts that when two gamma-ray peaks appear in the pulse profile, a dominant radio core peak appears in between the gamma-ray peaks. Our findings suggest that further improvements are required in describing both the radio and gamma- ray geometries.