It is proposed that sources of repeating fast radio bursts with detected periodicity are magnetars with extremely short initial spin periods at the protoneutron star stage, formed in binaries after tidal synchronization of their progenitor cores on late stages of thermonuclear burning in very close systems with orbital periods about a few days. This formation channel is in several respects different from evolution of progenitors of normal magnetars, and thus, it explains some differences between repeating and non-repeating sources of FRBs.

Fast radio bursts are bright, millisecond-scale radio flashes of yet unknown physical origin. Recently, their extragalactic nature has been demonstrated, and an increasing number of the sources have been found to repeat. Young, highly magnetized, isolated neutron stars - magnetars - have been suggested as the most promising candidates for fast radio burst progenitors owing to their energetics and high X-ray flaring activity. Here we report the detection with the Konus-Wind detector of a hard X-ray event of April 28, 2020, temporarily coincident with a bright, two-peak radio burst with properties remarkably similar to those of fast radio bursts. The source of the radio burst is located in the direction to the Galactic magnetar SGR 1935+2154, which recently entered an active state. We show that a separation between two peaks of the double-peaked X-ray burst is nearly the same as that of the radio peaks, confirming that the X-ray and radio emission most likely have a common origin. Thus, this is the first simultaneous detection of a fast radio burst from a Galactic magnetar and its high-energy counterpart. The total energy emitted in X-rays in this burst is typical of bright short magnetar bursts, but its spiky light curve and an unusual hardness of the energy spectrum distinguish the April 28 event among multiple 'ordinary' flares detected from SGR 1935+2154 previously. This, and a recent non-detection of radio emission from two typical bright soft bursts from the same magnetar, may imply the existence of a special class of magnetar flares capable of producing fast radio bursts.

Lyutikov (2002) predicted "radio emission from soft gamma-ray repeaters (SGRs) during their bursting activity". Detection of a Mega-Jansky radio burst in temporal coincidence with high energy bursts from a Galactic magnetar SGR 1935+2154 confirms that prediction. Similarity of this radio event with Fast Radio Bursts (FRBs) suggests that FRBs are produced within magnetar magnetospheres. We demonstrate that SGR 1935+2154 satisfies the previously derived constraints on the physical parameters at the FRBs' loci. Coherent radio emission is generated in the inner parts of the magnetosphere at $r< 100 R_{\rm NS}$. The radio emission is produced by the yet unidentified plasma emission process, occurring during the initial stages of reconnection events.

The possibility of compensating atmospheric influence in an experiment on precision measurement of gravitational redshift using the "RadioAstron" spacecraft (SC) is discussed. When a signal propagates from a ground-based tracking station to a spacecraft and back, interaction with the ionosphere and troposphere makes considerable contribution to the frequency shift. A brief overview of the physical effects determining this contribution is given, and the principles of calculation and compensation of the corresponding frequency distortions of radio signals are described. Then these approaches are used to reduce the atmospheric frequency shift of the "RadioAstron" spacecraft signal. The spacecraft hardware allows working in two communication modes: "one-way" and "two-way", in addition, two communication channels at different frequencies work simultaneously. "One-way" (SC - ground-based tracking station) communication mode, a signal is synchronized by the on board hydrogen frequency standard. The "two-way" (SC - ground-based tracking station - SC ) mode is synchronized by the ground hydrogen standard. The calculations performed allow us to compare the quality of compensation of atmospheric fluctuations performed by various methods and choose the optimal one.

Modernized version of the combined opto-acoustical gravitational wave detector OGRAN is presented. Located in the deep underground of the Baksan Neutrino Observatory this setup is aimed to work on the program of collapsing stars searching for in multi-channel manner with the neutrino telescope BUST. The both instruments have the sensitivity allowing a registration of collapses in our Galaxy as rare events with the estimated probability 0:03 year-1. The OGRAN narrow band sensitivity at the kilohertz frequency is limited by its acoustical mode thermal noise achieving 10^{-20} in term of metric perturbations. A possible algorithm of the joint data analysis for the both instruments is developed and resulting formulae of the right detection probability are given. A future increasing of the OGRAN sensitivity associated with the moderate cooling (nitrogen temperature) of the acoustical mode is also discussed.

Orbits of close-in planets can shrink significantly due to dissipation of tidal energy in a host star. This process can result in star-planet coalescence within the Galactic lifetime. In some cases, such events can be accompanied by an optical or/and UV/X-ray transient. Potentially, these outbursts can be observed in near future with new facilities such as LSST from distances about few Mpc. We use a population synthesis model to study this process and derive the rate of star-planet mergers of different types. Mostly, planets are absorbed by red giants. However, these events, happening with the rate about 3 per year, mostly do not produce detectable transients. The rate of mergers with main sequence stars depends on the effectiveness of tidal dissipation; for reasonable values of stellar tidal quality factor, such events happen in a Milky Way-like galaxy approximately once in 70 yrs or more rarely. This rate is dominated by planets with low masses. Such events do not produce bright transients having maximum luminosities $\lesssim 10^{36.5}$erg s$^{-1}$. Brighter events, related to massive planets, with maximum luminosity $\sim 10^{37.5}$--$10^{38}$erg s$^{-1}$, have the rate nearly five times smaller.

Coalescences of neutron stars and white dwarfs are relatively frequent phenomena, outnumbering other types of compact object mergers (neutron stars and black holes without involving white dwarfs) altogether. Such event potentially can produce not only optical, but also an X-ray burst. Transient source CDF-S XT2 [1] can be an example of this type of events as suggested by [2]. In this note we estimate the rate of these transients in the field of view of X-ray instruments on-board Spectrum-RG satellite. We demonstrate that during four years of the survey program several thousand of events related to neutron star - white dwarf mergers might appear in the field of view of eROSITA. Collimation of X-ray emission can reduce this number. Smaller, but comparable number of transients is expected in the case of ART-XC telescope. However, due to relatively short duration - $\lesssim 10^4$ s, - mostly such transients might be visible just in one scan of telescopes ($\sim 40$ s), and so only a few photons are expected to be detected which makes definite identification without additional information very problematic.

In this note we propose that recently discovered radio pulsar J0250+5854 with 23.5 sec spin period is presently at the Hall attractor stage. This can explain low temperature and absence of magnetar-like activity of this source together with its spin period and period derivative. We present results of calculations of the evolution of this source in a simple model of magnetic field decay. The neutron star could start its evolution as a magnetar with initial field $\sim 10^{14}-10^{15}$ G for realistic range of parameter $Q$ describing crust imperfections. Future measurements of surface temperature and age of this neutron star might help to probe this hypothesis.

In this note we discuss the possibility of detecting the accompanying X-ray emission from sources of fast radio bursts with the eROSITA telescope onboard the Spektr-RG observatory. It is shown that during four years of the survey program about 300 bursts are expected to appear in the field of view of eROSITA. About 1\% of them will be detected by ground-based radio telescopes. For a total energy release $\sim~10^{46}$~ergs depending on the spectral parameters and absorption in the interstellar and intergalactic media, an X-ray flare can be detected from distances from $\sim 1$~Mpc (thermal spectrum with $kT=200$~keV and strong absorption) up to $\sim1$~Gpc (power-law spectrum with photon index $\Gamma=2 $ and realistic absorption). Thus, eROSITA observations might help to provide important constraints on parameters of sources of fast radio bursts, or may even allow to identify the X-ray transient counterparts, which will help to constrain models of fast radio bursts generation.

Planet migration originally refers to protoplanetary disks, which are more massive and dense than typical accretion disks in binary systems. We study planet migration in an accretion disk in a binary system consisting of a solar-like star hosting a planet and a red giant donor star. The accretion disk is fed by a stellar wind. %, disk self-gravity is neglected. We use the $\alpha$-disk model and consider that the stellar wind is time-dependent. Assuming the disk is quasi-stationary we calculate its temperature and surface density profiles. In addition to the standard disk model, when matter is captured by the disk at its outer edge, we study the situation when the stellar wind delivers matter on the whole disc surface inside the accretion radius with the rate depending on distance from the central star. Implying that a planet experiences classical type I/II migration we calculate migration time for a planet on a circular orbit coplanar with the disk. Potentially, rapid inward planet migration can result in a planet-star merger which can be accompanied by an optical or/and UV/X-ray transient. We calculate timescale of migration for different parameters of planets and binaries. Our results demonstrate that planets can fall on their host stars within the lifetime of the late-type donor for realistic sets of parameters.

We report on the first X-ray observation of the 0.28 s isolated radio pulsar PSR J1154--6250 obtained with the XMM-Newton observatory in February 2018. A point-like source is firmly detected at a position consistent with that of PSR J1154--6250. The two closest stars are outside the 3$\sigma$ confidence limits of the source position and thus unlikely to be responsible for the observed X-ray emission. The energy spectrum of the source can be fitted equally well either with an absorbed power-law with a steep photon index $\Gamma\approx 3.3$ or with an absorbed blackbody with temperature $kT=0.21\pm 0.04$~keV and emitting radius $R_\mathrm{BB} \approx 80$ m (assuming a distance of 1.36~kpc). The X-ray luminosity of $4.4\times 10^{30}$ erg s$^{-1}$ derived with the power-law fit corresponds to an efficiency of $\eta_X = L^\mathrm{unabs}_X/\dot E = 4.5\times 10^{-3}$, similar to those of other old pulsars. The X-ray properties of PSR J1154--6250 are consistent with an old age and suggest that the spatial coincidence of this pulsar with the OB association Cru OB1 is due to a chance alignment.

The phenomenon of fast radio bursts (FRBs) was discovered in 2007. These are powerful (0.1-100 Jy) single radio pulses with durations of several milliseconds, large dispersion measures, and record high brightness temperatures suggesting coherent emission mechanism. As of time of writing, 32 FRBs were recorded. There is also one repeating source from which already hundreds of bursts have been detected. The FRB rate is estimated to amount up to several thousand per day over the sky, and their isotropic sky distribution likely suggests a cosmological origin. Since the discovery, different hypotheses on the possible FRB nature have been proposed, however up to now the origin of these transient events remains obscure. The most prospective models treat them as being related to a bursting emission from magnetars -- neutron stars whose activity is due to their magnetic field dissipation, -- or as being analogs of giant pulses observed from several radio pulsars . Future increase in the statistics of the observed FRBs and an improvement upon the FRB population characteristics will make possible to use them as a new tool to probe the intergalactic medium and to test fundamental physical theories.

We discuss the different signals, in gravitational and electromagnetic waves, emitted during the merger of two compact stars. We will focus in particular on the possible contraints that those signals can provide on the equation of state of dense matter. Indeed, the stiffness of the equation of state and the particle composition of the merging compact stars, strongly affect e.g. the life time of the post-merger remnant and its gravitational wave signal, the emission of the short gamma-ray-burst, the amount of ejected mass and the related kilonova. The first detection of gravitational waves from the merger of two compact stars in August 2017, GW170817, and the subsequent detections of its electromagnetic counterparts, GRB170817A and AT2017gfo, is the first example of the era of "multi-messenger astronomy": we discuss what we have learned from this detection on the equation of state of compact stars and we provide a tentative interpretation of this event, within the two families scenario, as due to the merger of a hadronic star with a quark star.

