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.

Novae are the observational manifestations of thermonuclear runaways on the surface of accreting white dwarfs (WDs). Although novae are an ubiquitous phenomenon, their properties at low metallicity are not well understood. Using the publicly-available stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we model the evolution of accreting carbon-oxygen WDs and consider models which accrete matter with metallicity Z=0.02 or $10^{-4}$. We consider both models without mixing and with matter enriched by CO-elements assuming that mixing occurs in the process of accretion (with mixing fraction 0.25). We present and contrast ignition mass, ejected mass, recurrence period and maximum luminosity of novae for different WD masses and accretion rates for these metallicities and mixing cases. We find that models with Z = 0.02 have ignition masses and recurrence periods smaller than models with low Z, while the ejected mass and maximum luminosity are larger. Retention efficiency during novae outbursts decreases with increasing metallicity. In our implementation, inclusion of mixing at the H/He interface reduces accreted mass, ejected mass and recurrence period as compared to the no-mixing case, while the maximum luminosity becomes larger. Retention efficiency is significantly reduced, becoming negative in most of our models. For ease of use, we provide a tabular summary of our results.

Close double degenerate binaries are one of the favoured progenitor channels for type Ia supernovae, but it is unclear how many suitable systems there are in the Galaxy. We report results of a large radial velocity survey for double degenerate (DD) binaries using the UVES spectrograph at the ESO VLT (ESO SN Ia Progenitor surveY - SPY). Exposures taken at different epochs are checked for radial velocity shifts indicating close binary systems. We observed 689 targets classified as DA (displaying hydrogen-rich atmospheres), of which 46 turned out to possess a cool companion. We measured radial velocities (RV) of the remaining 643 DA white dwarfs. We managed to secure observations at two or more epochs for 625 targets, supplemented by eleven objects meeting our selection criteria from literature. The data reduction and analysis methods applied to the survey data are described in detail. The sample contains 39 double degenerate binaries, only four of which were previously known. 20 are double-lined systems, in which features from both components are visible, the other 19 are single-lined binaries. We provide absolute RVs transformed to the heliocentric system suitable for kinematic studies. Our sample is large enough to sub-divide by mass: 16 out of 44 low mass targets (<= 0.45 Msun) are detected as DDs, while just 23 of the remaining 567 with multiple spectra and mass >0.45 Msun are double. Although the detected fraction amongst the low mass objects (36.4 +/- 7.3%) is significantly higher than for the higher-mass, carbon/oxygen-core dominated part of the sample (3.9 +/- 0.8%), it is lower than the detection efficiency based upon companion star masses >= 0.05 Msun. This suggests either companion stars of mass < 0.05 Msun, or that some of the low mass white dwarfs are single.

We consider initial stage of the evolution of AM CVn type stars with white dwarf donors, which is accompanied by thermonuclear explosions in the layer of accreted He. It is shown that the accretion never results in detonation of He and accretors in AM CVn stars finish their evolution as massive WDs. We found, for the first time, that in the outbursts the synthesis of n-rich isotopes, initiated by the ${\mathrm{^{22}{Ne}(\alpha,n)^{25}Mg}}$ reaction becomes possible.

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.

We investigate the evolutionary properties of AM~CVn stars with massive white dwarf donors and accretors. As a representative of them we consider a binary initially composed by a 0.30 M$_\odot$, He WD and a 1.02 M$_\odot$, CO WD. We evaluate the time-dependent mass transfer rate from the donor and compute the evolution of the accretor, accounting for the effects of mass exchange on the evolution of orbital parameters. We model the thermal response of the accreting CO WD with the FUN evolutionary code coupled to a full nuclear network, from H to Bi, including more than 700 isotopes linked by about 1000 nuclear processes. We find that accretors in these systems evolve through the stages of steady He-burning and mild and strong He-flashes and become at the end CO WDs capped by a massive ($\sim 0.1$ M$_\odot$) He-rich buffer. During He-flashes (both mild and strong) the temperature in the He-shell increases above $3\times 10^8$ K, so that the ${^{22}Ne}(\alpha,n){^{25}Mg}$ reaction becomes efficient and $n$-rich isotopes can be produced. During the RLOF episodes triggered by strong non-dynamical He-flashes matter enriched in $\alpha$-elements and $n$-rich isotopes is ejected, polluting the interstellar medium. Our results strongly suggest that massive AM CVn systems with WD donors do not experience a final very strong dynamical He-flash driving an explosive event like SN .Ia. Though the ejected matter is highly enriched in heavy isotopes, the relative contribution of massive AM CVn systems to the Galactic chemical evolution is, most probably, negligible due to their expected paucity.

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.

Using a hybrid binary population synthesis approach, we modelled the formation and evolution of populations of accreting WDs for differing star formation histories. We found that the delay time distribution of SNe Ia in the single degenerate scenario is inconsistent with observations. Additionally, we found that our predicted X-ray and UV emission of populations of accreting WDs are consistent with the X-ray luminosities of early-type galaxies observed by Chandra and the HeII 4686\AA/H\beta line ratio measured in stacked SDSS spectra of passively evolving galaxies. Moreover, we found that the majority of current novae in elliptical-like galaxies have low-mass WDs, long decay times, long recurrence periods and are relatively faint. In contrast, the majority of current novae in spiral-like galaxies have massive WDs, short decay times, short recurrence periods and are relatively bright. Our predicted distribution of mass-loss timescales in an M31-like galaxy is consistent with observations for Andromeda.

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.

We present results of Monte Carlo simulation aiming at the estimate of the frequency of semi-detached Algol-type binaries among the stars observed as single ones. When account is made for various detection biases (mostly due to inclination of orbits), the fraction of Algols among Galactic disk stars appears to be 0.1--0.2\%. However, this number should be regarded as a lower limit only, since there are still unaccounted selection effects and other types of photometrically unresolved binaries. Hidden binarity appears to be an important phenomenon that should be taken into account when considering stellar statistics and construction of fundamental relations between stellar parameters.

Statistics of orbital parameters of binary stars as well as statistics of their physical characteristics bear traces of star formation history. However, statistical investigations of binaries are complicated by lacking or incomplete observational data and by a number of observational selection effects. Visual binaries are the most numerous observed binaries, the number of pairs exceeds 130000. The most complete list of presently known visual binary stars was compiled by cross-matching objects and combining data of the three largest catalogues of visual binaries. It was supplemented by the data on parallaxes, multicolor photometry, spectral characteristics of the stars of this list taken from other catalogues. This allowed us to compensate partly for the lack of observational data for these objects. Combined data allowed us to check validity of observational values and to investigate statistics of the orbital and physical parameters of visual binaries. Corrections for incompleteness of observational data are discussed. Obtained datasets and modern distributions of binary parameters will be used to reconstruct the initial distributions and parameters of the function of star formation for binary systems.

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.

Visual binary stars are the most abundant class of observed binaries. The most comprehensive list of data on visual binaries compiled recently by cross-matching the largest catalogues of visual binaries allowed a statistical investigation of observational parameters of these systems. The dataset was cleaned by correcting uncertainties and misclassifications, and supplemented with available parallax data. The refined dataset is free from technical biases and contains 3676 presumably physical visual pairs of luminosity class V with known angular separations, magnitudes of the components, spectral types, and parallaxes. We also compiled a restricted sample of 998 pairs free from observational biases due to the probability of binary discovery. Certain distributions of observational and physical parameters of stars of our dataset are discussed.

