Theoretical population number studies of binaries with at least one black hole (BH) component are obviously depending on whether or not BHs receive a (natal) kick during their formation. Several observational facts seem to indicate that indeed BHs receive a kick during their formation. In the present paper we discuss additional evidence. The progenitors of wind fed high mass X-ray binaries (HMXB) with a BH component (BH HMXB) are WR+OB binaries where the Wolf-Rayet (WR) star will finally collapse and form the BH. Starting from the observed population of WR+OB binaries in the Solar Neighborhood we predict the population of wind fed BH HMXBs as a function of the BH-natal kick. The simulations reveal that when WR stars collapse into a BH with zero or low kick, we would expect 100 or more wind fed BH HMXBs in the Solar Neighborhood whereas only one is observed (Cyg X-1). We consider this as evidence that either WR components in binaries end their life as a neutron star or that they collapse into BHs both accompanied by a supernova explosion imparting significant (natal) kicks.

Despite 50 years of extensive binary research we have to conclude that the Roche lobe overflow/mass transfer process that governs close binary evolution is still poorly understood. It is the scope of the present paper to lift a tip of the veil by studying the spin-up and spin-down processes of the O-type components of WR+O binaries. We critically analyze the available observational data of rotation speeds of the O-type components in WR+O binaries. By combining a binary evolutionary code and a formalism that describes the effects of tides in massive stars with an envelope in radiative equilibrium, we compute the corresponding rotational velocities during the Roche lobe overflow of the progenitor binaries. In all the studied WR+O binaries, we find that the O-type stars were affected by accretion of matter during the RLOF of the progenitor. This means that common envelope evolution which excludes any accretion onto the secondary O-star, has not played an important role to explain the WR+O binaries. Moreover, although it is very likely that the O-type star progenitors were spun-up by the mass transfer, many ended the RLOF/mass transfer phase with a rotational velocity that is significantly smaller than the critical rotation speed. This may indicate that during the mass transfer phase there is a spin-down process which is of the same order as, although significantly less than that of the spin-up process. We propose a Spruit-Tayler type dynamo spin-down suggested in the past to explain the rotation speeds of the mass gainers in long-period Algols.

The role of mass and momentum exchanges in close massive binaries is very important in the subsequent evolution of the components. Such exchanges produce several observational signatures such as asynchronous rotation and altered chemical compositions, that remain after the stars detach again. We investigated these effects for the close O-star binary LSS 3074 (O4 f + O6-7 :(f):), which is a good candidate for a past Roche lobe overflow (RLOF) episode because of its very short orbital period, P = 2.185 days, and the luminosity classes of both components. We determined a new orbital solution for the system. We studied the photometric light curves to determine the inclination of the orbit and Roche lobe filling factors of both stars. Using phase-resolved spectroscopy, we performed the disentangling of the optical spectra of the two stars. We then analysed the reconstructed primary and secondary spectra with the CMFGEN model atmosphere code to determine stellar parameters, such as the effective temperatures and surface gravities, and to constrain the chemical composition of the components. We confirm the apparent low stellar masses and radii reported in previous studies. We also find a strong overabundance in nitrogen and a strong carbon and oxygen depletion in both primary and secondary atmospheres, together with a strong enrichment in helium of the primary star. We propose several possible evolutionary pathways through a RLOF process to explain the current parameters of the system. We confirm that the system is apparently in overcontact configuration and has lost a significant portion of its mass to its surroundings. We suggest that some of the discrepancies between the spectroscopic and photometric properties of LSS 3074 could stem from the impact of a strong radiation pressure of the primary.

We first discuss a large set of evolutionary calculations of close binaries with a B-type primary at birth and with a period such that the Roche lobe overflow starts during the core hydrogen burning phase of the primary (intermediate mass and massive case A binaries). The evolution of both components is followed simultaneously allowing us to check for the occurrence of contact binaries. We consider various models to treat a contact system and the influence of these models on the predicted Algol-type system population is investigated. We critically discuss the available observations of Algol-type binaries with a B-type primary at birth. Comparing these observations with the predictions allows us to put constraints on the contact star physics. We find that mass transfer in Algols is most probably not conservative, that contact during this phase does not necessarily lead to a merger, and that angular momentum loss must be moderate.

