Possible models for the generation of electromagnetic (EM) radiation during the coalescence of neutron star-black hole binaries are considered. The mass of the remnant disk around the black hole during the coalescence of neutron stars and black holes is calculated by taking into account the equation of state for neutron stars and the rotation of the binary components before the coalescence. The parameters of binary systems before the coalescence (the mass ratio, the component rotation, the neutron star magnetic field) are calculated by the population synthesis method. The derived mass of the remnant disk around the black hole after the coalescence is used to estimate the kinetic energy of the relativistic jet launched by the Blandford-Znajek mechanism. A disk mass of more than $\sim 0.05 M_\odot$ required for the formation of short gamma-ray bursts is shown to be obtained in no more than 1-10\% of the coalescences (depending on the equation of state). Less efficient common envelopes (a large parameter $\alpha_{CE}$) lead to a noticeably larger percentage of events with astrophysically interesting EM energy release. For binaries with a large mass ratio, in which a magnetized neutron star is not subjected to tidal disruption before the coalescence, the possibility of the formation of an electrically charged rotating black hole (Wald charge) is considered and estimates of the maximum EM power released by such a black hole after the coalescence are made. The conversion of the emitted gravitational waves into electromagnetic ones in the relativistic lepton plasma generated in coalescing pulsar-black hole binaries at the pre-coalescence stage is discussed.

The available data on the chirp mass distribution of the black holes in the coalescing binaries in O1-O3 LIGO/Virgo runs are analyzed and compared with theoretical expectations based on the hypothesis that these black holes are primordial with log-normal mass spectrum. The inferred best-fit mass spectrum parameters, $M_0=17 M_\odot$ and $\gamma=0.9$, fall within the theoretically expected range and shows excellent agreement with observations. On the opposite, binary black hole models based on massive binary star evolution require additional adjustments to reproduce the observed chirp mass distribution.

We model the Galactic population of detached binaries that harbor black holes with (0.5-1.7) solar mass companions -- the remnants of case B mass-exchange that rapidly cross Hertzsprung gap after the termination of the Roche-lobe overflow or as He-shell burning stars. Several such binaries can be currently present in the Galaxy. The range of black hole masses in them is about 4 to 10 solar ones, the orbital periods are tens to hundreds day. The unique BH-binary LB-1 fits well into this extremely rare class of double stars.

Modern astrophysical methods of determination of spins of rotating stellar-mass black hole in close binaries and of supermassive black holes in active galactic nuclei are briefly discussed. Effective spins of coalescing binary black holes derived from LIGO/Virgo gravitational wave observations are specially addressed. We consider three types of coalescing binaries: double black holes, black hole-neutron star binaries and primordial double black holes. The effective spins of coalescing astrophysical binary black holes and black holes with neutron stars are calculated for two plausible models of black hole formations from stellar core collapses (without or with additional fallback from the stellar envelope) taking into account the stellar metallicity and star formation rate evolution in the Universe. The calculated distributions do not contradict the reported LIGO/Virgo observations. The effective spins of primordial coalescing stellar-mass black holes can reach a few per cent due to accretion spin-up in the cold external medium.

We present a new model of the population of symbiotic X-ray binaries (SyXBs) that takes into account non-stationary character of quasi-spherical sub-sonic accretion of the red giant's stellar wind onto slowly rotating neutron stars. Updates of the earlier models are given, which include more strict criteria of slow NS rotation for plasma entry into the NS magnetosphere via Rayleigh-Taylor instability, as well as more strict conditions for settling accretion for slow stellar winds, with an account of variations in the specific angular momentum of captured stellar wind in eccentric binaries. These modifications enabled a more adequate description of the distributions of observed systems over binary orbital periods, NS spin periods and their X-ray luminosity in the $\sim 10^{32}-10^{36}$~erg s$^{-1}$ range and brought their model Galactic number into reasonable agreement with the observed one. Reconciliation of the model and observed orbital periods of SyXBs requires a low efficiency of matter expulsion from common envelopes during the evolution that results in the formation of NS-components of symbiotic X-ray systems.

The possible formation mechanisms of massive close binary black holes (BHs) that can merge in the Hubble time to produce powerful gravitational wave bursts detected during advanced LIGO O1 and O2 science runs include the evolution from field low-metallicity massive binaries, the dynamical formation in dense stellar clusters and primordial BHs. Different formation channels produce different source distributions of total masses ${M_\mathrm{tot}}$ and effective spins $\chi_\mathrm{eff}$ of coalescing binary BHs. Using a modified \textsc{bse} code, we carry out extensive population synthesis calculations of the expected effective spin and total mass distributions from the standard field massive binary formation channel for different metallicities of BH progenitors (from zero-metal Population III stars up to solar metal abundance), different initial rotations of the binary components, stellar wind mass loss prescription, different BH formation models and a range of common envelope efficiencies. The stellar rotation is treated in two-zone (core-envelope) approximation using the effective core-envelope coupling time and with an account of the tidal synchronization of stellar envelope rotation during the binary system evolution. The results of our simulations, convolved with the metallicity-dependent star-formation history, show that the total masses and effective spins of the merging binary black holes detected during LIGO O1-O2 runs but the heaviest one (GW170729) can be simultaneously reproduced by the adopted BH formation models. Noticeable effective spin of GW170729 requires additional fallback from the rotating stellar envelope.