Theoretical modelling of the evolution of classical and recurrent novae plays an important role in studies of binary evolution, nucleosynthesis and accretion physics. However, from a theoretical perspective the observed statistical properties of novae remain poorly understood. In this paper, we have produced model populations of novae using a hybrid binary population synthesis approach for differing star formation histories (SFHs): a starburst case (elliptical-like galaxies), a constant star formation rate case (spiral-like galaxies) and a composite case (in line with the inferred SFH for M31). We found that the nova rate at 10\;Gyr in an elliptical-like galaxy is $\sim 10-20$ times smaller than a spiral-like galaxy with the same mass. The majority of novae in elliptical-like galaxies at the present epoch are characterized by low mass white dwarfs (WDs), long decay times, relatively faint absolute magnitudes and long recurrence periods. In contrast, the majority of novae in spiral-like galaxies at 10\;Gyr have massive WDs, short decay times, are relatively bright and have short recurrence periods. The mass loss time distribution for novae in our M31-like galaxy is in agreement with observational data for Andromeda. However, it is possible that we underestimate the number of bright novae in our model. This may arise in part due to the present uncertainties in the appropriate bolometric correction for novae.

Accreting white dwarfs (WDs) with non-degenerate companions are expected to emit in soft X-rays and the UV, if accreted H-rich material burns stably. They are an important component of the unresolved emission of elliptical galaxies, and their combined ionizing luminosity may significantly influence the optical line emission from warm ISM. In an earlier paper we modeled populations of accreting WDs, first generating WD with main-sequence, Hertzsprung gap and red giant companions with the population synthesis code \textsc{BSE}, and then following their evolution with a grid of evolutionary tracks computed with \textsc{MESA}. Now we use these results to estimate the soft X-ray (0.3-0.7keV), H- and He II-ionizing luminosities of nuclear burning WDs and the number of super-soft X-ray sources for galaxies with different star formation histories. For the starburst case, these quantities peak at $\sim 1$ Gyr and decline by $\sim 1-3$ orders of magnitude by the age of 10 Gyr. For stellar ages of $\sim$~10 Gyr, predictions of our model are consistent with soft X-ray luminosities observed by Chandra in nearby elliptical galaxies and He II 4686$\AA/\rm{H}{\beta}$ line ratio measured in stacked SDSS spectra of retired galaxies, the latter characterising the strength and hardness of the UV radiation field. However, the soft X-ray luminosity and He~II~4686$\AA/\rm{H}{\beta}$ ratio are significantly overpredicted for stellar ages of $\lesssim 4-8$ Gyr. We discuss various possibilities to resolve this discrepancy and tentatively conclude that it may be resolved by a modification of the typically used criteria of dynamically unstable mass loss for giant stars.

We study the physical and evolutionary properties of the "WD family" of AM CVn stars by computing realistic models of IDD systems. We evaluate self-consistently both the mass transfer rate from the donor, as determined by GW emission and interaction with the binary companion, and the thermal response of the accretor to mass deposition. We find that, after the onset of mass transfer, all the considered systems undergo a strong non-dynamical He-flash. However, due to the compactness of these systems, the expanding accretors fill their Roche lobe very soon, thus preventing the efficient heating of the external layers of the accreted CO WDs. Moreover, due to the loss of matter from the systems, the orbital separations enlarge and mass transfer comes to a halt. The further evolution depends on the value of $\dot{M}$, after the donors fill again their lobe. On one hand, if the accretion rate, as determined by the actual value of (M$_{don}$,M$_{acc}$), is high enough, the accretors experience several He-flashes of decreasing strength and then quiescent He-burning sets in. Later on, since the mass transfer rate in IDD is a permanently decreasing function of time, accretors experience several recurrent strong flashes. On the other hand, for intermediate and low values of $\dot{M}$, the accretors enter the strong flashes accretion regime. As expected, in all the considered systems the last He-flash is the strongest one, even if a dynamical event never occurs. When the mass accretion rate decreases below (2-3)$\times 10^{-8} M_\odot yr^{-1}$, the compressional heating of the He-shell becomes less efficient than the neutrino cooling, so that all the accretors in the considered systems evolve into massive degenerate objects. Our results suggest that SNe .Ia or type Ia Supernovae due to Edge-Lit Detonation in the WD family of AM CVn stars should be much more rare than previously expected.

The behaviour of carbon-oxygen white dwarfs (WDs) subject to direct helium accretion is extensively studied. We aim to analyze the thermal response of the accreting WD to mass deposition at different time scales. The analysis has been performed for initial WDs masses and accretion rates in the range (0.60 - 1.02) Msun and 1.e-9 - 1.e-5 Msun/yr, respectively. Thermal regimes in the parameters space M_{WD} - dot{M}_{He}, leading to formation of red-giant-like structure, steady burning of He, mild, strong and dynamical flashes have been identified and the transition between those regimes has been studied in detail. In particular, the physical properties of WDs experiencing the He-flash accretion regime have been investigated in order to determine the mass retention efficiency as a function of the accretor total mass and accretion rate. We also discuss to what extent the building-up of a He-rich layer via H-burning could be described according to the behaviour of models accreting He-rich matter directly. Polynomial fits to the obtained results are provided for use in binary population synthesis computations. Several applications for close binary systems with He-rich donors and CO WD accretors are considered and the relevance of the results for the interpretation of He-novae is discussed.

Accreting, nuclear-burning white dwarfs have been deemed to be candidate progenitors of type Ia supernovae, and to account for supersoft X-ray sources, novae, etc. depending on their accretion rates. We have carried out a binary population synthesis study of their populations using two algorithms. In the first, we use the binary population synthesis code \textsf{BSE} as a baseline for the "rapid" approach commonly used in such studies. In the second, we employ a "hybrid" approach, in which we use \textsf{BSE} to generate a population of white dwarfs (WD) with non-degenerate companions on the verge of filling their Roche lobes. We then follow their mass transfer phase using the detailed stellar evolution code \textsf{MESA}. We investigate the evolution of the number of rapidly accreting white dwarfs (RAWDs) and stably nuclear-burning white dwarfs (SNBWDs), and estimate the type Ia supernovae (SNe Ia) rate produced by "single-degenerate" systems (SD). We find significant differences between the two algorithms in the predicted numbers of SNBWDs at early times, and also in the delay time distribution (DTD) of SD SNe Ia. Such differences in the treatment of mass transfer may partially account for differences in the SNe Ia rate and DTD found by different groups. Adopting 100\% efficiency for helium burning, the rate of SNe Ia produced by the SD-channel in a Milky-way-like galaxy in our calculations is $2.0\times10^{-4}\rm{yr}^{-1}$, more than an order of magnitude below the observationally inferred value. In agreement with previous studies, our calculated SD DTD is inconsistent with observations.

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Mergings of compact binary stars are expected to be the most important sources for the forthcoming gravitational-wave (GW) astronomy. In the first part of the review, we discuss observational manifestations of close binary stars with NS and/or black components and their merger rate, crucial points in the formation and evolution of compact stars in binary systems, including the treatment of the natal kicks which NSs and BHs acquire during the core collapse of massive stars and the common envelope phase of binary evolution, which are most relevant to the merging rates of NS-NS, NS-BH and BH-BH binaries. The second part of the review is devoted mainly to formation and evolution of binary WDs and their observational manifestations, including their role as progenitors of cosmologically important thermonuclear SN Ia. We also consider AM CVn-stars which are thought to be the best verification binary GW sources for future low-frequency GW space interferometers.

We conducted a multi-wavelength study to unveil the properties of the extremely long-period cataclysmic variable V479 And. We performed series of observations, including moderate to high spectral resolution optical spectrophotometry, X-ray observations with Swift, linear polarimetry and near-IR photometry. This binary system is a low-inclination ~ 17^o system with a 0.594093(4) day orbital period. The absorption line complex in the spectra indicate a G8--K0 spectral type for the donor star, which has departed from the zero-age main sequence. This implies a distance to the object of about 4 kpc. The primary is probably a massive 1.1-1.4 Msun magnetic white dwarf, accreting matter at a rate M(dot) > 10^-10 Msun/ yr. This rate can be achieved if the donor star fills its corresponding Roche lobe, but there is little observational evidence for a mass-transfer stream in this system. An alternative explanation is a stellar wind from the donor star, although such a high rate mass loss is not anticipated from a subgiant. If the strongly magnetic white dwarf in V479 And. is confirmed by future observations, the system the polar with the longest observed orbital period. We also discuss the evolutionary state of V479 And.

