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.

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.

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.