Different accretion regimes onto magnetized NSs in HMXBs are considered: wind-fed supersonic (Bondi) regime at high accretion rates $\dot M\gtrsim 4\times 10^{16}$~[g s$^{-1}$] , subsonic settling regime at lower $\dot M$ and supercritical disc accretion during Roche lobe overflow. In wind-fed stage, NSs in HMXBs reach equilibrium spin periods $P^*$ proportional to binary orbital period $P_b$. At supercritical accretion stage, the system may appear as a pulsating ULX. Population synthesis of Galactic HMXBs using standard assumptions on the binary evolution and NS formation is presented. Comparison of the model $P^*$ -- $P_b$ (the Corbet diagram), $P^*$ -- $L_x$ and $P_b$ -- $L_x$ distributions with those for the observed HMXBs (including Be X-ray binaries) and pulsating ULXs suggests the importance of the reduction of $P^*$ in non-circular orbits, explaining the location of Be X-ray binaries in the model Corbet diagram, and the universal parameters of pulsating ULXs depending only on the NS magnetic fields.

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.

Possible formation mechanisms of massive close binary black holes that can merge in the Hubble time to produce powerful gravitational wave bursts detected during advanced LIGO O1 science run are briefly discussed. The pathways include the evolution from field low-metallicity massive binaries, the dynamical formation in globular clusters and primordial black holes. Low effective black hole spins inferred for LIGO GW150914 and LTV151012 events are discussed. Population synthesis calculations of the expected spin and chirp mass distributions from the standard field massive binary formation channel are presented for different metallicities (from zero-metal Population III stars up to solar metal abundance). We conclude that that merging binary black holes can contain systems from different formation channels, discrimination between which can be made with increasing statistics of mass and spin measurements from ongoing and future gravitational wave observations.

Using population synthesis, we study a double-degenerate (DD) scenario for SNe Ia, aiming to estimate the maximum possible contribution to the rate of SNe from this scenario and the dependence of the delay-time distribution (DTD) on it. We make an extreme assumption that all mergers of super-Chandrasekhar pairs of CO white dwarfs (WDs) and mergers of CO WDs more massive than 0.47 $M_\odot$ with hybrid or helium WDs more massive than 0.37$M_\odot$ produce SNe Ia. The models are parametrized by the product of the common envelope efficiency and the parameter of binding energy of stellar envelopes $\alpha_{ce}\lambda$, which we vary between 0.25 and 2. The best agreement with observations is obtained for $\alpha_{ce}\lambda$=2. A substantial contribution to the rate of SNe Ia is provided by the pairs with a hybrid WD. The estimated Galactic rate of SNe Ia is $6.5 10^{-3}$ per yr (for the mass of the bulge and thin disk equal to $7.2 10^{10} M_\odot$), which is comparable to the observational estimate $(5.4\pm0.12) 10^{-3}$ per yr. The model DTD for 1 to 8 Gyr range satisfactorily fits DTD for SNe Ia in the field galaxies (Maoz et al. 2012). For this epoch, model DTD is a power low with an index -1.64. At earlier and later epochs our DTD has a deficit of events, as in other studies. Marginal agreement with observational DTD is achieved even if only CO+CO WD with $M_1\geq0.8\,M_\odot$ and $M_2\geq0.6\,M_\odot$ produce SNe Ia. A better agreement of observed and model DTD may be obtained if tidal effects are weaker than assumed and/or metallicity of population is much lower than solar.

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.

Symbiotic X-ray binaries (SyXBs) comprise a rare class of low-mass X-ray binaries. We study the Galactic SyXBs, which we consider as detached binaries composed of low-mass giants and wind-fed neutron star companions, by simulation of the interaction of a magnetized neutron star (NS) with its environment and utilizing a population synthesis code. We focus mainly on the parameters that influence observational appearance of the SyXB: the donor wind velocity (vw) and the angular momentum distribution in the shell of matter settling onto NS. We estimate the birthrate of SyXB as $\sim 4.1\times 10^{-5}$ yr$^{-1}$ to $ \sim 6.6\times 10^{-6}$ yr$^{-1}$ and their number in the Galaxy as $\sim$(100 -- 1000). Assumed stellar wind velocity from cool giants is the input parameter that influences the model SyXBs population most. Among known SyXBs or candidate systems, 4U 1954+31 and IGR J16358-4724 in which NS have very long spin periods may host quasi-spherically accreting NSs. GX 1+4 has a peculiar long-term spin behaviour and it may also be a quasi-spherical wind-accreting source. We cannot identify whether there are wind-fed accretion disks in 4U 1700+24, Sct X-1, IRXS J180431.1-273932 and 2XMM J174016.0-290337.

