  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.