We briefly review main observational properties of fast radio bursts (FRBs) and discuss two most popular hypothesis for the explanation of these enigmatic intense millisecond radio flashes. FRBs most probably originate on extragalactic distances, and their rate on the sky is about a few thousand per day with fluences above $\sim$~1~Jy~ms (or with fluxes larger than few tenths of Jy). Two leading scenarios describing these events include strong flares of magnetars and supergiant pulses of young radio pulsars with large rotational energy losses, correspondingly. At the moment, it is impossible to choose between these models. However, new telescopes can help to solve the puzzle of FRBs in near future.

SXP 1062 is a Be X-ray binary located in the Small Magellanic Cloud. It hosts a long-period X-ray pulsar and is likely associated with the supernova remnant MCSNR J0127-7332. In this work we present a multi-wavelength view on SXP 1062 in different luminosity regimes. We consider monitoring campaigns in optical (OGLE survey) and X-ray (SWIFT telescope). During these campaigns a tight coincidence of X-ray and optical outbursts is observed. We interpret this as typical Type I outbursts as often detected in Be X-ray binaries at periastron passage of the neutron star. To study different X-ray luminosity regimes in depth, during the source quiescence we observed it with XMM-Newton while Chandra observations followed an X-ray outburst. Nearly simultaneously with Chandra observations in X-rays, in optical the RSS/SALT telescope obtained spectra of SXP 1062. On the basis of our multi-wavelength campaign we propose a simple scenario where the disc of the Be star is observed face-on, while the orbit of the neutron star is inclined with respect to the disc. According to the model of quasi-spherical settling accretion our estimation of the magnetic field of the pulsar in SXP 1062 does not require an extremely strong magnetic field at the present time.

HD 49798/RX J0648.0--4418 is a peculiar X-ray binary with a hot subdwarf (sdO) mass donor. The nature of the accreting compact object is not known, but its spin period $P=13.2$~s and $\dot P =-2.15 \times 10^{-15}$s~s$^{-1}$, prove that it can be only either a white dwarf or a neutron star. The spin-up has been very stable for more than 20 years. We demonstrate that the continuous stable spin-up of the compact companion of HD 49798 can be best explained by contraction of a young white dwarf with an age $\sim 2$~Myrs. This allows us to interpret all the basic parameters of the system in the framework of an accreting white dwarf. We present examples of binary evolution which result in such systems. If correct, this is the first direct evidence for a white dwarf contraction on early evolutionary stages.

Recently, numerical calculations of the magnetic field evolution in neutron stars demonstrated the possible existence of a Hall attractor, a stage at which the evolution of the field driven by the Hall cascade ends. The existence of such a stage in neutron star magnetic evolution is very important, and can be potentially probed by observations. Here we discuss three types of objects which could have reached this stage. First, we briefly describe the evolution of normal radio pulsars with ages about a few hundred thousand years. Then we analyse in more detail observations of RX J1856.5-3754, one of the Magnificent Seven, focusing on the surface temperature distribution and comparing model predictions with the temperature map inferred from X-ray observations. Finally, we discuss the necessity of the Hall attractor stage to explain the hypothetical existence of accreting magnetars. We conclude that at the moment there is no direct confirmation of the Hall attractor stage in known sources. However, more detailed observations in the near future can demonstrate existence (or absence) of this stage of the crustal magnetic field evolution.

Several candidates for accreting magnetars have been proposed recently by different authors. Existence of such systems contradicts the standard magnetic field decay scenario where a large magnetic field of a neutron star reaches $\lesssim$ few$\times 10^{13}$G at ages $\gtrsim 1$ Myr. Among other sources, the high mass X-ray binary 4U0114+65 seems to have a strong magnetic field around $10^{14}$ G. We develop a new Bayesian estimate for the kinematic age and demonstrate that 4U0114+65 has kinematic age 2.4-5 Myr ($95\%$ credential interval) since the formation of the neutron star. We discuss which conditions are necessary to explain the potential existence of magnetars in accreting high-mass binaries with ages about few Myrs and larger. Three necessary ingredients are: the Hall attractor to prevent rapid decay of dipolar field, relatively rapid cooling of the crust in order to avoid Ohmic decay due to phonons, and finally, low values of the parameter $Q$ to obtain long Ohmic time scale due to impurities. If age and magnetic field estimates for proposed accreting magnetars are correct, then these systems set the strongest limit on the crust impurity for a selected sample of neutron stars and provide evidence in favour of the Hall attractor.

The existence of strange quark stars has been proposed many years ago. More recently, the possible co-existence of a first family composed of "normal" neutron stars with a second family of strange quark stars has been proposed as a solution of problems related to the maximum mass and to the minimal radius of these compact stellar objects. In this paper we study the mass distribution of compact objects formed in binary systems and the relative fractions of quark and neutron stars in different subpopulations. We incorporate the strange quark star formation model provided by the two-families scenario and we perform a large-scale population synthesis study in order to obtain the population characteristics. In our model, below a critical gravitational mass $M_\mathrm{max}^H- \Delta M$ only normal (hadron) neutron stars exist. Then in the mass range $(M_\mathrm{max}^H- \Delta M) \leqslant M \leqslant M_\mathrm{max}^H$ strange quark stars and neutron stars coexist. Finally, above $M_\mathrm{max}^H$ all compact objects are strange quark stars. We argue that $M_\mathrm{max}^H$ is in the range $\sim 1.5-1.6 M_\odot$. According to our results, the main channel for strange quark star formation in binary systems is accretion from a secondary companion on a neutron star.This opens the possibility of having explosive GRB-like phenomena not related to supernovae and not due to the merger of two neutron stars. The enhancement in the number of compact objects in the co-existence mass range $(M_\mathrm{max}^H- \Delta M) \leqslant M \leqslant M_\mathrm{max}^H$ is not very pronounced. The number of double strange quark star's systems is rather small with only a tiny fraction which merge within a Hubble time. This drastically limits the flux of strangelets produced by the merger, which turns out to be compatible with all limits stemming from Earth and lunar experiments.

The evolution of magnetic field in isolated neutron stars is one of the most important ingredients in the attempt to build a unified description of these objects. A prediction of field evolution models is the existence of an equilibrium configuration, in which the Hall cascade vanishes. Recent calculations have explored the field structure in this stage, called the Hall attractor. We use X-ray data of near-by, cooling neutron stars to probe this prediction, as these sources are surmised to be close to or at Hall attractor phase. We show that the source RX J1856.5-3754 might be closer to the attractor than other sources of its class. Our modelling indicates that the properties of surface thermal emission, assuming that the star is in the Hall attractor, are in contradiction with the spectral data of RX J1856.5-3754.

It has long been unclear if the small-scale magnetic structures on the neutron star (NS) surface could survive the fall-back episode. The study of the Hall cascade (Cumming, Arras and Zweibel 2004; Wareing and Hollerbach 2009) hinted that energy in small scales structures should dissipate on short timescales. Our new 2D magneto-thermal simulations suggest the opposite. For the first $\sim$10 kyrs after the fall-back episode with accreted mass $10^{-3} M_\odot$, the observed NS magnetic field appears dipolar, which is insensitive to the initial magnetic topology. In framework of the Ruderman & Sutherland (1975) vacuum gap model during this interval, non-thermal radiation is strongly suppressed. After this time the initial (i.e. multipolar) structure begins to re-emerge through the NS crust. We distinguish three evolutionary epochs for the re-emergence process: the growth of internal toroidal field, the advection of buried poloidal field, and slow Ohmic diffusion. The efficiency of the first two stages can be enhanced when small-scale magnetic structure is present. The efficient re-emergence of high order harmonics might significantly affect the curvature of the magnetospheric field lines in the emission zone. So, only after few $10^4$ yrs would the NS starts shining as a pulsar again, which is in correspondence with radio silence of central compact objects (CCOs). In addition, these results can explain the absence of good candidates for thermally emitting NSs with freshly re-emerged field among radio pulsars (Bogdanov, Ng and Kaspi 2014), as NSs have time to cool down, and supernova remnants can already dissipate.

Rotating proto-neutron stars can be important sources of gravitational waves to be searched for by present-day and future interferometric detectors. It was demonstrated by Imshennik that in extreme cases the rapid rotation of a collapsing stellar core may lead to fission and formation of a binary proto-neutron star which subsequently merges due to gravitational wave emission. In the present paper, we show that such dynamically unstable collapsing stellar cores may be the product of a former merger process of two stellar cores in a common envelope. We applied population synthesis calculations to assess the expected fraction of such rapidly rotating stellar cores which may lead to fission and formation of a pair of proto-neutron stars. We have used the BSE population synthesis code supplemented with a new treatment of stellar core rotation during the evolution via effective core-envelope coupling, characterized by the coupling time, $\tau_c$. The validity of this approach is checked by direct MESA calculations of the evolution of a rotating 15 $M_\odot$ star. From comparison of the calculated spin distribution of young neutron stars with the observed one, reported by Popov and Turolla, we infer the value $\tau_c \simeq 5 \times 10^5$ years. We show that merging of stellar cores in common envelopes can lead to collapses with dynamically unstable proto-neutron stars, with their formation rate being $\sim 0.1-1\%$ of the total core collapses, depending on the common envelope efficiency.

We discuss identification of possible counterparts and persistent sources related to Fast Radio Bursts (FRBs) in the framework of the model of supergiant pulses from young neutron stars with large spin-down luminosities. In particular, we demonstrate that at least some of sources of FRBs can be observed as ultraluminous X-ray sources (ULXs). At the moment no ULXs are known to be coincident with localization areas of FRBs. We searched for a correlation of FRB positions with galaxies in the 2MASS Redshift survey catalogue. Our analysis produced statistically insignificant overabundance ($p$-value $\approx 4\%$) of galaxies in error boxes of FRBs. In the very near future with even modestly increased statistics of FRBs and with the help of dedicated X-ray observations and all-sky X-ray surveys it will be possible to decisively prove or falsify the supergiant pulses model.

We discuss possible association of fast radio bursts (FRBs) with supergiant pulses emitted by young pulsars (ages $\sim$ tens to hundreds of years) born with regular magnetic field but very short -- few milliseconds -- spin periods. FRBs are extra-Galactic events coming from distances $d \lesssim 100$ Mpc. Most of the dispersion measure (DM) comes from the material in the freshly ejected SNR shell; for a given burst the DM should decrease with time. FRBs are not expected to be seen below $\sim 300 $ MHz due to free-free absorption in the expanding ejecta. A supernova might have been detected years before the burst; FRBs are mostly associated with star forming galaxies. The model requires that some pulsars are born with very fast spins, of the order of few milliseconds. The observed distribution of spin-down powers $\dot{E}$ in young energetic pulsars is consistent with equal birth rate per decade of $\dot{E}$. Accepting this injection spectrum and scaling the intrinsic brightness of FRBs with $\dot{E}$, we predict the following properties of a large sample of FRBs: (i) the brightest observed events come from a broad distribution in distances; (ii) for repeating bursts brightness either remains nearly constant (if the spin-down time is longer than the age of the pulsar) or decreases with time otherwise; in the latter case DM $\propto \dot{E}$.

