A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

We present the results of a population synthesis of radio and gamma-ray millisecond pulsars (MSPs) from the Galactic Disk (GD). Using 92 radio millisecond pulsars detected in 13 radio surveys and 54 Fermi MSPs detected as point sources in the first point source catalog, we establish six free parameters corresponding to the overall factor and the exponents of the period and period derivative dependence for each of the radio and gamma-ray empirical luminosity models. We test three high-energy emission models described by the two pole caustic Slot Gap, Outer Gap, and Pair Starved Polar Cap geometries. The simulated distributions of pulsar properties adequately describe the distributions of detected MSPs from the GD. We explore the gamma-ray emission from groups of MSPs in globular clusters and in the Galactic Bulge. The simulation predicts reasonable numbers of Fermi MSPs detected in the other point source catalogs and anticipates a bright future for Fermi observations of MSPs, expecting a total of ~ 170 MSP detections from the GD within ten years. Our numbers of simulated MSPs in globular clusters are in agreement with those derived from Fermi detections. The simulation predicts about 11,000 MSPs in the Galactic Bulge are required to explain the gamma-ray Galactic Center Excess.

Hard, non-thermal, persistent pulsed X-ray emission extending between 10 keV and $\sim 150$ keV has been observed in nearly ten magnetars. For inner-magnetospheric models of such emission, resonant inverse Compton scattering of soft thermal photons by ultra-relativistic charges is the most efficient production mechanism. We present angle-dependent upscattering spectra and pulsed intensity maps for uncooled, relativistic electrons injected in inner regions of magnetar magnetospheres, calculated using collisional integrals over field loops. Our computations employ a new formulation of the QED Compton scattering cross section in strong magnetic fields that is physically correct for treating important spin-dependent effects in the cyclotron resonance, thereby producing correct photon spectra. The spectral cut-off energies are sensitive to the choices of observer viewing geometry, electron Lorentz factor, and scattering kinematics. We find that electrons with energies $\lesssim 15$ MeV will emit most of their radiation below 250 keV, consistent with inferred turnovers for magnetar hard X-ray tails. More energetic electrons still emit mostly below 1 MeV, except for viewing perspectives sampling field line tangents. Pulse profiles may be singly- or doubly-peaked dependent upon viewing geometry, emission locale, and observed energy band. Magnetic pair production and photon splitting will attenuate spectra to hard X-ray energies, suppressing signals in the Fermi-LAT band. The resonant Compton spectra are strongly polarized, suggesting that hard X-ray polarimetry instruments such as X-Calibur, or a future Compton telescope, can prove central to constraining model geometry and physics.

Pulsed non-thermal quiescent emission between 10 keV and around 150 keV has been observed in $\sim10$ magnetars. For inner magnetospheric models of such hard X-ray signals, resonant Compton upscattering of soft thermal photons from the neutron star surface is the most efficient radiative process. We present angle-dependent hard X-ray upscattering model spectra for uncooled monoenergetic relativistic electrons. The spectral cut-off 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 observed turnovers in magnetar hard X-ray tails. Moreover, electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulses phases. Our spectral computations use new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.

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.

Nearly 1/3 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 |b| > 5 degrees. 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 |b| > 5 degrees to 106. 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 gamma-ray pulsar photons collected by the Fermi Large Area Telescope since June 2008, it is now possible to constrain gamma-ray geometrical models by comparing simulated and observed light-curve morphological characteristics. We assumed vacuum-retarded dipole pulsar magnetic field and tested simulated and observed morphological light-curve characteristics in the framework of two pole emission geometries, Polar Cap (PC), radio, and Slot Gap (SG), and Outer Gap (OG)/One Pole Caustic (OPC) emission geometries. We compared simulated and observed/estimated light-curve morphological parameters as a function of observable and non-observable pulsar parameters. The PC model gives the poorest description of the LAT pulsar light-curve morphology. The OPC best explains both the observed gamma-ray peak multiplicity and shape classes. The OPC and SG models describe the observed gamma-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 gamma-ray emission regions are likely to be located in the outer magnetosphere. 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 between observations and the OPC model 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.