In the present review we discuss the past and present status of the interacting OB-type binary frequency. We critically examine the popular idea that Be-stars and supergiant sgB[e] stars are binary evolutionary products. The effects of rotation on stellar evolution in general, stellar population studies in particular, and the link with binaries will be evaluated. Finally a discussion is presented of massive double compact star binary mergers as possible major sites of chemical enrichment of r-process elements and as the origin of recent aLIGO GW events.

In order to investigate the temporal evolution of binary populations in general, double compact star binaries and mergers in particular within a galactic evolution context, a most straightforward method is obviously the implementation of a detailed binary evolutionary model in a galactic chemical evolution code. To our knowledge, the Brussels galactic chemical evolution code is the only one that fully consistently accounts for the important effects of interacting binaries on the predictions of chemical evolution. With a galactic code that does not explicitly include binaries, the temporal evolution of the population of double compact star binaries and mergers can be estimated with reasonable accuracy if the delayed time distribution (DTD) for these mergers is available. The DTD for supernovae type Ia has been studied extensively the last decade. In the present paper we present the DTD for merging double neutron star binaries and mixed systems consisting of a neutron star and a black hole. The latter mergers are very promising sites for the production of r-process elements and the DTDs can be used to study the galactic evolution of these elements with a code that does not explicitly account for binaries.

We first give a short historical overview with some key facts of massive star population synthesis with binaries. We then discuss binary population codes and focus on two ingredients which are important for massive star population synthesis and which may be different in different codes. Population simulations with binaries is the third part where we consider the initial massive binary frequency, the RSG/WR and WC/WN and SNII/SNIbc number ratio's, the probable initial rotational velocity distribution of massive stars.

We investigate the contribution to the formation of type Ia supernovae of the single (a white dwarf accreting from a non-degenerate companion) and double (two merging white dwarfs) degenerate scenario, as well as various aspects of the binary evolution process leading to such a progenitor system. We aim to get a better insight into uncertainties and parameter spaces by means of a combined modelization, resulting in a reduction of the number of possible model approaches. This exclusion of parts of the parameter space is, as will be shown, often independent of uncertainties in the modelization. We use the combination of a population synthesis code with detailed binary evolution and a galactic chemical evolution model to predict the metallicity distribution of G-type dwarfs in the solar neighborhood. Because of the very long lifetime of these stars, this distribution is a good indicator of the entire chemical history of a region. By comparing the observed distribution with those predicted by assuming different type Ia supernova progenitors and evolutionary parameters (e.g. concerning mass and angular momentum loss and common envelope evolution), it is possible to constrain the possible combinations of assumptions. We find that in order to reproduce the observed G-dwarf metallicity distribution, it is absolutely necessary to include both the single and double degenerate scenario. The best match is obtained when all merging C-O white dwarfs contribute to the latter. The correspondence is also critically dependent on the assumptions about galaxy and star formation, e.g. the use of the two-infall model vs. a constant star formation rate. However, this does not affect the previous conclusion, which is consistent with the results obtained by investigating type Ia supernova delay time distributions in starburst galaxies.

With our galactic evolutionary code that contains a detailed intermediate mass and massive binary population model, we study the temporal evolution of the galactic population of double neutron star binaries, mixed systems with a neutron star and black hole component and double black hole binaries. We compute the merger rates of these relativistic binaries and we translate them into LIGO II detection rates. We demonstrate that accounting for the uncertainties in the relation 'initial mass-final mass' predicted by massive close binary evolution and due to the possible effect of large stellar wind mass loss during the luminous blue variable phase of a star with initial mass larger than 30-40 Mo and during the red supergiant phase of a star with initial mass smaller than 30-40 Mo when such a star is a binary component, the double black hole merger rate may be very small, contrary to predictions made by other groups. Hydrodynamic computations of r-process chemical yields ejected during the relativistic binary merger process have recently become available. With our galactic code that includes binaries it is then straightforward to calculate the temporal galactic evolution of the r-process elements ejected by these mergers. We conclude that except for the earliest evolutionary phase of the Galaxy (~the first 100 Myr) double compact star mergers may be the major production sites of r-process elements and it is probable that the mixed systems dominate this production over double neutron star binary mergers.