The isolated neutron star (INS) 2XMM J104608.7-594306 is one of the only two to be discovered through their thermal emission since the ROSAT era. In a first dedicated XMM-Newton observation of the source, we found intriguing evidence of a very fast spin period. We re-observed 2XMM J104608.7-594306 with XMM-Newton to better characterise the spectral energy distribution of the source, confirm the candidate spin period, and possibly constrain the pulsar spin-down. Statistically acceptable spectral fits and meaningful physical parameters for the source are only obtained when the purely thermal spectrum is modified by at least one line in absorption. The implied distance is consistent with a location in (or in front of) the Carina nebula, and radiation radii are compatible with emission originating on most of the surface. Non-thermal X-ray emission is ruled out at levels above 0.5% of the source luminosity. Unfortunately, the second XMM-Newton observation proved inconclusive in terms of confirming (discarding) the fast candidate spin, providing an upper limit on the pulsed fraction of the source that is very close to the limiting sensitivity for detecting the modulation found previously. In the absence of an unambiguous period determination and an estimate of the magnetic field, the nature of the source remains open to interpretation. Its likely association with the Carina cluster and its overall spectral properties (only shared by a handful of other peculiar INSs) disfavour a standard evolutionary path, or one in which the source was previously recycled by accretion in a binary system. The INS 2XMM J104608.7-594306 may be similar to Calvera (1RXS J141256.0+792204), a neutron star for which the scenario of an evolved `anti-magnetar' has been discussed. A better age estimate and deeper radio and gamma-ray limits are required to further constrain the evolutionary state of the neutron star.

I review several scenarios of magnetar formation in binary systems via spin-up of a progenitor due to interaction with its companion. Mostly, these evolutionary channels lead to formation of isolated magnetars, and indeed, all well-established sources of this class are single objects. However, some binaries can survive, and several candidates to accreting magnetars have been proposed. I discuss this issue, and conclude that new accretion models can explain properties of the proposed candidates without large magnetic field in correspondence with models of magnetic field evolution.

We analyse the origin of the magnetic field decay in normal radio pulsars found by us in a recent study. This decay has a typical time scale $\sim 4 \times 10^5$~yrs, and operates in the range $\sim 10^5$~--~few$\times 10^5$~yrs. We demonstrate that this field evolution may be either due to the Ohmic decay related to the scattering from phonons, or due to the Hall cascade which reaches the Hall attractor. According to our analysis the first possibility seems to be more reliable. So, we attribute the discovered field decay mainly to the Ohmic decay on phonons which is saturated at the age few$\times 10^5$~yrs, when a NS cools down to the critical temperature below which the phonon scattering does not contribute much to the resistivity of the crust. Some role of the Hall effect and attractor is not excluded, and will be analysed in our further studies.

After a brief review of population synthesis of close-by cooling neutron stars, I focus on the interpretation of dichotomy of spin periods of near-by coolers. The existence of two well separated groups -- short period ($\sim$0.1-0.3 s) radio pulsars and long period ($\sim$3-10 s) radio quiet sources, aka the Magnificent seven, -- can not be easily explained in unified models developed recently (Popov et al. 2010, Gull{\'o}n et al. 2014). I speculate that the most natural solution of the problem can be in bimodal initial magnetic field distribution related to the existence of an additional mechanism of field generation in magnetars.

We propose that observations of "hidden" magnetars in central compact objects can be used to probe crustal activity of neutron stars with large internal magnetic fields. Estimates based on calculations by Perna \& Pons (2011), Pons \& Rea (2012) and Kaminker et al. (2014) suggest that central compact objects, which are proposed to be "hidden" magnetars, must demonstrate flux variations on the time scale of months-years. However, the most prominent candidate for the "hidden" magnetars --- CXO J1852.6+0040 in Kes 79, -- shows constant (within error bars) flux. This can be interpreted by lower variable crustal activity than in typical magnetars. Alternatively, CXO J1852.6+0040 can be in a high state of variable activity during the whole period of observations. Then we consider the source 1E161348-5055 in RCW103 as another candidate. Employing a simple 2D-modeling we argue that properties of the source can be explained by the crustal activity of the magnetar type. Thus, this object may be supplemented for the three known candidates for the "hidden" magnetars among central compact objects discussed in literature.

We apply the model of subsonic settling accretion on to isolated neutron stars accreting from the interstellar medium (AINS). We show that in this regime the expected mean X-ray luminosity from AINS turns out to be 2-3 orders of magnitude as small as the maximum possible Bondi value, i.e. $10^{27}-10^{28}$~erg s$^{-1}$. The intrinsically unstable character of settling accretion due to long plasma cooling time leads to regular appearance of X-ray flares with a duration of about one hour and a maximum luminosity of about the Bondi value, $\sim10^{31}$~erg~s$^{-1}$. This feature can be used to distinguished AINS from other dim X-ray sources. With the sensitivity of the forthcoming all-sky X-ray surveys the expected number of the potentially detectable AINS can be from a few to ten.

Supernova fallback disks around neutron stars have been discussed to influence the evolution of the diverse neutron star populations. Slowly rotating neutron stars are most promising to find such disks. Searching for the cold and warm debris of old fallback disks, we carried out Herschel PACS (70 $\mu$m, 160 $\mu$m) and Spitzer IRAC (3.6 $\mu$m, 4.5 $\mu$m) observations of eight slowly rotating ($P\approx 3 - 11$ s) nearby ($<1$ kpc) isolated neutron stars. Herschel detected 160 $\mu$m emission ($>5\sigma$) at locations consistent with the positions of the neutron stars RX J0806.4-4123 and RX J2143.0+0654. No other significant infrared emission was detected from the eight neutron stars. We estimate probabilities of 63%, 33% and 3% that, respectively, none, one, or both Herschel PACS 160 $\mu$m detections are unrelated excess sources due to background source confusion or an interstellar cirrus. If the 160 $\mu$m emission is indeed related to cold (10 K to 22 K) dust around the neutron stars, this dust is absorbing and re-emitting $\sim 10$% to $\sim 20$% of the neutron stars' X-rays. Such high efficiencies would be at least three orders of magnitude larger than the efficiencies of debris disks around nondegenerate stars. While thin dusty disks around the neutron stars can be excluded as counterparts of the 160 $\mu$m emission, dusty asteroid belts constitute a viable option.

We use a modified pulsar current analysis to study magnetic field decay in radio pulsars. In our approach we analyse the flow, not along the spin period axis as has been performed in previous studies, but study the flow along the direction of growing characteristic age, $\tau=P/(2\dot P)$. We perform extensive tests of the method and find that in most of the cases it is able to uncover non-negligible magnetic field decay (more than a few tens of per cent during the studied range of ages) in normal radio pulsars for realistic initial properties of neutron stars. However, precise determination of the magnetic field decay timescale is not possible at present. The estimated timescale may differ by a factor of few for different sets of initial distributions of neutron star parameters. In addition, some combinations of initial distributions and/or selection effects can also mimic enhanced field decay. We apply our method to the observed sample of radio pulsars at distances $<10$ kpc in the range of characteristic ages $8 \times 10^4 < \tau < 10^6$ years where, according to our study, selection effects are minimized. By analysing pulsars in the Parkes Multibeam and Swinburne surveys we find that, in this range, the field decays roughly by a factor of two. With an exponential fit this corresponds to the decay time scale $\sim 4 \times 10^5$ yrs. With larger statistics and better knowledge of the initial distribution of spin periods and magnetic field strength, this method can be a powerful tool to probe magnetic field decay in neutron stars.

In this methodological note we discuss some details and peculiarities of the cosmic expansion as viewed by a realistic observer. We show that the velocity $v_\Theta$ related to a change (measured by observer's clock) of the angular distance, plays an important role in formation of a meaningful observed picture of the expansion of the universe. Usage of this velocity and the angular distance (in addition to the standard approach --- proper distance and corresponding velocity) allows to present the cosmic expansion in a more illustrative manner. These parameters play a key role in visualization of the expansion.

The variety of the observational appearance of young isolated neutron stars must find an explanation in the framework of some unifying approach. Nowadays it is believed that such scenario must include magnetic field decay, the possibility of magnetic field emergence on a time scale $\lesssim 10^4$--$10^5$ yrs, significant contribution of non-dipolar fields, and appropriate initial parameter distributions. We present our results on the initial spin period distribution, and suggest that inconsistences between distributions derived by different methods for samples with different average ages can uncover field decay or/and emerging field. We describe a new method to probe the magnetic field decay in normal pulsars. The method is a modified pulsar current approach, where we study pulsar flow along the line of increasing characteristic age for constant field. Our calculations, performed with this method, can be fitted with an exponential decay for ages in the range $8\times 10^4$--$3.5 \times 10^5$ yrs with a time scale $\sim 5 \times 10^5$ yrs. We discuss several issues related to the unifying scenario. At first, we note that the dichotomy, among local thermally emitting neutron stars, between normal pulsars and the Magnificent Seven remains unexplained. Then we discuss the role of high-mass X-ray binaries in the unification of neutron star evolution. We note, that such systems allow to check evolutionary effects on a time scale longer than what can be probed with normal pulsars alone. We conclude with a brief discussion of importance of discovering old neutron stars accreting from the interstellar medium.

Properties of the population of millisecond extragalactic radio bursts discovered by Thornton et al. (2013) are in good correspondence with the hypothesis that such events are related to hyperflares of magnetars, as was proposed by us after the first observation of an extragalactic millisecond radio burst by Lorimer et al. (2007). We also point that some of multiple millisecond radio bursts from M31 discovered by Rubio-Herrera et al. (2013) also can be related to weaker magnetar bursts.

We discuss the origin of strong crustal magnetic field in one of central compact objects (CCOs) --- a neutron star PSR J1852+0040 in the supernova remnant Kes 79. Taking into account its relatively long present day spin period we conclude that the field could not be generated via a dynamo mechanism. If this neutron star indeed is a magnetar with field submerged during a strong fall-back episode, then it argues against the dynamo field origin in magnetars. Otherwise, Kes 79 is not a close relative of normal magnetars. A discovery of an anti-magnetar with a millisecond period and strong crustal field identifiable, for example, due to large pulse fraction, would be the proof of the dynamo field origin. Existence of such sources is in correspondence with the present standard picture of neutron star unification. However, the fraction of magnetars with submerged fields can be small --- few percent of the total number of CCOs.

Recently, Lee et al. used Gaussian mixture models (GMM) to study the radio pulsar population. In the distribution of normal pulsars in the P-dotP plane, they found four clusters. We develop this approach further and apply it to different synthetic pulsar populations in order to determine whether the method can effectively select groups of sources that are physically different. We check several combinations of initial conditions as well as the models of pulsar evolution and the selection effects. We find that in the case of rapidly evolving objects, the GMM is oversensitive to parameter variations and does not produce stable results. We conclude that the method does not help much to identify the sub-populations with different initial parameters or/and evolutionary paths. For the same reason, the GMM does not discriminate effectively between theoretical population synthesis models of normal radio pulsars.

We analyze different possibilities to explain the wide initial spin period distribution of radio pulsars presented by Noutsos et al. (2013). With a population synthesis modeling we demonstrate that magnetic field decay can be used to interpret the difference between the recent results by Noutsos et al. (2013) and those by Popov & Turolla (2012), where a much younger population of NSs associated with supernova remnants with known ages has been studied. In particular, an exponential field decay with tau_mag = 5Myrs can produce a "tail" in the reconstructed initial spin period distribution (as obtained by Noutsos et al. 2013) up to P_0 > 1 s starting with a standard gaussian with P_0 = 0.3 s and sigma_P = 0.15 s. Another option to explain the difference between initial spin period distributions from Noutsos et al. (2013) and Popov & Turolla (2012) -- the emerging magnetic field -- is also briefly discussed.

The large number of isolated black holes (IBHs) in the Galaxy, estimated to be 10^8, implies a very high density of 10^-4 pc^-3 and an average distance between IBHs of 10 pc. Our study shows that the magnetic flux, accumulated on the horizon of an IBH because of accretion of interstellar matter, allows the Blandford-Znajeck mechanism to be activated. Thus, electron-positron jets can be launched. We have performed 2D numerical modelling which allowed the jet power to be estimated. Their inferred properties make such jets a feasible electron accelerator which, in molecular clouds, allows electron energy to be boosted up to 1 PeV. For the conditions expected in molecular clouds the radiative cooling time should be comparable to the escape time. Thus these sources can contribute both to the population of unidentified point-like sources and to the local cosmic ray (CR) electron spectrum. The impact of the generated electron CRs depends on the diffusion rate inside molecular clouds (MCs). If the diffusion regime in a MC is similar to Galactic diffusion, the produced electrons should rapidly escape the cloud and contribute to the Galactic CR population at very high energies >100 TeV. However, due to the modest jet luminosity (at the level of 10^35 erg s^-1) and low filling factor of MC, these sources cannot make a significant contribution to the spectrum of cosmic ray electrons at lower energies. On the other hand, if the diffusion within MCs operates at a rate close to the Bohm limit, the CR electrons escaping from the source should be confined in the cloud, significantly contributing to the local density of CRs. The IC emission of these locally-generated CRs may explain the variety of gamma ray spectra detected from nearby MCs.

We present a multi-wavelength (X-ray, $\gamma$-ray, optical and radio) study of HD 194816, a late O-type X-ray detected spectroscopic binary. X-ray spectra are analyzed and the X-ray photon arrival times are checked for pulsation. In addition, newly obtained optical spectroscopic monitoring data on HD 164816 are presented. They are complemented by available radio data from several large scale surveys as well as the \emph{FERMI} $\gamma$-ray data from its \emph{Large Area Telescope}. We report the detection of a low energy excess in the X-ray spectrum that can be described by a simple absorbed blackbody model with a temperature of $\sim$ 50 eV as well as a 9.78 s pulsation of the X-ray source. The soft X-ray excess, the X-ray pulsation, and the kinematical age would all be consistent with a compact object like a neutron star as companion to HD 164816. The size of the soft X-ray excess emitting area is consistent with a circular region with a radius of about 7 km, typical for neutron stars, while the emission measure of the remaining harder emission is typical for late O-type single or binary stars. If HD 164816 includes a neutron star born in a supernova, this supernova should have been very recent and should have given the system a kick, which is consistent with the observation that the star HD 164816 has a significantly different radial velocity than the cluster mean. In addition we confirm the binarity of HD 164816 itself by obtaining an orbital period of 3.82 d, projected masses $m_1 {\rm sin}^{3} i$ = 2.355(69) M$_\odot$, $m_2 {\rm sin}^{3} i$ = 2.103(62) M$_\odot$ apparently seen at low inclination angle, determined from high-resolution optical spectra.

SXP 1062, a newly discovered Be/X-ray binary in the Small Magellanic Cloud, provides the first example of a robust association with a supernova remnant (SNR). The short age estimated for the SNR qualifies SXP 1062 as the youngest known source in its class, ~ 10^4 yr. Here we discuss possible evolutionary scenarios for SXP 1062 in the attempt to reconcile its long spin period, P=1062 s, and short age. Our results indicate that SXP 1062 may host a neutron star born with a large initial magnetic field, typically in excess of ~ 10^14 G, which then decayed to ~ 10^13 G.

A subpopulation of neutron stars (NSs), known as central compact objects (CCOs) in supernova remnants, are suspected to be low-field objects basing on $P$-$\dot P$ measurements for three of them. The birth rate of low-field NSs is probably comparable with the birth rate of normal radio pulsars. However, among compact objects in High-Mass X-ray Binaries (HMXBs) we do not see robust candidates for low-field NSs. We propose that this contradiction can be solved if magnetic fields of CCOs was buried due to strong fall-back, and then the field emerges on the time scale $10^4$-$10^5$ yrs.

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) -- which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact objects in supernova remnants (CCOs) -- represent a key element in understanding the neutron star phenomenology. We report the results of an observational campaign to study the properties of the source 2XMM J104608.7-594306. Its evolutionary state is investigated by means of deep dedicated observations obtained with XMM-Newton, the ESO Very Large Telescope, as well as publicly available gamma-ray data from the Fermi and AGILE missions. The observations confirm previous expectations and reveal a unique type of object. The source, which is likely within the Carina Nebula, has a soft spectrum with absorption features and no magnetospheric emission. The optical counterpart is fainter than V=27 and no gamma-ray emission is significantly detected. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, P=18.6 ms. We interpret these new results in the light of the observed properties of the currently known neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of these scenarios can satisfactorily explain the collective properties of 2XMM J104608.7-594306, although it may be related to the still poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AO11), will help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin down.

We present estimates of initial spin periods, $P_0$, for radio pulsars associated with supernova remnants. By using the published data on 30 objects, we were able to derive a reliable estimate for the initial spin period, assuming standard magneto-dipole spin-down (braking index n=3), in many cases. Our set of estimates is still not sufficient to infer the exact shape of the initial period distribution. However, we show that a gaussian distribution with mean and deviation $\sim 0.1$ s is consistent with our results, while flat, wide distributions and very narrow ones are disfavored.

Some modern models of neutron star evolution predict that initially large magnetic fields rapidly decay down to some saturation value $\sim {\rm few}\times 0^{13}$ G and weaker magnetic fields do not decay significantly (Pons). It is difficult to check the predictions of this model for initially highly magnetized objects on the time scale of a few million years. We propose to use Be/X-ray binaries for this purpose. We apply several methods to estimate magnetic fields of neutron stars in these accreting systems using the data obtained by the RXTE satellite (Galache). Only using the most modern approach for estimating the magnetic field strengths of long period NSs as proposed by Shakura we are able to obtain a field distribution compatible with predictions of the theoretical model of field decay of Pons.

The newly discovered Be/X-ray binary in the Small Magellanic Cloud, SXP 1062, provides the first example of a robust association with a supernova remnant (SNR). The short age estimated for the SNR qualifies SXP 1062 as the youngest known source in its class, tau ~ 1e4 yr. As such, it allows to test current models of magneto-rotational evolution of neutron stars ina still unexplored regime. Here we discuss possible evolutionary scenarios for SXP 1062 in the attempt to reconcile its long spin period, P=1062 s, and short age. Although several options can be considered, like an anomalously long initial period or the presence of a fossil disc, our results indicate that SXP 1062 may host a neutron star born with a large initial magnetic field, typically in excess of ~ 1e14 G, which then decayed to ~ 1e13 G.

In these lectures presented at Baikal summer school on physics of elementary particles and astrophysics 2011, I present a wide view of neutron star astrophysics with special attention paid to young isolated compact objects and studies of the properties of neutron star interiors using astronomical methods.

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.

We have carried out a search for radio emission at 820 MHz from six X-ray dim isolated neutron stars with the Robert C. Byrd Green Bank Radio Telescope. No transient or pulsed emission was found using fast folding, fast Fourier transform, and single-pulse searches. The corresponding flux limits are about 0.01 mJy for pulsed emission, depending on the integration time for the particular source and assuming a duty cycle of 2%, and 20 mJy for single dispersed pulses. These are the most sensitive limits to date on radio emission from X-ray dim isolated neutron stars. There is no evidence for isolated radio pulses, as seen in a class of neutron stars known as rotating radio transients. Our results imply that either the radio luminosities of these objects are lower than those of any known radio pulsars, or they could simply be long-period nearby radio pulsars with high magnetic fields beaming away from the Earth. To test the latter possibility, we would need around 40 similar sources to provide a 1 sigma probability of at least one of them beaming toward us. We also give a detailed description of our implementation of the Fast Folding Algorithm.

We consider the possible existence of a common channel of evolution of binary systems, which results in a gamma-ray burst during the formation of a black hole or the birth of a magnetar during the formation of a neutron star. We assume that the rapid rotation of the core of a collapsing star can be explained by tidal synchronization in a very close binary. The calculated rate of formation of rapidly rotating neutron stars is qualitatively consistent with estimates of the formation rate of magnetars. However, our analysis of the binarity of newly-born compact objects with short rotational periods indicates that the fraction of binaries among them substantially exceeds the observational estimates. To bring this fraction into agreement with the statistics for magnetars, the additional velocity acquired by a magnetar during its formation must be primarily perpendicular to the orbital plane before the supernova explosion, and be large.

Over the last decade, X-ray observations unveiled the existence of several classes of isolated neutron stars (INSs) which are radio-quiet or exhibit radio emission with properties much at variance with those of ordinary radio pulsars. The identification of new sources is crucial in order to understand the relations among the different classes and to compare observational constraints with theoretical expectations. A recent analysis of the 2XMMp catalogue provided less than 30 new thermally emitting INS candidates. Among these, the source 2XMM J104608.7-594306 appears particularly interesting because of the softness of its X-ray spectrum and of the present upper limits in the optical, which imply a logarithmic X-ray-to-optical flux ratio greater than 3.1, corrected for absorption. We present the X-ray and optical properties of 2XMM J104608.7-594306 and discuss its nature in the light of two possible scenarios invoked to explain the X-ray thermal emission from INSs: the release of residual heat in a cooling neutron star, as in the seven radio-quiet ROSAT-discovered INSs, and accretion from the interstellar medium. We find that the present observational picture of 2XMM J104608.7-594306 is consistent with a distant cooling INS with properties in agreement with the most up-to-date expectations of population synthesis models: it is fainter, hotter and more absorbed than the seven ROSAT sources and possibly located in the Carina Nebula, a region likely to harbour unidentified cooling neutron stars. The accretion scenario, although not entirely ruled out by observations, would require a very slow (~10 km/s) INS accreting at the Bondi-Hoyle rate.

The role of binary progenitors of neutron stars in the apparent distribution of space velocities and spin-velocity alignment observed in young pulsars is studied. A Monte-Carlo synthesis of pulsar population from single and binary stars with different assumptions about the NS natal kick model (direction distribution, amplitude, and kick reduction in binary progenitors which experienced mass exchange due to Roche lobe overflow with initial masses on the main sequence from the range 8-11 $M_\odot$) is performed. The calculated spin-velocity alignment distributions are compared with observational data obtained from radio polarization measurements. The observed space velocity of pulsars is found to be mostly shaped by the natal kick velocity form and its amplitude; the fraction of binaries is not important here for reasonably large kicks. The distribution of kick direction relative to the spin axis during the formation of a NS is found to affect strongly the spin-velocity correlation of pulsars. Comparison with observed pulsar spin-velocity angles favours a sizeable fraction of binary progenitors and the kick-spin angle $\sim 5-20^\circ$. The form of the initial binary mass ratio distribution does not affect our results.

I discuss several scenarios to explain properties of the radio transient source GCRT J1745-3009. Namely, a highly magnetized neutron star on the propeller or georotator stage, a transient propeller, and an ejector in a binary system are discussed. Simple populational estimates favor the transient propeller model.

We propose a new method to detect observational appearance of Dark Matter axions. The method utilizes observations of neutron stars (NSs) in radio. It is based on the conversion of axions to photons in strong magnetic fields of NSs (Primakoff effect). Whether the conversion takes place, the radio spectrum of the object would have a very distinctive feature -- a narrow spike at a frequency corresponding to the rest mass of the axion. For example, if the coupling constant of the photon-axion interaction is $M=10^{10}$ GeV, the density of Dark Matter axions is $\rho=10^{-24} {\rm g cm^{-3}}$, and the axion mass is $5 {\rm \mu eV}$, then a flux from a strongly magnetized ($10^{14}$ G) NS at the distance 300 pc from the Sun is expected to be about few tenths of mJy at the frequency $\approx 1200$ MHz in the bandwidth $\approx 3$ MHz. Close-by X-ray dim isolated neutron stars are proposed as good candidates to look for such radio emission.

Here I discuss possible relations between free precession of neutron stars, Tkachenko waves inside them and glitches. I note that the proposed precession period of the isolated neutron star RX J0720.4-3125 (Haberl et al. 2006) is consistent with the period of Tkachenko waves for the spin period 8.4s. Based on a possible observation of a glitch in RX J0720.4-3125 (van Kerkwijk et al. 2007), I propose a simple model, in which long period precession is powered by Tkachenko waves generated by a glitch. The period of free precession, determined by a NS oblateness, should be equal to the standing Tkachenko wave period for effective energy transfer from the standing wave to the precession motion. A similar scenario can be applicable also in the case of the PSR B1828-11.

In this short note I discuss the hypothesis that bursting activity of magnetars evolves in time analogously to the glitching activity of normal radio pulsars (i.e. sources are more active at smaller ages), and that the increase of the burst rate follows one of the laws established for glitching radio pulsars. If the activity of soft gamma repeaters decreases in time in the way similar to the evolution of core-quake glitches ($\propto t^{5/2}$), then it is more probable to find the youngest soft gamma repeaters, but the energy of giant flares from these sources should be smaller than observed $10^{44}$ --$10^{46}$ ergs as the total energy stored in a magnetar's magnetic field is not enough to support thousands of bursts similar to the prototype 5 March 1979 flare.

We propose that the strong millisecond extragalactic radio burst (mERB) discovered by Lorimer et al. (2007) may be related to a hyperflare from an extragalactic soft gamma-ray repeater. The expected rate of such hyperflares, $\sim$ 20 - 100 d$^{-1}$ Gpc$^{-3}$, is in good correspondence with the value estimated by Lorimer et al. The possible mechanism of radio emission can be related to the tearing mode instability in the magnetar magnetosphere as discussed by Lyutikov (2002), and can produce the radio flux corresponding to the observed $\sim$ 30 Jy from the mERB using a simple scaling of the burst energy.

We have carried out a search for radio emission from six X-ray dim isolated neutron stars (XDINSs) observed with the Robert C. Byrd Green Bank Radio Telescope (GBT) at 820 MHz. No bursty or pulsed radio emission was found down to a 4sigma significance level. The corresponding flux limit is 0.01-0.04 mJy depending on the integration time for the particular source and pulse duty cycle of 2%. These are the most sensitive limits yet on radio emission from these objects.

We present our new advanced model for population synthesis of close-by cooling NSs. Detailed treatment of the initial spatial distribution of NS progenitors and a detailed ISM structure up to 3 kpc give us an opportunity to discuss the strategy to look for new isolated cooling NSs. Our main results in this respect are the following: new candidates are expected to be identified behind the Gould Belt, in directions to rich OB associations, in particular in the Cygnus-Cepheus region; new candidates, on average, are expected to be hotter than the known population of cooling NS. Besides the usual approach (looking for soft X-ray sources), the search in 'empty' $\gamma$-ray error boxes or among run-away OB stars may yield new X-ray thermally emitting NS candidates.

We compute the probability that intermediate mass black holes (IMBHs) capture companions due to dynamical interactions and become accreting sources, and explore the possibility that the accreting IMBHs would appear as ultra-luminous X-ray sources (ULXs). We focus on IMBHs originating from low-metallicity Population III stars. Two channels of IMBH formation are considered: from primordial halos in the framework of hierarchical clustering, and from non-mixed, zero-metallicity primeval gas in galactic discs. IMBHs can form binary systems due to tidal captures of single stars and exchange interactions with existing binary systems in galactic discs. We find that neither formation mechanism of the accreting IMBH binary is able to provide enough sources to explain the observed population of ULXs. Even at sub-ULX luminosity, the total number of accreting IMBHs with $L> 10^{36}$ erg s$^{-1}$ with dynamically captured companions is found to be $< 0.01$ per galaxy.

In these lecture notes I briefly discuss the present day situation and new discoveries in astrophysics of neutron stars focusing on isolated objects. The latter include soft gamma repeaters, anomalous X-ray pulsars, central compact objects in supernova remnants, the Magnificent seven, and rotating radio transients. In the last part of the paper I describe available tests of cooling curves of neutron stars and discuss different additional constraints which can help to confront theoretical calculations of cooling with observational data.

We revisit the case of magneto-rotational evolution of neutron stars (NSs) with accounting for changes of the angle \chi between spin and magnetic axes. This element of the evolution of NSs is very important for age estimates and population modeling, but usually it is neglected. In the framework of two models of energy losses we demonstrate that unless NSs are born with the inclination angle \chi_0 very close (within 1 degree) to the position of maximal losses, pulsars with short initial periods quickly reach regions of small energy losses (aligned rotators in the case of magneto-dipole losses and orthogonal rotators -- in the case of current losses) without significant spin-down. This means that either most of known NSs should be born with relatively long periods close to presently observed values, or the initial inclination angles should be large (close to 90) for the magnetodipole model and very small (close to 0) for the current losses model, or that both models cannot be applied for the whole evolutionary track of a typical NS. In particular, magnetar candidates should be born with periods close to the observed ones within the scope of both models of energy losses. We discuss how these considerations can influence population synthesis models and age estimates for different types of isolated NSs. Focusing on the model of current losses we illustrate our conclusion with evolutionary track reconstructions on the P--\sin \chi plane. We conclude, that most probably both existing models of energy losses -- magneto-dipole and longitudinal current -- require serious modifications, on the other hand, estimates obtained under the standard assumption of \sin \chi=1=const does not have solid theoretical ground and should be taken with care.

We discuss a set of tests which confront observations of cooling compact objects and theories of their thermal evolution. As an example we apply the recently developed $\mathrm{Log N}$-$\mathrm{Log S}$ test of compact star cooling theories to hybrid stars with a color superconducting quark matter core, we also apply and discuss other existing tests. While there is not yet a microscopically founded superconducting quark matter phase which would fulfill constraints from cooling phenomenology, we explore the hypothetical 2SC+X phase and show that the magnitude and density-dependence of the X-gap can be chosen to satisfy a set of tests: the temperature -- age ($\mathrm{T}$-$\mathrm{t}$) test, the $\mathrm{Log N}$-$\mathrm{Log S}$ test, the brightness constraint, and the mass spectrum constraint. Some recent modifications of the population synthesis model used to obtain the $\mathrm{Log N}$-$\mathrm{Log S}$ distribution are briefly discussed. In addition, we propose to use the age-distance diagram as a new tool to study the local population of young isolated neutron stars.

The Magnificent Seven have all been discovered by their exceptional soft X-ray spectra and high ratios of X-ray to optical flux. They all are considered to be nearby sources. Searching for similar neutron stars with larger distances, one expects larger interstellar absorption resulting in harder X-ray counterparts. Current interstellar absorption treatment depends on chosen abundances and scattering cross-sections of the elements as well as on the 3D distribution of the interstellar medium. After a discussion of these factors we use the comprehensive 3D measurements of the Local Bubble by Lallement et al. 2003 to construct two simple models of the 3D distribution of the hydrogen column density. We test these models by using a set of soft X-ray sources with known distances. Finally, we discuss possible applications for distance estimations and population synthesis studies.

We discuss three topics related to the neutron star (NS) mass spectrum. At first we discuss the possibility to form low-mass ($ M \stackrel{<}{\sim} 1 M_{\odot}$) and suggest this is possible only due to fragmentation of rapidly rotating proto-NSs. Such low-mass NSs should have very high spatial velocities which could allow identification. A critical assessment of this scenario is given. Secondly, we discuss mass growth due to accretion for NSs in close binary systems. With the help of numerical population synthesis calculations we derive the mass spectrum of massive ($M > 1.8 M_{\odot}$) NSs. Finally, we discuss the role of the mass spectrum in population studies of young cooling NSs. We formulate a kind of {\it mass constraint} which can be helpful, in our opinion, in discussing different competive models of the thermal evolution of NSs.

We highlight similarities between recently discovered Rotating Radio Transients and X-ray Dim Isolated Neutron Stars. In particular, it is shown that X-ray Dim Isolated Neutron Stars have a birthrate comparable to that of Rotating Radio Transients. On the contrary, magnetars have too low a formation rate to account for the bulk of the radio transient population. The consequences of the recent detection of a thermal X-ray source associated with one of the Rotating Radio Transients on the proposed scenarios for these sources are also discussed.

We apply the recently developed Log N - Log S test of compact star cooling theories for the first time to hybrid stars with a color superconducting quark matter core. While there is not yet a microscopically founded superconducting quark matter phase which would fulfill constraints from cooling phenomenology, we explore the hypothetical 2SC+X phase and show that the magnitude and density-dependence of the X-gap can be chosen to satisfy a set of tests: temperature - age (T - t), the brightness constraint, Log N - Log S, and the mass spectrum constraint. The latter test appears as a new conjecture from the present investigation.

We present a compilative catalogue of Be/X-ray binaries and candidates in the Galaxy and in the Large and Small Magellanic Clouds. This catalogue contains 130 sources and provides information on names and spectral types of optical components, distances, spin characteristics of neutron stars, and on orbital and X-ray properties of binary systems. We give brief comments on each object and provide necessary references to original data.

We consider nearby young black holes formed after supernova explosions in close binaries whose secondary components are currently observed as the so-called runaway stars. Using data on runaway stars and making reasonable assumptions about the mechanisms of supernova explosions and binary breakup, we estimate the present position of nearby young black holes. For two objects, we obtained relatively small error regions ($\sim 50$-100 deg$^2$). The possibility of detecting these nearby young black holes is discussed.

The study of thermal emission from isolated, cooling neutron stars plays a key role in probing the physical conditions of both the star crust and the core. The comparison of theoretical models for the star thermal evolution with the surface temperature derived from X-ray observations of sources of different age is one of the main tools to investigate the properties of the interior and constrain the equation of state. Here we propose to use population synthesis studies as an independent approach to test the physics governing the star cooling. Theoretical Log N - Log S distributions depend on the assumed neutron star thermal evolution. We have computed distributions for several different cooling scenarios and found that comparison with the observed Log N - Log S of isolated neutron stars is effective in discriminating among cooling models. Among the eleven cooling models considered in this paper, all of which may reproduce the observed temperature vs. age diagram, only at most three can explain the Log N - Log S distribution of close-by cooling neutron stars. The Log N - Log S test, being a ``global'' one and despite some limitations, appears capable of complementing the standard temperature vs. age test used up to now.

Among a dozen known magnetar candidates there are no binary objects. As an estimate of a fraction of binary neutron stars is about 10% it is reasonable to address the question of solitarity of magnetars, to estimate theoretically the fraction of binary objects among them, and to mark o probable companions. We present population synthesis calculations of binary systems. Our goal is to estimate the number of neutron stars originated from progenitors with enhanced rotation, as such compact objects can be expected to have large magnetic fields, ie. they can be magnetars. The fraction of such neutron stars in our calculations is about 13-16%. Most of these objects are isolated due to coalescences of components prior to a neutron star formation, or due to a system disruption after a supernova explosion. The fraction of such neutron stars in survived binaries is about 1% or lower. Their most numerous companions are black holes.

We propose that the best sites to search for SGRs outside the Local group are galaxies with active massive star formation. Different possibilities to observe SGR activity from these sites are discussed. In particular we searched for giant flares from nearby galaxies ($\sim 2$ -- 4 Mpc) M82, M83, NGC 253, and NGC 4945 in the BATSE data. No candidates alike giant SGR flares were found. The absence of such detections implies that the rate of giant flares with energy release in the initial spike above $0.5 \cdot 10^{44}$ erg is less then 1/25 yr$^{-1}$ in our Galaxy. However, hyperflares similar to the one of 27 December 2004 can be observed from larger distances. Nevertheless, we do not see any significant excess of short GRBs from the Virgo galaxy cluster and from galaxies Arp 299 and NGC 3256 with extremely high star formation rate. This implies that the galactic rate of hyperflares with energy release $\sim 10^{46}$ erg is less than $\sim 10^{-3}$ yr$^{-1}$. With this constraint the fraction of possible extragalactic SGR hyperflares among BATSE short GRBs should not exceed few percents. We present a list of short GRBs coincident with galaxies mentioned above, and discuss the possibility that some of them are SGR giant flares. We propose that the best target for observations of extragalactic SGR flares by {\it Swift} is the Virgo cluster.

We propose that the best sites to search for SGRs outside the Local group are galaxies with active star formation. Different possibilities to observe SGR activity from these sites are discussed. In particular we searched for giant flares from M82, M83, NGC 253, and NGC 4945 in the BATSE data. We present a list of potential candidates, however, in our opinion no good candidates alike giant SGR flares were found. Still, hyperflares similar to the one of 27 December 2004 can be observed from larger distances. From the BATSE data we select 5 candidates coincident with the galaxies Arp 299 and NGC 3256 which have very high rate of starformation, and propose that they can be examples of hyperflares; but this result has low statistical significance.

We follow the approach used by Davies and Pringle (1981) and discuss the trans-sonic substage of the propeller regime. This substage is intermediate between the supersonic and subsonic propeller substages. In the trans-sonic regime an envelope around a magnetosphere of a neutron star passes through a kind of a reorganization process. The envelope in this regime consists of two parts. In the bottom one turbulent motions are subsonic. Then at some distance $r_\mathrm{s}$ the turbulent velocity becomes equal to the sound velocity. During this substage the boundary $r_\mathrm{s}$ propagates outwards till it reaches the outer boundary, and so the subsonic regime starts. We found that the trans-sonic substage is unstable, so the transition between supersonic and subsonic substages proceeds on the dynamical time scale. For realistic parameters this time is in the range from weeks to years.

As it is well known for stiff equations of state an existence of neutron stars with masses $\ga 2 M_{\sun}$ is possible. Especially interesting possibility is opened if the equation of state based on the Skyrme theory is realized in nature. This equation of state was proposed recently by Ouyed and Butler. We discuss different channels of formation of massive rapidly rotating neutron stars. We use a population synthesis code to estimate numbers of massive neutron stars on different evolutionary stages. A neutron star increases its mass by accretion from a secondary companion. Significant growth of a neutron star mass due to accretion is possible only for certain values of initial parameters of the binary. In this paper we show that significant part of massive neutron stars with $M\ga 2 M_{\sun}$ can be observed as millisecond radio pulsars, as X-ray sources in pair with white dwarfs, and as accreting neutron stars with very low magnetic fields.

We briefly describe different astrophysical applications of a population synthesis method. In some details we discuss the population synthesis of close binary systems and of isolated neutron stars.

Using modern high-resolution observations of extragalactic compact radio sources we obtain an estimate of the upper bound on a photon electric charge at the level $e_{\gamma} \lesssim 3 \cdot 10^{-33}$ of elementary charge (assuming the photon charge to be energy independent). This is three orders of magnitude better than the limit obtained with radio pulsar timing. Also we set a limit on a photon charge in the gamma-ray band (energies about 0.1 MeV). In future the estimate made for extragalactic sources can be significantly improved.

I discuss young close-by compact objects, recent supernovae in the solar neighbourhood, and point to their connection with cosmic ray studies. Especially the role of the Gould Belt is underlined.

An age--distance diagram for close-by young isolated neutron stars of different types is introduced and discussed. It is shown that such a diagram can be a useful tool for studying of populations of sources if their detectability strongly depends on age and distance.

Here we present a compilation of data on Be/X-ray binaries. On the whole we include 90 objects into our catalogue. Brief comments on each object are provided.

We present new population synthesis calculations of close young neutron stars. In comparison with our previous investigation we use a different neutron star mass spectrum and different initial spatial and velocity distributions. The results confirm that most of ROSAT dim radioquiet isolated neutron stars had their origin in the Gould Belt. Several tens of young neutron stars can be identified in future in ROSAT data at low galactic latitudes and some of them also can be EGRET unidentified sources.

We propose that the bimodal nature of the kick velocity distribution of radio pulsars is connected with the dichotomy between hadronic stars ({\it i.e.} neutron stars with no quark matter content) and quark stars. Bimodality can appear due to different mechanisms of explosion which leads to the formation of two types of compact stars or due to two different sets of parameters mastering a particular kick mechanism. The low velocity maximum (at $\sim 100$ km s$^{-1}$) is connected with hadronic star formation, whereas the second peak corresponds to quark stars. In the model of delayed collapse of hadronic stars to quark stars (Berezhiani et al. 2003\nocite{bbd2003}) quark deconfinement leads to a second energy release, and to a second kick, in addition to the kick imparted to the newly formed hadronic star during the supernova explosion. If the electromagnetic rocket mechanism can give a significant contribution to pulsar kicks, then the high velocity peak can be connected with the shorter initial spin periods of quark stars with respect to hadronic stars. We discuss {\it pro et contra} of these scenarios.

We suggest that low-mass compact objects (hadron stars, quark stars) with $M<1 M_{\odot}$ can appear only due to fragmentation of rapidly rotating proto-neutron stars. Such low-mass stars receive large kicks due to an explosion of a lighter companion in a pair of fragments, or due to dynamical ejection of one of the lighter components in the case when three bodies are formed. As far as low-mass compact objects are expected to be slowly cooling in all popular models of thermal evolution possible candidates are expected to be found among hot high velocity sources.

We present a simple model for the evolutionary states of neutron stars in pre-low-mass X-ray binaries during the main sequence stage of donor stars. It is shown that for typical parameters some of the neutron stars in these systems can not accrete matter from winds of their companions in a stable way. Accretors are found if neutron stars have magnetic fields about $10^{13}$ G and higher and/or in close systems with orbital periods about a few days and shorter.

We briefly review selected results in astrophysics of neutron stars (NSs) obtained during the last two years, focusing on isolated radioquiet objects. We discuss in some details population synthesis of close-by isolated NSs, spectra of INSs, detection of spectral features in these sources (including cyclotron features), recent results on velocity distribution of NSs and accretion onto INSs.

We briefly discuss the evolutionary path and observational appearance of isolated neutron stars (INSs) focusing on radioquiet objects. There are many reasons to believe that these sources are extremely elusive once the star surface has cooled down: their high spatial velocities, the long propeller stage and/or the very low accretion efficiency. We describe recent population synthesis models of close-by young INSs, highlighting the major difficulties encountered in the past by these simulations in reproducing the observed properties of known sources. As we show, a likely possibility is that most of the INSs in the Solar proximity are young (less than few Myrs) neutron stars born in the Gould Belt. To stay hot enough and sustain X-ray emission for a time $\approx 1$ Myr, they probably need to be low- to medium-massive, with $M$ less than$\sim 1.35 M_\odot$.

We model populations of isolated neutron stars in the Galaxy following their orbital and magneto-rotational evolution. It is shown that accretors become more abundant than coolers at fluxes below $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1},$ and one can predict that about one accreting neutron star per square degree should be observed at the {\it Chandra} and {\it Newton} flux limits of $\sim 10^{-16}$ erg cm$^{-2}$ s$^{-1}.$ The soft ROSAT sources associated with isolated neutron stars can be relatively young cooling objects only if the neutron star birth rate in the Solar vicinity during the last $\sim 10^6$ yr is higher than that inferred from radiopulsar observations.

The origin of the local population of young, cooling neutron stars is investigated with a population synthesis model taking into account the contribution of neutron stars born in the Gould Belt, in addition to those originating in the Galactic disk. We estimate their emission in the soft X-ray band as a function of distance and age and construct the Log N -- Log S distribution. It is shown that the inclusion of neutron stars from the Gould Belt provides a good fit to the observed Log N -- Log S distribution. As the Sun is situated inside the Gould Belt, one can naturally explain the comparative local overabundance of massive progenitors and can remove the difficulty of the deficit of relatively bright ($\ga 0.1$ ROSAT PSPC cts s$^{-1}$) cooling neutron stars previously reported from models where only neutron stars from the Galactic disk were accounted for.

We present Log N -- Log S distribution for close-by young isolated neutron stars. On the basis of this distribution it is shown that the seven ROSAT isolated neutron stars (if they are young cooling objects) are genetically related to the Gould Belt. We predict, that there are about few tens unidentified close-by young isolated neutron stars in the ROSAT All-Sky Survey. The possibility that these seven peculiar sources contain a neutron star less massive and more magnetized then in ordinary radiopulsars is also discussed. In the aftermath of relatively close recent supernova explosions (1 kpc around the Sun, a few Myrs ago), a few black holes might have been formed, according to the local initial mass function. We thus discuss the possibility of determining approximate positions of close-by isolated black holes using data on runaway stars and simple calculations of binary evolution and disruption.

In this paper we review present state of theoretical and observational exploration of isolated neutron stars. We mainly discuss objects which do not show usual radiopulsar activity. We include into our review the following types of objects: dim X-ray sources in the Galactic disk and in globular clusters; compact X-ray sources in supernova remnants; anomalous X-ray pulsars; soft gamma repeaters; Geminga. Also we briefly discuss related types of objects: isolated black holes, quark stars and supernovae. We review basic physical processes connected with evolution of isolated neutron stars: cooling; magnetic field decay; period evolution; magnetospheric processes, accretion. We discuss populational synthesis of isolated neutron stars of different types.

We suggest that the seven radio-quiet isolated neutron stars observed with ROSAT are young cooling objects associated to recent near-by supernova explosions which formed runaway stars and the Local Bubble, affecting the topology of the interstellar medium in the vicinity of the Sun (within a few hundred parsecs). In the aftermath of these explosions, a few black holes might have been formed, according to the local initial mass function. We thus discuss the possibility of determining approximate positions of close-by isolated black holes using data on runaway stars and simple calculations of binary evolution and disruption.

We model Galactic populations of accreting and cooling isolated neutron stars in the attempt to explore their link with a new class of dim soft X-ray sources revealed by ROSAT. For accretors we follow the magneto-rotational and dynamical evolution in the Galactic potential and a realistic large scale distribution of the interstellar medium is used. Under standard assumptions old neutron stars enter the accretor stage only if their magnetic field exceeds $\approx 10^{11}$--$10^{12}$ G. We predict about 1 source per square degree for fluxes $\approx 10^{-15}$--$10^{-16}$ erg cm$^{-2}$s$^{-1}$ in the energy range 0.5-2 keV. Cooling neutron stars are explored within a simpler model of local sources, including however interstellar absorption. They are found to be significantly less abundant at low fluxes, $<0.1$ sources per square degree, but dominate over accretors at higher flux levels ($\approx 10^{-12}$--$10^{-11}$ erg cm$^{-2}$s$^{-1}$). We suggest that the faint sources observed by ROSAT may be young cooling neutron stars with typical age $\la 10^6$ yrs, if the total number of young neutron stars in the Solar proximity is $\sim 10$ times higher than inferred from radiopulsars statistics.

In this paper we briefly review our recent results on evolution and properties of isolated neutron stars (INSs) in the Galaxy. As the first step we discuss stochastic period evolution of INSs. We briefly discuss how an INS's spin period evolves under influence of interaction with turbulized interstellar medium. To investigate statistical properties of the INS population we calculate a {\it census} of INSs in our Galaxy. Then we show that for exponential field decay the range of minimum value of magnetic moment, $\mu_b$: $\sim 10^{29.5}\ge \mu_b \ge 10^{28} {\rm G} {\rm cm}^3$, and the characteristic decay time, $t_d$: $\sim 10^8\ge t_d \ge 10^7 {\rm yrs}$, can be excluded assuming the standard initial magnetic momentum, $\mu_0=10^{30} {\rm G} {\rm cm}^3$, if accreting INSs are observed. It is shown that old accreting INSs become more abundant than young cooling INSs at X-ray fluxes below $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. The weak {\it ROSAT} sources, associated with INSs, can be young cooling objects, if the NSs birth rate in the solar vicinity during the last $\sim 10^6$ yr was much higher than inferred from radiopulsar observations.

We show, that different types of isolated neutron stars (INSs) show evidence in favor of magnetic field decay in these objects, and discuss how observations of INSs can help to constrain models of field decay.

In this paper we present calculations of period distribution for old accreting isolated neutron stars (INSs). At the age about a few billions years low velocity INSs come to the stage of accretion. At that stage their period evolution is governed by magnetic breaking and accreted angular momentum. Due to turbulence of the interstellar medium (ISM) accreted momentum can both accelerate and decelerate rotation of an INS and spin evolution has chaotic character. Calculations show that for constant magnetic fields INSs have relatively long spin periods, $\ge 10^4$--$10^5$ s, depending on parameters of INSs and an ISM. Due to long periods INSs have high spin up/spin down rates, which should fluctuate on a time scale about a few years.

In this paper we present calculations of period distribution for old accreting isolated neutron stars (INSs). After few billion years of evolution low velocity INSs come to the stage of accretion. At this stage INS's period evolution is governed by magnetic braking and angular momentum accreted. Since the interstellar medium is turbulized accreted momentum can either accelerate or decelerate spin of an INS, therefore the evolution of period has chaotic character. Our calculations show that in the case of constant magnetic field accreting INSs have relatively long spin periods (some hours and more, depending on INS's spatial velocity, its magnetic field and density of the surrounding medium). Such periods are much longer than the values measured by {\emph ROSAT} for 3 radio-silent isolated neutron stars. Due to long periods INSs should have high spin up/down rates, $\dot p$, which should fluctuate on a time scale of $\sim 1$ yr.

With a simple model from the point of view of population synthesis we try to verify an interesting suggestion made by Pfahl & Rappaport (2001) that dim sources in globular clusters (GCs) can be isolated accreting neutron stars (NSs). Simple estimates show, that we can expect about 0.5-1 accreting isolated NS per typical GC with $M=10^5 M_{\odot}$ in correspondence with observations. Properties of old accreting isolated NSs in GCs are briefly discussed. We suggest that accreting NSs in GCs experienced significant magnetic field decay.

The thermal, spin and magnetic evolution of neutron stars in the old low mass binaries is first explored. Recycled to very short periods via accretion torques, the neutron stars lose their magnetism progressively. If accretion proceeds undisturbed for 100 Myrs these stars can rotate close to break up with periods far below the minimum observed of 1.558 ms. We investigate their histories using population synthesis models to show that a tail should exist in the period distribution below 1.558 ms. The search of these ultrafastly spinning neutron stars as pulsars can help discriminating among the various equations of state for nuclear matter, and can shed light into the physics of binary evolution. The evolution of isolated neutron stars in the Galaxy is explored beyond the pulsar phase. Moving through the tenuous interstellar medium, these old solitary neutron stars lose their rotational energy. Whether also their magnetism fades is still a mystery. A population synthesis model has revealed that only a tiny fraction of them is able to accrete from the interstellar medium, shining in the X-rays. There is the hope that these solitary stars will eventually appear as faint sources in the Chandra sky survey. This might give insight on the long term evolution of the magnetic field in isolated objects.

In this article we briefly review our recent results on evolution and properties of isolated neutron stars (INSs) in the Galaxy. As the first step we calculate a {\it census} of INSs in our Galaxy. We infer a lower bound for the mean kick velocity of NSs, $ <V> \sim $(200-300) ${\rm km s^{-1}}$. The same conclusion is reached for both a constant magnetic field ($B\sim 10^{12}$ G) and for a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. These results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to $\sim$few percents. Then we show that the range of minimum value of magnetic moment, $\mu_b$: $\sim 10^{29.5}\ge \mu_b \ge 10^{28} {\rm G} {\rm cm}^3$, and the characteristic decay time, $t_d$: $\sim 10^8\ge t_d \ge 10^7 {\rm yrs}$, can be excluded assuming the standard initial magnetic momentum, $\mu_0=10^{30} {\rm G} {\rm cm}^3$, if accreting INSs are observed. For these parameters an INS would never reach the stage of accretion from the interstellar medium even for a low space velocity of the star and high density of the ambient plasma. The range of excluded parameters increases for lower values of $\mu_0$. It is shown that old accreting INSs become more abundant than young cooling INSs at X-ray fluxes below $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. We can predict that about one accreting INS per square degree should be observed at the {\it Chandra} and {\it Newton} flux limits of $\sim 10^{-16}$ erg cm$^{-2}$ s$^{-1}.$ The weak {\it ROSAT} sources, associated with INSs, can be young cooling objects, if the NSs birth rate in the solar vicinity during the last $\sim 10^6$ yr was much higher than inferred from radiopulsar observations.

In some models of $\gamma$-ray bursts super-strong electric fields ($E\sim 10^{14} {\rm statvolt cm}^{-1}$) have been surmized. Such large fields may provide copious pair production through vacuum polarization. Here we examine various astrophysical scenarios where huge electric field can be present. It is shown that when the conditions for quantistic vacuum breakdown are met, pair creation through conventional processes is likely to be always more important. Although of great physical interest, quantum pair production seems therefore of little astrophysical relevance.

The influence of exponential magnetic field decay (MFD) on the spin evolution of isolated neutron stars is studied. The ROSAT observations of several X-ray sources, which can be accreting old isolated neutron stars, are used to constrain the exponential and power-law decay parameters. We show that for the exponential decay the ranges of minimum value of magnetic moment, $\mu_b$, and the characteristic decay time, $t_d$, $\sim 10^{29.5}\ge \mu_b \ge 10^{28} {\rm G} {\rm cm}^3$, $\sim 10^8\ge t_d \ge 10^7 {\rm yrs}$ are excluded assuming the standard initial magnetic moment, $\mu_0=10^{30} {\rm G} {\rm cm}^3$. For these parameters, neutron stars would never reach the stage of accretion from the interstellar medium even for a low space velocity of the stars and a high density of the ambient plasma. The range of excluded parameters increases for lower values of $\mu_0$. We also show, that, contrary to exponential MFD, no significant restrictions can be made for the parameters of power-law decay from the statistics of isolated neutron star candidates in ROSAT observations. Isolated neutron stars with constant magnetic fields and initial values of them less than $\mu_0 \sim 10^{29} {\rm G} {\rm cm}^3$ never come to the stage of accretion. We briefly discuss the fate of old magnetars with and without MFD, and describe parameters of old accreting magnetars.

In this short note we discuss the influence of power-law magnetic field decay on the evolution of old accreting isolated neutron stars. We show, that, contrary to exponential field decay (Popov & Prokhorov 2000), no additional restrictions can be made for the parameters of power-law decay from the statistics of isolated neutron star candidates in ROSAT observations. We also briefly discuss the fate of old magnetars with and without field decay, and describe parameters of old accreting magnetars.

This paper is a continuation and development of our previous articles (Popov et al., 1997, 1998). We use `Scenario Machine'' (Lipunov et al., 1996b) -- the population synthesis simulator (for single binary systems calculations the program is available in WWW: http://xray.sai.msu.ru/ (Nazin et al., 1998)) -- to calculate evolution of populations of several types of X-ray sources during the first 20 Myrs after a starformation burst. We examined the evolution of 12 types of X-ray sources in close binary systems (both with neutron stars and with black holes) for different parameters of the IMF -- slopes: $\alpha=1$, $\alpha=1.35$ and $\alpha=2.35$ and upper mass limits, $M_{up}$: 120 $M_{\odot}$, 60 $M_{\odot}$ and 40 $M_{\odot}$. Results, especially for sources with black holes, are very sensitive to variations of the IMF, and it should be taken into account when fitting parameters of starformation bursts. Results are applied to several regions of recent starformation in different galaxies: Tol 89, NGC 5253, NGC 3125, He 2-10, NGC 3049. Using known ages and total masses of starformation bursts (Shaerer at al., 1998) we calculate expected numbers of X-ray sources in close binaries for different parameters of the IMF. Usually, X-ray transient sources consisting of a neutron star and a main sequence star are most abundant, but for very small ages of bursts (less than $\approx 4$ Myrs) sources with black holes can become more abundant.

I discuss the nature of the compact X-ray source in the center of the supernova remnant RCW 103. Several models, based on the accretion onto a compact object are analyzed. I show that it is more likely that the central X-ray source is an accreting neutron star than an accreting black hole. I also argue that models of a disrupted binary system consisting of an old accreting neutron star and a new one observed as a 69-ms pulsar are most favored.

Estimates of the magnetic field of neutron stars in X-ray pulsars are obtained using the hypothesis of the equilibrium period for disk and wind accretion and also from the BATSE data on timing of X-ray pulsars using the observed maximum spin-down rate. Cyclotron lines at energies $\ge 100$ keV in several Be-transient are predicted for future observations. We suggest a new method of estimating distances to X-ray pulsars and their magnetic fields. Using observations of fluxes and period variations in the model of disk accretion one can estimate the magnetic momentum of a neutron star and the distance to X-ray pulsar. As an illustration the method is applied to the system GROJ1008-57.

The paucity of old isolated accreting neutron stars in ROSAT observations is used to derive a lower limit on the mean velocity of neutron stars at birth. The secular evolution of the population is simulated following the paths of a statistical sample of stars for different values of the initial kick velocity, drawn from an isotropic Gaussian distribution with mean velocity $0\leq < V>\leq 550$ ${\rm km s^{-1}}$. The spin-down, induced by dipole losses and the interaction with the ambient medium, is tracked together with the dynamical evolution in the Galactic potential, allowing for the determination of the fraction of stars which are, at present, in each of the four possible stages: Ejector, Propeller, Accretor, and Georotator. Taking from the ROSAT All Sky Survey an upper limit of $\sim 10$ accreting neutron stars within $\sim 140$ pc from the Sun, we infer a lower bound for the mean kick velocity, $ < V>\ga 200-300$ ${\rm km s^{-1}}$. The same conclusion is reached for both a constant ($B\sim 10^{12}$ G) and a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. Present results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to $\la 1%$.

In this poster I briefly review several articles on astrophysics of old isolated neutron stars, which were published in 1994-99 by my co-authors and myself.

The paucity of old isolated accreting neutron stars in ROSAT observations is used to derive a lower limit on the mean velocity of neutron stars at birth. The secular evolution of the population is simulated following the paths of a statistical sample of stars for different values of the initial kick velocity, drawn from an isotropic Gaussian distribution with mean velocity $0\leq < V>\leq 550$ ${\rm km s^{-1}}$. The spin--down, induced by dipole losses and the interaction with the ambient medium, is tracked together with the dynamical evolution in the Galactic potential, allowing for the determination of the fraction of stars which are, at present, in each of the four possible stages: Ejector, Propeller, Accretor, and Georotator. Taking from the ROSAT All Sky Survey an upper limit of $\sim 10$ accreting neutron stars within $\sim 140$ pc from the Sun, we infer a lower bound for the mean kick velocity, $ < V>\gtrsim 200-300$ ${\rm km s^{-1}},$ corresponding to a velocity dispersion $\sigma_V\gtrsim 125-190$ km s$^{-1}$. The same conclusion is reached for both a constant magnetic field ($B\sim 10^{12}$ G) and a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. Such high velocities are consistent with those derived from radio pulsar observations. Present results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to less than 1%.

Estimates of the magnetic field of neutron stars in X-ray pulsars are obtained using the hypothesis of the equilibrium period for disk and wind accretion and also from the BATSE data on timing of X-ray pulsars using the observed maximum spin-down rate. Cyclotron lines at energies $\ge 100$ keV in several Be-transient are predicted for future observations.

In this paper we semianalyticaly evaluate influence of the exponential decay of magnetic field on the fate of isolated neutron stars. The fact of ROSAT observations of several X-ray sources, which can be accreting old isolated neutron stars gives us an opportunity to put some limits on the parameters of the exponential decay. We argue, that, if most part of neutron stars have approximately the same decay and initial parameters, then the combinations of the bottom magnetic momentum, $\mu_b$, in the range $\sim 10^{28}-10^{29.5} {\rm G} {\rm cm}^3$ and characteristic time scale, $t_d$, in the range $\sim 10^7-10^8 {\rm yrs}$ for standard initial magnetic momentum, $\mu_0=10^{30} {\rm G} {\rm cm}^3$, can be excluded, because for that sets of parameters neutron stars never come to the stage when accretion of the interstellar medium on their surfaces is possible even for low velocity of neutron stars and relatively high density of the interstellar medium. The region of excluded parameters increases with $\mu_0$ decreasing.

We suggest a new method of estimating distances to X-ray pulsars and their magnetic fields. Using observations of fluxes and period variations in the model of disk accretion one can estimate the magnetic momentum of a neutron star and the distance to X-ray pulsar. As an illustration the method is applied to the system GROJ1008-57. Estimates of the distance: approximately 6 kpc, and the magnetic momentum: approximately $4\cdot10^{31} E\cdot cm^3$, are obtained.

This paper is a continuation and development of our previous articles (Popov et al., 1997, 1998). We use ``Scenario Machine'' (Lipunov et al., 1996b) -- the population synthesis simulator (for single binary systems calculations the program is available in WWW at: http://xray.sai.msu.su/sciwork/scenario.html (Nazin et al., 1998) -- to calculate evolution of populations of several types of X-ray sources during the first 20 Myrs after a starformation burst. We examined the evolution of 12 types of X-ray sources in close binary systems (both with neutron stars and with black holes) for different parameters of the IMF -- slopes: $\alpha=1, \alpha=1.35$ and $\alpha=2.35$ and upper mass limits: 120 M_{\odot}, 60 M_{\odot} and 40 M_{\odot}. Results, especially for sources with black holes, are very sensitive to variations of the IMF, and it should be taken into account when fitting parameters of starformation bursts. Results are applied to several regions of recent starformation in different galaxies: Tol 89, NGC 5253, NGC 3125, He 2-10, NGC 3049. Using known ages and total masses of starformation bursts (Shaerer at al., 1998) we calculate expected numbers of X-ray sources in close binaries for different parameters of the IMF. Usially, X-ray transient sources consisting of a neutron star and a main sequence star are most abundant, but for very small ages of bursts (less than $\approx 4$ Myrs) sources with black holes can become more abundant.

In this article we use "Scenario Machine" - the population synthesis simulator- to calculate the evolution of the populations of the selected types of X-ray sources after the starformation burst with the total mass in binaries $1.5\cdot10^6 M_{\odot}$ during the first 20 Myr after the burst. Sources of the four types were calculated: transient sources- accreting neutron stars with Be- stars; accreting neutron stars in pair with supergiants; Cyg X-1-like sources- black holes with supergiants; superaccreting black holes. We used two values of the $\alpha - coefficient$ in the mass-function: 2.35 (Salpeter's function) and 1.01 ("flat spectrum"). The calculations were made for two values of the upper limit of the mass-function: 120 and 30 $M_{\odot}$. For the flat spectrum the number of sources of all types significantly increased. Decreasing of the upper mass limit below the critical mass of a black hole formation increase (for the "flat spectrum") the number of transient sources with neutron stars up to $\approx 300$. We give approximating formulae for the time dependence of source numbers.

We studied possible evolution of the rotational period and the magnetic field of the X-ray source RX J0720.4-3125 assuming this source to be an isolated neutron star accreting from the interstellar medium. Magnetic field of the source is estimated to be $10^6 - 10^9$ Gs (most probably $\approx 2\cdot 10^8$ Gs), and it is difficult to explain the observable rotational period 8.38 s without invoking hypothesis of the magnetic field decay. We used the model of ohmic decay of the crustal magnetic field. The estimates of the accretion rate ($10^{-14} - 10^{-16} M_\odot/yr$), velocity of the source relative to the interstellar medium ($10 - 50 $ km/s), the neutron star age ($2\cdot 10^9 - 10^{10}$ yrs) are obtained.

In this poster I discuss the nature of the compact X-ray source in the center of the supernova remnant RCW 103. Several models, based on the accretion onto a compact object such as a neutron star or a black hole (isolated or binary), are analyzed. I show that it is more likely that the central X-ray source is an accreting neutron star than an accreting black hole. I also argue that models of a disrupted binary system consisting of an old accreting neutron star and a new one observed as a 69-ms X-ray and radio pulsar are most favored.

I discuss the nature of the compact X-ray source inside the supernova remnant RCW 103. Several models, based on the accretion onto a compact object such as a neutron star or a black hole (isolated or binary), are analyzed. I show that it is more likely that the X-ray source is an accreting neutron star than an accreting black hole. I also argue that models of a binary system with an old accreting neutron star are most favored.

We use "Scenario Machine" -- the population synthesis simulator -- to calculate the evolution of populations of selected types of X-ray sources after a starformation burst with the total mass in binaries (1--1.5) \cdot 10^6 M_{\odot}$ during the first 20 Myr after a burst. We present here the results of two sets of runs of the program. In the first set we examined the following types of close binaries: transient sources-- neutron stars with Be- stars; ``X-ray pulsars''-- neutron stars in pairs with supergiants; Cyg X-1-like sources-- black holes with supergiants; ``SS443-like sources''-- superaccreting black holes. We used two values of the exponent $\alpha$ in the initial mass function: 2.35 (Salpeter's function) and 1 (``flat spectrum''). The calculations were made for the folowing values of the upper limit of the mass function: 120 and 30 $M_{\odot}$. For the ``flat spectrum'', suggested in (Contini et al, 1995), the number of sources of all types significantly increased. With the ``flat spectrum'' and with the upper mass limit 120 $M_{\odot}$ we obtained hundreds of sources of all calculated types. Decreasing of the upper mass limit below the critical mass of a black hole formation increase the number of transient sources with neutron stars up to $\approx 300$. In the second set we examined the evolution of of 12 other types of X-ray sources for $\alpha=1$, $\alpha=1.35$ and $\alpha=2.35$ and for three upper mass limits: 120 $M_{\odot}$, 60 $M_{\odot}$, 40 $M_{\odot}$ (see Perez-Olia & Colina 1995 for the reasons for such upper limits) on the same time scale 20 Myr after a star formation burst.

We studied possible evolution of the rotational period and the magnetic field of the X-ray source RX J0720.4-3125 assuming this source to be an isolated neutron star accreting interstellar medium. Magnetic field of the source is estimated to be $10^6 - 10^9$ G, and it is difficult to explain observed rotational period 8.38 s without invoking hypothesis of the magnetic field decay. We used the model of ohmic decay of the crustal magnetic field. The estimates of accretion rate ($10^{-14} - 10^{-16} M_\odot/yr$), velocity of the source relative to interstellar medium ($10 - 50 $ km/s), neutron star age ($2\cdot 10^9 - 10^{10}$ yrs) are obtained.

We present here the computer model of the distribution of the luminosity, produced by old isolated neutron stars and black holes accreting from the interstellar medium in the Galaxy in R-Z plane. We show, that the luminosity distributions have a torus-like structure, with the maximum at $\approx 5-6 kpc$.

In this work we give detail consideration of the possible scenario of evolution of isolated neutron stars and determine some characteristics of X-ray pulsars from their spin period evolution (the paper containes more detailed abstract).

We present here the computer model of the distribution of the luminosity, produced by old isolated neutron stars (OINSs) accreting from the interstellar medium (ISM). We show, that for different mean velocities of OINSs the distribution of the luminosity has a torus-like structure, with the maximum at $\approx 5 kpc$.

We present a WWW-version of the {\it Scenario Machine} - a computer code designed to calculate the evolution of close binary stellar systems. The Internet users can directly access to the code and calculate binary evolutionary tracks with parameters at the user's will. The program is running on the {\it Pentium} server of the Division of the Relativistic Astrophysics of the Sternberg Astronimical Institute (http://xray.sai.msu.su/ ). The results are presented both in the form of tables and graphic diagrams. The work is always in progress. More possibilities for Internet users are intended to become available in the near future.

We model the evolution of the selected types of X-ray sources (both with BHs and NSs) for the star formation burst at the Galactic center. For the currently assumed starburst age of 4-7 millions years, our results are consistent with the recent {\it GRANAT} X-ray observations of the Galactic center and predict a substantial number of X-ray binaries with black holes. In addition we made a very simple model of the spatial distribution of the X-ray sources, which is also consistent with the {\it GRANAT} observations. The results can be applied to other starformation regions.

X-ray binary population at the Galactic center is proposed here to form as a result of a starburst, which is synthesized by using the ``Scenario Machine'' code of binary star evolution. For the currently assumed starburst age of 4-7 millions years, our results are consistent with the recent {\it GRANAT} X-ray observations of the Galactic center and predict a substantial number of X-ray binaries with black holes.

The distribution of neutron stars (NS's) is determined by stationary solution of the Fokker-Planck equation. In this work using the observed period changes for four systems: Vela X-1, GX 301-2, Her X-1 and Cen X-3 we determined D, the 'diffusion coefficient',-parameter from the Fokker-Planck equation. Using strong dependence of D on the velocity for Vela X-1 and GX 301-2, systems accreting from a stellar wind, we determined the stellar wind velocity. For different assumptions for a turbulent velocity we obtained $V=(660-1440) km s ^{-1}$. It is in good agreement with the stellar wind velocity determined by other methods. We also determined the specific characteristic time scales for the 'diffusion processes' in X-ray pulsars. It is of order of 200 sec for wind-fed pulsars and 1000-10000 sec for the disk accreting systems.