Pair cascades from millisecond pulsars (MSPs) may be a primary source of Galactic electrons and positrons that contribute to the increase in positron flux above 10 GeV as observed by PAMELA and AMS-02. The Fermi Large Area Telescope (LAT) has increased the number of detected gamma-ray MSPs tremendously. Light curve modelling furthermore favours abundant pair production in MSP magnetospheres, so that models of primary cosmic-ray positrons from pulsars should include the contribution from the larger numbers of MSPs and their potentially higher positron output per source. We model the contribution of Galactic MSPs to the terrestrial cosmic-ray electron / positron flux by using a population synthesis code to predict the source properties of present-day MSPs. We simulate pair spectra assuming an offset-dipole magnetic field which boosts pair creation rates. We also consider positrons and electrons that have additionally been accelerated to very high energies in the strong intrabinary shocks in black widow (BW) and redback (RB) binary systems. We transport these particles to Earth by calculating their diffusion and the radiative energy losses they suffer in the Galaxy using a model. Our model particle flux increases for non-zero offsets of the magnetic polar caps. We find that pair cascades from MSP magnetospheres contribute only modestly around a few tens of GeV to the measured fluxes. BW and RB fluxes may reach a few tens of percent of the observed flux up to a few TeV. Future observations should constrain the source properties in this case.

Observations by the Fermi Large Area Telescope of gamma-ray millisecond pulsar 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 millisecond pulsars, impacting Galactic stellar population models. We investigate the contribution of Galactic millisecond pulsars to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day millisecond pulsars. 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 and redback systems. Since millisecond pulsars 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 millisecond pulsars are only modest contributors around a few tens of GeV to the lepton fluxes measured by AMS-02, PAMELA, and Fermi, after which this component cuts off. The contribution by black widows and redbacks may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.

Pulsars are believed to be sources of relativistic electrons and positrons. The abundance of detections of gamma-ray millisecond pulsars by Fermi Large Area Telescope coupled with their light curve characteristics that imply copious pair production in their magnetospheres, motivated us to investigate this old pulsar population as a source of Galactic electrons and positrons and their contribution to the enhancement in cosmic-ray positron flux at GeV energies. We use a population synthesis code to predict the source properties (number, position, and power) of the present-day Galactic millisecond pulsars, taking into account the latest Fermi and radio observations to calibrate the model output. Next, we simulate pair cascade spectra from these pulsars using a model that invokes an offset-dipole magnetic field. We assume free escape of the pairs from the pulsar environment. We then compute the cumulative spectrum of transported electrons and positrons at Earth, following their diffusion and energy losses as they propagate through the Galaxy. Our results indicate that the predicted particle flux increases for non-zero offsets of the magnetic polar caps. Comparing our predicted spectrum and positron fraction to measurements by AMS-02, PAMELA, and Fermi, we find that millisecond pulsars are only modest contributors at a few tens of GeV, after which this leptonic spectral component cuts off. The positron fraction is therefore only slightly enhanced above 10 GeV relative to a background flux model. This implies that alternative sources such as young, nearby pulsars and supernova remnants should contribute additional primary positrons within the astrophysical scenario.

In more than four years of observation the Large Area Telescope on board the Fermi satellite has identified pulsed $\gamma$-ray emission from more than 80 young pulsars, providing light curves with high statistics. Fitting the observations with geometrical models can provide estimates of the magnetic obliquity $\alpha$ and aspect angle $\zeta$, yielding estimates of the radiation beaming factor and luminosity. Using $\gamma$-ray emission geometries (Polar Cap, Slot Gap, Outer Gap, One Pole Caustic) and radio emission geometry, we fit $\gamma$-ray light curves for 76 young pulsars and we jointly fit their $\gamma$-ray plus radio light curves when possible. We find that a joint radio plus $\gamma$-ray fit strategy is important to obtain ($\alpha$, $\zeta$) estimates that can explain simultaneous radio and $\gamma$-ray emission. The intermediate-to-high altitude magnetosphere models, Slot Gap, Outer Gap, and One pole Caustic, are favoured in explaining the observations. We find no evolution of $\alpha$ on a time scale of a million years. For all emission geometries our derived $\gamma$-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. For all models, the correlation between $\gamma$-ray luminosity and spin-down power is consistent with a square root dependence. The $\gamma$-ray luminosities obtained by using our beaming factors 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 that is consistent with the theoretical prediction for the Slot Gap model.