[abridged] Aims. We model the number and properties of ultracompact X-ray binaries (UCXBs) in the Galactic Bulge. The objective is to learn about the formation of UCXBs and their evolution, such as the onset of mass transfer and late-time behavior. Methods. The binary population synthesis code SeBa and stellar evolutionary tracks are used to model the UCXB population. The luminosity behavior of UCXBs is predicted using long-term X-ray observations of the known UCXBs and the thermal-viscous disk instability model. Results. In our model, the majority of UCXBs initially have a helium burning star donor. In the absence of a mechanism that destroys old UCXBs, we predict (0.2 - 1.9) x 10^5 UCXBs in the Galactic Bulge, mostly at orbital periods longer than 60 min (a large number of long-period systems also follows from the observed short-period UCXB population). About 5 - 50 UCXBs should be brighter than 10^35 erg/s, mostly persistent sources with orbital periods shorter than 30 min and with degenerate helium and carbon-oxygen donors. This is one order of magnitude more than the observed number of (probably) three. Conclusions. This overprediction of short-period UCXBs implies that fewer systems are formed, or that a super-Eddington mass transfer rate is more difficult to survive than assumed. The very small number of observed long-period UCXBs with respect to short-period UCXBs, the surprisingly high luminosity of the observed UCXBs with orbital periods around 50 min, and the properties of the PSR J1719-1438 system all point to much faster UCXB evolution than expected from angular momentum loss via gravitational wave radiation alone. Old UCXBs, if they still exist, probably have orbital periods longer than 2 h and have become very faint due to either reduced accretion or quiescence, or have become detached. UCXBs are promising candidate progenitors of isolated millisecond radio pulsars.

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.

The properties of wind accretion in symbiotic X-ray binaries (SyXBs) consisting of red-giant and magnetized neutron star (NS) are discussed. The spin-up/spin-down torques applied to NS are derived based on a hydrodynamic theory of quasi-spherical accretion onto magnetized NSs. In this model, a settling subsonic accretion proceeds through a hot shell formed around the NS magnetosphere. The accretion rate onto the NS is determined by the ability of the plasma to enter the magnetosphere.Due to large Reynolds numbers in the shell, the interaction of the rotating magnetosphere with plasma initiates a subsonic turbulence. The convective motions are capable of carrying the angular momentum through the shell. We carry out a population synthesis of SyXBs in the Galaxy with account for the spin evolution of magnetized NS. The Galactic number of SyXBs with bright (M_v<1) low-mass red-giant companion is found to be from \sim 40 to 120, and their birthrate is \sim 5\times 10^{-5}-10^{-4} per year. According to our model, among known SyXBs, Sct X-1 and IRXS J180431.1-273932 are wind-fed accretors. GX 1+4 lies in the transition from the wind-fed SyXBs to SyXBs in which the giants overflow their Roche lobe. The model successfully reproduces very long NS spins (such as in IGR J16358-4724 and 4U 1954+31) without the need to invoke very strong magnetic fields.

We investigate the thermal response of CO WDs accreting He-rich matter directly from their companions in binary systems. Our results suggest that the He-channel cannot provide progenitors for majority of "normal" SNe Ia.

The rates of SN Ia for double-degenerate and single-degenerate scenario are computed for the models of spiral and elliptical galaxies. The number of nuclearly burning white dwarfs (NBWD) is traced. The data favours double-degenerate scenario and suggests lower number of NBWD per unit mass in ellipticals and their lower average mass as one of the reasons for the difference in the number of supersoft X-ray sources observed in the galaxies of different types.

We analyze temporal evolution of the number of accreting white dwarfs with shell hydrogen burning in semidetached and detached binaries. We consider a stellar system in which star formation lasts for 10 Gyr with a constant rate, as well as a system in which the same amount of stars is formed in a single burst lasting for 1 Gyr. Evolution of the number of white dwarfs is confronted to the evolution of occurrence rate of events that usually are identified with SN Ia or accretion-induced collapses, i.e. with accumulation of Chandrasekhar mass by a white dwarf or a merger of a pair of CO white dwarfs with total mass not lower than the Chandrasekhar one. In the systems with a burst of star formation, at $t=$10 Gyr observed supersoft X-ray sources, most probably, are not precursors of SN Ia. The same is true for an overwhelming majority of the sources in the systems with constant star formation rate. In the systems of both kinds mergers of white dwarfs is the dominant SN Ia scenario. In symbiotic binaries, accreting CO-dwarfs do not accumulate enough mass for SN Ia explosion, while ONeMg-dwarfs finish their evolution by an accretion-induced collapse with formation of a neutron star.

We present a detailed investigation of SBS1150+599A, a close binary star hosted by the planetary nebula PN G135.9+55.9 (TS01, Stasinska et al, 2009). The nebula, located in the Galactic halo, is the most oxygen-poor one known to date and is the only one known to harbor a double degenerate core. We present XMM-Newton observations of this object, which allowed the detection of the previously invisible component of the binary core, whose existence was inferred so far only from radial velocity and photometric variations. The parameters of the binary system were deduced from a wealth of information via three independent routes using the spectral energy distribution (from the infrared to X-rays), the light and radial velocity curves, and a detailed model atmosphere fitting of the stellar absorption features of the optical/UV component. We find that the cool component must have a mass of 0.54+/-0.2 Msun, an average effective temperature, Teff, of 58000+/-3000 K, a mean radius of 0.43+/-0.3 Rsun, a gravity log g=5.0+/-0.3, and that it nearly fills its Roche lobe. Its surface elemental abundances are found to be: 12 + log He/H = 10.95+/-0.04 dex, 12 + log C/H = 7.20+/-0.3 dex, 12 + log N/H < 6.92 and 12 + log O/H < 6.80, in overall agreement with the chemical composition of the planetary nebula. The hot component has Teff = 160-180 kK, a luminosity of about ~10e4 Lsun and a radius slightly larger than that of a white dwarf. It is probably bloated and heated as a result of intense accretion and nuclear burning on its surface in the past. The total mass of the binary system is very close to Chandrasekhar limit. This makes TS01 one of the best type Ia supernova progenitor candidates. We propose two possible scenarios for the evolution of the system up to its present stage.

The planetary nebula TS 01 (also called PN G 135.9+55.9 or SBS 1150+599A), with its record-holding low oxygen abundance and its double degenerate close binary core (period 3.9 h), is an exceptional object located in the Galactic halo. We have secured observational data in a complete wavelength range in order to pin down the abundances of half a dozen elements in the nebula. The abundances are obtained via detailed photoionization modelling taking into account all the observational constraints (including geometry and aperture effects) using the pseudo-3D photoionization code Cloudy_3D. The spectral energy distribution of the ionizing radiation is taken from appropriate model atmospheres. Both stellar components contribute to the ionization: the ``cool'' one provides the bulk of hydrogen ionization, and the ``hot'' one is responsible for the presence of the most highly charged ions, which explains why previous attempts to model the nebula experienced difficulties. The nebular abundances of C, N, O, and Ne are found to be respectively, 1/3.5, 1/4.2, 1/70, and 1/11 of the Solar value, with uncertainties of a factor 2. Thus the extreme O deficiency of this object is confirmed. The abundances of S and Ar are less than 1/30 of Solar. Standard models of stellar evolution and nucleosynthesis cannot explain the abundance pattern observed in the nebula. To obtain an extreme oxygen deficiency in a star whose progenitor has an initial mass of about 1 msun requires an additional mixing process, which can be induced by stellar rotation and/or by the presence of the close companion. We have computed a stellar model with initial mass of 1 msun, appropriate metallicity, and initial rotation of 100 kms, and find that rotation greatly improves the agreement between the predicted and observed abundances.