Binary population synthesis (BPS) modelling is a very effective tool to study the evolution and properties of close binary systems. The uncertainty in the parameters of the model and their effect on a population can be tested in a statistical way, which then leads to a deeper understanding of the underlying physical processes involved. To understand the predictive power of BPS codes, we study the similarities and differences in the predicted populations of four different BPS codes for low- and intermediate-mass binaries. We investigate whether the differences are caused by different assumptions made in the BPS codes or by numerical effects. To simplify the complex problem of comparing BPS codes, we equalise the inherent assumptions as much as possible. We find that the simulated populations are similar between the codes. Regarding the population of binaries with one WD, there is very good agreement between the physical characteristics, the evolutionary channels that lead to the birth of these systems, and their birthrates. Regarding the double WD population, there is a good agreement on which evolutionary channels exist to create double WDs and a rough agreement on the characteristics of the double WD population. Regarding which progenitor systems lead to a single and double WD system and which systems do not, the four codes agree well. Most importantly, we find that for these two populations, the differences in the predictions from the four codes are not due to numerical differences, but because of different inherent assumptions. We identify critical assumptions for BPS studies that need to be studied in more detail.

From the eclipsing binary VW Cep, 1963 epochs of primary and secondary minima have been registered since 1926. Because VW Cep is a shallow contact binary the direction of the flow from one component to the other could possibly be reversed. We aim to fit the observed eclipses into a scenario in which reversals may occur. Orbital period changes can not be fitted into a unique quadratic ephemeris, because of discontinuities in the (O-C) vs E graph. We try to obtain better agreement through the allowance of RLOF reversal in this binary system. From 1926 until 1943 the orbital period increase fits into a model of mass transfer from the less massive component. From 1946 on the direction of the mass transfer changed so that the most massive component became the donor star, as shown by a systematic period decrease. A sudden period increase in 1999 introduces a currently underway episode of mass transfer from the less massive component. The ephemeris of VW Cep can be fitted into a model in which the direction of the mass transfer was reversed twice since 1926.

In the present paper we show that within all the uncertainties that govern the process of Roche lobe overflow in Case Br type massive binaries, it can not be excluded that a significant fraction of them merge and become single stars. We demonstrate that at least some of them will spend most of their core helium burning phase as hydrogen rich blue stars, populating the massive blue supergiant region and/or the massive Be type star population. The evolutionary simulations let us suspect that these mergers will explode as luminous hydrogen rich stars and it is tempting to link them to at least some super luminous supernovae.

Despite the significance of Type Ia supernovae (SNeIa) in many fields in astrophysics, SNeIa lack a theoretical explanation. The standard scenarios involve thermonuclear explosions of carbon/oxygen white dwarfs approaching the Chandrasekhar mass; either by accretion from a companion or by a merger of two white dwarfs. We investigate the contribution from both channels to the SNIa rate with the binary population synthesis (BPS) code SeBa in order to constrain binary processes such as the mass retention efficiency of WD accretion and common envelope evolution. We determine the theoretical rates and delay time distribution of SNIa progenitors and in particular study how assumptions affect the predicted rates.

The chemical processes during the Asymptotic Giant Branch (AGB) evolution of intermediate mass single stars predict most of the observations of the different populations in Globular Clusters although some important issues still need to be further clarified. In particular, to reproduce the observed anticorrelations of Na-O and Al-Mg, chemically enriched gas lost during the AGB phase of intermediate mass single stars must be mixed with matter with a pristine chemical composition. The source of this matter is still a matter of debate. Furthermore, observations reveal that a significant fraction of the intermediate mass and massive stars are born as components of close binaries. We will investigate the effects of binaries on the chemical evolution of Globular Clusters and on the origin of matter with a pristine chemical composition that is needed for the single star AGB scenario to work. We use a population synthesis code that accounts for binary physics in order to estimate the amount and the composition of the matter returned to the interstellar medium of a population of binaries. We demonstrate that the mass lost by a significant population of intermediate mass close binaries in combination with the single star AGB pollution scenario may help to explain the chemical properties of the different populations of stars in Globular Clusters.

We compare the results of three binary population synthesis codes to understand the differences in their results. As a first result we find that when equalizing the assumptions the results are similar. The main differences arise from deviating physical input.

We analyze multicolor light curves and high resolution optical spectroscopy of the eclipsing binary and Double Periodic Variable OGLE 05155332-6925581. According to Mennickent et al., this system shows a significant change in the long non-orbital photometric cycle, a loop in the color-magnitude diagram during this cycle and discrete spectral absorption components that were interpreted as evidence of systemic mass loss. We find that the best fit to the multi-band light curves requires a circumprimary optically thick disc with a radius about twice the radius of the more massive star. The spectroscopy indicates a mass ratio of 0.21+-0.02 and masses for the hot and cool stars of 9.1+-0.5 and 1.9+-0.2 M_sun, respectively. A comparison with synthetic binary-star evolutionary models indicates that the system has an age of 4.76E7 years, is in the phase of rapid mass transfer, the second one in the life of this binary, in a Case-B mass-exchange stage. Donor-subtracted H_alpha profiles show the presence of double emission formed probably in an optically thin circumstellar medium, while the variable HeI profile and the H_beta absorption wings are probably formed in the optically thick circumprimary disc. The model that best fit the observations shows the system with a relatively large mass transfer rate of dM/dt = 3.1E-6 M_sun/yr. However, the orbital period remains relatively stable during almost 15 years. This observation suggests that the hot-spot mass-loss model proposed by other authors is not adequate in this case, and that some other mechanism is efficiently removing angular momentum from the binary. Furthermore, our observations suggest that the DPV phenomenon could have an important effect in the balance of mass and angular momentum in the system.

Context. The chemical processes during the Asymptotic Giant Branch (AGB) evolution of intermediate mass single stars predict most of the observations of the different populations in globular clusters although some important issues still need to be further clarified. In particular, to reproduce the observed anticorrelations of Na-O and Al-Mg, chemically enriched gas lost during the AGB phase of intermediate mass single stars must be mixed with matter with a pristine chemical composition. The source of this matter is still a matter of debate. Furthermore, observations reveal that a significant fraction of the intermediate mass and massive stars are born as components of close binaries. Aims. We will investigate the effects of binaries on the chemical evolution of Globular Clusters and on the origin of matter with a pristine chemical composition that is needed for the single star AGB scenario to work Methods. We use a population synthesis code that accounts for binary physics in order to estimate the amount and the composition of the matter returned to the interstellar medium of a population of binaries. Results. We demonstrate in the present paper that the mass lost by a significant population of intermediate mass close binaries in combination with the single star AGB pollution scenario may help to explain the chemical properties of the different populations of stars in Globular Clusters.

Using a population number synthesis code with detailed binary evolution, we calculate the distribution of the number of type Ia supernovae as a function of time after starburst. This is done for both main progenitor scenarios (single degenerate and double degenerate), but also with various evolutionary assumptions (such as mass transfer efficiency, angular momentum loss, and common envelope description). The comparison of these theoretically predicted delay time distributions with observations in elliptical galaxies then allows to constrain the evolutionary scenarios and parameters. From the morphological shape of the distributions, we conclude that all supernovae Ia cannot be produced through the single degenerate scenario alone, with the best match being obtained when both scenarios contribute. Within the double degenerate scenario, most systems go through a phase of quasi-conservative, stable Roche lobe overflow. We propose stellar rotation as a possible solution for the underestimation of the observed absolute number of events, as is the case in many theoretical population synthesis studies. A brief comparison with these other studies is made, showing good correspondence under the nontrivial condition of equivalent assumptions. We also investigate the influence of different supernova Ia progenitors and evolutionary parameters on the theoretical distribution of the iron abundance of G-type dwarfs in the Galactic disk. These stars are good indicators of the entire chemical history of the Galaxy, and their predicted metallicity distribution can also be compared to the observational ones. This again limits the number of acceptable combinations of assumptions. Supporting previous results, the best correspondence is found in the case where both the single and double degenerate scenario contribute.