The quantum electrodynamical (QED) process of Compton scattering in strong magnetic fields is commonly invoked in atmospheric and inner magnetospheric models of x-ray and soft gamma-ray emission in high-field pulsars and magnetars. A major influence of the field is to introduce resonances at the cyclotron frequency and its harmonics, where the incoming photon accesses thresholds for the creation of virtual electrons or positrons in intermediate states with excited Landau levels. At these resonances, the effective cross section typically exceeds the classical Thomson value by over 2 orders of magnitude. Near and above the quantum critical magnetic field of 44.13 TeraGauss, relativistic corrections must be incorporated when computing this cross section. This paper presents formalism for the QED magnetic Compton differential cross section valid for both subcritical and supercritical fields, yet restricted to scattered photons that are below pair creation threshold. Calculations are developed for the particular case of photons initially propagating along the field, mathematically simple specializations that are germane to interactions involving relativistic electrons frequently found in neutron star magnetospheres. This exposition of relativistic, quantum, magnetic Compton cross sections treats electron spin dependence fully, since this is a critical feature for describing the finite decay lifetimes of the intermediate states. The formalism employs both the Johnson and Lippmann (JL) wave functions and the Sokolov and Ternov (ST) electron eigenfunctions of the magnetic Dirac equation. The ST states are formally correct for self-consistently treating spin-dependent effects that are so important in the resonances. Relatively compact analytic forms for the cross sections are presented that will prove useful for astrophysical modelers.

With the large sample of young gamma-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 gamma-ray gap models (Polar Cap, Slot Gap, Outer Gap, and One Pole Caustic) normalizing to the number of detected radio pulsars in select group of 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 beams requires an increase by a factor of ~10 of the predicted luminosity to produce a reasonable number of gamma-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. Standard distributions in birth location and pulsar spin-down power (Edot) fail to reproduce the LAT findings: all models under-predict the number of LAT pulsars with high Edot, and they cannot explain the high probability of detecting both the radio and gamma-ray beams at high Edot. The beaming factor remains close to 1 over 4 decades in Edot evolution for the slot gap whereas it significantly decreases with increasing age for the outer gaps. The evolution of the slot-gap luminosity with Edot is compatible with the large dispersion of gamma-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 gamma-ray pulsars therefore provides a fresh perspective on the early evolution of the luminosity and beam width of the gamma-ray emission from young pulsars, calling for thin and more luminous gaps.

We present preliminary results of a pulsar population synthesis of normal pulsars from the Galactic disk using a Markov Chain Monte Carlo method to better understand the parameter space of the assumed model. We use the Kuiper test, similar to the Kolmogorov-Smirnov test, to compare the cumulative distributions of chosen observables of detected radio pulsars with those simulated for various parameters. Our code simulates pulsars at birth using Monte Carlo techniques and evolves them to the present assuming initial spatial, kick velocity, magnetic field, and period distributions. Pulsars are spun down to the present, given radio and gamma-ray emission characteristics, filtered through ten selected radio surveys, and a {\it Fermi} all-sky threshold map. Each chain begins with a different random seed and searches a ten-dimensional parameter space for regions of high probability for a total of one thousand different simulations before ending. The code investigates both the "large world" as well as the "small world" of the parameter space. We apply the K-means clustering algorithm to verify if the chains reveal a single or multiple regions of significance. The outcome of the combined set of chains is the weighted average and deviation of each of the ten parameters describing the model. While the model reproduces reasonably well the detected distributions of normal radio pulsars, it does not replicate the predicted detected $\dot P - P$ distribution of {\it Fermi} pulsars. The simulations do not produce sufficient numbers of young, high-$\dot E$ pulsars in the Galactic plane.

The extremely efficient process of resonant Compton upscattering by relativistic electrons in high magnetic fields is believed to be a leading emission mechanism of high field pulsars and magnetars in the production of intense X-ray radiation. New analytic developments for the Compton scattering cross section using Sokolov & Ternov (S&T) states with spin-dependent resonant widths are presented. These new results display significant numerical departures from both the traditional cross section using spin-averaged widths, and also from the spin-dependent cross section that employs the Johnson & Lippmann (J&L) basis states, thereby motivating the astrophysical deployment of this updated resonant Compton formulation. Useful approximate analytic forms for the cross section in the cyclotron resonance are developed for S&T basis states. These calculations are applied to an inner magnetospheric model of the hard X-ray spectral tails in magnetars, recently detected by RXTE and INTEGRAL. Relativistic electrons cool rapidly near the stellar surface in the presence of intense baths of thermal X-ray photons. We present resonant Compton cooling rates for electrons, and the resulting photon spectra at various magnetospheric locales, for magnetic fields above the quantum critical value. These demonstrate how this scattering mechanism has the potential to produce the characteristically flat spectral tails observed in magnetars.

