An overview of the results of observations for the transient X-ray pulsar 4U 0115+63, a member of a binary system with a Be star, since its discovery to the present day (~40 years) based on data from more than dozen observatories and instruments is presented. A overall light curve and the history of change in the spin frequency of the neutron star over the entire history of its observations, which also includes the results of recent measurements made by the INTEGRAL observatory during the 2004, 2008, and 2011 outbursts, are provided. The source's energy spectra have also been constructed from the INTEGRAL data obtained during the 2011 outburst for a dynamic range of its luminosities (~10^{37} - 7 x 10^{37} erg/s). We show that apart from the fundamental harmonic of the cyclotron absorption line at energy ~11 keV, its four higher harmonics at energies ~24, 35.6, 48.8, and 60.7 keV are detected in the spectrum. We have performed a detailed analysis of the source's spectra in the 4-28 keV energy band based on all of the available RXTE archival data obtained during bright outbursts in 1995-2011. We have confirmed that modifying the source's continuum model can lead to the disappearance of the observed anticorrelation between the energy of the fundamental harmonic of the cyclotron absorption line and the source's luminosity. Thus, the question about the evolution of the cyclotron absorption line energy with the luminosity of the X-ray pulsar 4U 0115+63 remains open and a physically justified radiation model for X-ray pulsars is needed to answer it.

We present a new web-tool - Net-PSICoNS - for population synthesis of isolated near-by cooling neutron stars (NSs). The main aim is to provide an easy test of models of the NS thermal evolution which can be used by groups studying this subject. A user can upload cooling curves for a set of masses, modify the mass spectrum if necessary, change radii to fit the EoS used for cooling curve calculations, and then a population synthesis of close-by isolated cooling NSs is performed. The output includes the Log N -- Log S distribution confronted with the ROSAT observations and several other sets of data. In this paper, we summarize the idea of the test proposed by Popov et al. (2006), and present a user's manual for the web-tool.

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 evolution 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 WD's birth rate and upper limit of the X-ray spectrum temperature: DA birth rate $= 0.61\times 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 perform population synthesis studies of different types of neutron stars taking into account the magnetic field decay. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G]) > \sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible. We apply the obtained field distribution and the model of decay to study long-term evolution of neuton stars till the stage of accretion from the interstellar medium. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0 > \sim 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.

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 taking 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_0\ga 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.

Context: Isolated cooling neutron stars with thermal X-ray emission remain rarely detected objects despite many searches investigating the ROSAT data. Aims: We simulate the population of close-by young cooling neutron stars to explain the current observational results. Given the inhomogeneity of the neutron star distribution on the sky it is particularly interesting to identify promising sky regions with respect to on-going and future searches. Methods: Applying a population synthesis model the inhomogeneity of the progenitor distribution and the inhomogeneity of the X-ray absorbing interstellar medium are considered for the first time. The total number of observable neutron stars is derived with respect to ROSAT count rates. In addition, we present sky maps of neutron star locations and discuss age and distance distributions of the simulated neutron stars. Implications for future searches are discussed. Results: With our advanced model we can successfully explain the observed logN - logS distribution of close-by neutron stars. Cooling neutron stars will be most abundant in the directions of rich OB associations. New candidates are expected to be identified behind the Gould Belt, in particular in the Cygnus-Cepheus region. They are expected to be on average younger and then hotter than the known population of isolated cooling neutron stars. In addition, we propose to use data on runaway stars to search for more radio-quiet cooling neutron stars.

We perform population synthesis studies of different types of neutron stars (thermally emitting isolated neutron stars, normal radio pulsars, magnetars) taking into account the magnetic field decay and using results from the most recent advances in neutron star cooling theory. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. For this purpose, we fix a baseline neutron star model (all microphysics input), and other relevant parameters to standard values (velocity distribution, mass spectrum, birth rates ...), allowing to vary the initial magnetic field strength. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G])>\sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life (only about a factor of 2 for low field neutron stars but more than an order of magnitude for magnetars), and a mass distribution function dominated by low mass objects. This model explains satisfactorily all known populations. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible.