We report the results of a long-term spectroscopic monitoring of the A-type supergiant with the B[e] phenomenon 3 Pup = HD 62623. We confirm earlier findings that it is a binary system. The orbital parameters were derived using cross-correlation of the spectra in a range of 4460-4632 A, which contains over 30 absorption lines. The orbit was found circular with a period of $137.4\pm0.1$ days, radial velocity semi-amplitude $K_{1} = 5.0\pm0.8$ km s$^{-1}$, systemic radial velocity $\gamma = +26.4\pm2.0$ km s$^{-1}$, and the mass function $f(m) = (1.81^{+0.97}_{-0.76})\times10^{-3}$ M$_{\odot}$. The object may have evolved from a pair with initial masses of $\sim$6.0 M$_{\odot}$ and $\sim$3.6 M$_{\odot}$ with an initial orbital period of $\sim$5 days. Based on the fundamental parameters of the A-supergiant (luminosity $\log$ L/L$_{\odot} = 4.1\pm$0.1 and effective temperature T$_{\rm eff} = 8500\pm$500 K) and evolutionary tracks of mass-transferring binaries, we found current masses of the gainer M$_{2} = 8.8\pm$0.5 M$_{\odot}$ and donor M$_{1} = 0.75\pm0.25$ M$_{\odot}$. We also modeled the object's IR-excess and derived a dust mass of $\sim 5\,\times10^{-5}$ M$_{\odot}$ in the optically-thin dusty disk. The orbital parameters and properties of the H$\alpha$ line profile suggest that the circumstellar gaseous disk is predominantly circumbinary. The relatively low mass of the gainer led us to a suggestion that 3 Pup should be excluded from the B[e] supergiant group and moved to the FS CMa group. Overall these results further support our original suggestion that FS CMa objects are binary systems, where an earlier mass-transfer caused formation of the circumstellar envelope.

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.

Some massive, merging black holes (BH) may be descendants of binary O stars. The evolution and mass transfer between these O stars determines the spins of their progeny BH. These will be measurable with future gravitational wave detectors, incentivizing the measurement of the spins of O stars in binaries. We previously measured the spins of O stars in Galactic Wolf-Rayet (WR) + O binaries. Here we measure the vsini of four LMC and two SMC O stars in WR + O binaries to determine whether lower metallicity might affect the spin rates. We find that the O stars in Galactic and Magellanic WR + O binaries display average vsini = 258 +/- 18 km/s and 270 +/- 15 km/s, respectively. Two LMC O stars measured on successive nights show significant line width variability, possibly due to differing orbital phases exhibiting different parts of the O stars illuminated differently by their WR companions. Despite this variability, the vsini are highly super-synchronous but distinctly subcritical for the O stars in all these binaries; thus we conclude that an efficient mechanism for shedding angular momentum from O stars in WR + O binaries must exist. This mechanism, probably related to Roche lobe overflow-created dynamo magnetic fields, prevents nearly 100% breakup spin rates, as expected when RLOF operates, as it must, in these stars. A Spruit-Tayler dynamo and O star wind might be that mechanism.

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.

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.

From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

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.

The black holes (BH) in merging BH-BH binaries are likely progeny of binary O stars. Their properties, including their spins, will be strongly influenced by the evolution of their progenitor O stars. The remarkable observation that many single O stars spin very rapidly can be explained if they accreted angular momentum from a mass-transferring, O-type or Wolf-Rayet companion before that star blew up as a supernova. To test this prediction, we have measured the spin rates of eight O stars in Wolf-Rayet (WR) + O binaries, increasing the total sample size of such O stars' measured spins from two to ten. Polarimetric and other determinations of these systems' sin i allow us to determine an average equatorial rotation velocity from HeI (HeII) lines of $v_e$ = 348 (173) km/s for these O stars, with individual star's $v_e$ from HeI (HeII) lines ranging from 482 (237) to 290 (91) km/s. We argue that the $\sim$ 100\% difference between HeI and HeII speeds is due to gravity darkening. Super-synchronous spins, now observed in all 10 O stars in WR + O binaries where it has been measured, are strong observational evidence that Roche lobe overflow mass transfer from a WR progenitor companion has played a critical role in the evolution of WR+OB binaries. While theory predicts that this mass transfer rapidly spins-up the O-type mass gainer to a nearly break-up rotational velocity $v_e \sim 530 $ km/s, the observed average $v_e$ of the O-type stars in our sample is 65\% that speed. This demonstrates that, even over the relatively short WR-phase timescale, tidal and/or other effects causing rotational spin-down must be efficient. A challenge to tidal synchronization theory is that the two longest-period binaries in our sample (with periods of 29.7 and 78.5 days) unexpectedly display super-synchronous rotation.

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.

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.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

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.

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.

Context: Starbursts, and particularly their high-mass stars, play an essential role in the evolution of galaxies. The winds of massive stars not only significantly influence their surroundings, but the mass loss also profoundly affects the evolution of the stars themselves. In addition to the evolution of each star, the evolution of the dense cores of massive starburst clusters is affected by N-body interactions, and the formation of very massive stars via mergers may be decisive for the evolution of the cluster. Aims: To introduce an advanced diagnostic method of O-type stellar atmospheres with winds, including an assessment of the accuracy of the determinations of abundances, stellar and wind parameters. Methods: We combine consistent models of expanding atmospheres with detailed stellar evolutionary calculations of massive and very massive single stars with regard to the evolution of dense stellar clusters. Accurate predictions of the mass loss rates of very massive stars requires a highly consistent treatment of the statistical equilibrium and the hydrodynamic and radiative processes in the expanding atmospheres. Results: We present computed mass loss rates, terminal wind velocities, and spectral energy distributions of massive and very massive stars of different metallicities, calculated from atmospheric models with an improved level of consistency. Conclusions: Stellar evolutionary calculations using our computed mass loss rates show that low-metallicity very massive stars lose only a very small amount of their mass, making it unlikely that very massive population III stars cause a significant helium enrichment of the interstellar medium. Solar-metallicity stars have higher mass-loss rates, but these are not so high to exclude very massive stars formed by mergers in dense clusters from ending their life massive enough to form intermediate-mass black holes.

We first discuss the stellar and wind parameters of zeta Pup resulting from detailed UV diagnostics. These parameters together with the runaway nature of the star can most easily be explained by dynamical binary-binary or binary-single star interactions in dense stellar clusters. In this case zeta Pup is most likely the merger of at least two massive stars.

Massive stars are among the most important objects in the Universe and many (most?) of them are formed in binaries. A selection of observational and theoretical facts that illustrate the importance of binaries and the evolution of massive and very massive stars in clusters with special emphasis on massive binaries have been summarized in two recent review papers (Vanbeveren, 2009, 2010). The present paper can be considered as an addendum of both reviews.

Before binary components interact, they evolve as single stars do. We therefore first critically discuss massive single star processes which affect their evolution, stellar wind mass loss and rotation in particular. Next we consider binary processes and focus on the effect of rotation on binary evolution and on the mass transfer during Roche lobe overflow. The third part highlights the importance of close pairs on the comprehension of the evolution of stellar populations in starburst regions.

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 this review I first summarize why binaries are key objects in the study of stellar populations, key objects to understand the evolution of star clusters, key objects to understand galaxies and thus the universe. I then focus on 4 specific topics: 1. the formation (via binaries) and evolution of very massive stars in dense clusters and the importance of stellar wind mass loss. I discuss preliminary computations of wind mass loss rates of very massive stars performed with the Munich hydrodynamical code, and the influence of these new rates on the possible formation of an intermediate mass black hole in the cluster MGG 11 in M82 2. the evolution of intermediate mass binaries in a starburst with special emphasis on the variation of the SN Ia rate (the delayed time distribution of SN Ia). A comparison with SN ia rates in elliptical galaxies may provide important clues on the SN Ia model as well as on the evolution of the SN Ia progenitors 3. the evolution of the double neutron stars mergers in a starburst (the delayed time distribution of these mergers) and what this tells us about the suggestion that these mergers may be important production sites of r-process elements 4. the possible effect of massive binaries on the self-enrichment of globular clusters.