For inner magnetospheric models of hard X-ray and gamma-ray emission in high-field pulsars and magnetars, resonant Compton upscattering is anticipated to be the most efficient process for generating continuum radiation. This is due in part to the proximity of a hot soft photon bath from the stellar surface to putative radiation dissipation regions in the inner magnetosphere. Moreover, because the scattering process becomes resonant at the cyclotron frequency, the effective cross section exceeds the classical Thomson value by over two orders of magnitude, thereby enhancing the efficiency of continuum production and the cooling of relativistic electrons. This paper presents computations of the electron cooling rates for this process, which are needed for resonant Compton models of non-thermal radiation from such highly-magnetized pulsars. The computed rates extend previous calculations of magnetic Thomson cooling to the domain of relativistic quantum effects, sampled near and above the quantum critical magnetic field of 44.13 TeraGauss. This is the first exposition of fully relativistic, quantum magnetic Compton cooling rates for electrons, and it employs both the traditional Johnson and Lippman cross section, and a newer Sokolov and Ternov (ST) formulation of Compton scattering in strong magnetic fields. Such ST formalism is formally correct for treating spin-dependent effects that are important in the cyclotron resonance, and has not been addressed before in the context of cooling by Compton scattering. The QED effects are observed to profoundly lower the rates below extrapolations of the familiar magnetic Thomson results, as expected, when recoil and Klein-Nishina reductions become important.

We compare population synthesis results for inner and outer magnetosphere emission models with the various characteristics measured in the first LAT pulsar catalogue for both the radio-loud and radio-weak or radio-quiet gamma-ray pulsars. We show that all models fail to reproduce the observations: for each model there is a lack of luminous and energetic objects that suggest a non dipolar magnetic field structure or spin-down evolution. The large dispersion that we find in the simulated gamma-ray luminosity versus spin-down power relation does not allow to use the present trend seen in the Fermi data to distinguish among models. For each model and each Fermi detected pulsar, we have generated light curves as a function of obliquity and inclination angles. The theoretical curves were fitted to the observed one, using a maximum-likelihood approach, to derive the best-fit orientations and to compare how well each model can reproduce the data. Including the radio light-curve gives an additional key constraint to restrict the orientation space

We present results of our pulsar population synthesis of normal pulsars from the Galactic disk using our previously developed computer code. On the same footing, we use slot gap and outer gap models for gamma-ray emission from normal pulsars to obtain statistics of radio-loud and radio-quiet gamma-ray pulsars. From recently improved understanding of HII and star forming regions in the Galaxy, we develop a new surface density model of the birth location of neutron stars. We explore models of neutron star evolution with magnetic field-decay, and with different initial period and magnetic field distributions. We present preliminary results including simulated population statistics that are compared with recent detections by Fermi of normal, isolated pulsars.

The Fermi Gamma-Ray Space Telescope has discovered many gamma-ray pulsars, both as radio-loud objects and radio-quiet or radio-weak pulsars that have been identified through blind period searches. The latter presumably have gamma-ray beams originating from high altitudes in the magnetosphere, resulting in little or no overlap with the radio beam. The exponential cut-off of the emission at high energy also points to a medium- or high-altitude origin of the gamma rays. Population synthesis offers means to study two complementary aspects of the gamma-ray pulsar population: energetics and light curve morphology. The latter (peak multiplicity, position, separation, asymmetries, ...) can be used to constrain the beam(s) geometry (their origin within the open magnetosphere, evolution of the radiating region with pulsar power and age, ...). It can also be used to constrain the obliquity and orientation of a pulsar in the framework of a given accelerator model. We present preliminary results from new simulations and light-curve analysis that indicate that the sample of Fermi pulsars retains much of the geometrical information of the parent gamma-ray pulsar population.

We present results of our pulsar population synthesis of normal pulsars from the Galactic disk using our previously developed Monte-Carlo code. From our studies of observed 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 luminosity and ratio of radio core-to-cone peak fluxes. In this model, short period pulsars are more cone-dominated. We explore models of neutron star evolution with and without magnetic field decay, and with different initial period distributions. We present preliminary results including simulated population statistics that are compared with the observed radio pulsar population. The evolved neutron star populations resulting from this simulation can be used to model distributions of gamma-ray pulsars for comparison to Fermi results.

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 sweepback of the open field lines. While these effects are essential in 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.