Liberal evolution of interacting binaries has been proposed previously by several authors in order to meet various observed binary characteristics better than conservative evolution does. Since Algols are eclipsing binaries the distribution of their orbital periods is precisely known. The distribution of their mass ratios contains however more uncertainties. We try to reproduce these two distributions theoretically using a liberal scenario in which the gainer star can lose mass into interstellar space as a consequence of its rapid rotation and the energy of a hot spot. In a recent paper (Van Rensbergen et al. 2010, A&A) we calculated the liberal evolution of binaries with a B-type primary at birth where mass transfer starts during core hydrogen burning of the donor. In this paper we include the cases where mass transfer starts during hydrogen shell burning and it is our aim to reproduce the observed distributions of the system parameters of Algol-type semi-detached systems. Our calculations reveal the amount of time that an Algol binary lives with a well defined value of mass ratio and orbital period. We use these data to simulate the distribution of mass ratios and orbital periods of Algols. Binaries with a late B-type initial primary hardly lose any mass whereas those with an early B primary evolve in a non-conservative way. Conservative binary evolution predicts only ~ 12 % of Algols with a mass ratio q above 0.4. This value is raised up to ~ 17 % using our scenario of liberal evolution, which is still far below the ~ 45 % that is observed. Observed orbital periods of Algol binaries larger than one day are faithfully reproduced by our liberal scenario. Mass ratios are reproduced better than with conservative evolution, but the resemblance is still poor.

Two formation scenarios are investigated for type Ia supernovae in elliptical galaxies: the single degenerate scenario (a white dwarf reaching the Chandrasekhar limit through accretion of matter transferred from its companion star in a binary) and the double degenerate scenario (the inspiraling and merging of two white dwarfs in a binary as a result of the emission of gravitational wave radiation). A population number synthesis code is used, which includes the latest physical results in binary evolution and allows to differentiate between certain physical scenarios (such as the description of common envelope evolution) and evolutionary parameters (such as the mass transfer efficiency during Roche lobe overflow). The thus obtained theoretical distributions of type Ia supernova delay times are compared to those that are observed, both in morphological shape and absolute number of events. The critical influence of certain parameters on these distributions is used to constrain their values. The single degenerate scenario alone is found to be unable in reproducing the morphological shape of the observational delay time distribution, while use of the double degenerate one (or a combination of both) does result in fair agreement. Most double degenerate type Ia supernovae are formed through a normal, quasi-conservative Roche lobe overflow followed by a common envelope phase, not through two successive common envelope phases as is often assumed. This may cast doubt on the determination of delay times by using analytical formalisms, as is sometimes done in other studies. The theoretical absolute number of events in old elliptical galaxies lies a factor of at least three below the rates that are observed. While this may simply be the result of observational uncertainties, a better treatment of the effects of rotation on stellar structure could mitigate the discrepancy.

We investigate the contribution of different formation scenarios for type Ia supernovae in elliptical galaxies. The single degenerate scenario (a white dwarf accreting from a late main sequence or red giant companion) is tested against the double degenerate scenario (the spiral-in and merging of two white dwarfs through the emission of gravitational wave radiation). We use a population number synthesis code incorporating the latest physical results in binary evolution, and allowing to differentiate between certain physical scenarios (e.g. description of common envelope evolution) and evolutionary parameters (e.g. mass transfer efficiency during Roche lobe overflow). The obtained theoretical distributions of the delay times of type Ia supernovae are compared to those which are observed, both in morphological shape and absolute number of events. The critical dependency of these distributions on certain parameters is used to constrain the values of the latter. We find that the single degenerate scenario alone cannot explain the morphological shape of the observational delay time distribution, while the double degenerate scenario (or a combination of both) can. Most of these double degenerate type Ia supernovae are created through a normal quasi-conservative Roche lobe overflow followed by a common envelope phase, not through two successive common envelope phases. This may cast doubt on the use in other studies of analytical formalisms to determine delay times. In terms of absolute number, theoretical supernova Ia rates in old elliptical galaxies lie a factor of at least three below the observed ones. We propose a solution involving the effect of rotation on the evolution of intermediate mass binaries.

