The principal goal of this paper is to specify conditions of detection of unresolved binaries by multicolor photometry. We have developed a method for estimating the critical distance at which an unresolved binary of given mass and age can be detected. The method is applied to the photometric system of the planned Lyra-B spaceborne experiment. We have shown that some types of unresolved binary stars can be discovered and distinguished from single stars solely by means of photometric observations.

We propose an improved method for the atmospheric extinction reduction within optical photometry. Our method is based on the simultaneous multicolor observations of photometric standards. Such data are now available within the modern wide-field sky surveys and contain a large amount of information about instant atmospheric conditions. So, it became possible to estimate the extinction parameters on the basis of a quite short observational dataset and, hence, to trace the rapid stars twinkling accurately. Having been developed for a new MiniMegaTORTORA observational system, the proposed method can be adopted for a wide range of modern observational programs.

The space photometric project "Lyra" is developed now in Russia.The project purpose is determination of the photometric information and coordinates of the natural and artificial space objects, from the brightest ones to 16^m in visual lights. It is supposed to obtain the data for about 40-400 million objects from board of International Space Station, using an astronomical telescope with diameter of the main mirror of 0.5 m. The observations will be carried out in a scanning mode. Photometry will be obtained in 10 spectral bands. The expected uncertainty of magnitudes for objects of 16^m in the V-band is 0.001^m. The main results of experiment should become: 1) creation of spatial model of the Galaxy on distances to 3 kpc from the Sun; 2) specification of physical parameters of stars and models of star evolution; 3) discovering of a huge number (to several millions) of variable stars and determination of their variability parameters. The central wavelengths of the 10 bands of "Lyra" photometric system will be at 195, 218, 270, 350, 440, 550, 700, 825, 930 and 1000 nm. It is shown that combinations of various colour indices will allow to determine confidently both effective temperature and metallicity of stars. The presence of a 218 nm band allows to determine confidently interstellar extinction on stars of O-F spectral classes. The photometric system will make it possible to separate halo stars from disk stars and to define physical parameters of their atmospheres.

We present here an updated description of the "Scenario Machine" code. This tool is used to carry out a population synthesis of binary stars. Previous version of the description can be found at http://xray.sai.msu.ru/~mystery//articles/review/contents.html

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 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.

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 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.

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 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.

Maximum gravitational wave emission from binary stars in eccentric orbits occurs near the periastron passage. We show that for a stationary distribution of binary neutron stars in the Galaxy, several high-eccentricity systems with orbital periods in the range from tens of minutes to several days should exist that emit broad gravitational-wave pulses in the frequency range 1-100 mHz. The space interferometer LISA could register the pulsed signal from these system at a signal-to-noise ratio level $S/N>5\sqrt{5}$ in the frequency range $\sim 10^{-3}-10^{-1}$ Hz during one-year observational time. Some detection algorithms for such a signal are discussed.

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.

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.

The observed mass distribution for the compact remnants of massive stars (neutron stars and black holes) and its relationship to possible mechanisms for the ejection of the envelopes of type II and Ib/c supernovae is analyzed. The conclusion is drawn that this distribution can be obtained only by a magneto-rotational mechanism for the supernovae with sufficiently long time of the field amplification, and a soft equation of state for neutron stars with limiting masses \sim1.5-1.6M_\odot. Some consequences of this hypothesis are discussed.

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.

Unresolved coalescing binary neutron stars in the Galaxy and beyond produce a stochastic gravitational wave background within the LISA frequency band, which can be potentially dangerous for possible detection of fundamentally more interesting relic cosmological backgrounds. Here we address the question what effects the unavoidable eccentricity of orbits of these systems should have on the properties of this backgrounds. In particularly, we examine starting from which frequency one-year observations with LISA can be secure from contamination by the noise produced by galactic binary neutron stars in eccentric orbits. We come to the conclusion that harmonics from such binaries do not contribute significantly above $\sim 10^{-3}$ Hz.

Stochastic gravitational wave background produced by a stationary coalescing population of binary neutron stars in the Galaxy is calculated. This background is found to constitute a confusion limit within the LISA frequency band up to a limiting frequency $\NUlim{}\sim 10^{-3}$ Hz, leaving the frequency window $\sim 10^{-3}$--$10^{-2}$ Hz open for the potential detection of cosmological stochastic gravitational waves and of signals involving massive black holes out to cosmological distances.

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 first generation of long-baseline laser interferometric detectors of gravitational waves will start collecting data in 2001-2003. We carefully analyse their planned performance and compare it with the expected strengths of astrophysical sources. The scientific importance of the anticipated discovery of various gravitatinal wave signals and the reliability of theoretical predictions are taken into account in our analysis. We try to be conservative both in evaluating the theoretical uncertainties about a source and the prospects of its detection. After having considered many possible sources, we place our emphasis on (1) inspiraling binaries consisting of stellar mass black holes and (2) relic gravitational waves. We draw the conclusion that inspiraling binary black holes are likely to be detected first by the initial ground-based interferometers. We estimate that the initial interferometers will see 2-3 events per year from black hole binaries with component masses 10-15M_\odot, with a signal-to-noise ratio of around 2-3, in each of a network of detectors consisting of GEO, VIRGO and the two LIGOs. It appears that other possible sources, including coalescing neutron stars, are unlikely to be detected by the initial instruments. We also argue that relic gravitational waves may be discovered by the space-based interferometers in the frequency interval 2x10^{-3}-10^{-2} Hz, at the signal-to-noise ratio level around 3.