Radio-quiet gamma-ray pulsars like Geminga may account for a number of the unidentified EGRET sources in the Galaxy. The number of Geminga-like pulsars is very sensitive to the geometry of both the gamma-ray and radio beams. Recent studies of the shape and polarization of pulse profiles of young radio pulsars have provided evidence that their radio emission originates in wide cone beams at altitudes that are a significant fraction (1 -10%) of their light cylinder radius. Such wide radio emission beams will be visible at a much larger range of observer angles than the narrow core components thought to originate at lower altitude. Using 3D geometrical modeling that includes relativistic effects from pulsar rotation, we study the visibility of such radio cone beams as well as that of the gamma-ray beams predicted by slot gap and outer gap models. From the results of this study one can obtain revised predictions for the fraction of Geminga-like, radio quiet pulsars present in the gamma-ray pulsar population.

We present results of our pulsar population synthesis of ordinary isolated 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 $\gamma$-ray emission from both millisecond and ordinary 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.

Our recent studies of pulsar population statistics suggest that improvements of radio and gamma-ray beam geometry and luminosity models require further refinement. The goal of this project is to constrain the viewing geometry for some radio pulsars, especially three-peaked pulse profiles, in order to limit the uncertainty of the magnetic inclination and impact angles. We perform fits of the pulse profile and position angle sweep of radio pulsars for the available frequencies. We assume a single core and conal beams described by Gaussians. We incorporate three different size cones with frequency dependence from the work of Mitra & Deshpande (1999). We obtain separate spectral indices for the core and cone beams and explore the trends of the ratio of core to cone peak fluxes. This ratio is observed to have some dependence with period. However, we cannot establish the suggested functional form of this ratio as indicated by the work of Arzoumanian, Chernoff & Cordes (2002).

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 expediently used in astrophysical models. The rates are derived using two popular eigenstate formalisms, namely that due to Sokolov and Ternov, and that due to Johnson and 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 and Ternov formulation, which is another reason why it is preferable when electron spin information must be explicitly retained.

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 multi-beam pulsar survey in our select group of nine radio surveys, doubling our sample of radio pulsars. We adopted with some modifications the radio beam geometry of Arzoumanian, Chernoff & Cordes (2002). For the $\gamma$-ray beam, we have assumed the slot gap geometry described in the work of Muslimov & Harding (2003). To account for the shape of the distribution of radio pulsars in the $\dot P-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 7 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 multi-beam 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.

The EGRET telescope on the Compton Gamma-Ray Observatory detected over 200 sources and the majority of these are still unidentified. At least three subpopulations of EGRET sources have been associated with the Galaxy: bright sources lying along the Galactic plane, weaker sources spatially correlated with the Gould Belt and a high-latitude, halo population. Many of these sources may be pulsars and there are more than two-dozen radio pulsars in or near EGRET source error boxes, most of them recently discovered in the Parkes Multibeam Survey. We present results from several population synthesis studies of pulsars in the Galaxy, which predict the number of pulsars detected by gamma-ray and radio surveys assuming models for the high-energy emission beam and its relation to the radio beam. Future gamma-ray pulsar detections by AGILE and GLAST together with the recent large rise in the radio pulsar population will give greatly improved statistics. The relative numbers of radio and gamma-ray pulsars detected in the plane and in the Gould Belt populations will provide an important discriminator between models of high-energy emission and will profoundly affect the numbers and birth rate of pulsars in the Galaxy.

We simulate the characteristics of the Galactic population of radio and $\gamma$-ray pulsars using Monte Carlo techniques. At birth, neutron stars are spatially distributed in the Galactic disk, with supernova-kick velocities, and randomly dispersed in age back to $10^9$ years. They are evolved in the Galactic gravitational potential to the present time. From a radio luminosity model, the radio flux is filtered through a selected set of radio-survey parameters. $\gamma$-ray luminosities are assigned using the features of recent polar cap acceleration models invoking space-charge-limited flow, and a pulsar death valley further attenuates the population of radio-loud pulsars. Assuming a simple emission geometry with aligned radio and $\gamma$-ray beams of 1 steradian solid angle, our model predicts that EGRET should have seen 7 radio-loud and 1 radio-quiet, $\gamma$-ray pulsars. With much improved sensitivity, GLAST, on the other hand, is expected to observe 76 radio-loud and 74 radio-quiet, $\gamma$-ray pulsars of which 7 would be identified as pulsed sources. We also explore the effect of magnetic field decay on the characteristics of the radio and $\gamma$-ray pulsar populations. Including magnetic field decay on a timescale of 5 Myr improves agreement with the radio pulsar population and increases the predicted number of GLAST detected pulsars to 90 radio-loud and 101 radio-quiet (9 pulsed) $\gamma$-ray pulsars. The lower flux threshold allows GLAST to detect $\gamma$-ray pulsars at larger distances than those observed by the radio surveys used in this study.