We study the formation of ultra-compact binaries (AM CVn stars and ultra-compact X-ray binaries) with emphasis on the surface chemical abundances of the donors in these systems. Hydrogen is not convincingly detected in the spectra of these systems. Three different proposed formation scenarios involve different donor stars, white dwarfs, helium stars or evolved main-sequence stars. Using detailed evolutionary calculations we show that the abundances of helium WD donors and evolved main-sequence stars are close to equilibrium CNO-processed material, and the detailed abundances correlate with the core temperature and thus mass of the MS progenitors. Evolved MS donors typically have traces of H left. For hybrid or CO white dwarf donors, the carbon and oxygen abundances depend on the temperature of the helium burning and thus on the helium core mass of the progenitors. For helium star donors in addition to their mass, the abundances depend strongly on the amount of helium burnt before mass transfer starts and can range from unprocessed and thus almost equal to CNO-processed matter, to strongly processed and thus C/O rich and N-deficient. We briefly discuss the relative frequency of these cases for helium star donors, based on population synthesis results. Finally we give diagnostics for applying our results to observed systems and find that the most important test is the N/C ratio, which can indicate the formation scenario as well as, in some cases, the mass of the progenitor of the donor. In addition, if observed, the N/O, O/He and O/C ratios can distinguish between helium star and WD donors. Applied to the known systems we find evidence for WD donors in the AM CVn systems GP Com, CE 315 and SDSS J0804+16 and evidence for hybrid WD or very evolved helium star donors in the UCXBs 4U 1626-67 and 4U 0614+09. [Abridged]

This white paper briefly describes the astrophysics of ultra-compact binaries, with emphasis of the challenges and opportunities in the next decade.

We consider evolutionary models for the population of short-period (<10 hr) low-mass black-hole binaries (LMBHBs) and compare them with observations of soft X-ray transients (SXTs). We show that assuming strongly reduced magnetic braking (as suggested by us before for low-mass semidetached binaries) the calculated masses and effective temperatures of secondaries are encouragingly close to the observed masses and effective temperatures (as inferred from their spectra) of donor stars in short-period LMBHBs. Theoretical mass-transfer rates in SXTs are consistent with the observed ones if one assumes that accretion discs in these systems are truncated (``leaky''). We find that the population of short-period SXTs is formed mainly by systems which had unevolved or slightly evolved main-sequence donors ($M_2 \lesssim 1.2M_\odot$) with a hydrogen abundance in the center $X_c > 0.35$ at the Roche-lobe overflow (RLOF). Longer period (0.5 - 1) day) SXTs might descend from systems with initial donor masses of about 1 $M_\odot$ and $X_c < 0.35$. Thus, one can explain the origin of short period LMBHB without invoking donors with cores almost totally depleted of hydrogen. Our models suggest that, unless the currently accepted empirical estimates of mass-loss rates by winds for massive O-stars and Wolf-Rayet stars are significantly over-evaluated, a very high efficiency of common-envelope ejection is necessary to form short-period LMBHBs.

We present results of a systematic investigation of the evolution of low-mass (0.35, 0.40, and 0.65 solar mass) helium donors in semidetached binaries with accretors - white dwarfs. In the initial models of evolutionary sequences abundance of helium in the center is between $\simeq 0.1$ and 0.98. Results of computations may be applied to the study of the origin and evolutionary state of AM CVn stars. It is shown that the minimum orbital periods of the systems only weakly depend on the total mass of the system and evolutionary state of the donor at RLOF and are equal to 9-11 min. The scatter in the mass-exchange rates at a given orbital period in the range between period minimum and about 40 min. does not exceed $\sim 2.5$. At orbital periods exceeding about 20 min. mass-losing stars are weakly degenerate homogeneous cooling objects and abundances of He, C, N, O, Ne in the matter lost by them depend on the extent of He-depletion at RLOF. For the systems which are currently considered as the most probable model candidates for AM CVn stars with helium donors these abundances are, approximately, Y$>$0.4, X_C$<$0.3, X_O$<$0.25, X_N$<$0.005. At orbital periods greater than about 40 min. the timescale of mass-loss begins to exceed thermal time-scale of the donors, the latter begin to contract, they become more degenerate and, apparently, "white-dwarf" and "helium-star" populations of AM CVn stars merge.

We discuss the formation, evolution and observational parameters of the population of short-period ($<10$ hr) low-mass black-hole binaries (LMBHB). Their evolution is determined by the orbital angular momentum loss and/or nuclear evolution of the donors. All observed semidetached LMBHB are observed as soft X-ray transients (SXTs). The absence of observed short-period stable luminous X-ray sources with black holes and low-mass optical components suggests that upon RLOF by the donor, the angular-momentum losses are substantially reduced. The model with reduced angular-momentum loss reasonably well reproduces the masses and effective temperatures of the observed secondaries of SXTs. Theoretical mass-transfer rates in SXTs are consistent with those deduced from observations only if the accretion discs in LMBHB are truncated. The population of short-period LMBHB is formed mainly by systems which at RLOF had unevolved or slightly evolved donors (abundance of hydrogen in the center $> 0.35$). Our models suggest that a very high efficiency of common envelopes ejection is necessary to form LMBHB.

We consider evolutionary models for the population of short-period (<10 hr) low-mass black-hole binaries (LMBHB) and compare them with observations of soft X-ray transients (SXT). Evolution of LMBHB is determined by nuclear evolution of the donors and/or orbital angular momentum loss due to magnetic braking by the stellar wind of the donors and gravitational wave radiation. We show that the absence of observed stable luminous LMBHB implies that upon RLOF by the low-mass donor angular momentum losses are substantially reduced with respect to the Verbunt and Zwaan "standard" prescription for magnetic braking. Under this assumption masses and effective temperatures of the model secondaries of LMBHB are in a satisfactory agreement with the masses and effective temperatures (as inferred from their spectra) of the observed donors in LMBHB. Theoretical mass-transfer rates in SXTs are consistent with the observed ones if one assumes that accretion discs in these systems are truncated ("leaky"). We find that the population of short-period SXT is formed mainly by systems which had unevolved or slightly evolved (X_c > 0.35) donors at the RLOF. Longer period (0.5 - 1 day) SXT might descend from systems with initial donor mass about 1 solar and X_c < 0.35. It is unnecessary to invoke donors with almost hydrogen-depleted cores to explain the origin of LMBHB. Our models suggest that a very high efficiency of common-envelopes ejection is necessary to form LMBHB, unless currently commonly accepted empirical estimates of mass-loss rates by winds for pre-WR and WR-stars are significantly over-evaluated.

Aims. Determination of the present evolutionary state and the restrictions on the initial mass ratios of RZ Cas, KO Aql and S Equ. Methods. Comparison of mass gaining stars with evolutionary models of single stars with the same mass and subsequent comparison with accretion tracks from simultaneous conservative binary evolution. Results. The gainers are in an early main sequence stage (Xc greater than 0.5), with KO Aql being almost unevolved (assuming quasi-thermal equilibrium). The initial donor/gainer mass ratios Mdi/Mgi must be larger than three to obtain the present mass and luminosity of the gainers.

