Extending the population synthesis method to isolated young cooling white dwarfs we are able to confront our model assumptions with observations made in ROSAT All-Sky Survey (Fleming et al., 1996). This allows us to check model parameters such as evo
lution of spectra and separation of heavy elements in DA WD envelopes. It seems like X-ray spectrum temperature of these objects is given by the formula T_{X-ray} = min(T_eff, T_max). We have obtained DA WDs birth rate and upper limit of the X-ray spectrum temperature: DA birth rate $= 0.61times 10^{-12}$ in cubic parsec per year and T_max = 41000 K. These values are in good correspondence with values obtained by other authors (Liebert et al., 2004; Wolff et al., 1996). From this fact we also conclude that our population synthesis method is applicable to the population of close-by isolated cooling white dwarfs as well as to the population of the isolated cooling neutron stars.
We study evolution of isolated neutron stars on long time scale and calculate distribution of these sources in the main evolutionary stages: Ejector, Propeller, Accretor, and Georotator. We compare different initial magnetic field distributions takin
g into account a possibility of magnetic field decay, and include in our calculations the stage of subsonic Propeller. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0ga 10^{13}$ G) reach the accretion regime within the Galactic lifetime if their kick velocities are not too large. The fact that in previous studies made $>$10 years ago, such objects were not considered results in a slight increase of the Accretor fraction in comparison with earlier conclusions. Most of the neutron stars similar to the Magnificent seven are expected to become accreting from the interstellar medium after few billion years of their evolution. They are the main predecestors of accreting isolated neutron stars.