This paper explores the effects of strong magnetic fields on the Compton scattering of relativistic electrons. Recent studies of upscattering and energy loss by relativistic electrons that have used the non-relativistic, magnetic Thomson cross section for resonant scattering or the Klein-Nishina cross section for non-resonant scattering do not account for the relativistic quantum effects of strong fields ($ > 4 \times 10^{12}$ G). We have derived a simplified expression for the exact QED scattering cross section for the broadly-applicable case where relativistic electrons move along the magnetic field. To facilitate applications to astrophysical models, we have also developed compact approximate expressions for both the differential and total polarization-dependent cross sections, with the latter representing well the exact total QED cross section even at the high fields believed to be present in environments near the stellar surfaces of Soft Gamma-Ray Repeaters and Anomalous X-Ray Pulsars. We find that strong magnetic fields significantly lower the Compton scattering cross section below and at the resonance, when the incident photon energy exceeds $m_ec^2$ in the electron rest frame. The cross section is strongly dependent on the polarization of the final scattered photon. Below the cyclotron fundamental, mostly photons of perpendicular polarization are produced in scatterings, a situation that also arises above this resonance for sub-critical fields. However, an interesting discovery is that for super-critical fields, a preponderance of photons of parallel polarization results from scatterings above the cyclotron fundamental. This characteristic is both a relativistic and magnetic effect not present in the Thomson or Klein-Nishina limits.

Gamma-rays from pulsars can be efficiently attenuated in their magnetospheres via the mechanism of single-photon pair production and also the exotic QED process of photon splitting, which become prolific in fields approaching the quantum critical value of $B_{cr}=4.41\times 10^{13}$ Gauss. Recently we have published results of our modelling of strongly-magnetized $\gamma$-ray pulsars, which focused on the escape or attenuation of photons emitted near the pole at the neutron star surface in dipole fields, in a Schwarzschild metric. We found that pair production and splitting totally inhibit emission above around 10--30 MeV in PSR1509-58, whose surface field is inferred to be as high as $0.7B_{cr}$. Our model pulsar spectra are consistent with the EGRET upper limits for PSR1509-58 for a wide range of polar cap sizes. Here we review the principal predictions of our attenuation analysis, and identify how its powerful observational diagnostic capabilities relate to current and future gamma-ray experiments. Diagnostics include the energy of the gamma-ray turnover and the spectral polarization, which constrain the estimated polar cap size and field strength, and can determine the relative strength of splitting and pair creation.

Magnetic photon splitting, a QED process that becomes important only in magnetic fields approaching the quantum critical value, B_cr = 4.41 X 10^13 Gauss, is investigated as a mechanism for attenuation of gamma-rays emitted near the surface of strongly-magnetized pulsars. We model photon splitting attenuation and subsequent splitting cascades in gamma-ray pulsars, including the dipole field and curved spacetime geometry of the neutron star magnetosphere. We focus specifically on PSR1509-58, which has the highest surface magnetic field of all the gamma-ray pulsars (B_0 = 3 X 10^13 Gauss). We find that splitting will not be important for most gamma-ray pulsars, i.e. those with B_0 <~ 0.2 B_cr, but will be important for gamma-ray pulsars having B_0 >~ 0.3 B_cr, where the splitting attenuation lengths and escape energies become comparable to or less than those for pair production. We compute Monte Carlo spectral models for PSR1509-58. We find that photon splitting, or combined splitting and pair production, can explain the unusually low cutoff energy (between 2 and 30 MeV) of PSR1509-58, and that the model cascade spectra, which display strong polarization, are consistent with the observed spectral points and upper limits for polar cap emission at a range of magnetic colatitudes up to ~ 25 degrees.