We study the evolution and fate of solar composition supermassive stars in the mass range 60 - 1000 Msun. Our study is relevant both for very massive objects observed in young stellar complexes as well as supermassive stars that may potentially form through runaway stellar collisions. We predict the outcomes of stellar evolution employing a mass-loss prescription that is consistent with the observed Hertzsprung-Russell Diagram location of the most massive stars. We compute a series of stellar models with an appropriately modified version of the Eggleton evolutionary code. We find that super-massive stars with initial masses up to 1000 Msun end their lives as objects less massive than ~150 Msun. These objects are expected to collapse into black holes (with M < 70 Msun) or explode as pair-instability supernovae. We argue that if ultraluminous X-ray sources (ULXs) contain intermediate-mass black holes, these are unlikely to be the result of runaway stellar collisions in the cores of young clusters.

We follow the evolution of helium stars of initial mass $(2.2 - 2.5) M_\odot$, and show that they undergo off-center carbon burning, which leaves behind ${\mathbf \sim 0.01 M_\odot}$ of unburnt carbon in the inner part of the core. When the carbon-oxygen core grows to Chandrasekhar mass, the amount of left-over carbon is sufficient to ignite thermonuclear runaway. At the moment of explosion, the star will possess an envelope of several $0.1 M_{\odot}$, consisting of He, C, and possibly some H, perhaps producing a kind of peculiar SN. Based on the results of Waldman and Barkat (2007) for accreting white dwarfs, we expect to get thermonuclear runaway at a broad range of $\rho_c \approx (1 - 6) \times 10^9 \mathrm{g cm^{-3}}$, depending on the amount of residual carbon. We verified the feasibility of this scenario by showing that in a close binary system with initial masses $(8.5 + 7.7) M_{\odot}$ and initial period of 150 day the primary produces a helium remnant of $2.3 M_{\odot}$ that evolves further like the model we considered.

We present X-ray observations of the close binary nucleus of the planetary nebula PN G135.9+55.9 obtained with the XMM satellite. The nebula is the most oxygen-poor PN known to date and is located in the Galactic halo. It is known to harbor a close binary nucleus of which only one component can be observed in optical-UV range. New X-ray observations show that the invisible component is a very hot compact star. This finding allows us to reconstruct the immediate past of the object and predict its future. The parameters of the binary components we determine strongly suggest that the precursor was a symbiotic supersoft X-ray source that finished its life by Roche lobe overflow. PN G135.9+55.9 is an excelent candidate for a future type Ia supernova.

We present X-ray observations of the close binary nucleus of the planetary nebula SBS 1150+599A obtained with the XMM-Newton satellite. Only one component of the binary can be observed in optical-UV. New X-ray observations show that the previously invisible component is a very hot compact star. This finding allows us to deduce rough values for the basic parameters of the binary. With a high probability the total mass of the system exceeds Chandrasekhar limit and makes the SBS 1150+599A one of the best candidate for a supernova type Ia progenitor.

Context. NGC 300 X-1 is the second extragalactic candidate, after IC 10 X-1, in the rare class of Wolf-Rayet/compact object X-ray binary systems exemplified in the Galaxy by Cyg X-3. From a theoretical point of view, accretion onto a black hole in a detached system is possible for large orbital periods only if the mass of the relativistic object is high or the velocity of the accreted wind is low. Aims. We analysed a 2 week SWIFT XRT light curve of NGC 300 X-1 and searched for periodicities. Methods. Period searches were made using Lomb-Scargle periodogram analysis. We evaluated the confidence level using Monte Carlo simulations. Results. A period of 32.8+-0.4h (3 sigma error) was found for NGC 300 X-1 with a confidence level >99%. Furthermore, we confirm the high irregular variability during the high flux level, as already observed in the XMM-Newton observations of the source. A folded XMM-Newton light curve is shown, with a profile that is in agreement with SWIFT. The mean absorbed X-ray luminosity in the SWIFT observations was 1.5x10^38 erg/s, close to the value derived from the XMM-Newton data. Conclusions. While Cyg X-3 has a short period of 4.8 h, the period of NGC 300 X-1 is very close to that of IC 10 X-1 (34.8+-0.9 h). These are likely orbital periods. Possibility of formation of accretion disk for such high orbital periods strongly depends on the terminal velocity of the Wolf-Rayet star wind and black-hole mass. While low masses are possible for wind velocities < 1000 km/s, these increase to several tens of solar masses for velocities > 1600 km/s and no accretion disk may form for terminal velocities larger than 1900 km/s.

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Binary NSs and BHs are thought to be the primary astrophysical sources of gravitational waves (GWs) within the frequency band of ground-based detectors, while compact binaries of WDs are important sources of GWs at lower frequencies to be covered by space interferometers (LISA). Major uncertainties in the current understanding of properties of NSs and BHs most relevant to the GW studies are discussed, including the treatment of the natal kicks which compact stellar remnants acquire during the core collapse of massive stars and the common envelope phase of binary evolution. We discuss the coalescence rates of binary NSs and BHs and prospects for their detections, the formation and evolution of binary WDs and their observational manifestations. Special attention is given to AM CVn-stars -- compact binaries in which the Roche lobe is filled by another WD or a low-mass partially degenerate helium-star, as these stars are thought to be the best LISA verification binary GW sources.

We show that even most massive initial solar composition stars hardly form black holes with mass exceeding about 100 solar mass.

We have carried out a detailed study of symbiotic stars with white dwarf accretors by means of a population synthesis code. We estimate the total number of symbiotic stars with white dwarf accretors in the Galaxy as 1,200 - 15,000. This range is compatible with observational estimates. Two crucial physical parameters that define the birthrate and number of symbiotic stars are the efficiency of accretion by white dwarfs (which greatly depends on the separation of components after common envelope stage and stellar wind velocity) and the mass of the hydrogen layer which the white dwarf can accumulate prior to the hydrogen ignition. The theoretical estimate of the Galactic occurrence rate of symbiotic novae ranges from about 1.3 to about 13.5 per year out of which weak symbiotic novae comprise about 0.5 to 6.0 per year, depending on the model assumptions. We simulate the distributions of symbiotic stars over orbital periods, masses of components, mass-loss rates of cool components, mass-accretion rates of hot components and luminosity of components. Agreement with observations is reasonable.

We present results of a population synthesis study for semidetached short orbital period binaries which contain low-mass(<1.5 Msun) donors and black hole (>4 Msun) accretors. Evolution of these binaries is determined by nuclear evolution of the donors and/or orbital angular momentum loss due to magnetic braking by the stellar wind of the donors and gravitational wave radiation. According to our model, the estimated total number of this type of black-hole binaries in the Galaxy is about 10000. If the magnetic braking is described by the Verbunt & Zwaan formula, the model predicts around 3000 transient systems with periods >2 hours and around 300 luminous stable systems with periods between 3 and 8 hours. Several dozens of these bright systems should be above the RXTE ASM sensitivity limit. The absence of such systems implies that angular momentum losses are reduced by a factor more than 2 with respect to the Verbunt & Zwaan prescription. We show that it is unlikely that the transient behaviour of black-hole short-period X-ray binaries is explained by the evolved nature of the stellar companion. A substantial fraction of black-hole binaries with periods >3 hours could be faint with truncated, stable cold accretion discs as proposed by Menou et al. Most of the semidetached black-hole binaries are expected to have periods shorter than ~2 hours. Properties of such, still to be observed, very small mass-ratio (q<0.02) binaries are different from those of their longer period cousins.