In a liberal evolutionary scenario, mass can escape from a binary during eras of fast mass transfer. We calculate the mass lost by binaries with a B-type primary at birth where mass transfer starts during hydrogen core burning of the donor. We simulate the distribution of mass-ratios and orbital periods for those interacting binaries. The amount of time the binary shows Algol characteristics within different values of mass-ratio and orbital period has been fixed from conservative and liberal evolutionary calculations. We use these data to simulate the distribution of mass-ratios and orbital periods of Algols with the conservative as well as the liberal model. We compare mass-ratios and orbital periods of Algols obtained by conservative evolution with those obtained by our liberal model. Since binaries with a late B-type primary evolve almost conservatively, the overall distribution of mass-ratios will only yield a few Algols more with high mass-ratios than conservative calculations do. Whereas the simulated distribution of orbital periods of Algols fits the observations well, the simulated distribution of mass-ratios produces always too few systems with large values.

The observed distribution of periods and mass ratios of Algols with a B type primary at birth was updated. Conservative evolution fails to produce the large fraction with a high mass ratio: i.e. q in [0.4-0.6]. Interacting binaries thus have to lose mass before or during Algolism. During RLOF mass is transferred continuously from donor to gainer. The gainer spins up; sometimes up to critical velocity. Equatorial material on the gainer is therefore less bound to the system. The material coming from the donor through the first Langrangian point impinges violently on the surface of the gainer or the edge of the accretion disc, creating a hot spot in the area of impact. The sum of rotational energy (fast rotation) and radiative energy (hot spot) depends on the mass-loss rate. The sum of both energies on a test mass located in the impact area equals exactly its binding energy at some critical value. As long as the mass transfer rate is smaller than this value the gainer accepts all the mass coming from the donor: RLOF happens conservatively. But as soon as the critical rate is exceeded the gainer will acquire no more than the critical value and RLOF runs into its liberal era. Low mass binaries never achieve mass-loss rates larger than the critical value. Intermediate mass binaries evolve mainly conservatively but mass will be blown away from the system during the short era of rapid mass transfer soon after RLOF-ignition. Binaries with 9+5.4 solar masses (P in [2-4] d) evolve almost always conservatively. Only during some 20,000 years the gainer is not capable of grasping all the material that comes from the donor. During this short lapse of time a significant fraction of the mass of the system is blown into interstellar space. The mass ratio bin [0.4-0.6] is now much better represented.

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 first summarize work that has been done on the effects of binaries on theoretical population synthesis of stars and stellar phenomena. Next, we highlight the influence of stellar dynamics in young clusters by discussing a few candidate UFOs (unconventionally formed objects) like intermediate mass black holes, Eta Carinae, Zeta Puppis, Gamma Velorum and WR 140.

In the present paper we combine an N-body code that simulates the dynamics of young dense stellar systems with a massive star evolution handler that accounts in a realistic way for the effects of stellar wind mass loss. We discuss two topics: 1. The formation and the evolution of very massive stars (with a mass >120 Mo) is followed in detail. These very massive stars are formed in the cluster core as a consequence of the successive (physical) collison of 10-20 most massive stars of the cluster (the process is known as runaway merging). The further evolution is governed by stellar wind mass loss during core hydrogen burning and during core helium burning (the WR phase of very massive stars). Our simulations reveal that as a consequence of runaway merging in clusters with solar and supersolar values, massive black holes can be formed but with a maximum mass of 70 Mo. In small metallicity clusters however, it cannot be excluded that the runaway merging process is responsible for pair instability supernovae or for the formation of intermediate mass black holes with a mass of several 100 Mo. 2. Massive runaways can be formed via the supernova explosion of one of the components in a binary (the Blaauw scenario) or via dynamical interaction of a single star and a binary or between two binaries in a star cluster. We explore the possibility that the most massive runaways (e.g., zeta Pup, lambda Cep, BD+433654) are the product of the collision and merger of 2 or 3 massive stars.