The present paper reviews massive star (initial mass smaller than 120 M0) and very massive star (initial mass larger than 120 M0) evolution. I will focus on evolutionary facts and questions that may critically affect predictions of population and spectral synthesis of starburst regions. We discuss the ever-lasting factor 2 or more uncertainty in the stellar wind mass loss rates. We may ask ourselves if stellar rotation is one of the keys to understand the universe, why so many massive stars are binary components and why binaries are ignored or are considered as the poor cousins by some people? And finally, do ultra luminous X-ray sources harbor an intermediate mass black hole with a mass of the order of 1000 M0?

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.

The early evolution of dense stellar systems is governed by massive single star and binary evolution. Core collapse of dense massive star clusters can lead to the formation of very massive objects through stellar collisions ($M\geq$ 1000 \msun). Stellar wind mass loss determines the evolution and final fate of these objects, and decides upon whether they form black holes (with stellar or intermediate mass) or explode as pair instability supernovae, leaving no remnant. We present a computationaly inexpensive evolutionary scheme for very massive stars that can readily be implemented in an N-body code. Using our new N-body code 'Youngbody' which includes a detailed treatment of massive stars as well as this new scheme for very massive stars, we discuss the formation of intermediate mass and stellar mass black holes in young starburst regions. A more detailed account of these results can be found in Belkus et al. 2007.

We discuss differences between massive single star and massive close binary population number synthesis predictions of WR stars. We show that the WC/WN number ratio as function of metallicity depends significantly on whether or not binaries are included. Furthermore, the observed WC(+OB)/WN(+OB) number ratio in the Solar neighborhood seems to indicate that the WR mass loss rates are lower by another factor two compared to recently proposed clumping corrected formalisms. We then demonstrate that the observed lower luminosity distribution of single WN stars can be explained in a satisfactory way by massive single star evolutionary computations where the red supergiant phase is calculated using a stellar wind mass loss rate formalism that is based on recent observations.

Core collapse of dense massive star clusters is unavoidable and this leads to the formation of massive objects, with a mass up to 1000 $\msun$ and even larger. When these objects become stars, stellar wind mass loss determines their evolution and final fate, and decides upon whether they form black holes (with normal mass or with intermediate mass) or explode as a pair instability supernova. In the present paper, we discuss the evolution of very massive stars and we present a convenient evolution recipe that can be implemented in a gravitational N-body code to study the dynamics of dense massive clusters.

Understanding the galaxy in which we live is one of the great intellectual challenges facing modern science. With the advent of high quality observational data, the chemical evolution modeling of our galaxy has been the subject of numerous studies in the last years. However, all these studies have one missing element which is the evolution of close binaries. Reason: their evolution is very complex and single stars only perhaps can do the job. (Un)Fortunately at present we know that a significant fraction of the observed intermediate mass and massive stars are members of a binary or multiple system and that certain objects can only be formed through binary evolution. Therefore galactic studies that do not account for close binaries may be far from realistic. We implemented a detailed binary population in a galactic chemical evolutionary model. Notice that this is not something simple like replacing chemical yields. Here we discuss three topics: the effect of binaries on the evolution of 14N, the evolution of the type Ia supernova rate and the effects on the G-dwarf distribution, the link between the evolution of the r-process elements and double neutron star mergers (candidates of short gamma-ray burst objects).

Low metallicity very massive stars with an initial mass between 140 Mo and 260 Mo can be subdivided into two groups: those between 140 Mo and 200 Mo which produce a relatively small amount of Fe, and those with a mass between 200 Mo and 260 Mo where the Fe-yield ejected during the supernova explosion is enormous. We first demonstrate that the inclusion of the second group into a chemical evolutionary model for the Solar Neighbourhood predicts an early temporal evolution of Fe which is at variance with observations whereas it can not be excluded that the first group could have been present. We then show that a low metallicity binary with very massive components (with a mass corresponding to the first group) can be an efficient site of primary 14N production through the explosion of a binary component that has been polluted by the pair instability supernova ejecta of its companion. When we implement these massive binary 14N yields in a chemical evolution model, we conclude that very massive close binaries may be important sites of 14N enrichment during the early evolution of the Galaxy.