Cygnus X-3 is a strong X-ray source (L_X about 10^38 erg/s) which is thought to consist of a compact object, accreting matter from a helium star. We find analytically that the estimated ranges of mass-loss rate and orbital-period derivative for Cyg X-3 are consistent with two models: i) the system is detached and the mass loss from the system comes from the stellar wind of a massive helium star, of which only a fraction that allows for the observed X-ray luminosity is accreted, or ii) the system is semidetached and a Roche-lobe-overflowing low- or moderate-mass helium donor transfers mass to the compact object, followed by ejection of its excess over the Eddington rate from the system. These analytical results appear to be consistent with evolutionary calculations. By means of population synthesis we find that currently in the Galaxy there may exist ~1 X-ray binary with a black hole that accretes from a >~ 7 MSun Wolf-Rayet star and ~1 X-ray binary in which a neutron star accretes matter from a Roche-lobe-overflowing helium star with mass <~ 1.5 MSun. Cyg X-3 is probably one of these systems.

By means of population synthesis we model the Galactic ensemble of helium stars. It is assumed that all helium stars are formed in binaries. Under this assumption, single helium stars are produced by the mergers of helium remnants of components of close binaries (mainly, by merging helium white dwarfs) and by disruption of binaries with helium components in supernovae explosions. The estimate of the total birthrate of helium stars in the Galaxy is 0.043 yr$^{-1}$, their total number is estimated as $4 \times 10^6$. The rate of binarity in the total sample is 76%. We construct a subsample of low-mass ($M_{\rm He} \lesssim 2 \ms$) helium stars limited by observational selection effects: stellar magnitude ($V_{\rm He} \leq 16$), ratio of stellar magnitudes of components in binaries ($V_{\rm He}\leq V_{\rm comp}$), lower limit of the semiamplitude of radial velocity that is necessary for discovery of binarity ($K_{min}=30$ km/s). The parameters of this ``observable'' sample are in satisfactory agreement with the parameters of the observed ensemble of sdB stars. In particular, in the selection-limited sample binarity rate is 58%. We analyze the relations between orbital periods and masses of helium stars and their companions in systems with different combinations of components. We expect that overwhelming majority ($\sim 90$%) of unobserved components in binary sdB stars are white dwarfs, predominantly, carbon-oxygen ones.

A review of the basic principles of the population synthesis for binary stars is presented. We discuss the break-up of low and intermediate mass close binaries over different evolutionary scenarios and, as an example, briefly consider results of the population synthesis for SN Ia.

Thanks to the rapid increase of observations of Supernovae Ia, we may now claim that the Universe is accelerating. SN Ia are believed to be good standard candles, and after correcting the observations by various methods a precise Hubble diagram results. One of the most important problems to solve in astrophysics today is to find the progenitors of SN Ia. Candidates are to be found in close binary systems where one of the components may accumulate Chandrasekhar mass and then explode. We show that the AM CVn systems may contribute to the SNe population, but will not be the dominant contributor.

We discuss application of population synthesis for binary stars to progenitors of SN Ia. We show that the only candidate systems able to support the rate of SNe Ia $\sim 10^{-3} yr^{-1}$ both in old and young populations are merging white dwarfs. In young populations ($\sim 1$ Gyr) edge-lit detonations in semidetached systems with nondegenerate helium star donors are also able to support a similar rate. The estimated current Galactic rate of SN Ia with single-degenerate progenitors is $\sim 10^{-4} yr^{-1}$.

We report on systematic radial velocity surveys for white dwarf - white dwarf binaries (double degenerates - DDs) including SPY (ESO Supernovae Ia progenitor survey) recently carried out at the VLT. A large sample of DD will allow us to put strong constrains on the phases of close binary evolution of the progenitor systems and to perform an observational test of the DD scenario for supernovae of type Ia. We explain how parameters of the binaries can be derived from various methods. Results for a sample of DDs are presented and discussed.

We discuss the model for the Galactic sample of the AM CVn systems with orbital periods P less than 1500 s, that can be detected in the optical and/or X-ray bands and may be resolved by the gravitational waves detector LISA. At P=3 to 10 min all detectable systems are X-ray sources. At P greater than 10 min most systems are only detectable in the V-band. About 30% of the X-ray sources is also detectable in the V-band. About 10,000 AM CVn systems might be resolved by LISA; this is comparable to the number of detached double white dwarfs that can be resolved. Several hundreds of AM CVn LISA sources might be also detectable in the V- and/or X-ray bands.

We model the population of AM CVn systems in the Galaxy and discuss the detectability of these systems with optical, X-ray and gravitational wave detectors. We concentrate on the short period (P < 1500 s) systems, some of which are expected to be in a phase of direct impact accretion. Using a self-consistent model for the star formation history and radial distribution of stars in the Galaxy plus simple models for the emission of optical and X-ray radiation from the AM CVn systems and interstellar absorption, we derive the sample of short-period AM CVn systems that can be detected in the optical and/or X-ray bands. At the shortest periods the detectable systems are all X-ray sources, some with periods as short as three minutes. At periods above 10 minutes most detectable systems are optical sources. About one third of the X-ray sources are also detectable in the optical band. We also calculate the gravitational wave signal of the short-period AM CVn systems. We find that potentially several thousand AM CVn systems can be resolved by the gravitational wave detector LISA, comparable to the expected number of detached double white dwarfs that can be resolved. We estimate that several hundreds of the AM CVn systems resolvable by LISA, are also detectable in the optical and/or X-ray bands.

We have considered scenarios for the evolution of close binaries resulting in the formation of semi-detached systems in which a white dwarf can accumulate Chandrasekhar mass by accretion from a main-sequence or subgiant companion with $M\sim 2 M_\odot$. These white dwarfs, probably, explode as type Ia supernovae or collapse with formation of neutron stars. We have carried out a population synthesis study for these systems and have estimated the occurrence rate of such events in the Galaxy, depending on the parameter of common envelopes, mass-exchange rate in the binary, reaction of the main-sequence component on accretion of helium in the intervening phase of evolution. We have found that the model occurrence rate of SNe Ia in semi-detached systems is $\simeq 0.2 \times 10^{-3}$ yr$^{-1}$, i. e., it does not exceed $\sim 10%$ of the observational estimate of the Galactic occurrence rate of SNe Ia.

We propose a simple model for the bolometric light curve of type Ia supernova exploding in a dense circumstellar (CS) envelope to describe the light curves of supernovae 2002ic and 1997cy. The modeling shows that at the radius $\sim7\times10^{15}$ cm the density of CS envelopes around both supernovae is similar. The mass of the CS envelope around SN 1997cy is close to $5 M_{\odot}$, while the characteristic time of the ejection of this envelope does not exceed 600 yr. We analyze two possible evolutionary scenarios which might lead to the explosion of type Ia supernova inside a dense CS envelope: accretion on CO white dwarf in the symbiotic binary and evolution of a single star with the initial mass of about $8 M_{\odot}$. If the conjecture about the explosion of type Ia supernova in a dense CS envelope is correct in the case of SN 2002ic and SN 1997cy then the rapid loss of the red supergiant envelope and the subsequent explosion of the CO white dwarf are synchronized by certain mechanism. This mechanism might be related to the contraction of the white dwarf as it approaches the Chandrasekhar limit. We show that formation of a (super)Chandrasekhar mass due to the merger of a CO white dwarf and CO core of a red supergiant with consequent explosion is unlikely, since it does not provide the required synchronization of the rapid mass loss and explosion.

We report on a large survey for double degenerate (DD) binaries as potential progenitors of type Ia supernovae with the UVES spectrograph at the ESO VLT (SN Ia Progenitor surveY - SPY). About 560 white dwarfs were checked for radial velocity variations until now. Ninety new DDs have been discovered, including short period systems with masses close to the Chandrasekhar mass.

We report on a large survey for double degenerate (DD) binaries as potential progenitors of type Ia supernovae with the UVES spectrograph at the ESO VLT (ESO SN Ia Progenitor surveY - SPY).

