Magnetic pair creation $\gamma \to e^+e^-$ has been at the core of radio pulsar paradigms and central to polar cap models of gamma-ray pulsars for over three decades. The Fermi gamma-ray pulsar population now exceeds 140 sources and has defined an important part of Fermi's science legacy. Among the population characteristics well established is the common occurrence of exponential turnovers in their spectra in the 1--10 GeV range. These turnovers are too gradual to arise from magnetic pair creation in the strong magnetic fields of pulsar inner magnetospheres. By demanding insignificant photon attenuation precipitated by such single-photon pair creation, the energies of these turnovers for Fermi pulsars can be used to compute lower bounds for the typical altitude of GeV band emission. This paper explores such pair transparency constraints below the turnover energy, and updates earlier altitude bound determinations of that have been deployed in various Fermi pulsar papers. For low altitude emission locales, general relativistic influences are found to be important. Rotational aberration influences are also explored, and are found to be small at low altitudes, except near the magnetic pole. The altitude bounds are typically in the range of 2-7 stellar radii for the young Fermi pulsar population, and provide key information on the emission altitude in radio quiet pulsars that do not possess double-peaked pulse profiles. The bound for the Crab pulsar is at a much higher altitude, with the putative detection by MAGIC out to 350-400 GeV implying a lower bound of 310km to the emission region, i.e., approximately 20% of the light cylinder radius. These results are also extended to the super-critical field domain, where it is found that emission in magnetars originating below around 10 stellar radii will not appear in the Fermi-LAT band.

We present new photometric and spectroscopic observations for 2M 1533+3759 (= NSVS 07826147). It has an orbital period of 0.16177042 day, significantly longer than the 2.3--3.0 hour periods of the other known eclipsing sdB+dM systems. Spectroscopic analysis of the hot primary yields Teff = 29230 +/- 125 K, log g = 5.58 +/- 0.03 and log N(He)/N(H) = -2.37 +/- 0.05. The sdB velocity amplitude is K1 = 71.1 +/- 1.0 km/s. The only detectable light contribution from the secondary is due to the surprisingly strong reflection effect. Light curve modeling produced several solutions corresponding to different values of the system mass ratio, q(M2/M1), but only one is consistent with a core helium burning star, q=0.301. The orbital inclination is 86.6 degree. The sdB primary mass is M1 = 0.376 +/- 0.055 Msun and its radius is R1 = 0.166 +/- 0.007 Rsun. 2M1533+3759 joins PG0911+456 (and possibly also HS2333+3927) in having an unusually low mass for an sdB star. SdB stars with masses significantly lower than the canonical value of 0.48 Msun, down to as low as 0.30 Msun, were theoretically predicted by Han et al. (2002, 2003), but observational evidence has only recently begun to confirm the existence of such stars. The existence of core helium burning stars with masses lower than 0.40--0.43 Msun implies that at least some sdB progenitors have initial main sequence masses of 1.8--2.0 Msun or more, i.e. they are at least main sequence A stars. The secondary is a main sequence M5 star.

The principal aim of this project is to determine the structural parameters of the rapidly pulsating subdwarf B star PG 0911+456 from asteroseismology. Our work forms part of an ongoing programme to constrain the internal characteristics of hot B subdwarfs with the long-term goal of differentiating between the various formation scenarios proposed for these objects. First comparisons of asteroseismic values with evolutionary theory look promising, however it is clear that more targets are needed for meaningful statistics to be derived. The observational pulsation periods of PG 0911+456 were extracted from rapid time-series photometry using standard Fourier analysis techniques. Supplemented by spectroscopic estimates of the star's mean atmospheric parameters, they were used as a basis for the "forward modelling" approach in asteroseismology. The latter culminates in the identification of one or more "optimal" models that can accurately reproduce the observed period spectrum. This naturally leads to an identification of the oscillations detected in terms of degree l and radial order k, and infers the structural parameters of the target. From the photometry it was possible to extract 7 independent pulsation periods in the 150-200 s range with amplitudes between 0.05 and 0.8 % of the star's mean brightness. An asteroseismic search of parameter space identified several models that matched the observed properties of PG 0911+456 well, one of which was isolated as the "optimal" model on the basis of spectroscopic and mode identification considerations. All the observed pulsations are identified with low-order acoustic modes with degree indices l=0,1,2 and 4, and match the computed periods with a dispersion of only ~0.26 %.

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.

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