The effects of binaries on population studies of stars and stellar phenomena have been investigated over the past 3 decades by many research groups. Here we will focus mainly on the work that has been done recently in Brussels and we will consider the following topics: the effect of binaries on overall galactic chemical evolutionary models and on the rates of different types of supernova, the population of point-like X-ray sources where we distinguish the standard high mass X-ray binaries and the ULXs, a UFO-scenario for the formation of WR+OB binaries in dense star systems. Finally we critically discuss the possible effect of rotation on population studies.

We present a set of evolutionary computations for binaries with a B-type primary at birth. Some liberal computations including loss of mass and angular momentum during binary evolution are added to an extensive grid of conservative calculations. Our computations are compared statistically to the observed distributions of orbital periods and mass ratios of Algols. Conservative Roche Lobe Over Flow (RLOF) reproduces the observed distribution of orbital periods decently but fails to explain the observed mass ratios in the range 0.4-1. In order to obtain a better fit the binaries have to lose a significant amount of matter, without transferring too much angular momentum.

In the present paper we report on first results of a project in Brussels where we study the effects of stellar dynamics on the evolution of young dense stellar systems using the 3 decades expertise in massive star evolution and our population (number and spectral) synthesis code. We highlight an unconventionally formed object scenario (UFO-scenario) for Wolf Rayet binaries and study the effects of a luminous blue variable-type instability wind mass loss formalism on the formation of intermediate mass black holes.

Understanding the galaxy in which we live is one of the great intellectual challenges facing modern science. With the advent of high quality observational data, the chemical evolution modeling of our galaxy has been the subject of numerous studies in the last years. However, all these studies have one missing element which is the evolution of close binaries. Reason: their evolution is very complex and single stars only perhaps can do the job. (Un)Fortunately at present we know that the majority of the observed stars are members of a binary or multiple system and that certain objects can only be formed through binary evolution. Therefore galactic studies that do not account for close binary evolution may be far from realistic. Because of the large expertise developed through the years in stellar evolution in general and binary evolution in particular at the Brussels Astrophysical Institute, we found ourselves in a privileged position to be the first to do chemical evolutionary simulations with the inclusion of detailed binary evolution. The complexity of close binary evolution has kept many astronomers from including binary stars into their studies. However, it is not always the easiest way of living that gives you the most excitement and satisfaction.

In the present paper we discuss a selection of facts and questions related to observations and evolutionary calculations of massive single stars and massive stars in interacting binaries. We focus on the surface chemical abundances, the role of stellar winds, the early Be-stars, the high mass X-ray binaries, the effects of rotation on stellar evolution. Finally, we present an unconventionally formed object scenario (a UFO-scenario) of WR binaries in dense stellar environments.

A comparison between the observed O-type star and the WR star populations and the theoretically predicted ones depends on the effects of stellar wind mass loss during various phases and rotation on stellar evolution and, last not least, on the effects of binaries. Obviously, both populations depend on the massive star formation rate. In the present paper we show that the rate in the Solar Neighbourhood is fluctuating on a timescale which is smaller than the evolutionary timescale of a massive star. We demonstrate that this fluctuating behaviour affects both population enormously, in such a way that a comparison between observations and theoretical prediction to test stellar evolutionary models is ambiguous.