New diffraction-limited speckle images of the Red Rectangle in the wavelength range 2.1--3.3 microns with angular resolutions of 44--68 mas and previous speckle images at 0.7--2.2 microns revealed well-resolved bright bipolar outflow lobes and long X-shaped spikes originating deep inside the outflow cavities. This set of high-resolution images stimulated us to reanalyze all infrared observations of the Red Rectangle using our two-dimensional radiative transfer code. The new detailed modeling, together with estimates of the interstellar extinction in the direction of the Red Rectangle enabled us to more accurately determine one of the key parameters, the distance D=710 pc with model uncertainties of 70 pc, which is twice as far as the commonly used estimate of 330 pc. The central binary is surrounded by a compact, massive (M=1.2 Msun), very dense dusty torus with hydrogen densities reaching n_H=2.5x10^12 cm^-3 (dust-to-gas mass ratio rho_d/rho~0.01). The bright component of the spectroscopic binary HD 44179 is a post-AGB star with mass M*=0.57 Msun, luminosity L*=6000 Lsun, and effective temperature T*=7750 K. Based on the orbital elements of the binary, we identify its invisible component with a helium white dwarf with Mwd~0.35 Msun, Lwd~100 Lsun, and Twd~6x10^4 K. The hot white dwarf ionizes the low-density bipolar outflow cavities inside the dense torus, producing a small HII region observed at radio wavelengths. We propose an evolutionary scenario for the formation of the Red Rectangle nebula, in which the binary initially had 2.3 and 1.9 Msun components at a separation of 130 Rsun. The nebula was formed in the ejection of a common envelope after Roche lobe overflow by the present post-AGB star.

We study the formation of neon-enriched donor stars in ultracompact X-ray binaries (orbital periods P<80 min) and show that their progenitors have to be low-mass (0.3 - 0.4 solar mass) ``hybrid'' white dwarfs (with CO cores and thick helium mantles). Stable mass transfer is possible if in the initial stages of mass exchange mass is lost from the system, taking away the specific orbital angular momentum of the accretor (``isotropic re-emission''). The excess of neon in the transferred matter is due to chemical fractionation of the white dwarf which has to occur prior to the Roche lobe overflow by the donor. The estimated lower limit of the orbital periods of the systems with neon-enriched donors is close to 10 min. We show that the X-ray pulsar 4U 1626-67, which likely also has a neon-enriched companion, may have been formed via accretion induced collapse of an oxygen-neon white dwarf accretor if the donor was a hybrid white dwarf.

In the course of our search for double degenerate (DD) binaries as potential progenitors of type Ia supernovae with the UVES spectrograph at the ESO VLT (ESO SN Ia Progenitor surveY - SPY) we discovered that the white dwarf HE 1414-0848 is a double-lined DA+DA binary with an orbital period of P = 12h 25m 39s. Semi-amplitudes of 128km/s and 100km/s are derived for the individual components. The amplitude ratio and the measured difference in gravitational redshift is used to estimate the masses of the individual components: 0.55Msol and 0.71Msol. Hence the total mass of the HE 1414-0848 system is 1.26Msol, only 10% below the Chandrasekhar limit. The results of a model atmosphere analysis are consistent with our mass estimated from the orbit. Temperatures of the individual components are also determined. Possible scenarios for the formation of this system are discussed. The system will merge due to loss of angular momentum via gravitational wave radiation after two Hubble times. HE 1414-0848 does not qualify as a SN Ia progenitor, but it is the most massive close DD known today.

We have started a large survey for double degenerate (DD) binaries as potential progenitors of type Ia supernovae with the UVES spectrograph at the ESO VLT (ESO SN Ia Progenitor surveY - SPY). About 400 white dwarfs were checked for radial velocity variations during the first 15 months of this project, twice the number of white dwarfs investigated during the last 20 years. We give an overview of the SPY project and present first results Fifty four new DDs have been discovered, seven of them double lined (only 18 and 6 objects of these groups were known before, respectively). The final sample is expected to contain 150 to 200 DDs. Eight new pre-cataclysmic binaries were also detected. SPY is the first DD survey which encompasses also non-DA white dwarfs. SPY produces an immense, unique sample of very high resolution white dwarf spectra, which provides a lot of spin-off opportunities. We describe our projects to exploit the SPY sample for the determination of basic parameters, kinematics, and rotational velocities of white dwarfs. A catalogue with a first subset of our white dwarf data has already been published by Koester et al. 2001.

We investigate the evolution of a binary system which initially has an orbital period of 2^h 17^m and contains a 0.5 M_sun helium star with a white dwarf companion of 0.97 M_sun, similarly to suggested SN Ia candidate progenitor KPD 1930+2752. We show that the helium star completes core helium burning and becomes a white dwarf before components merge. The most probable outcome of the merger of components is formation of a massive white dwarf, despite initially the total mass of the system is above the Chandrasekhar mass.

We review the properties of Galactic binaries containing two compact objects, as derived by means of population synthesis. Using this information we calculate the gravitational wave signal of these binaries. At frequencies below f ~< 2 mHz the double white dwarf population forms an unresolved background for the low-frequency gravitational wave detector LISA. Above this limit some few thousand double white dwarfs and few tens of binaries containing neutron stars will be resolved. Of the resolved double white dwarfs ~500 have a total mass above the Chandrasekhar limit. About ~95 of these have a measurable frequency change allowing a determination of their chirp mass. We discuss the properties of the resolved systems.

We study two models for AM CVn stars: white dwarfs accreting (i) from a helium white dwarf companion and (ii) from a helium-star donor. We show that in the first model possibly no accretion disk forms at the onset of mass transfer. The stability and the rate of mass transfer then depend on the tidal coupling between the accretor and the orbital motion. In the second model the formation of AM CVn stars may be prevented by detonation of the CO white dwarf accretor and the disruption of the system. With the most favourable conditions for the formation of AM CVn stars we find a current Galactic birth rate of 6.8 10^{-3} yr^{-1}. Unfavourable conditions give 1.1 10^{-3} yr^{-1}. The expected total number of the systems in the Galaxy is 9.4 10^{7} and 1.6 10^{7}, respectively. We model very simple selection effects to get some idea about the currently expected observable population and discuss the (quite good) agreement with the observed systems.

We present a new model for the Galactic population of close double white dwarfs. The model accounts for the suggestion of the avoidance of a substantial spiral-in during mass transfer between a giant and a main-sequence star of comparable mass and for detailed cooling models. It agrees well with the observations of the local sample of white dwarfs if the initial binary fraction is close to 50% and an ad hoc assumption is made that white dwarfs with mass less than about 0.3 solar mass cool faster than the models suggest. About 1000 white dwarfs brighter than V=15 have to be surveyed for detection of a pair which has total mass greater than the Chandrasekhar mass and will merge within 10 Gyr.

We model the population of double white dwarfs in the Galaxy and find a better agreement with observations compared to earlier studies, due to two modifications. The first is the treatment of the first phase of unstable mass transfer and the second the modelling of the cooling of the white dwarfs. A satisfactory agreement with observations of the local sample of white dwarfs is achieved if we assume that the initial binary fraction is ~ 50% and that the lowest mass white dwarfs (M < 0.3 Msun) cool faster than the most recently published cooling models predict. With this model we find a Galactic birth rate of close double white dwarfs of 0.05 yr^{-1}, a birth rate of AM CVn systems of 0.005 yr^{-1}, a merger rate of pairs with a combined mass exceeding the Chandrasekhar limit (which may be progenitors of SNe Ia) of 0.003 yr^{-1} and a formation rate of planetary nebulae of 1 yr^{-1}. We estimate the total number of double white dwarfs in the Galaxy as 2.5 10^8. In an observable sample with a limiting magnitude V_lim = 15 we predict the presence of ~855 white dwarfs of which ~220 are close pairs. Of these 10 are double CO white dwarfs of which one has a combined mass exceeding the Chandrasekhar limit and will merge within a Hubble time.