A careful analysis of X-ray light curves and pulse profiles of Her X-1 obtained over more than 20 years strongly evidences for free precession of a magnetized neutron star with rotational axis inclined to the orbital plane as a central clock underlying the observed 35-day period. Strong asymmetric X-ray illumination of the optical star atmosphere leads to the formation of gaseous streams coming out of the orbital plane and forming a tilted accretion disk around the neutron star. Such a disk precesses due to tidal forces and dynamical action of gaseous streams from the secondary companion. The locking of these torques with neutron star precession makes the net disk precession period to be very close to that of the neutron star free precession.

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.

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.

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%$.

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%.

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 show that the distribution of observed energies of GRB with known redshifts can be explained by the hypothesis of the standard energy release E_0=5 10^51 ergs. Two situations are possible, either the beaming angle differs from burst to burst, or there is a universal emssion diagram in each burst, the observed difference being due to different viewing angles.

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.

The formation and evolution of binary black holes (BH) is studied using the modern evolutionary scenario for very massive stars with high mass loss (Vanbeveren et al. 1998). Main sequence stars with masses $M>35 M_\odot$ are assumed to form a BH in the end of their nuclear evolution. The mass of BH formed is parametrized as $M_{bh}=k_{bh}\times M_{SN}$, where $M_{SN}$ is the mass of the pre-supernova star taken from evolutionary calculations, $k_{bh}\le 1$. The possibility is explored that a newly formed BH acquires a kick velocity 0-250 km/s. Binary BH are found to merge within the Hubble time at an appreciable rate only for non-zero kick velocities. We calculate the galactic merging rates of binary BH systems, their detection rate by the initial laser interferometers, and the distributions of merging binary BH over orbital eccentricities at different frequencies. The distribution of angles between BH spins and the orbital angular momentum is also presented.

A strong X-ray illumination of the optical star atmosphere in Her X-1, asymmetric because of a partial shadowing by the tilted twisted accretion disk around central neutron star, leads to the formation of matter flows coming out of the orbital plane and crossing the line of sight before entering the disk. We suggest that the absorption of X-ray emission by this flow leads to the formation of pre-eclipse and anomalous dips of type I. These dips are observed during several orbits after turn-on both in the main-on and short-on state. Almost coherent action of tidal torques and matter streams enhances the disk wobbling which causes the disk edge to shield the X-ray source after the turn-on. Anomalous dips of type II and post-eclipse recovery appear due to this process only on the first orbit after turn-on.

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 consider the dependence of the internal structure of a neutron star in a close binary system on the semi-major axis of the binary orbit, focusing on the case when the Roche lobes of the components are nearly filled. We adopt a polytropic equation of state. The temporal evolution of the semi-major axis and its dependence on the mass ratio of the binary components and the polytropic index are determined. The calculation are carried out right up to the moment of contact, when quasi-stationary model becomes invalid. We analyze differences in the shapes of the pulses of gravitational radiation emitted by a pair of point masses and by a binary neutron star, taking into account its internal structure and tidal deformations.

An analysis of the publically available RXTE archive on Her X-1 including data on 23 34.85-day cycles is performed. The turn-on times for these cycles are determined. The number of cycles with a duration of 20.5 orbits has been found exceedingly larger than of shorter (20 orbits) or longer (21 orbits). A correlation between the duration of a cycle and its mean X-ray flux is noted. The mean X-ray light curve shows a very distinct short-on state. The anomalous X-ray absorption dip is found during the first orbit after the turn-on in the main-on state for the cycles starting near the binary phase 0.25, while is present during two successive orbits in the low-on state. The post-eclipse recovery feature have not been found in the main-on state but appears at least for two orbits during the low-on state. The pre-eclipse dips are present both in main-on and low-on state and demonstrate the behaviour as in early observations. The comparison of durations of the main-on and short-on states enabled us to constrain the accretion disk semi-thickness and its inclination to the orbital plane.

We show that the free precession of a triaxial body can naturally explain the anomalously rapid change of the X-ray pulse profile of Her X-1 observed by the HEAO-1 in September 1978 without requiring a large change in the moment of inertia.

