We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars $\gtrsim10^6$ and the gas accretion rate $\gtrsim10^5\,{\rm M}_\odot$/Myr, reminiscent of GC formation in high gas-density environments, such as -- but not restricted to -- the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This `conveyor-belt' production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests.

Observations of globular clusters (GCs) and field stars in the halos of the giant elliptical galaxy Cen A and the spiral galaxy M31 show a large range of cluster-to-star number ratios ('specific frequencies'). The cluster-to-star ratio decreases with increasing metallicity by a factor of 100-1000, at all galactocentric radii and with a slope that does not seem to depend on radius. In dwarf galaxies, the GCs are also more metal-poor than the field stars on average. These observations indicate a strong dependence of either the cluster formation efficiency or the cluster destruction rate on metallicity and environment. We aim to explain these trends by considering various effects that may influence the observed cluster-to-star ratio as a function of metallicity, environment and cosmological history. We show that both the cluster formation efficiency and the maximum cluster mass increase with metallicity, so they cannot explain the observed trend. Destruction of GCs by tidal stripping and dynamical friction destroy clusters over too small a range of galactocentric radii. We show that cluster destruction by tidal shocks from giant molecular clouds in the high-density formation environments of GCs becomes increasingly efficient towards high galaxy masses and, hence, towards high metallicities. The predicted cluster-to-star ratio decreases by a factor 100-1000 towards high metallicities and should only weakly depend on galactocentric radius due to orbital mixing during hierarchical galaxy merging, consistent with the observations. The observed, strong dependence of the cluster-to-star ratio on metallicity and the independence of its slope on galactocentric radius can be explained by cluster destruction and hierarchical galaxy growth. As a result, we find that the metallicity-dependence of the cluster-to-star ratio does not reflect a GC formation efficiency, but a survival fraction. (Abridged)

We present a new version of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). While previous versions of EMACSS reproduced clusters of single-mass stars, this version models clusters with an evolving stellar content. Stellar evolution dominates early evolution, and leads to: (1) reduction of the mean mass of stars due to the mass loss of high-mass stars; (2) expansion of the half-mass radius; (3) for (nearly) Roche Volume filling clusters, the induced escape of stars. Once sufficient relaxation has occurred (~ 10 relaxation times-scales), clusters reach a second, 'balanced' state whereby the core releases energy as required by the cluster as a whole. In this state: (1) stars escape due to tidal effects faster than before balanced evolution; (2) the half-mass radius expands or contracts depending on the Roche volume filling factor; and (3) the mean mass of stars increases due to the preferential ejection of low-mass stars. We compare the EMACSS results of several cluster properties against N-body simulations of clusters spanning a range of initial number of stars, mass, half-mass radius, and tidal environments, and show that our prescription accurately predicts cluster evolution for this database. Finally, we consider applications for EMACSS, such as studies of galactic globular cluster populations in cosmological simulations.

We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the rate of energy generation in the core of the cluster equals the diffusion rate of energy by two-body relaxation, which makes the code suitable for modelling clusters in weak tidal fields. In this new version we extend the model to include an unbalanced phase of evolution to describe the pre-collapse evolution and the accompanying escape rate such that clusters in strong tidal fields can also be modelled. We also add a prescription for the evolution of the core radius and density and a related cluster concentration parameter. The model simultaneously solves a series of first-order ordinary differential equations for the rate of change of the core radius, half-mass radius and the number of member stars N. About two thousand integration steps in time are required to solve for the entire evolution of a star cluster and this number is approximately independent of N. We compare the model to the variation of these parameters following from a series of direct N-body calculations of single-mass clusters and find good agreement in the evolution of all parameters. Relevant time-scales, such as the total lifetimes and core collapse times, are reproduced with an accuracy of about 10% for clusters with various initial half-mass radii (relative to their Jacobi radii) and a range of different initial N up to N = 65536. We intend to extend this framework to include more realistic initial conditions, such as a stellar mass spectrum and mass loss from stars. The EMACSS code can be used in star cluster population studies and in models that consider the co-evolution of (globular) star clusters and large scale structures.

Globular clusters (GCs), once thought to be well approximated as simple stellar populations (i.e. all stars having the same age and chemical abundance), are now known to host a variety of anomalies, such as multiple discrete (or spreads in) populations in colour-magnitude diagrams and abundance variations in light elements (e.g., Na, O, Al). Multiple models have been put forward to explain the observed anomalies, although all have serious shortcomings (e.g., requiring a non-standard initial mass function of stars and GCs to have been initially 10-100 times more massive than observed today). These models also do not agree with observations of massive stellar clusters forming today, which do not display significant age spreads nor have gas/dust within the cluster. Here we present a model for the formation of GCs, where low mass pre-main sequence (PMS) stars accrete enriched material released from interacting massive binary and rapidly rotating stars onto their circumstellar discs, and ultimately onto the young stars. As was shown in previous studies, the accreted material matches the unusual abundances and patterns observed in GCs. The proposed model does not require multiple generations of star-formation, conforms to known properties of massive clusters forming today, and solves the "mass budget problem" without requiring GCs to have been significantly more massive at birth. Potential caveats to the model as well as model predictions are discussed.

The evolution of the global stellar mass function (MF) of star clusters is studied based on a large set of N-body simulations of clusters with a range of initial masses, initial concentrations, in circular or elliptical orbits in different tidal environments. Models with and without initial mass segregation are included. The depletion of low mass stars in initially Roche-volume (tidal) filling clusters starts typically on a time scale of the order of the core collapse time. In clusters that are initially underfilling their Roche-volume it takes longer because the clusters have to expand to their tidal radii before dynamical mass loss becomes important. We introduce the concept of the differential mass function (DMF), which describes the changes with respect to the initial mass function (IMF). We show that the evolution of the DMF can be described by a set of very simple analytic expressions that are valid for a wide range of initial cluster parameters and for different IMFs. The agreement between this description and the models is very good, except for initially Roche-volume underfilling clusters that are severely mass segregated.

(Abridged) Interacting galaxies are well-known for their high star formation rates and rich star cluster populations, but the rapidly changing tidal field can also efficiently destroy clusters. We use numerical simulations of merging disc galaxies to investigate which mechanism dominates. The simulations include a model for the formation and dynamical disruption of the entire star cluster population. We find that the dynamical heating of clusters by tidal shocks is about an order of magnitude higher in interacting galaxies than in isolated galaxies. This is driven by the increased gas density, and is sufficient to destroy star clusters at a higher rate than new clusters are formed: the total number of clusters in the merger remnant is 2-50% of the amount in the progenitor discs, with low-mass clusters being disrupted preferentially. By adopting observationally motivated selection criteria, we find that the observed surplus of star clusters in nearby merging galaxies is caused by the bias to detect young, massive clusters. We provide a general expression for the survival fraction of clusters, which increases with the gas depletion time-scale. Due to the preferential disruption of low-mass clusters, the mass distribution of the surviving star clusters in a merger remnant develops a peak at a mass of about 10^3 Msun, which evolves to higher masses at a rate of 0.3-0.4 dex per Gyr. The peak mass initially depends weakly on the galactocentric radius, but this correlation disappears as the system ages. We discuss the similarities between the cluster populations of the simulated merger remnants and (young) globular cluster systems. Our results suggest that the combination of cluster formation and destruction should be widespread in the dense star-forming environments at high redshifts, which could provide a natural origin to present-day globular cluster systems.

We study the stellar cluster population in two adjacent fields in the nearby, face-on spiral galaxy, M83, using WFC3/HST imaging. The clusters are selected through visual inspection to be centrally concentrated, symmetric, and resolved on the images, which allows us to differentiate between clusters and likely unbound associations. We compare our sample with previous studies and show that the differences between the catalogues are largely due to the inclusion of large numbers of diffuse associations within previous catalogues. The luminosity function of the clusters is well approximated by a power-law with index, -2, over most of the observed range, however a steepening is seen at M_V = -9.3 and -8.8 in the inner and outer fields, respectively. Additionally, we show that the cluster population is inconsistent with a pure power-law mass distribution, but instead exhibits a truncation at the high mass end. If described as a Schechter function, the characteristic mass is 1.6 and 0.5 * 10^5 Msun, for the inner and outer fields, respectively, in agreement with previous estimates of other cluster populations in spiral galaxies. Comparing the predictions of the mass independent disruption (MID) and mass dependent disruption (MDD) scenarios with the observed distributions, we find that both models can accurately fit the data. However, for the MID case, the fraction of clusters destroyed (or mass lost) per decade in age is dependent on the environment, hence, the age/mass distributions of clusters are not universal. In the MDD case, the disruption timescale scales with galactocentric distance (being longer in the outer regions of the galaxy) in agreement with analytic and numerical predictions. Finally, we discuss the implications of our results on other extragalactic surveys, focussing on the fraction of stars that form in clusters and the need (or lack thereof) for infant mortality.

Using multi-wavelength imaging from the Wide Field Camera 3 on the Hubble Space Telescope we study the stellar cluster populations of two adjacent fields in the nearby face-on spiral galaxy, M83. The observations cover the galactic centre and reach out to ~6 kpc, thereby spanning a large range of environmental conditions, ideal for testing empirical laws of cluster disruption. The clusters are selected by visual inspection to be centrally concentrated, symmetric, and resolved on the images. We find that a large fraction of objects detected by automated algorithms (e.g. SExtractor or Daofind) are not clusters, but rather are associations. These are likely to disperse into the field on timescales of tens of Myr due to their lower stellar densities and not due to gas expulsion (i.e. they were never gravitationally bound). We split the sample into two discrete fields (inner and outer regions of the galaxy) and search for evidence of environmentally dependent cluster disruption. Colour-colour diagrams of the clusters, when compared to simple stellar population models, already indicate that a much larger fraction of the clusters in the outer field are older by tens of Myr than in the inner field. This impression is quantified by estimating each cluster's properties (age, mass, and extinction) and comparing the age/mass distributions between the two fields. Our results are inconsistent with "universal" age and mass distributions of clusters, and instead show that the ambient environment strongly affects the observed populations.

(Abridged) The formation and evolution of star cluster populations are related to the galactic environment. Cluster formation is governed by processes acting on galactic scales, and star cluster disruption is driven by the tidal field. In this paper, we present a self-consistent model for the formation and evolution of star cluster populations, for which we combine an N-body/SPH galaxy evolution code with semi-analytic models for star cluster evolution. The model includes star formation, feedback, stellar evolution, and star cluster disruption by two-body relaxation and tidal shocks. We apply the model by simulating a suite of 9 isolated disc galaxies and 24 galaxy mergers. The evolutionary histories of individual clusters in these simulations are discussed to illustrate how the environment of clusters changes in time and space. The resulting variability of the disruption rate with time and space affects the properties of star cluster populations. The combined effect of clusters escaping their dense formation sites (`cluster migration') and the preferential disruption of clusters residing in dense environments (`natural selection') implies that the mean disruption rate of the population decreases with cluster age. This affects the slope of the cluster age distribution, which becomes a function of the star formation rate density. The evolutionary histories of clusters in a galaxy merger vary widely and determine which clusters survive the merger. This impacts the age distributions and the locations of the surviving clusters at all times during a merger. We conclude that accounting for the interplay between the formation, disruption, and orbital histories of clusters enables a more sophisticated interpretation of observed properties of cluster populations, thereby extending the role of cluster populations as tracers of galaxy evolution.