Supernova rates (hypernova, type II, type Ib/c and type Ia) in a particular galaxy depend on the metallicity (i.e. on the galaxy age), on the physics of star formation and on the binary population. In order to study the time evolution of the galactic supernova rates, we use our chemical evolutionary model that accounts in detail for the evolution of single stars and binaries. In particular, supernovae of type Ia are considered to arise from exploding white dwarfs in interacting binaries and we adopt the two most plausible physical models: the single degenerate model and the double degenerate model. Comparison between theoretical prediction and observations of supernova rates in different types of galaxies allows to put constraints on the population of intermediate mass and massive close binaries. The temporal evolution of the absolute galactic rates of different types of SNe (including the SN Ia rate) is presented in such a way that the results can be directly implemented into a galactic chemical evolutionary model. Particularly for SNIa the inclusion of binary evolution leads to results considerably different from those in earlier population synthesis approaches, in which binary evolution was not included in detail.

A population number synthesis code follows in detail the evolution of a population of single stars and of close binaries. We use our code to simulate the population of neutron star - neutron star and black hole - neutron star binaries. We then combine our population number synthesis code with a galactic chemical evolutionary model in order to follow the time evolution of the formation and merger rate of these double compact star binaries and the resulting chemical enrichment of r-process elements, over the whole Galactic lifetime. It can be concluded that the neutron star/black hole merger process is able to reproduce the observed r-process enrichment of the Galaxy. However, we show that the latter conclusion depends critically on the physics of case BB Roche lobe overflow in binaries with a neutron star component and a hydrogen deficient core helium/helium shell burning star with a mass between 2.6 Mo and 6 Mo.

In the present review we summarize direct and indirect evidence that the massive close binary frequency is very large. We then discuss the binary evolutionary processes and we present a general massive close binary evolutionary scheme. Finally, we highlight the importance of massive close binaries for population number synthesis and the chemical evolution of galaxies.

In the present paper we investigate in how far stars with an initial mass larger than 40Mo affect the chemical enrichment of the Galaxy. We illustrate the importance for chemical yields of a most up to date treatment of the various stellar wind mass loss episodes in stellar evolutionary codes and we discuss the effects of a possible supernova-like outburst prior to massive black hole formation.

Recent models of galactic chemical evolution account for updated evolutionary models of massive stars (with special emphasis on stellar winds) and for the effects of intermediate mass and massive binaries. The results are summarised. We also present a critical discussion on possible effects of stellar rotation on overall galactic chemical evolutionary simulations.

The evolution of massive stars in general, massive close binaries in particular depend on processes where, despite many efforts, the physics are still uncertain. Here we discuss the effects of stellar wind as function of metallicity during different evolutionary phases and the effects of rotation and we highlight the importance for population (number) synthesis. Models are proposed for the X-ray binaries with a black hole component. We then present an overall scheme for a massive star PNS code. This code, in combination with an appropriate set of single star and close binary evolutionary computations predicts the O-type star and Wolf-Rayet star population, the population of double compact star binaries and the supernova rates, in regions where star formation is continuous in time.

In this paper we compute the time evolution of the elements (4He, 12C, 14N, 16O, 20Ne, 24Mg, 28Si, 32S, 40Ca and 56Fe) and of the supernova rates in the solar neighbourhood by means of a galactic chemical evolutionary code that includes in detail the evolution of both single and binary stars. Special attention is payed to the formation of black holes. Our main conclusions: in order to predict the galactic time evolution of the different types of supernovae, it is essential to compute in detail the evolution of the binary population, the observed time evolution of carbon is better reproduced by a galactic model where the effect is included of a significant fraction of intermediate mass binaries, massive binary mass exchange provides a possible solution for the production of primary nitrogen during the very early phases of galactic evolution, chemical evolutionary models with binaries or without binaries but with a detailed treatment of the SN Ia progenitors predict very similar age-metallicity relations and very similar G-dwarf distributions whereas the evolution of the yields as function of time of the elements 4He, 16O, 20Ne, 24Mg, 28Si, 32S and 40Ca differ by no more than a factor of two or three, the observed time evolution of oxygen is best reproduced when most of the oxygen produced during core helium burning in ALL massive stars serves to enrich the interstellar medium. This can be used as indirect evidence that (massive) black hole formation in single stars and binary components is always preceded by a supernova explosion.