Gravitational wave noise associated with unresolved binary stars in the Galaxy is studied with the special aim of determining the upper frequency at which it stops to contribute at the rms noise level of the proposed space-born interferometer (LISA). The upper limit to this background is derived from the statistics of SN Ia explosions, part of which can be triggered by binary white dwarf coalescences. The upper limiting frequency at which binary stochastic noise crosses LISA rms sensitivity is found to lie within the range 0.03-0.07 Hz, depending on the galactic binary white dwarf coalescence rate. To be reliably detectable by LISA, the energy density of relic cosmological background per logarithmic frequency interval should be Omega_{GW}h_{100}^2>10^{-8} at f>0.03 Hz.

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 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$.

Assuming binary neutron star mergings as a standard-candle model for GRBs, an independent estimate is obtained for the redshift of GRB970228 and GRB970508 with optical counterparts, using mean statistical properties of GRBs observed by BATSE. We derive z=0.7\pm 0.1 and z=1.9\pm 0.1 for GRB970228 and GRB970508 respectively, depending on the power-law index of the GRB spectrum s=-1.1\pm 0.3 and the value of the redshift of the initial star formation z_*=3-10 in a flat \Omega=1 Universe with a cosmological term \Lambda=0.7.

Merging rates of compact binaries (double neutron stars or black holes) are calculated based on the modern concept of binary stellar evolution. It is found that the initial laser interferometers with an rms-sensitivity of $10^{-21}$ at the frequency 100 Hz can detect 10-700 black holes and only \sim 1 neutron star coalescences in a 1-year integration time. Implications of the evolutionary effects to the cosmological origin of GRB are also discussed.

Galactic binary pulsars with different companion types (OB-star, white dwarf (WD), neutron star (NS), black hole (BH), or planet) assuming various phenomenological distributions for kick velocities of newborn NS are studied. The merging rates of relativistic binary stars (NS+NS, NS+BH, BH+BH) dependence on the kick velocity is specially adresed. The BH+BH merging is shown to occur, depending on parameters of BH formation, at a rate of one per 200,000 -- 500,000 years in a Milky Way type galaxy. The NS+NS merging rate $R_{ns}$ is found to be 1 per \sim 3,000 years for zero recoil, and decrease to one per 10,000 years even for the highest kick velocities of 400 km/s. That the merging rates derived from evolutionary calculations are by two order of magnitude higher than those based on binary pulsar statistics only, is suggested to be due to the fact that the observable binary pulsars in pairs with NS form only a fraction of the total number of binary NS systems. The expected detection rate of binary BH by a LIGO-type gravitational wave detector comparable with and even higher than the binary NS merging rate for a wide range of parameters.

Study of gravitational-radiation induced merging rates of relativistic binary stars (double neutron stars; neutron star + black hole; double black holes) shows that the first-generation gravitational wave interferometers with an rms-sensitivity of $10^{-21}$ at frequency 100 Hz can detect 10-700 black hole and only $\sim 1$ neutron star coalescences in a 1-year integration time in a wide range of stellar evolution parameters. It is notable that modern concepts of stellar evolution predict that the first detection of gravitational wave will independently discover black holes.

Based on evolutionary scenarios for binary stellar evolution we study the merging rates of relativistic binary stars (NS+NS, NS+BH, BH+BH) under different assumptions of BH formation. We find the BH+BH merging rate in the range one per 200,000 -- 500,000 year in a Milky-Way type galaxy, while the NS+NS merging rate $R_{ns}$ is approximately 10 times as high, which means that the expected event rate even for high mean kick velocities of NS up to 400 km/s is at least 30-50 binary NS mergings per year from within a distance of 200 Mpc. As typical BH is formed with a mass 3-10 times the NS mass (assumed 1.4 M$_\odot$), the rates obtained imply that the expected detection rate of binary BH by a LIGO-type gravitational wave detector is 10-100 times higher than the binary NS merging rate for a wide range of parameters.

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.

Evolution of the coalescence rate of double neutron stars (NS) and neutron star -- black hole (BH) binaries are computed for model galaxies with different star formation rates. Assuming gamma-ray bursts (GRB) to originate from NS+NS or NS+BH merging in distant galaxies, theoretical logN--logS distributions and <V/V_max> tests of gamma-ray bursts (GRB) are calculated for the first time taking the computed merging rates into account. We use a flat cosmological model (Omega=1) with different values of the cosmological constant Lambda and under various assumptions about the star formation history in galaxies. The calculated source evolution predicts a 5-10 times increase of the source statistics at count rates 3-10 times lower than the exising BATSE sensitivity limit. The most important parameter in fitting the 2nd BATSE catalogue is the initial redshift of star formation, which is found to be z_*=2-5 depending on a poorly determined average spectral index of GRB.