Most globular clusters have half-mass radii of a few pc with no apparent correlation with their masses. This is different from elliptical galaxies, for which the Faber-Jackson relation suggests a strong positive correlation between mass and radius. Objects that are somewhat in between globular clusters and low-mass galaxies, such as ultra-compact dwarf galaxies, have a mass-radius relation consistent with the extension of the relation for bright ellipticals. Here we show that at an age of 10 Gyr a break in the mass-radius relation at ~10^6 Msun is established because objects below this mass, i.e. globular clusters, have undergone expansion driven by stellar evolution and hard binaries. From numerical simulations we find that the combined energy production of these two effects in the core comes into balance with the flux of energy that is conducted across the half-mass radius by relaxation. An important property of this `balanced' evolution is that the cluster half-mass radius is independent of its initial value and is a function of the number of bound stars and the age only. It is therefore not possible to infer the initial mass-radius relation of globular clusters and we can only conclude that the present day properties are consistent with the hypothesis that all hot stellar systems formed with the same mass-radius relation and that globular clusters have moved away from this relation because of a Hubble time of stellar and dynamical evolution.

We describe the interplay between stellar evolution and dynamical mass loss of evolving star clusters, based on the principles of stellar evolution and cluster dynamics and on a grid of N-body simulations of cluster models. The cluster models have different initial masses, different orbits, including elliptical ones, and different initial density profiles. We use two sets of cluster models: initially Roche-lobe filling and Roche-lobe underfilling. We identify four distinct mass loss effects: (1) mass loss by stellar evolution, (2) loss of stars induced by stellar evolution and (3) relaxation-driven mass loss before and (4) after core collapse. Both the evolution-induced loss of stars and the relaxation-driven mass loss need time to build up. This is described by a delay-function of a few crossing times for Roche-lobe filling clusters and a few half mass relaxation times for Roche-lobe underfilling clusters. The relaxation-driven mass loss can be described by a simple power law dependence of the mass dM/dt =-M^{1-gamma}/t0, (with M in Msun) where t0 depends on the orbit and environment of the cluster. Gamma is 0.65 for clusters with a King-parameter W0=5 and 0.80 for more concentrated clusters with W0=7. For initially Roche-lobe underfilling clusters the dissolution is described by the same gamma=0.80. The values of the constant t0 are described by simple formulae that depend on the orbit of the cluster. The mass loss rate increases by about a factor two at core collapse and the mass dependence of the relaxation-driven mass loss changes to gamma=0.70 after core collapse. We also present a simple recipe for predicting the mass evolution of individual star clusters with various metallicities and in different environments, with an accuracy of a few percent in most cases. This can be used to predict the mass evolution of cluster systems.

Star clusters are found in all sorts of environments and their formation and evolution is inextricably linked to the star formation process. Their eventual destruction can result from a number of factors at different times, but the process can be investigated as a whole through the study of the cluster age distribution. Observations of populous cluster samples reveal a distribution following a power law of index approximately -1. In this work we use M33 as a test case to examine the age distribution of an archetypal cluster population and show that it is in fact the evolving shape of the mass detection limit that defines this trend. That is to say, any magnitude-limited sample will appear to follow a dN/dt=1/t, while cutting the sample according to mass gives rise to a composite structure, perhaps implying a dependence of the cluster disruption process on mass. In the context of this framework, we examine different models of cluster disruption from both theoretical and observational standpoints.

Evolutionary synthesis models are the prime method to construct models of stellar populations, and to derive physical parameters from observations. One of the assumptions for such models so far has been the time-independence of the stellar mass function. However, dynamical simulations of star clusters in tidal fields have shown the mass function to change due to the preferential removal of low-mass stars from clusters. Here we combine the results from dynamical simulations of star clusters in tidal fields with our evolutionary synthesis code GALEV to extend the models by a new dimension: the total cluster disruption time. We reanalyse the mass function evolution found in N-body simulations of star clusters in tidal fields, parametrise it as a function of age and total cluster disruption time and use this parametrisation to compute GALEV models as a function of age, metallicity and the total cluster disruption time. We study the impact of cluster dissolution on the colour (generally, they become redder) and magnitude (they become fainter) evolution of star clusters, their mass-to-light ratios (off by a factor of ~2 -- 4 from standard predictions), and quantify the effect on the cluster age determination from integrated photometry (in most cases, clusters appear to be older than they are, between 20 and 200%). By comparing our model results with observed M/L ratios for old compact objects in the mass range 10^4.5 -- 10^8 Msun, we find a strong discrepancy for objects more massive than 10^7 Msun (higher M/L). This could be either caused by differences in the underlying stellar mass function or be an indication for the presence of dark matter in these objects. Less massive objects are well represented by the models. The models for a range of total cluster disruption times are available online. (shortened)

Differences between the inferred star formation histories (SFHs) of star clusters and field stars seem to suggest distinct star formation processes for the two. The Small Magellanic Cloud (SMC) is an example of a galaxy where such a discrepancy is observed. We model the observed age distributions of the SMC clusters and field stars using a new population synthesis code, SPACE, that includes stellar evolution, infant mortality and cluster dissolution. We find that the two observed age distributions can be explained by a single SFH, thus eliminating the need to assume two separate mechanisms for star formation.

Aims. We study the connection between spatially resolved star formation and young star clusters across the disc of M51. Methods. We combine star cluster data based on B, V, and I-band Hubble Space Telescope ACS imaging, together with new WFPC2 U-band photometry to derive ages, masses, and extinctions of 1580 resolved star clusters using SSP models. This data is combined with data on the spatially resolved star formation rates and gas surface densities, as well as Halpha and 20cm radio-continuum (RC) emission, which allows us to study the spatial correlations between star formation and star clusters. Two-point autocorrelation functions are used to study the clustering of star clusters as a function of spatial scale and age. Results. We find that the clustering of star clusters among themselves decreases both with spatial scale and age, consistent with hierarchical star formation. The slope of the autocorrelation functions are consistent with projected fractal dimensions in the range of 1.2-1.6, which is similar to other galaxies, therefore suggesting that the fractal dimension of hierarchical star formation is universal. Both star and cluster formation peak at a galactocentric radius of 2.5 and 5 kpc, which we tentatively attribute to the presence of the 4:1 resonance and the co-rotation radius. The positions of the youngest (<10 Myr) star clusters show the strongest correlation with the spiral arms, Halpha, and the RC emission, and these correlations decrease with age. The azimuthal distribution of clusters in terms of kinematic age away from the spiral arms indicates that the majority of the clusters formed 5-20 Myr before their parental gas cloud reached the centre of the spiral arm.

To obtain an accurate description of broad-band photometric star cluster evolution, certain effects should be accounted for, such as the preferential loss of low-mass stars and the retain of stellar remnants. Moreover, the IMF and metallicity affect photometry as well. Due to their advanced evolutionary state, globular clusters (GCs) are interesting test cases for cluster models describing these effects. In this paper we describe cluster models in which the above effects are included. The photometric evolution of clusters is predicted, and the results are applied to Galactic GCs. The presented cluster models represent an analytical description of the evolution of the underlying stellar mass function due to stellar evolution and dynamical cluster dissolution. Stellar remnants are included by using initial-remnant mass relations, while cluster photometry is computed from the Padova 1999 isochrones. Our study shows that the preferential loss of low-mass stars, the retain of stellar remnants, as well as IMF and metallicity variations strongly influence the magnitude, colour and mass-to-light ratio evolution of clusters. The different effects can be clearly separated with our models. The models are applied to the Galactic GC population, of which the magnitude, colour and mass-to-light ratio ranges are well reproduced. We also show that the mass-to-light ratios of clusters of similar ages and metallicities cannot be assumed to be constant for all cluster luminosities. Instead, mass-to-light ratio increases with cluster luminosity and mass. These models underline the importance of detailed cluster models when considering cluster photometry. With the analytic framework provided in this paper, observed cluster properties can be interpreted in a more complete perspective. (Abridged)

Whether or not there exists a physical upper mass limit for star clusters is still unclear. HST/ACS data for the rich cluster population in the interacting galaxy M51 enables us to investigate this in more detail. We investigate whether the cluster luminosity function (LF) in M51 shows evidence for an upper limit to the mass function. The variations of the LF parameters with position on the disk are addressed. We determine the cluster LF for all clusters in M51 falling within our selection criteria, as well as for several subsets of the sample. In that way we can determine the properties of the cluster population as a function of galactocentric distance and background intensity. By comparing observed and simulated LFs we can constrain the underlying cluster initial mass function and/or cluster disruption parameters. A physical upper mass limit for star clusters will appear as a bend dividing two power law parts in the LF, if the cluster sample is large enough to sample the full range of cluster masses. The location of the bend in the LF is indicative of the value of the upper mass limit. The slopes of the power laws are an interplay between upper masses, disruption and fading. The LF of the cluster population of M51 is better described by a double power law than by a single power law. We show that the cluster initial mass function is likely to be truncated at the high mass end. We conclude from the variation of the LF parameters with galactocentric distance that both the upper mass limit and the cluster disruption parameters are likely to be a function of position in the galactic disk. At higher galactocentric distances the maximum mass is lower, cluster disruption slower, or both.

The analysis of the age distributions of star cluster samples of different galaxies has resulted in two very different empirical models for the dissolution of star clusters: the Baltimore model and the Utrecht model. I describe these two models and their differences. The Baltimore model implies that the dissolution of star clusters is mass independent and that about 90% of the clusters are destroyed each age dex, up to an age of about a Gyr, after which point mass-dependent dissolution from two-body relaxation becomes the dominant mechanism. In the Utrecht model, cluster dissolution occurs in three stages: (i) mass-independent infant mortality due to the expulsion of gas up to about 10 Myr; (ii) a phase of slow dynamical evolution with strong evolutionary fading of the clusters lasting up to about a Gyr; and (iii) a phase dominated by mass dependent-dissolution, as predicted by dynamical models. I describe the cluster age distributions for mass-limited and magnitude-limited cluster samples for both models. I refrain from judging the correctness of these models.

Aims: We study a peculiar object with a projected position close to the nucleus of M51. It is unusually large for a star cluster in M51 and we therefore investigate the three most likely options to explain this object: (a) a background galaxy, (b) a cluster in the disk of M51 and (c) a cluster in M51, but in front of the disk. Methods: We use HST/ACS and HST/NICMOS broad-band photometry to study the properties of this object. Assuming the object is a star cluster, we fit the metallicity, age, mass and extinction using simple stellar population models. Assuming the object is a background galaxy, we estimate the extinction from the colour of the background around the object. We study the structural parameters of the object by fitting the spatial profile with analytical models. Results: We find de-reddened colours of the object which are bluer than expected for a typical elliptical galaxy, and the central surface brightness is brighter than the typical surface brightness of a disc galaxy. It is therefore not likely that the object is a background galaxy. Assuming the object is a star cluster in the disc of M51, we estimate an age and mass of 0.7 Gyr and 2.2 x 10^5 \msun, respectively (with the extinction fixed to E(B-V) = 0.2). Considering the large size of the object, we argue that in this scenario we observe the cluster just prior to final dissolution. If we fit for the extinction as a free parameter, a younger age is allowed and the object is not close to final dissolution. Alternatively, the object could be a star cluster in M51, but in front of the disc, with an age of 1.4 Gyr and mass M = 1.7 x 10^5 \msun. Its effective radius is between ~12-25 pc. This makes the object a "fuzzy star cluster", raising the issue of how an object of this age would end up outside the disc.

The dissolution time (t_dis) of clusters in a tidal field does not scale with the ``classical'' expression for the relaxation time. First, the scaling with N, and hence cluster mass, is shallower due to the finite escape time of stars. Secondly, the cluster half-mass radius is of little importance. This is due to a balance between the relative tidal field strength and internal relaxation, which have an opposite effect on t_dis, but of similar magnitude. When external perturbations, such as encounters with giant molecular clouds (GMC) are important, t_dis for an individual cluster depends strongly on radius. The mean dissolution time for a population of clusters, however, scales in the same way with mass as for the tidal field, due to the weak dependence of radius on mass. The environmental parameters that determine t_dis are the tidal field strength and the density of molecular gas. We compare the empirically derived t_dis of clusters in six galaxies to theoretical predictions and argue that encounters with GMCs are the dominant destruction mechanism. Finally, we discuss a number of pitfalls in the derivations of t_dis from observations, such as incompleteness, with the cluster system of the SMC as particular example.

Two aspects of our recent N-body studies of star clusters are presented: (1) What impact does mass segregation and selective mass loss have on integrated photometry? (2) How well compare results from N-body simulations using NBODY4 and STARLAB/KIRA?

We report on neon abundances derived from {\it Spitzer} high resolution spectral data of eight Wolf-Rayet (WR) stars using the forbidden line of [\ion{Ne}{3}] 15.56 microns. Our targets include four WN stars of subtypes 4--7, and four WC stars of subtypes 4--7. We derive ion fraction abundances $\gamma$ of Ne^{2+} for the winds of each star. The ion fraction abundance is a product of the ionization fraction $Q_{\rm i}$ in stage i and the abundance by number ${\cal A}_E$ of element E relative to all nuclei. Values generally consistent with solar are obtained for the WN stars, and values in excess of solar are obtained for the WC stars.

We re-analyze the age distribution (dN/dt) of star clusters in the Small Magellanic Cloud (SMC) using age determinations based on the Magellanic Cloud Photometric Survey. For ages younger than 3x10^9 yr the dN/dt distribution can be approximated by a power-law distribution, dN/dt propto t^-beta, with -beta=-0.70+/-0.05 or -beta=-0.84+/-0.04, depending on the model used to derive the ages. Predictions for a cluster population without dissolution limited by a V-band detection result in a power-law dN/dt distribution with an index of ~-0.7. This is because the limiting cluster mass increases with age, due to evolutionary fading of clusters, reducing the number of observed clusters at old ages. When a mass cut well above the limiting cluster mass is applied, the dN/dt distribution is flat up to 1 Gyr. We conclude that cluster dissolution is of small importance in shaping the dN/dt distribution and incompleteness causes dN/dt to decline. The reason that no (mass independent) infant mortality of star clusters in the first ~10-20 Myr is found is explained by a detection bias towards clusters without nebular emission, i.e. cluster that have survived the infant mortality phase. The reason we find no evidence for tidal (mass dependent) cluster dissolution in the first Gyr is explained by the weak tidal field of the SMC. Our results are in sharp contrast to the interpretation of Chandar et al. (2006), who interpret the declining dN/dt distribution as rapid cluster dissolution. This is due to their erroneous assumption that the sample is limited by cluster mass, rather than luminosity.

We use HST/ACS observations of the spiral galaxy M51 in F435W, F555W and F814W to select a large sample of star clusters with accurate effective radius measurements in an area covering the complete disc of M51. We present the dataset and study the radius distribution and relations between radius, colour, arm/interarm region, galactocentric distance, mass and age. We select a sample of 7698 (F435W), 6846 (F555W) and 5024 (F814W) slightly resolved clusters and derive their effective radii by fitting the spatial profiles with analytical models convolved with the point spread function. The radii of 1284 clusters are studied in detail. We find cluster radii between 0.5 and ~10 pc, and one exceptionally large cluster candidate with a radius of 21.6 pc. The median radius is 2.1 pc. We find 70 clusters in our sample which have colours consistent with being old GC candidates and we find 6 new "faint fuzzy" clusters in, or projected onto, the disc of M51. The radius distribution can not be fitted with a power law, but a log-normal distribution provides a reasonable fit to the data. This indicates that shortly after the formation of the clusters from a fractal gas, their radii have changed in a non-uniform way. We find an increase in radius with colour as well as a higher fraction of redder clusters in the interarm regions, suggesting that clusters in spiral arms are more compact. We find a correlation between radius and galactocentric distance which is considerably weaker than the observed correlation for old Milky Way GCs. We find weak relations between cluster luminosity and radius, but we do not observe a correlation between cluster mass and radius.

The development and progress of the studies of winds and mass loss from hot stars, from about 1965 up to now, is discussed in a personal historical perspective. The present state of knowledge about stellar winds, based on papers presented at this workshop, is described. About ten years ago the mechanisms of the winds were reasonably well understood, the mass loss rates were known, and the predictions of stellar evolution theory with mass loss agreed with observations. However, recent studies especially those based on FUSE observations, have resulted in a significant reduction of the mass loss rates, that disagrees with predictions from radiation driven wind models. The situation is discussed and future studies that can clarify the situation are suggested. I also discuss what is known about the dissolution of star clusters in different environments. The dissolution time can be derived from the mass and age distributions of cluster samples. The resulting dissolution times of clusters in the solar neighborhood (SN) and in interacting galaxies are shorter than predicted by two-body relaxation of clusters in a tidal field. Encounters with giant molecular clouds can explain the fate of clusters in the SN and are the most likely cause of the short lifetime of clusters in interacting galaxies.

In 1958 Jan Oort remarked that the lack of old clusters in the solar neighborhood (SN) implies that clusters are destroyed on a timescale of less than a Gyr. This is much shorter than the predicted dissolution time of clusters due to stellar evolution and two-body relaxation in the tidal field of the Galaxy. So, other (external) effects must play a dominant role in the destruction of star clusters in the solar neighborhood. We recalculated the survival time of initially bound star clusters in the solar neighborhood taking into account: (1) stellar evolution, (2) tidal stripping, (3) perturbations by spiral arms and (4) encounters with giant molecular clouds (GMCs). We find that encounters with GMCs are the most damaging to clusters. The resulting predicted dissolution time of these combined effects, t_dis=1.7 (Mi/10^4 M_sun)^0.67 Gyr for clusters in the mass range of 10^2 < M < 10^5 M_sun, is very similar to the disruption time of t_dis=1.3+/-0.5 (M/10^4 M_sun)^0.62 Gyr that was derived empirically from a mass limited sample of clusters in the solar neighborhood within 600 pc. The predicted shape of the age distribution of clusters agrees very well with the observed one. The comparison between observations and theory implies a surface star formation rate (SFR) near the sun of 3.5x10^-10 M_sun yr^-1 pc^-2 for stars in bound clusters with an initial mass in the range of 10^2 to 3x10^4 M_sun. This can be compared to a total SFR of 7-10x10^-10 M_sun yr^-1 pc^-2 derived from embedded clusters or 3-7x10^-9 M_sun yr^-1 pc^-2 derived from field stars. This implies an infant mortality rate of clusters in the solar neighborhood between 50% and 95%, in agreement with the results of a study of embedded clusters.

We study the luminosity function (LF) of the star clusters in M51. Comparing the observed LF with the LF resulting from artificial cluster populations suggests that there exists an upper mass limit for clusters and that this limit and/or the cluster disruption varies with galactocentric distance.

We predict the survival time of initially bound star clusters in the solar neighbourhood taking into account: (1) stellar evolution, (2) tidal stripping, (3) shocking by spiral arms and (4) encounters with giant molecular clouds. We find that the predicted dissolution time is t_dis= 1.7 (M_i/10^4 M_sun)^0.67 Gyr for clusters in the mass range of 10^2 < M_i < 10^5 M_sun, where M_i is the initial mass of the cluster.. The resulting predicted shape of the logarithmic age distribution agrees very well with the empirical one, derived from a complete sample of clusters in the solar neighbourhood within 600 pc. The required scaling factor implies a star formation rate of 400 M_sun/Myr within 600 pc from the Sun or a surface formation rate of 3.5 10^-10 M_sun/(yr pc^2) for stars in bound clusters with an initial mass in the range of 10^2 to 3 10^4 M_sun.

We investigate encounters between giant molecular clouds (GMCs) and star clusters. We propose a single expression for the energy gain of a cluster due to an encounter with a GMC, valid for all encounter distances and GMC properties. This relation is verified with N-body simulations of cluster-GMC encounters and excellent agreement is found. The fractional mass loss from the cluster is 0.25 times the fractional energy gain. This is because 75% of the injected energy goes to the velocities of escaping stars, that are higher than the escape velocity. We derive an expression for the cluster disruption time (t_dis) based on the mass loss from the simulations, taking into account the effect of gravitational focusing by the GMC. The disruption time depends on the cluster mass (M_c) and half-mass radius (r_h) as t_dis=2.0 S (M_c/10^4 M_sun)(3.75 pc/r_h)^3 Gyr, with S=1 for the solar neighbourhood and inversely proportional to the GMC density. The observed shallow relation between cluster radius and mass gives t_dis a power-law dependence on the mass with index 0.7, similar to that found from observations and from simulations of clusters dissolving in tidal fields (0.62). The constant of 2.0 Gyr is about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving under the combined effect of galactic tidal field and stellar evolution. It is somewhat higher than the observationally determined value of 1.3 Gyr. It suggests, however, that the combined effect of tidal field and encounters with GMCs can explain the lack of old open clusters in the solar neighbourhood. GMC encounters can also explain the (very) short disruption time that was observed for star clusters in the central region of M51, since there rho_n is an order of magnitude higher than in the solar neighbourhood.

We present an analysis of optical lightcurves of Small Magellanic Cloud Be-type stars. Observations show that (1) optical excess flux is correlated with near-IR excess flux indicating a similar mechanism and (2) the lightcurves can trace out ``loops'' in a colour-magnitude diagram. A simple model for the time dependent bound-free and free-free emission produced by an outflowing circumstellar disk gives reasonable fits to the observations.

(Abridged) Some 240 blue stars in the Small Magellanic Cloud are investigated on their fantastic irregular continuum variability. We report here two results regarding these stars. First, their optical flux excess is correlated to their near-IR flux excess, as determined from optical EROSII light curves and 2MASS measurements. Second, the relation between optical colour and magnitude is observed to be bi-valued in 40% of the cases, resulting in a "loop" when the light curve is presented in a colour-magnitude diagram. We argue that optical variability for a large fraction of the variable stars is due to variations in the amount of bound-free and free-free radiation. We do simple model calculations that allow us to interpret the observed colour-magnitude variability as due to an outflowing ionized circumstellar disk. The mass loss of the central star is variable, i.e. on or off. Once the star stops losing mass, the disk evolves naturally into a ring. The observed bi-valued colour-magnitude relation is the transition of a partially optically thick to a fully optically thin disk. Significantly, the loop is traversed clockwise by outflowing matter, but anti-clockwise by infalling matter. It is observed that the material is generally outflowing, but few cases of inflow are also observed.

We exploit the superb resolution of the new HST/ACS mosaic image of M51 to select a large sample of young (< 1 Gyr) star clusters in the spiral disk, based on their sizes. The image covers the entire spiral disk in B, V, I and H_alpha, at a resolution of 2 pc per pixel. The surface density distribution of 4357 resolved clusters shows that the clusters are more correlated with clouds than with stars, and we find a hint of enhanced cluster formation at the corotation radius. The radius distribution of a sample of 769 clusters with more accurate radii suggests that young star clusters have a preferred effective radius of ~3 pc, which is similar to the preferred radius of the much older GCs. However, in contrast to the GCs, the young clusters in M51 do not show a relation between radius and galactocentric distance. This means that the clusters did not form in tidal equilibrium with their host galaxy, nor that their radius is related to the ambient pressure.

We calculated the broad-band photometric evolution of unresolved star clusters, including the preferential loss of low-mass stars due to mass segregation. The stellar mass function of a cluster evolves due to three effects: (a) the evolution of massive stars; (b) early tidal effects reduce the mass function independently of the stellar mass; (c) after mass segregation has completed, tidal effects preferentially remove the lowest-mass stars from the cluster. Results: (1) During the first ~40% of the lifetime of a cluster the cluster simply gets fainter due to the loss of stars by tidal effects. (2) Between ~40 and ~80% of its lifetime the cluster gets bluer due to the loss of low-mass stars. This will result in an underestimate of the age of clusters if standard cluster evolution models are used (0.15 -- 0.5 dex). (3) After ~80% of the total lifetime of a cluster it will rapidly get redder. This is because stars at the low-mass end of the main sequence, which are preferentially lost, are bluer than the AGB stars that dominate the light at long wavelengths, resulting in an age overestimate. (4) Clusters with mass segregation and the preferential loss of low-mass stars evolve along almost the same tracks in colour-colour diagrams as clusters without mass segregation. Therefore it will be difficult to distinguish this effect from that due to the cluster age for unresolved clusters, unless the total lifetime of the clusters can be estimated. (5) The changes in the colour evolution of unresolved clusters due to the preferential loss of low-mass stars will affect the determination of the SFHs. (6) The preferential loss of low-mass stars might explain the presence of old (~13 Gyr) clusters in NGC 4365 which are photometrically disguised as intermediate-age clusters (2 - 5 Gyr). [Abridged]

We present the luminosity function (LF) of star clusters in M51 based on HST/ACS observations taken as part of the Hubble Heritage project. The clusters are selected based on their size and with the resulting 5990 clusters we present one of the largest cluster samples of a single galaxy. We find that the LF can be approximated with a double power-law distribution with a break around M_V = -8.9. On the bright side the index of the power-law distribution is steeper (a = 2.75) than on the faint-side (a = 1.93), similar to what was found earlier for the ``Antennae'' galaxies. The location of the bend, however, occurs about 1.6 mag fainter in M51. We confront the observed LF with the model for the evolution of integrated properties of cluster populations of Gieles et al., which predicts that a truncated cluster initial mass function would result in a bend in, and a double power-law behaviour of, the integrated LF. The combination of the large field-of view and the high star cluster formation rate of M51 make it possible to detect such a bend in the LF. Hence, we conclude that there exists a fundamental upper limit to the mass of star clusters in M51. Assuming a power-law cluster initial mass function with exponentional cut-off of the form NdM ~ M^-b * exp(-M/M_C)dM, we find that M_C = 10^5 M_sun. A direct comparison with the LF of the ``Antennae'' suggests that there M_C = 4*10^5 M_sun.

We present a study of the optical spectrum of the fascinating B[e] star Hen 2-90 based on new high-resolution observations taken with FEROS at the ESO 1.52m telescope in La Silla (Chile). The recent HST image of Hen 2-90 (Sahai et al. 2002) reveals a bipolar, highly ionized region, a neutral disk-like structure seen almost perfectly edge-on, and an intermediate region of moderate ionization. The slits of our observations cover the same innermost region of Hen 2-90 as the HST image, which allows us to combine the observations. Our spectra contain a huge amount of permitted and forbidden emission lines of atoms in different stages of ionization. In addition, the line wings deliver velocity information of the emitting region. We find correlation between the different ionization states of the elements and the velocities derived from the line profiles: the highly ionized atoms have the highest outflow velocity, while the neutral lines have the lowest. When combining the velocity information with the HST image of Hen 2-90, it seems that a non-spherical stellar wind model is a good option to explain the ionization and spatial distribution of the circumstellar material. Our modeling of the forbidden emission lines results in strong evidence for Hen 2-90 being a compact planetary nebulae that has undergone a superwind phase of high, non-spherical mass loss, most probably triggered by a central star that was rotating with about 80% of the critical velocity. We find a total mass loss rate during this superwind phase on the order of 3 x 10^(-5) M_sun/yr.

We investigate the ionization structure in the non-spherical winds and disks of B[e] stars. Especially the luminous B[e] supergiants seem to have outflowing disks which are neutral in hydrogen already close to the stellar surface. The existence of neutral material so close to the central star is surprising and needs to be investigated in detail. We perform our model calculations mainly in the equatorial plane, trying to find a plausible scenario that leads to recombination in the vicinity of the hot stars. Two different approaches are presented that both result in a hydrogen neutral equatorial region. We especially focus on the influence of stellar rotation which is known to play a significant role in shaping non-spherical winds. We show that a rotating star can have a neutral equatorial wind, even without the need of a density enhancement due to bi-stability or wind compression, but simply due to graviational darkening in combination with the known high mass-loss rates of the B[e] stars and especially the B[e] supergiants.

We explore whether we can constrain the shape of the INITIAL mass distribution of the star cluster population in M82's ~1 Gyr-old post-starburst region "B", in which the present-day cluster mass function (CMF) is closely approximated by a log-normal distribution. We conclude that the M82 B initial CMF must have had a mean mass very close to that of the "equilibrium" CMF of Vesperini (1998). Consequently, if the presently observed M82 B CMF has remained approximately constant since its formation, as predicted, then the INITIAL CMF must have been characterized by a mean mass that was only slightly larger than the present mean mass. From our detailed analysis of the expected evolution of CMFs, we conclude that our observations of the M82 B CMF are inconsistent with a scenario in which the 1 Gyr-old cluster population originated from an initial power-law mass distribution. Our conclusion is supported by arguments related to the initial density in M82 B, which would have been unphysically high if the present cluster population were the remains of an initial power-law distribution.

We report on a study of young star cluster complexes in the spiral galaxy M51. Recent studies have confirmed that star clusters do not form in isolation, but instead tend to form in larger groupings or complexes. We use {\it HST} broad and narrow band images (from both {\it WFPC2} and {\it ACS}), along with {\it BIMA}-CO observations to study the properties and investigate the origin of the e complexes. We find that the complexes are all young ($< 10$ Myr), have sizes between $\sim$85 and $\sim$240 pc, and have masses between 3-30 $\times 10^{4} M_{\odot}$. Unlike that found for isolated young star clusters, we find a strong correlation between the complex mass and radius, namely $M\propto R^{2.33 \pm 0.19}$. This is similar to that found for giant molecular clouds (GMCs). By comparing the mass-radius relation of GMCs in M51 to that of the complexes we can estimate the star formation efficiency within the complexes, although this value is heavily dependent on the assumed CO-to-H$_2$ conversion factor. The complexes studied here have the same surface density distribution as individual young star clusters and GMCs. If star formation within the complexes is proportional to the gas density at that point, then the shared mass-radius relation of GMCs and complexes is a natural consequence of their shared density profiles. We briefly discuss possibilities for the lack of a mass-radius relation for young star clusters. We note that many of the complexes show evidence of merging of star clusters in their centres, suggesting that larger star clusters can be produced through the build up of smaller clusters.

In this work we concentrate on the evolution of the cluster population of the interacting galaxy M51 (NGC 5194), namely the timescale of cluster disruption and possible variations in the cluster formation rate. We present a method to compare observed age vs. mass number density diagrams with predicted populations including various physical input parameters like the cluster initial mass function, cluster disruption, cluster formation rate and star bursts. If we assume that the cluster formation rate increases at the moments of the encounters with NGC 5195, we find an increase in the cluster formation rate of a factor of 3, combined with a disruption timescale which is slightly higher then when assuming a constant formation rate (t_4 = 200 Myr vs. 100 Myr). The measured cluster disruption time is a factor of 5 shorter than expected on theoretical grounds. This implies that the disk of M51 is not a preferred location for survival of young globular clusters, since even clusters with masses of the order of 10^6 M_sun will be destroyed within a few Gyr.

We present a simple analytical description of the disruption of star clusters in a tidal field, which agrees excellently with detailed N-body simulations. The analytic expression can be used to predict the mass and age histograms of surviving clusters for any cluster initial mass function and any cluster formation history. The method is applied to open clusters in the solar neighbourhood, based on the new cluster sample of Kharchenko et al. From a comparison between the observed and predicted age distributions in the age range between 10 Myr to 3 Gyr we find the following results: (1) The disruption time of a 10^4 M_sun cluster in the solar neighbourhood is about 1.3+/-0.5 Gyr. This is a factor 5 shorter than derived from N-body simulations of clusters in the tidal field of the galaxy. (2) The present starformation rate in bound clusters within 600 pc from the Sun is 5.9+/-0.8 * 10^2 M_sun / Myr, which corresponds to a surface star formation rate in bound clusters of 5.2+/-0.7 10^(-10) M_sun/yr/pc^2. (3) The age distribution of open clusters shows a bump between 0.26 and 0.6 Gyr when the cluster formation rate was 2.5 times higher than before and after. (4) The present star formation rate in bound clusters is half as small as that derived from the study of embedded clusters. The difference suggests that half of the clusters in the solar neighbourhood become unbound within 10 Myr. (5) The most massive clusters within 600 pc had an initial mass of 3*10^4 M_sun. This is in agreement with the statistically expected value based on a cluster initial mass function with a slope of -2, even if the physical upper mass limit is as high as 10^6 M_sun.

(abridged) We study the optical spectrum of the exciting B[e] star Hen 2-90 based on new high-resolution observations that cover the innermost 2". Our investigation is splitted in two parts, a qualitative study of the presence of the numerous emission lines and the classification of their line profiles which indicate a circumstellar environment of high complexity, and a quantitative analysis of numerous forbidden lines, e.g. [OI], [OII], [OIII], [SII], [SIII], [ArIII], [ClII], [ClIII] and [NII]. We find a correlation between the different ionization states of the elements and the velocities derived from the line profiles: the highly ionized atoms have the highest outflow velocity while the neutral lines have the lowest outflow velocity. The recent HST image of Hen 2-90 reveals a bipolar, highly ionized region, a neutral disk-like structure and an intermediate region of moderate ionization. It seems that a non-spherical stellar wind model is a good option to explain the ionization and spatial distribution of the circumstellar material. We modelled the forbidden lines under the assumption of a non-spherically symmetric wind based on the HST image. We find that in order to fit the observed line luminosities, the mass flux, surface temperature, and terminal wind velocities need to be latitude dependent, which might be explained in terms of a rapidly rotating central star. A rotation speed of 75-80 % of the critical velocity has been derived. The total mass loss rate of the star was determined to be of order 3 10^{-5} M_sun/yr. Such a wind scenario and the fact that compared to solar abundances C, O, and N seem to be underabundant while S, Ar and Cl have solar abundances, might be explained in terms of a rapidly rotating post-AGB star.

(ABRIDGED) We assess the systematic uncertainties in (young) cluster age, mass, and - to a lesser extent - extinction and metallicity determinations, based on broad-band imaging observations with the Hubble Space Telescope. Our aim here is to intercompare the results obtained using a variety of commonly used modelling techniques, specifically with respect to our own extensively tested multi-dimensional approach. Any significant differences among the resulting parameters are due to the details of the various, independently developed modelling techniques used, rather than to the stellar population models themselves. Despite the model uncertainties and the selection effects inherent to most methods used, we find that the PEAKS in the relative age and mass distributions of a given young (<= 10^9 yr) cluster system can be derived relatively robustly and consistently, to accuracies of sigma_t = Delta <log(Age / yr)> <= 0.35 and sigma_M = Delta <log(M_cl / Msun)> <= 0.14, respectively, assuming Gaussian distributions in cluster ages and masses for reasons of simplicity. The peaks in the relative mass distributions can be obtained with a higher degree of confidence than those in the relative age distributions, as exemplified by the smaller spread among the peak values of the mass distributions derived. This implies that mass determinations are mostly insensitive to the approach adopted. We reiterate that as extensive a wavelength coverage as possible is required to obtain robust and internally consistent age and mass estimates for the individual objects, with reasonable uncertainties. Finally, we conclude that the ACTUAL filter systems used for the observations should be used for constructing model colours, instead of using conversion equations, to achieve more accurate derivations of ages and masses.

We use archival {\it Hubble Space Telescope} observations of broad-band images from the ultraviolet (F255W-filter) through the near infrared (NICMOS F160W-filter) to study the star cluster population of the interacting spiral galaxy M51. We obtain age, mass, extinction, and effective radius estimates for 1152 star clusters in a region of $\sim 7.3 \times 8.1$ kpc centered on the nucleus and extending into the outer spiral arms. In this paper we present the data set and exploit it to determine the age distribution and relationships among the fundamental parameters (i.e. age, mass, effective radius). Using this dataset we find: {\it i}) that the cluster formation rate seems to have had a large increase $\sim$ 50-70 Myr ago, which is coincident with the suggested {\it second passage} of its companion, NGC 5195, {\it ii}) a large number of extremely young ($<$ 10 Myr) star clusters, which we interpret as a population of unbound clusters of which a large majority will disrupt within the next $\sim$10 Myr, and {\it iii)} that the distribution of cluster sizescan be well approximated by a power-law with exponent, $ -\eta = -2.2 \pm 0.2$, which is very similar to that of Galactic globular clusters, indicating that cluster disruption is largely independent of cluster radius. In addition, we have used this dataset to search for correlations among the derived parameters. In particular, we do not find any strong trends between the age and mass, mass and effective radius, nor between the galactocentric distance and effective radius. There is, however, a strong correlation between the age of a cluster and its extinction, with younger clusters being more heavily reddened than older clusters.

We compare the absolute visual magnitude of the majority of bright O stars in the sky as predicted from their spectral type with the absolute magnitude calculated from their apparent magnitude and the Hipparcos parallax. We find that many stars appear to be much fainter than expected, up to five magnitudes. We find no evidence for a correlation between magnitude differences and the stellar rotational velocity as suggested for OB stars by Lamers et al. (1997), whose small sample of stars is partly included in ours. Instead, by means of a simulation we show how these differences arise naturally from the large distances at which O stars are located, and the level of precision of the parallax measurements achieved by Hipparcos. Straightforwardly deriving a distance from the Hipparcos parallax yields reliable results for one or two O stars only. We discuss several types of bias reported in the literature in connection with parallax samples (Lutz-Kelker, Malmquist) and investigate how they affect the O star sample. In addition, we test three absolute magnitude calibrations from the literature (Schmidt-Kaler et al. 1982; Howarth & Prinja 1989; Vacca et al. 1996) and find that they are consistent with the Hipparcos measurements. Although O stars conform nicely to the simulation, we notice that some B stars in the sample of Lamers et al. (1997) have a magnitude difference larger than expected.

The observed average lifetime of the population of star clusters in the Solar Neighbourhood, the Small Magellanic Cloud and in selected regions of M51 and M33 is compared with simple theoretical predictions and with the results of N-body simulations. The empirically derived lifetimes (or disruption times) of star clusters depend on their initial mass as t_dis ~ Mcl^0.60 in all four galaxies. N-body simulations have shown that the predicted disruption time of clusters in a tidal field scales as t_dis^pred ~ t_rh^0.75 t_cr^0.25, where t_rh is the initial half-mass relaxation time and t_cr is the crossing time for a cluster in equilibrium. We show that this can be approximated accurately by t_dis^pred ~ M_cl^0.62 for clusters in the mass range of about 10^3 to 10^6 M_sun, in excellent agreement with the observations. Observations of clusters in different extragalactic environments show that t_dis also depends on the ambient density in the galaxies where the clusters reside. Linear analysis predicts that the disruption time will depend on the ambient density of the cluster environment as t_dis ~ rho_amb^-0.5. This relation is consistent with N-body simulations.

We have carried out a detailed analysis of the IUE archival high resolution spectra of the classical nova V1974 Cyg 1992. In addition to the P Cygni and emission lines, two shortward shifted absorption systems are present. Evidence is given that these absorptions originate in two separate expanding shells, outside the wind layers where the emission lines are formed. The outer main shell, containing most of the matter ejected at outburst, produces the so-called ``principal absorption line system'', and the inner faster moving second shell produces the so-called ``diffuse--enhanced absorption line system''. The outflow velocity of the two shells increases exponentially with time reaching a value of about 1750 km/s and 2900 km/s, respectively. We suggest that the acceleration of the shells is the result of increasing line-radiation pressure due to the UV-brightening of the star as the effective radius decreases. Around day 60 the second shell has overtaken the slower moving principal system shell, and merged with it. This explains: the sudden disappearance of the diffuse line system near that date, the upward jump of 240 km/s in velocity of the principal system and the first detection of hard X-ray emission on day 63. This velocity jump indicates that the main shell is about 4 times more massive than the second shell. The deceleration suffered by the diffuse-enhanced system after the shock provides a shock temperature of about 1.6 KeV, in fairly good agreement with the temperature of the observed hard X-ray emission.

Observations and N-body simulations both support a simple relation for the disruption time of a cluster as a function of its mass of the form: t_dis = t_4 * (M/10^4 Msun)^gamma. The scaling factor t_4 seems to depend strongly on the environment. Predictions and observations show that gamma ~ 0.64 +/- 0.06. Assuming that t_dis ~ M^0.64 is caused by evaporation and shocking implies a relation between the radius and the mass of a cluster of the form: r_h ~ M^0.07, which has been observed in a few galaxies. The suggested relation for the disruption time implies that the lower mass end of the cluster initial mass function will be disrupted faster than the higher mass end, which is needed to evolve a young power law shaped mass function into the log-normal mass function of old (globular) clusters.

We present the age and mass distribution of star clusters in M51. The structural parameters are found by fitting cluster evolution models to the spectral energy distribution consisting of 8 HST-WFPC2 pass bands. There is evidence for a burst of cluster formation at the moment of the second encounter with the companion NGC5195 (50-100 Myr ago) and a hint for an earlier burst (400-500 Myr ago). The cluster IMF has a power law slope of -2.1. The disruption time of clusters is extremely short (< 100 Myr for a 10^4 Msun cluster).

From optical spectropolarimetry of Cygnus A we conclude that the scattering medium in the ionization cones in Cygnus A is moving outward at a speed of 170+-34 km/s, and that the required momentum can be supplied by the radiation pressure of an average quasar. Such a process could produce a structure resembling the observed ionization cones, which are thought to result from shadowing by a circumnuclear dust torus. We detect a polarized red wing in the [O III] emission lines arising from the central kiloparsec of Cygnus A. This wing is consistent with line emission created close to the boundary of the broad-line region.

We report the discovery of 5 new Herbig Ae/Be candidate stars in the Small Magellanic Cloud in addition to the 2 reported in Beaulieu et al. (2001). We discuss these 7 HAeBe candidate stars in terms of (1) their irregular photometric variability, (2) their near infrared emission, (3) their Halpha emission and (4) their spectral type. One star has the typical photometric behaviour that is observed only among Pre-Main Sequence UX Orionis type stars. The objects are more luminous than Galactic HAeBe stars and Large Magellanic Cloud HAeBe candidates of the same spectral type. The stars were discovered in a systematic search for variable stars in a subset of the EROS2 database consisting of 115,612 stars in a field of 24x24 arcmin in the Small Magellanic Cloud. In total we discovered 504 variable stars. After classifying the different objects according to their type of variability, we concentrate on 7 blue objects with irregular photometric behaviour. We cross-identified these objects with emission line catalogues from Simbad and JHK photometry from 2MASS. The analysis is supplemented with obtained narrow and broad band imaging. We discuss their variability in terms of dust obscuration and bound-free and free-free emission. We estimate the influence of metallicity on the circumstellar dust emission from pre-main sequence stars.

We use the ages, masses and metallicities of the rich young star cluster systems in the nearby starburst galaxies NGC 3310 and NGC 6745 to derive their cluster formation histories and subsequent evolution. We further expand our analysis of the systematic uncertainties involved in the use of broad-band observations to derive these parameters by examining the effects of a priori assumptions on the individual cluster metallicities. The age (and metallicity) distributions of both the clusters in the circumnuclear ring in NGC 3310 and of those outside the ring are statistically indistinguishable, but there is a clear and significant excess of higher-mass clusters IN the ring compared to the non-ring cluster sample; it is likely that the physical conditions in the starburst ring may be conducive for the formation of higher-mass star clusters, on average, than in the relatively more quiescent environment of the main galactic disc. For the NGC 6745 cluster system we derive a median age of ~10 Myr. NGC 6745 contains a significant population of high-mass "super star clusters", with masses in the range 6.5 <= log(M_cl/M_sun) <= 8.0. This detection supports the scenario that such objects form preferentially in the extreme environments of interacting galaxies. The age of the cluster populations in both NGC 3310 and NGC 6745 is significantly lower than their respective characteristic cluster disruption time-scales. This allows us to obtain an independent estimate of the INITIAL cluster mass function slope, alpha = 2.04(+- 0.23)(+0.13)(-0.43) for NGC 3310, and 1.96(+- 0.15)(+- 0.19) for NGC 6745, respectively, for masses M_cl >= 10^5 M_sun and M_cl >= 4 x 10^5 M_sun. These mass function slopes are consistent with those of other young star cluster systems in interacting and starburst galaxies.

The non-spherically symmetric winds of B[e] supergiants are investigated. An empirical density distribution is chosen that accounts for the density concentrations and ratios derived from observations, and our model winds are assumed to contain only hydrogen and helium. We first calculate the approximate ionization radii for H and He and compare the results with the ionization fractions calculated from the more accurate ionization balance equations. We find that winds with a r^-2 density distribution turn out to reach a constant ionization fraction as long as the wind density is low, i.e. in polar direction. For the high density equatorial regions, however, we find that the winds become neutral just above the stellar surface of the hot and massive B[e] supergiants forming a disk-like neutral region. In such a disk molecules and dust can form even very near the hot central star.

The spectra of stars with the B[e] phenomenon are dominated by features that are related to physical conditions of circumstellar material around these objects and are not intrinsic to the stars. Because of this, they form a very heterogeneous group. This group contains objects with different evolutionary stages. Lamers et al. (1998) have suggested a new designation with five sub-groups, which indicate the evolutionary stage. They are: supergiants, pre-main sequence or Herbig Ae/Be, compact planetary nebulae, symbiotic and unclassified. The unclassified group has many objects that need a better study to resolve their evolutionary status. Forbidden lines can be a useful tool to solve this problem. They can give informations about chemical composition, ionization and density of the circunstellar medium and probably the evolutionary phase of these objects. We analize spectra of some galactic objects, obtained with FEROS and B&C spectrograph at 1.52 telescope in ESO (La Silla-Chile), with a special focus on the forbidden lines. We have studied the spectra of 5 B[e] stars of uncertain evolutionary stage. We find that one of them is a pre-WN star, the other four are supergiant B[e] stars.

The spectra of stars with the B[e] phenomenon are dominated by features that are related to physical conditions of circumstellar material around these objects and are not intrinsic to the stars. Previous studies have shown that emission lines present in the optical spectra of these objects are formed in an equatorial rotating disk. The analysis of high and low resolution spectra, obtained by us with 1.52 telescope in ESO for some Southern Galactic B[e] Stars, can give us information about the structure and velocity of the disk. We will describe the analysis of the unclassified B[e] star Hen 2-90.

The currently most popular models for the dynamical evolution of star clusters predict that the power-law Cluster Luminosity Functions (CLFs) of young star cluster systems will be transformed rapidly into the universal Gaussian CLFs of old Milky Way-type "globular" cluster systems. Here, we provide the first evidence for a turn-over in the intermediate-age, approximately 1 Gyr-old CLF in the center of the nearby starburst galaxy M82, which very closely matches the universal CLFs of old Milky Way-type globular cluster systems. This provides an important test of both cluster disruption theories and hierarchical galaxy formation models. It also lends strong support to the scenario that these young cluster systems may eventually evolve into old Milky Way-type globular cluster systems. M82's proximity, its shortest known cluster disruption time-scale of any galaxy, and its well-defined peak of cluster formation make it an ideal candidate to probe the evolution of its star cluster system to fainter luminosities, and thus lower masses, than has been possible for any galaxy before.

ABRIDGED: We obtain new age and mass estimates for the star clusters in M82's fossil starburst region B, based on improved fitting methods. Our new age estimates confirm the peak in the age histogram attributed to the last tidal encounter with M81; we find a peak formation epoch at slightly older ages than previously published, log(t_peak / yr) = 9.04, with a Gaussian sigma of Delta log(t_width) = 0.273. Cluster disruption has removed a large fraction of the older clusters. Adopting the expression for the cluster disruption time-scale of t_dis(M)= t_dis^4 (M/10^4 Msun)^gamma with gamma = 0.62 (Paper I), we find that the ratios between the real cluster formation rates in the pre-burst phase (log(t/yr) <= 9.4), the burst-phase (8.4 < log(t/yr) < 9.4) and the post-burst phase (log(t/yr) <= 8.4) are about 1:2:1/40. The mass distribution of the clusters formed during the burst shows a turnover at log(M_cl/Msun) ~ 5.3 which is not caused by selection effects. This distribution can be explained by cluster formation with an initial power-law mass function of slope alpha=2 up to a maximum cluster mass of M_max = 3 x 10^6 Msun, and cluster disruption with a normalisation time-scale t_dis^4 / t_burst = (3.0 +/- 0.3) x 10^{-2}. For a burst age of 1 x 10^9 yr, we find that the disruption time-scale of a cluster of 10^4 Msun is t_dis^4 ~ 3 x 10^7 years, with an uncertainty of approximately a factor of two. This is the shortest disruption time-scale known in any galaxy.

A new point of view on the dynamics of the circumstellar nebula M1-67 around the run-away Wolf-Rayet (WR) star WR 124 is presented. We found that it has been interacting with the surrounding ISM and has formed a bow shock due to its high velocity of about 180 km/s relative to the local ISM. The star is about 1.3 parsec away from the front of this bow shock. The outbursts that are responsible for the nebula are assumed to be discrete outbursts that occurred inside this bow shock. The ejecta collide with this bow shock shortly after the outburst. After the collision, they are dragged away by the pressure of the ISM, along the surface of the bow shock. The bow shock is oriented in such way that we are looking from the rear into this paraboloid, almost along the main axis. Evidence for this is given firstly by the fact that the far hemisphere is much brighter than the near hemisphere, secondly by the fact that there is hardly any emission found with radial velocities higher than the star's radial velocity, thirdly by the fact that the star looks to be in the centre of the nebula, as seen from Earth, and finally by the asymmetric overall velocity distribution of the nebula, which indicates higher radial velocities in the centre of the nebula, and lower velocities near the edges. We find evidence for at least two discrete outbursts that occurred inside this bow shock. For these outbursts, we find expansion velocities of about 150 km/s and dynamical timescales of about 8 and 20 kyr, which are typical values for LBV outbursts. We therefore conclude that M1-67 originates from several outbursts that occurred inside the bow shock around WR 124, during an LBV phase that preceded the current WR phase of the star.

If the cluster formation rate is constant and the disruption time of clusters depends on their initial mass as t_dis= t_4 x (M_i/10^4 Msun)^gamma, where M_i is the initial mass of a cluster, t_4 is the disruption time of a cluster of M_i=10^4 Msun and gamma is an exponent, then the values of t_4 and gamma can be derived in a very simple way from the age and mass histograms of large homogeneous samples of clusters with reliable age and/or mass determinations. We demonstrate the method and derive the values of t_4 and gamma from cluster samples in selected fields of four galaxies: M51, M33, SMC and the solar neighbourhood. The values of the exponent gamma are the same within their uncertainties in the four galaxies with gamma = 0.62 +- 0.06. However, the disruption time t_4 of clusters of 10^4 Msun is very different in the different locations. It is shortest in the inner region of M51 (t_4 = 40 Myr) and highest in the SMC (t_4=8 Gyr).

We present a new method to derive the cluster disruption time in selected regions of galaxies from the mass or age distribution of magnitude-limited cluster samples. If the disruption time of clusters in a region of a galaxy depends on their initial mass as t_4 x (M_cluster/10^4 M_sun)^gamma and if the cluster formation rate is constant, then the mass and age distributions of the observed clusters will each show two powerlaw relations. The values of t_4 and gamma can be derived from these relations. We used this method to derive the cluster disruption time in specific regions in four galaxies: the inner region of M51, a region of M33, the SMC and the solar neighbourhood. The values of gamma are the same in the four galaxies within the uncertainty and the mean value is gamma= 0.62 +- 0.06. However the disruption time t_4 of a cluster of 10^4 M_sun is very different in the different galaxies. The clusters in the SMC have the longest disruption time, t_4 = 8 Gyr, and the clusters at 1 to 3 kpc from the nucleus of M51 have the shortest disruption time of t_4 = 0.04 Gyr. The disruption time of clusters 1 to 5 kpc from the nucleus of M33 is t_4 = 0.13 Gyr and for clusters within 1 kpc from the Sun we find t_4 = 1.0 Gyr.

We study the cluster population in a region of 3.2x3.2 kpc^2 in the inner spiral arms of the intergacting galaxy M51, at a distance of about 1 to 3 kpc from the nucleus, based on HST--WFPC2 images taken through five broadband and two narrowband filters. We found 877 cluster candidates and we derived their ages, initial masses and extinctions by comparing their energy distribution with the Starburst99 cluster models. We describe the 3 and 2-dimensional least-square energy fitting method that was used (3DEF, 2DEF). The lack of [OIII] emission in even the youngest clusters with strong H-alpha emission, indicates the absence of the most massive stars and suggests a mass upper limit of about 25 to 30 solar masses. The mass versus age distribution of the clusters shows a drastic decrease in the number of clusters with age, which indicates that cluster disruption is occurring on a timescale of about 10 Myr for low mass clusters. The cluster initial mass function for clusters younger than 10 Myr has an exponent of alpha = 2.0 (+- 0.05) We derived the cluster formation history from clusters with an initial mass larger than 10^4 solar masses. There is no evidence for a peak in the cluster formation rate within a factor two at about 200 to 400 Myr ago, i.e. at the time of the interaction with the companion galaxy NGC 5194.

We present the results of an analysis of stellar clusters in the interacting galaxy, M51, using HST-WFPC2 observations. The observations were made in five broad band filters; f336w (U), f439w (B), f555w (V), f675w (R), f814w (I). By comparing the spectral energy distribution (SED) of each source with those those of evolutionary spectral synthesis models, we have derived an age, mass, and extinction for each source. The comparison was done using the 3DEF method, that is based on a least $\chi^{2}$ test. We find 533 sources that are resolved and well fit with cluster models. Using this sample we study the cluster formation history, the cluster initial mass function (cIMF), and search for evidence of cluster disruption. We find that the cIMF is well represented by a power law with index of -2.05. In addition, we confirm the cluster disruption results by Boutloukos and Lamers (2002).

The stellar winds of massive stars show large changes in mass-loss rates and terminal velocities during their evolution from O-star through the Luminous Blue Variable phase to the Wolf-Rayet phase. The luminosity remains approximately unchanged during these phases. These large changes in wind properties are explained in the context of the radiation driven wind theory, of which we consider four different models. They are due to the evolutionary changes in radius, gravity and surface composition and to the change from optically thin (in continuum) line driven winds to optically thick radiation driven winds.

We present the first detection of 12 CO J=2->1 and 12 CO J=1->0 emission from the LBV AG Carinae. AG Carinae resides in a region which is very rich in molecular gas with complex motions. We find evidence of a slow outflow of molecular gas, expanding at ~ 7 km/s. This emission appears spatially unresolved. We argue that it is spatially localised, rather than extended, and possibly associated with the immediate circumstellar region of AG Carinae. Does it originate from a circumstellar envelope, similar to carbon stars, or from a circumstellar disk? The option of the circumstellar disk is preferable because it is consistent with additional independent indications for the existence of wind asymmetries in close proximity to the central star, found from spectropolarimetry and analysis of the UV and optical line profiles, and it provides the conditions of density and shielding necessary for the survival of the CO molecules in proximity to such a hot star (Teff ~ 14000 K - 20000 K). In the assumption that the CO emission originated when AG Carinae was in an evolved state, we derive a lower limit to the mass of molecular gas of 2.8 solar masses. This is smaller, but still comparable with the mass of ionized gas present in the circumstellar environment (4.2 solar masses), with the implication that the molecular gas fraction can contribute significantly to the overall mass lost from the central star in its post main sequence evolution.

We studied the HST-WFPC2 observations of the inner kpc of the interacting galaxy M51 in six bands from 2550 to 8140 Angstrom. The images show an oval shaped area ("bulge") of 11x16 arcsec or 450x650 pc around the nucleus, dominated by a smooth population of old stars with overimposed dust lanes. We found 30 bright point-like sources in the bulge of M51; many of these are bright in the UV. They are located in elongated "strings" which follow the general pattern of the dust lanes. The spectral energy distributions of the pointsources are compared with those of clusters or single stars. We present arguments that the pointsources are isolated massive stars or small groups of a few massive stars. The energy distributions are best fitted with those of single stars with 4000<Teff<50000 K, 4.2<log L/Lsun<7.2, and 12<Mstar<200 Msun. We have derived upper limits to the total mass of low mass stars that could be "hiding" within the point sources. For the bluest sources the upper limit is only a few hundred Msun. We conclude that massive stars are formed outside clusters (or in very low mass clusters) in the bulge of M51. The estimated star formation rate in the bulge is 1 to 2x10^(-3) Msun/yr. This star formation rate could be sustained for about 200 to 400 Myrs, suggesting that the ongoing massive star formation in the bulge of M51 is fed/triggered by the interaction with its companion about 400 Mrs ago. Theoretical cloud-models show that CO is destroyed in low extinction clouds near the bright nucleus of M51, so that H_2, [OI] and [CII] are the dominants coolants. This condition favors the formation of massive stars, in agreement with the observations. This mode of massive-star formation in the bulge of M51 may resemble the star formation in the early Universe, when the CO and dust contents were low due to the low metallicity.

Using data from the EROS2 microlensing survey, we report the discovery of two blue objects with irregular photometric behaviour of $\Delta V \sim 0.1$-0.4 mag on time scales of 20 to 200 days. They show a bluer when fainter behaviour. Subsequent spectra taken with the ESO 3.6m telescope show spectral type B4eIII and B2eIV-V with strong $H \alpha$ emission. These objects resemble the Herbig AeBe but also classical Be stars. At this stage, it is not possible to distinguish unambiguously between pre-main sequence and classical Be nature. If we favour the pre-main sequence interpretation, they are more luminous than the luminosity upper limit for Galactic HAeBe stars. The same was found for the HAeBe candidates in the LMC. This might be due to a shorter accretion time scale ($\tau = M_*/\dot{M}$), or the smaller dust content during the pre-main sequence evolution of SMC and LMC stars.

We use the analysis of the heavy element abundances (C, N, O, S) in circumstellar nebulae around Luminous Blue Variables to infer the evolutionary phase in which the material has been ejected. (1) We discuss the different effects that may have changed the gas composition of the nebula since it was ejected (2) We calculate the expected abundance changes at the stellar surface due to envelope convection in the red supergiant phase. If the observed LBV nebulae are ejected during the RSG phase, the abundances of the LBV nebulae require a significantly smaller amount of mass to be lost than assumed in evolutionary models. (3) We calculate the changes in the surface composition during the main sequence phase by rotation induced mixing. If the nebulae are ejected at the end of the MS-phase, the abundances in LBV nebulae are compatible with mixing times between 5 x 10^6 and 1 x 10^7 years. The existence of ON stars supports this scenario. (4) The predicted He/H ratio in the nebulae are significantly smaller than the current observed photospheric values of their central stars. Combining various arguments we show that the LBV nebulae are ejected during the blue SG phase and that the stars have not gone through a RSG phase. The chemical enhancements are due to rotation induced mixing, and the ejection is possibly triggered by near-critical rotation. During the ejection, the outflow was optically thick, which resulted in a large effective radius and a low effective temperature. This also explains the observed properties of LBV dust.

We have calculated a grid of massive star wind models and mass-loss rates for a wide range of metal abundances between 1/100 and 10 Z/Zsun. The calculation of this grid completes the Vink et al. (2000) mass-loss recipe with an additional parameter Z. We have found that the exponent of the power law dependence of mass loss vs. metallicity is constant in the range between 1/30 and 3 Z/Zsun. The mass-loss rate scales as Mdot \propto Z^0.85 Vinf^p with p = -1.23 for stars with Teff \ga 25000 K, and p = -1.60 for the B supergiants with Teff \la 25000 K. Taking also into account the metallicity dependence of Vinf, using the power law dependence Vinf \propto Z^0.13 from Leitherer et al. (1992), the overall result of mass loss as a function of metallicity can be represented by Mdot \propto Z^0.69 for stars with Teff \ga 25000 K, and Mdot \propto Z^0.64 for B supergiants with Teff \la 25000 K. Our mass-loss predictions are successful in explaining the observed mass-loss rates for Galactic and Small Magellanic Cloud O-type stars, as well as in predicting the observed Galactic bi-stability jump. Hence, we believe that our predictions are reliable and suggest that our mass-loss recipe be used in future evolutionary calculations of massive stars at different metal abundance. A computer routine to calculate mass loss is publicly available.

We report optical and near-infrared spectra of the central star of the radio source G79.29+0.46, a candidate luminous blue variable. The spectra contain numerous narrow (FWHM < 100 km/s emission lines of which the low-lying hydrogen lines are the strongest, and resemble spectra of other LBVc's and B[e] supergiants. A few prominent infrared lines are unidentified. The terminal wind speed is determined from H alpha to be 110 km/s. The strength of H alpha implies the presence of a very dense wind. Extended emission from H alpha and [N II] was detected but appears to be associated with the Cygnus X region rather than the radio source. Both diffuse interstellar bands and interstellar absorption lines are present in the optical spectrum of the central star, suggesting that there are both diffuse and molecular cloud components to the extinction and implying a minimum distance of 1 kpc and minimum luminosity of \~10^5 L(sun) for the star. The new spectra and their analysis indicate a low excitation, dense, and slowly expanding wind and support the LBVc classification.

We have calculated mass-loss rates for a grid of wind models covering a wide range of stellar parameters and have derived a mass-loss recipe for two ranges of effective temperature at either side of the bi-stability jump around spectral type B1. For a large sample of O stars, it is shown that there is now good agreement between these new theoretical mass-loss rates that take multiple scattering into account and observations. Agreement between the observed and new theoretical wind momenta increases confidence in the possibility to derive distances to luminous stars in distant stellar systems using the Wind momentum Luminosity Relation. For the winds of the B stars there is an inconsistency in the literature between various mass-loss rate determinations from observations by different methods. One group of the determinations for B stars does follow the new theoretical relation, while another group does not. The lack of agreement between the observed mass-loss rates derived by different methods may point to systematic errors in mass-loss determinations from observations for B stars. We show that our theoretical mass-loss recipe is reliable and recommend it be used in evolutionary calculations.

We have formulated the momentum equation for sectorial line driven winds from rotating stars including: (a) the oblateness of the star, (b) gravity darkening (von Zeipel effect), (c) conservation of angular momentum, (d) line driving specified by the force multiplier parameters (k, alpha, delta), (e) finite disk correction factors for an oblate star with gravity darkening for both the continuum and the line driving. The equations are solved numerically. We calculated the distribution of the mass flux and the wind velocity from the pole to the equator for the winds of B[e]-supergiants. Rotation decreases the terminal velocity in the equatorial region but hardly affects the wind velocity from the poles; it enhances the mass flux from the poles while the mass flux from the equator remains nearly the same. These effects increase with increasing rotation rates. We also calculated models with a bi-stability jump around 25 000 K, using force multipliers recently calculated with a Monte Carlo technique. In this case the mass flux increases drastically from the pole to the equator and the terminal velocity decreases drastically from pole to equator. This produces a density contrast in the wind rho(equator/pole) of about a factor 10 independent of the rotation rate of the star. We suggest that the observed density contrast of a factor 100 of the disks of B[e] stars may be reached by taking into account the wind compression due to the trajectories of the gas above the critical point towards the equatorial plane.

Observations of the B0.2V star tau Scorpii have revealed unusual stellar wind characteristics: red-shifted absorption in the far-ultraviolet O VI resonance doublet up to +250 km/s, and extremely hard X-ray emission implying gas at temperatures in excess of 10^7 K. We describe a phenomenological model to explain these properties. We assume the wind of tau Sco consists of two components: ambient gas in which denser clumps are embedded. The clumps are optically thick in the UV resonance lines primarily responsible for accelerating the ambient wind. The reduced acceleration causes the clumps to slow and even infall, all the while being confined by the ram pressure of the outflowing ambient wind. We calculate detailed trajectories of the clumps in the ambient stellar wind, accounting for a line radiation driving force and the momentum deposited by the ambient wind in the form of drag. We show these clumps will fall back towards the star with velocities of several hundred km/sec for a broad range of initial conditions. The infalling clumps produce X-ray emitting plasmas with temperatures in excess of (1-6)x10^7 K in bow shocks at their leading edge. The infalling material explains the peculiar red-shifted absorption wings seen in the O VI doublet. The required mass loss in clumps is 3% - 30% ofthe total mass loss rate. The model developed here can be generally applied to line-driven outflows with clumps or density irregularities. (Abstract Abridged)

We study the origin of the observed bi-stability jump in the terminal velocity of the winds of supergiants near spectral type B1. To this purpose, we have calculated a grid of wind models and mass-loss rates for these stars. The models show that the mass-loss rate 'jumps' by a factor of five around spectral type B1. Up to now, a theoretical explanation of the observed bi-stability jump was not yet provided by radiation driven wind theory. The models demonstrate that the subsonic part of the wind is dominated by the line acceleration due to Fe. The elements C, N and O are important line drivers in the supersonic part of the wind. We demonstrate that the mass-loss rate 'jumps' due to an increase in the line acceleration of Fe III below the sonic point. Finally, we discuss the possible role of the bi-stability jump on the mass loss during typical variations of Luminous Blue Variable stars.

We will review some results based on the EROS and the on going EROS-2 microlensing surveys to search for dark matter in the Galactic halo via microlensing effects on LMC/SMC stars. Microlensing surveys provide systematic observations for millions of stars over a long period in both Magellanic Clouds and give birth to very powerful database. We will quickly review the results obtained on pulsating stars and we will detail the similarity and differences of the Herbig Ae/Be stars in the Magellanic Clouds compared to their Galactic counterpart.

We present high-quality HST/GHRS spectra in the Hydrogen L alpha spectral region of Vega and Sirius-A. Thanks to the signal-to-noise ratio achieved in these observations and to the similarity of the two spectra, we found clear evidence of emission features in the low flux region, 1190-1222 A. These emission lines can be attributed unambiguously to Fe II and Cr II transitions. In this spectral range, silicon lines are observed in absorption. We built a series of non-LTE model atmospheres with different, prescribed temperature stratification in the upper atmosphere and treating Fe II with various degrees of sophistication in non-LTE. Emission lines are produced by the combined effect of the Schuster mechanism and radiative interlocking, and can be explained without the presence of a chromosphere. Silicon absorption lines and the L alpha profile set constraints on the presence of a chromosphere, excluding a strong temperature rise in layers deeper than tau_R ~ 10^{-4}.

We have searched for periodicity in the photometric and spectroscopic variations of the enigmatic star HD 101584 and found a long-term variability on a typical time scale of 1700 days (4.7 years) and a highly significant period of 218+-0.7 day. This period is most prominently present in the photometric indices which are a measure for the Balmer jump (Geneva d and Stromgren c). The Doppler velocities of the high-excitation photospheric absorption lines (HeI and CII) seem to be variable with the photometric period. Our data favor the 218 day period for the Doppler velocities with a small probability that the true period is 436 days. We argue that HD 101584 is a close (highly) eccentric 218 day binary system with a low-mass unseen secondary. The photometric and Doppler variations are attributed to changes in the velocity law and mass-loss rate of the stellar wind which lead to asymmetric line profiles and a phase dependent Balmer discontinuity. Binary interaction is responsible for the changes in velocity law and mass-loss rate leading to the observed phenomena.

High-resolution optical spectra (UES/WHT) of the extreme metal-poor post-AGB star HR4049 were obtained at four different orbital phases. The spectra cover the wavelength region from 3650ang to 10850ang at a resolution of R=5.2x10^4. These observations are supplemented with four high-resolution spectra of the NaI D1 & D2 and CaII K lines at R=10^5 (CAT/CES). The optical spectrum shows 217 spectral lines: the Balmer series (H alpha - H35), the Paschen series (P9 - P23), NI, OI and numerous CI lines. We show that the lines of H alpha, H beta, H gamma and NaI D show significant changes in profile between different observation dates. Nine components were identified in the profile of the NaI D lines of which three are circumstellar and six interstellar. The stronger CI lines are asymmetric and we derive a post-AGB mass-loss of Mdot= 6x10^-7Msol/yr from the asymmetry. The [OI] 6300ang line has been detected in emission at the system velocity and we argue that the emission is from an almost edge-on disk with a radius of about 20 Rstar.

We present the detection of molecular absorption lines in the optical spectrum of the post-AGB star HD56126. The C2 Phillips A^{1}\Pi_{u}-X^{1}\Sigma^{+}_{g} (1,0), (2,0), and (3,0); Swan d^{3}\Pi_{g}-a^{3}\Pi_{u} (0,0) and (1,0); and CN Red system A^{2}\Pi-X^{2}\Sigma^{+} (1,0), (2,0), (3,0), and (4,0) bands have been identified. From the identification of the molecular bands we find an expansion velocity of 8.5+-0.6 kms independent of excitation condition or molecular specie. On the basis of the expansion velocity, rotational temperatures, and molecular column densities we argue that the line-forming region is the AGB remnant. This is in agreement with the expansion velocity derived from the CO lines. We find column densities of log N_C2=15.3+-0.3 cm-2 and log N_CN=15.5+-0.3 cm-2, and rotational temperatures of T_rot=242+-20 K and T_rot=24+-5 K respectively for C2 and CN. By studying molecular line absorption in optical spectra of post-AGB stars we have found a new tracer of the AGB remnant. From comparison with the results of CO and IR observations it is possible to obtain information on non-spherical behavior of the AGB remnant. Using different molecules with different excitation conditions it should be possible to study the AGB remnant as a function of the distance to the star, and thus as a function of the evolutionary status of the star on the AGB.

We present a study of low- and high-resolution ultraviolet, high-resolution optical CAT/CES spectra and ultraviolet, optical and infrared photometry of the peculiar supergiant HD101584. From the photometry we learn that the ultraviolet and optical energy distribution cannot be fitted in a consistent way and we need a model in which the UV and optical energy distribution are formed by different gas. The Geneva photometry is best fitted to a B9II Kurucz model, Teff=12000+-1000K and log g=3.0 +-1.0, with an extinction of E(B-V)=0.49+-0.05. The observed spectral features in the spectrum of HD101584 are classified in eight different categories based on the velocity, shape of profile and the identification. The high-excitation HeI(chi=20.87eV), NII(chi=18.40eV), CII (chi=14.39eV) and NI (chi=10.29eV) optical absorption lines are formed in the photosphere of a late B-star (e.g. B8-9I-II). These absorption lines show radial velocity variations which are attributed to binary motion, with the secondary being a white dwarf or a low-mass main sequence object. The low-excitation P-Cygni lines in the optical and UV are formed in the wind. The number density of absorption lines in the UV is so large that the wind spectrum acts as an iron curtain in front of the B-star. The terminal velocity of the wind of $v_{\infty}=100\pm30$ km~s$^{-1}~$ is consistent with the star being a low-mass post-AGB star and the low effective gravity is attributed to the presence of a nearby, unseen, secondary. We estimate a mass-loss rate of \.{M}$\approx 10^{-8}$~M$_{\odot}$~yr$^{-1}$. Narrow absorption and emission lines are observed which are formed in a circumsystem disk with a typical radius of $10^{2} R_{\ast}$.

We present a list of post-AGB stars showing molecular line absorption and emission in the optical spectrum. Two objects show CH+, one in emission and one in absorption, and ten stars show C2 and CN in absorption. The Doppler velocities of the C2 lines and the rotational temperatures indicate that the line forming region is the AGB remnant. An analysis of the post-AGB stars of which CO millimeter data is available suggests that the C2} expansion velocity is of the same order as the CO expansion velocity. HD56126 has been studied in detail and we find a mass-loss rate of Mdot=2.8e-4 Msol/yr, fC2=2.4e-8 and fCN=1.3e-8. The mass loss derived from C2 is significantly larger than Mdot=1.2e-5 Msol/yr derived from CO. We find that all objects with the 21mu feature in emission show C2 and CN absorption, but not all objects with C2 and CN detections show a 21mu feature.

Wray 977, the B supergiant companion of the X-ray pulsar GX301-2, should be classified as B1~Ia+, based on a comparison of its optical spectrum to that of zeta1 Sco, a well-known B hypergiant and ONE OF THE BRIGHTEST STARS IN THE GALAXY. The classification of Wray 977 as a hypergiant results in a new distance determination of the binary system, i.e. 5.3 kpc (previously 1.8 +/- 0.4 kpc). The ``average'' X-ray luminosity of the pulsar is then ~10^37 erg/s, in good agreement with the predicted X-ray luminosity resulting from accretion of a dense, low-velocity (v_infty = 400 km/s) stellar wind. A mass-loss rate of <= 10^-5 M_sun/yr is estimated from the Halpha profile. A new upper limit for the inclination of the system is derived which provides a lower limit to the (present) mass of Wray 977 (48 M_sun). Regarding current binary evolution scenarios, the empirical lower mass limit for BLACK-HOLE formation in a binary increases to M_zams >= 50 M_sun.