We analyze a flux-limited sample of persistent and bright (with 2-10 keV fluxes exceeding 1.4e-10 erg/s/cm2) low-mass X-ray binaries (LMXBs) in our Galaxy. It is demonstrated that the majority of binary systems with X-ray luminosities below logL(erg/sec)~37.3 have unevolved secondary companions (except for those with white dwarf donors), while systems with higher X-ray luminosity predominantly harbor giant donors. Mass transfer in binary systems with giants significantly shortens their life time thus steepening the X-ray luminosity function of LMXBs at high luminosity. We argue that this is the reason why the LMXB luminosity function constructed in the last years from observations of sources in our and distant galaxies demonstrates a break at logL(erg/sec)~37.3.

Coalescing neutron star binaries are believed to be the most reliable sources for ground-based detectors of gravitational waves and likely progenitors of short gamma-ray bursts. In the process of coalescence, magnetic fields of neutron stars can induce interesting observational manifestations and affect the form of gravitational wave signal. In this papaer we use the population synthesis method to model the expected distribution of neutron star magnetic fields during the coalescence under different assumptions on the initial parameters of neutron stars and their magnetic field evolution. We discuss possible elecotrmagnetic phenomena preceding the coalescence of magnetized neutron star binaries and the effect of magnetic field on the gravitational wave signal. We find that a log-normal (Gaussian in logarithms) distribution of the initial magnetic fields of neutron stars, which agrees with observed properties of radio pulsars, produces the distribution of the magnetic field energy during the coalescence that adequately describes the observed luminosity function of short gamma-ray bursts under different assumptions on the field evolution and initial parameters of neutron stars. This agreement lends further support to the model of coalescing neutron star binaries as progenitors of gamma-ray bursts.

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.

We study the effect of the neutron star spin -- kick velocity alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. Two scenarios of the neutron star formation are considered: when the kick is always present and when it is small or absent if a neutron star is formed in a binary system due to electron-capture degenerate core collapse. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the expected galactic rate of binary neutron star coalescences compared to calculations with randomly directed kicks. The spin-kick correlation also leads to a much narrower NS spin-orbit misalignment.

We study the effect of the kick velocity -- neutron star spin alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the galactic rate of binary neutron star coalescences up to an order of magnitude with respect to the rates calculated for random kicks. The spin-kick correlation also leads to much narrower NS spin-orbit misalignments compared to random kicks.

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.

Population synthesis is used to model the number of neutron stars in globular clusters that are observed as LMXBs and millisecond PSRs. The dynamical interaction between binary and single stars in a GC are assumed to take place with a permanently replenished "background" of single stars whose density distribution keeps track with the cluster evolution as a whole and evolution of single stars. We use the hypothesis (Podsiadlowski et al) that NS forming in binary systems from components with initial masses \sim 8-12 M_\odot during the electron-capture collapse of the degenerate O-Ne-Mg core do not acquire a high space velocities (kicks). The remaining NSs (i.e. from single stars with M>8 M_\odot or binary comonents with M>12 M_\odot) are assumed to be born with high kicks, as found from obsrevations of single pulsars (Hobbs et al. 2005). Under this assumption, a sizeable fraction of NSs remain in GCs (about 1000 NSs in a GC with a mass of 5\times 10^5 M_\odot). The number of ms PSRs formed in the cluster via accretion spin-up in binaries is then about 10, which is consistent with observations. Our modelling reproduces the observed shape of the X-ray luminosity function for accreting NSs in binaries with normal and degenerate components and the distribution of spin periods of ms PSRs in GCs under the assumption of accretion-driven magnetic field decay of NSs up to a bottom value of 10^8 G. The number of LMXBs and ms PSRs dynamically expelling from GCs is also calculated.

X-ray luminosity function derived from observations of X-ray sources in galactic bulges can be explained by prinicipal evolutionary relations for mass accretion rate onto the compact object. The observed mean distribution of individual X-ray luminosities of galactic LMXB is satisfactorily described by a symmetric quasi-Lorentzian curve. The flux variance for bright sources is found to be proportional to the mean luminosity. Such a distribution does not change the slope of the power law luminosity function of the source population, which is expected from the dependence of the mass transfer rate on the mass of the Roche-lobe filling non-degenerate optical component.

The population synthesis method is used to study the possibility of explaining the appreciable fraction (20^+12_15%) of the intergalactic (no-host) type Ia supernovae observed in galaxy clusters (Gal-Yam ete al. 2003) by binary whote dwarf merginngs in the cores of globular clusters. In a typical globular cluster, the number of merging double white dwarfs is fount to be smaller than 10^{-13} per year per average cluster star during the entire evolution of the cluster, which is a factor of 3 higher than in a Milky-Way-type galaxy. From 5 to 30% of the merging white dwarfs are dynamically expelled from the cluster with barycenter velocities up to 150 km/s. SN Ia explosions during the mergers of binary white dwarfs in dense star clusters may account for \sim 1% of the total rate of SN Ia in the central parts of galaxy clusters if the baryon mass fraction in such star clusters is \sim 0.3%.

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.

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.