We performed populations synthesis calculations of single stars and binaries and show that binary evolution is extremely important for Galactic astronomy. We review several binary evolution models and conclude that they give quite different results. These differences can be understood from the assumptions related to how mass is transfered in the binary systems. Most important are 1) the fraction of mass that is accreted by the companion star during mass transfer, 2) the amount of specific angular momentum which is carried away with the mass that leaves the binary system.

The unique core-mass - radius relation for giants with degenerate helium cores enables us to reconstruct the evolution of three observed double helium white dwarfs with known masses of both components. The last mass transfer phase in their evolution must have been a spiral-in. In the formalism proposed by Webbink (1984), we can constrain the efficiency of the deposition of orbital energy into the envelope to be $1 \la \alpha \la 6$, for an envelope structure parameter $\lambda=0.5$. We find that the two standard mass transfer types (stable mass transfer and spiral-in) are both unable to explain the first phase of mass transfer for these three binaries. We use a parametric approach to describe mass transfer in low-mass binaries, where both stars are of comparable mass and find that the orbital characteristics of the observed double helium white dwarfs can be well reproduced if the envelope of the primary is lost with ~1.5 times the specific angular momentum of the initial binary. In this case no substantial spiral-in occurs.

We suggest a new mechanism to account for the loss of angular momentum in binaries with non-conservative mass exchange. It is shown that in some cases the loss of matter can result in increase of the orbital angular momentum of a binary. If included into consideration in evolutionary calculations, this mechanism appreciably extends the range of mass ratios of components for which mass exchange in binaries is stable. It becomes possible to explain the existence of some observed cataclysmic binaries with high donor/accretor mass ratio, which was prohibited in conservative evolution models.

We discuss the cosmic history of supernovae on the basis of various assumptions and recent data on the star formation history. We show that supernova rates as a function of redshift can be used to place significant constraints on progenitor models, on the star formation history, and on the importance of dust obscuration. We demonstrate that it is unlikely that the current observational indications for the existence of a cosmological constant are merely an artifact of the dominance of different progenitor classes at different redshift intervals.

We discuss the formation of pulsars with massive companions in eccentric orbits. We demonstrate that the probability for a non-recycled radio pulsar to have a white dwarf as a companion is comparable to that of having an old neutron star as a companion. Special emphasis is given to PSR B1820-11 and PSR B2303+46. Based on population synthesis calculations we argue that PSR B1820-11 and PSR B2303+46 could very well be accompanied by white dwarfs with mass > 1.1Msun. For PSR B1820-11, however, we can not exclude the possibility that its companion is a main-sequence star with a mass between about 0.7Msun and 5Msun.

We describe radial velocity observations of a large sample of apparently single white dwarfs (WDs), obtained in a long-term effort to discover close, double-degenerate (DD) pairs which might comprise viable Type Ia Supernova (SN Ia) progenitors. We augment the WD sample with a previously observed sample of apparently single subdwarf B (sdB) stars, which are believed to evolve directly to the WD cooling sequence after the cessation of core helium burning. We have identified 18 new radial velocity variables, including five confirmed sdB+WD short-period pairs. Our observations are in general agreement with the predictions of the theory of binary star evolution. We describe a numerical method to evaluate the detection efficiency of the survey and estimate the number of binary systems not detected due to the effects of varying orbital inclination, orbital phase at the epoch of the first observation, and the actual temporal sampling of each object in the sample. Follow-up observations are in progress to solve for the orbital parameters of the candidate velocity variables.

The evolution of the comoving cosmic merger-rate density of neutron star binaries n_c(z) is calculated using a distribution of their merging times provided by population-synthesis computations of binary stars. We adopt an exponential law for the star formation rate with various timescales for different morphological types of galaxies. For elliptical galaxies also an initial burst of star formation, lasting one Gyr, is considered. The resulting n_c(z) of most models agree with the form n_c(z) \propto (1+z)^{1.5-2} for z < 2, which has been proposed for the source population of gamma-ray bursts. Assuming a standard candle luminosity, the computed brightness distribution is consistent with the BATSE results if bursts at the peak flux threshold, P = 0.4 photons cm^{-2}/s, are located at a limiting redshift of 1.9 to 3.3. Progenitors of the systems producing gamma-ray bursts at small redshift (bright) are likely to host in spiral galaxies and star forming regions whereas these at high redshift (dim) reside mainly in elliptical galaxies. The location of a burst may be up to a Mpc away from the host galaxy.

We discuss the formation and evolution of binaries which contain neutron stars or black holes. It is shown that in a stellar system which for $10^{10}$ yr had star formation rate similar to the current one in the Galactic disc, the rate of neutron star mergers is $\sim 2 \times 10^{-5} yr{-1}$, consistent with the observational estimates. The rate of black hole and neutron star mergers is by an order of magnitude lower. The pairs of black holes form but do not merge because heavy mass loss by their progenitors efficiently moves the binary components apart.

We examine the possibility that supernovae type Ia (SN Ia) are produced by white dwarfs accreting from Roche-lobe filling evolved companions, under the assumption that a strong optically thick stellar wind from accretor is able to stabilize the mass transfer. We show that if a mass transfer phase on a thermal timescale precedes a nuclear burning driven phase, then such systems (of which the supersoft X-ray sources are a subgroup) can account for about 10% of the inferred SN Ia rate. In addition, we examine the cosmic history of the supernova rate, and we show that the ratio of the rate of SN Ia to the rate of supernovae produced by massive stars (supernovae of types II, Ib, Ic) should increase from about z = 1 towards lower redshifts.

By means of population synthesis we find that the expected Galactic number of black holes with massive helium star companions is \sim 100 and depends on an assumed threshold for M_{pre-BH}. The overwhelming majority of these systems has orbital periods in excess of 10 hr, with a maximum at \sim 100 hr, while under the Illarionov & Sunyaev (1975) disk formation criteria for accretion from the strong stellar wind of Wolf-Rayet star disk accretion is possible only for orbital periods below \sim 10 hr. However, the number of such short-period systems is vanishingly small. If the accretor in Cyg X-3 is a 10 M_sun black hole, then the accretion rate will be super-Eddington. Super-Eddington accretion may be responsible for the formation of jets in Cyg X-3 and may also support an X-ray luminosity as high as \sim 10^{39} erg/s. From the orbital period distribution for neutron stars with massive helium companions we find that if during the common envelope phase a neutron star accretes at Eddington rate and spins-up to the equilibrium period, then in most systems it acts as a ``propeller'' and accretion from the WR star wind is impossible. For the model with two massive helium stars as an immediate progenitor of Cyg X-3, requirement of accomodation of two WR stars in the post-common-envelope orbit combined with severe mass loss by them prevents formation of BH+WR systems with orbital periods less than several days.

The formation and evolution of binaries which contain two neutron stars or a neutron star with a black hole are discussed in detail. The evolution of the distributions in orbital period and eccentricity for neutron star binaries are studied as a function of time. In the model which fits the observations of high mass binary pulsars best the deposition of orbital energy into common envelopes has to be very efficient and a kick velocity distribution has to contain a significant contribution of low velocity kicks. The estimated age of the population has to be between several 100 Myr and 1 Gyr. The birthrate of binary neutron stars is about 3.4 10^{-5} per year (assuming 100% binarity) and their merger rate is about 2 x 10^{-5} per year. The merger rate of neutron star binaries is consistent with the estimated rate of gamma-ray bursts, if the latter are beamed into an opening angle of a few degrees. We argue that PSR B2303+46 is possibly formed in a scenario in which the common envelope is avoided while for the other three known high-mass binary pulsars a common envelope is required to explain their orbital period.