2007ApJ...671..713S

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 γ-ray emission. We present the preliminary results of our recent study and the implications for observing millisecond pulsars with GLAST and AGILE.