The influence of multiple stars on the high-mass stellar initial mass function and age dating of young massive star clusters

The study of young stellar populations has revealed that most stars are in binary or higher order multiple systems. In this study, the influence on the stellar initial mass function (IMF) of large quantities of unresolved multiple massive stars is investigated by taking into account the stellar evolution and photometrically determined system masses. The models, where initial masses are derived from the luminosity and colour of unresolved multiple systems, show that even under extreme circumstances (100 per cent binaries or higher order multiples), the difference between the power-law index of the mass function (MF) of all stars and the observed MF is small (≲0.1). Thus, if the observed IMF has the Salpeter index  α= 2.35  , then the true stellar IMF has an index not flatter than  α= 2.25  . Additionally, unresolved multiple systems may hide between 15 and 60 per cent of the underlying true mass of a star cluster. While already a known result, it is important to point out that the presence of a large number of unresolved binaries amongst pre-main-sequence stars induces a significant spread in the measured ages of these stars even if there is none. Also, lower mass stars in a single-age binary-rich cluster appear older than the massive stars by about 0.6 Myr.     

The origin of the initial mass function and its dependence on the mean Jeans mass in molecular clouds

We investigate the dependence of stellar properties on the mean thermal Jeans mass in molecular clouds. We compare the results from the two largest hydrodynamical simulations of star formation to resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical except for the radii of the clouds, which are chosen so that the mean densities and mean thermal Jeans masses of the clouds differ by factors of 9 and 3, respectively.
We find that the denser cloud, with the lower mean thermal Jeans mass, produces a higher proportion of brown dwarfs and has a lower characteristic (median) mass of the stars and brown dwarfs. This dependence of the initial mass function (IMF) on the density of the cloud may explain the observation that the Taurus star-forming region appears to be deficient in brown dwarfs when compared with the Orion Trapezium cluster. The new calculation also produces wide binaries (separations >20 au), one of which is a wide binary brown dwarf system.
Based on the hydrodynamical calculations, we develop a simple accretion/ejection model for the origin of the IMF. In the model, all stars and brown dwarfs begin with the same mass (set by the opacity limit for fragmentation) and grow in mass until their accretion is terminated stochastically by their ejection from the cloud through dynamically interactions. The model predicts that the main variation of the IMF in different star-forming environments should be in the location of the peak (due to variations in the mean thermal Jeans mass of the cloud) and in the substellar regime. However, the slope of the IMF at high masses may depend on the dispersion in the accretion rates of protostars.     

The Effect of the Initial Mass Function of Stars on the Burst Rate of GRBs

For the mechanism of production of γ-ray bursts (GRBs) it is rather generally recognized that the long-term γ-ray burst (LGRB) originates from the deaths of massive stars while the short-term γ-ray burst (SGRB) originates from the merging of close binaries. Therefore the speculation naturally follows that the number of LGRBs is directly proportional to the star formation rate (SFR). However, it is indicated from recent data analyses that this speculation does not fit the observations very well. It is considered that only massive stars with masses greater than a certain critical mass can produce the LGRB, so the initial mass function (IMF) of stars can significantly affect the production rate of LGRBs. In this paper it is considered that the IMF of stars can be used to explain the observed number distribution of the LGRBs with the redshift, and this has led to some good results.     

Determination of the IMF on the basis of a recently derived SFR history in the solar neighbourhood

The initial mass function (IMF) in the solar neighbourhood is determined on the basis of a recently derived history of the star formation rate (SFR) which shows the presence of a star formation burst about 8 Gyr ago. The observed present-day mass function (PDMF) is considered, and the average past distribution of stars of a given mass is estimated. Two cases are considered, namely (i) constant SFR, and (ii) variable SFR as derived from the new metallicity distribution of G dwarfs. The resulting IMF is compared with previous determinations by Scalo and Kroupa et al., and the variation with stellar mass of the slope of the IMF is compared with reference determinations in the literature.     

The proper-motion signal of unresolved binaries in the Hipparcos catalogue

We present an investigation of the differences between quasi-instantaneous stellar proper motions from the Hipparcos catalogue and long-term proper motions determined by combining Hipparcos and the Astrographic Catalogue. Our study is based on a sample of about 12000 stars of visual magnitude from 7 to 10 in two declination zones on the northern and equatorial sky. The distribution of the proper-motion differences shows an excess of large deviations. This is caused by the influence of orbital motion of unresolved binary systems. The proper-motion deviations provide statistical evidence for 360 astrometric binaries in the investigated zones, corresponding to about 2400 such binaries in the entire Hipparcos catalogue, in addition to those already known. In order to check whether the observed deviations are compatible with standard assumptions on the basic parameters of binary stars, we model the impact of orbital motion on the observed proper motions in a Monte Carlo simulation. We show that the simulation yields an acceptable approximation of the observations, if a binary frequency between 70% and 100% is assumed, i.e.if most of the stars in the sample are assumed to have a companion. Thus Hipparcos astrometric binaries confirm that the frequency of non-single stars among field stars is very high. We also investigate the influence of the mass function for the secondary component on the result of the simulation. A constant mass function and mass functions with moderate increase towards low masses lead to results, which are compatible with the observed proper-motion effects. A high preponderance of very-low-mass or substellar companions as produced, for example, by a M^—1 power law is not in agreement with the frequency of proper-motion deviations in our sample of stars.     

Accretion and dynamical interactions in small-N star-forming clusters: N= 5

We present results from high-resolution hydrodynamical simulations that explore the effects of small-scale clustering in star-forming regions. A large ensemble of small- N clusters with five stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data.
Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses that is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed initial mass function, in both the stellar and substellar regimes. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities ≥2 km s⁻¹, which might deprive low-mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long-period (uncircularized) binaries.
The results obtained support a mechanism in which a significant fraction of brown dwarfs form under similar circumstances as those of normal stars but are ejected from the common envelope of unstable multiple systems before their masses exceed the hydrogen burning limit. We predict that many close binary stars should have wide brown dwarf companions. Brown dwarfs, and, in general, very low-mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or substellar primaries being scarce. Where such binaries exist, they are expected either to be close enough (semimajor axis ∼10 au) to survive strong interactions with more massive binaries or to be born in very small molecular cloud cores.     

Early star formation and the evolution of the stellar initial mass function in galaxies

It has frequently been suggested in the literature that the stellar IMF in galaxies was top-heavy at early times. This would be plausible physically if the IMF depended on a mass-scale such as the Jeans mass that was higher at earlier times because of the generally higher temperatures that were present then. In this paper it is suggested, on the basis of current evidence and theory, that the IMF has a universal Salpeter-like form at the upper end, but flattens below a characteristic stellar mass that may vary with time. Much of the evidence that has been attributed to a top-heavy early IMF, including the ubiquitous G-dwarf problem, the high abundance of heavy elements in clusters of galaxies, and the high rate of formation of massive stars in high-redshift galaxies, can be accounted for with such an IMF if the characteristic stellar mass was several times higher during the early stages of galaxy evolution. However, significant variations in the mass-to-light ratios of galaxies and large amounts of dark matter in stellar remnants are not as easily explained in this way, because they require more extreme and less plausible assumptions about the form and variability of the IMF. Metal-free 'population III' stars are predicted to have an IMF that consists exclusively of massive stars, and they could help to account for some of the evidence that has been attributed to a top-heavy early IMF, as well as contributing importantly to the energetics and chemical enrichment of the early Universe.     

On the mass of dense star clusters in starburst galaxies from spectrophotometry

The mass of unresolved young star clusters derived from spectrophotometric data may well be off by a factor of 2 or more once the migration of massive stars driven by mass segregation is accounted for. We quantify this effect for a large set of cluster parameters, including variations in the stellar initial mass function (IMF), the intrinsic cluster mass, and mean mass density. Gas-dynamical models coupled with the Cambridge stellar evolution tracks allow us to derive a scheme to recover the real cluster mass given measured half-light radius, one-dimensional velocity dispersion and age. We monitor the evolution with time of the ratio of real to apparent mass through the parameter η. When we compute η for rich star clusters, we find non-monotonic evolution in time when the IMF stretches beyond a critical cut-off mass of  25.5 M_⊙  . We also monitor the rise of colour gradients between the inner and outer volume of clusters: we find trends in time of the stellar IMF power indices overlapping well with those derived for the Large Magellanic Cloud cluster NGC 1818 at an age of 30 Myr. We argue that the core region of massive Antennae clusters should have suffered from much segregation despite their low ages. We apply these results to a cluster mass function, and find that the peak of the mass distribution would appear to observers shifted to lower masses by as much as 0.2 dex. The star formation rate derived for the cluster population is then underestimated by from 20 to 50 per cent.     

A top-heavy stellar initial mass function in starbursts as an explanation for the high mass-to-light ratios of ultra-compact dwarf galaxies

It has been recently shown that the dynamical V -band mass-to-light ratios of compact stellar systems with masses from 10⁶ to  10⁸ M_⊙  are not consistent with the predictions from simple stellar population models. Top-heavy stellar initial mass functions (IMFs) in these so-called ultra-compact dwarf galaxies (UCDs) offer an attractive explanation for this finding, the stellar remnants and retained stellar envelopes providing the unseen mass. We therefore construct a model which quantifies by how much the IMFs of UCDs would have to deviate in the intermediate- and high-mass range from the canonical IMF in order to account for the enhanced   M / L_V   ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of the UCDs and the number of faint products of stellar evolution retained by them. Assuming that the IMF in the UCDs is a three-part power law equal to the canonical IMF in the low-mass range and taking 20 per cent as a plausible choice for the fraction of the remnants of high-mass stars retained by UCDs, the model suggests the exponent of the high-mass IMF to be ≈1.6 if the UCDs are  13 Gyr  old (i.e. almost as old as the Universe) or ≈1.0 if the UCDs are  7 Gyr  old, in contrast to 2.3 for the Salpeter–Massey IMF. If the IMF was as top heavy as suggested here, the stability of the UCDs might have been threatened by heavy mass loss induced by the radiation and evolution of massive stars. The central densities of UCDs must have been in the range  10⁶ to 10⁷ M_⊙ pc⁻³  when they formed with star formation rates of  10 to 100 M_⊙ yr⁻¹  .     

On the problem of secular variability in the stellar initial mass function

The hypothesis of secular variations in the stellar initial mass function (IMF) is studied. It is found that neither the present-day mass function of the nearby main sequence field stars nor the velocity distribution of these stars are contradictory with this hypothesis. The luminosity functions of unbiased kinematically defined age-groups of the nearby stars also provide no strong constraints. Simple evolution models with time-dependent IMF and star formation rate enable to fit the data of blue irregular galaxies. Some problems with an universal IMF are pointed out.     

The role of tidal interactions in star formation

Nearly all of the initial angular momentum of the matter that goes into each forming star must somehow be removed or redistributed during the formation process. The possible transport mechanisms and the possible fates of the excess angular momentum are discussed, and it is argued that transport processes in discs are probably not sufficient by themselves to solve the angular momentum problem, while tidal interactions with other stars in forming binary or multiple systems are likely to be of very general importance in redistributing angular momentum during the star formation process. Most, if not all, stars probably form in binary or multiple systems, and tidal torques in these systems can transfer much of the angular momentum from the gas around each forming star to the orbital motions of the companion stars. Tidally generated waves in circumstellar discs may contribute to the overall redistribution of angular momentum. Stars may gain much of their mass by tidally triggered bursts of rapid accretion, and these bursts could account for some of the most energetic phenomena of the earliest stages of stellar evolution, such as jet-like outflows. If tidal interactions are indeed of general importance, planet-forming discs may often have a more chaotic and violent early evolution than in standard models, and shock heating events may be common. Interactions in a hierarchy of subgroups may play a role in building up massive stars in clusters and in determining the form of the upper initial mass function (IMF) . Many of the processes discussed here have analogues on galactic scales, and there may be similarities between the formation of massive stars by interaction-driven accretion processes in clusters and the buildup of massive black holes in galactic nuclei.     

A highly abnormal massive star mass function in the Orion Nebula cluster and the dynamical decay of trapezium systems

The Orion Nebula cluster (ONC) appears to be unusual on two grounds: the observed constellation of the OB stars of the entire ONC and its Trapezium at its centre implies a time-scale problem given the age of the Trapezium, and an initial mass function (IMF) problem for the whole OB star population in the ONC. Given the estimated crossing time of the Trapezium, it ought to have totally dynamically decayed by now. Furthermore, by combining the lower limit of the ONC mass with a standard IMF it emerges that the ONC should have formed at least about 40 stars heavier than  5 M_⊙  while only 10 are observed. Using the N -body experiments we (i) confirm the expected instability of the Trapezium and (ii) show that beginning with a compact OB-star configuration of about 40 stars both the number of observed OB stars after 1 Myr within 1 pc radius and a compact trapezium configuration can be reproduced. These two empirical constraints thus support our estimate of 40 initial OB stars in the cluster. Interestingly, a more-evolved version of the ONC resembles the Upper Scorpius OB association. The N -body experiments are performed with the new C-code catena by integrating the equations of motion using the chain-multiple-regularization method. In addition, we present a new numerical formulation of the IMF.     

The Jeans mass and the origin of the knee in the IMF

We use numerical simulations of the fragmentation of a  1000 M_⊙  molecular cloud and the formation of a stellar cluster to study how the initial conditions for star formation affect the resulting initial mass function (IMF). In particular, we are interested in the relation between the thermal Jeans mass in a cloud and the knee of the IMF, i.e. the mass separating the region with a flat IMF slope from that typified by a steeper, Salpeter-like, slope. In three isothermal simulations with   M _Jeans= 1, 2  and  5 M_⊙  , the number of stars formed, at comparable dynamical times, scales roughly with the number of initial Jeans masses in the cloud. The mean stellar mass also increases (though less than linearly) with the initial Jeans mass in the cloud. It is found that the IMF in each case displays a prominent knee, located roughly at the mass scale of the initial Jeans mass. Thus clouds with higher initial Jeans masses produce IMFs which are shallow to higher masses. This implies that a universal IMF requires a physical mechanism that sets the Jeans mass to be near  1 M_⊙  . Simulations including a barotropic equation of state as suggested by Larson, with cooling at low densities followed by gentle heating at higher densities, are able to produce realistic IMFs with the knee located at  ≈1 M_⊙  , even with an initial   M _Jeans= 5 M_⊙  . We therefore suggest that the observed universality of the IMF in the local Universe does not require any fine tuning of the initial conditions in star forming clouds but is instead imprinted by details of the cooling physics of the collapsing gas.     

Unresolved Binaries and the Initial Mass Function

In the area of binary stars the important contribution of GAIA will be towards statistical studies of the binary star distribution functions. On the other hand, from the GAIA observations, more than 10 000 masses accurate to 1% can be expected. This improves accuracy of the stellar initial mass function (IMF) still poorly known, particularly for the low-mass stars. A goal of the present study is to simulate a solar vicinity model with 100% binaries for predicting observational distributions. Under the assumption that star formation rate is constant, a random pairing of objects drawn from a pre-assumed single star power-law IMF is generated. Stellar evolution and selection effects are taken into account. Orbital parameters are assumed to be spread according to the common distributions. It is shown that under given assumptions the power-law IMF does lead to quasi-lognormal turnovers for the resulting observational mass function. The described model can be used to estimate, under reasonable assumptions on the IMF shape, the number of binaries of the different types that can be observed with GAIA. This revised version was published online in July 2006 with corrections to the Cover Date.     

The massive star binary fraction in young open clusters – I. NGC 6231 revisited

We present the results of a long-term high-resolution spectroscopy campaign on the O-type stars in NGC 6231. We revise the spectral classification and multiplicity of these objects and we constrain the fundamental properties of the O-star population. Almost three quarters of the O-type stars in the cluster are members of a binary system. The minimum binary fraction is 0.63, with half the O-type binaries having an orbital period of the order of a few days. The eccentricities of all the short-period binaries are revised downward, and henceforth match a normal period–eccentricity distribution. The mass ratio distribution shows a large preference for O + OB binaries, ruling out the possibility that, in NGC 6231, the companion of an O-type star is randomly drawn from a standard initial mass function. Obtained from a complete and homogeneous population of O-type stars, our conclusions provide interesting observational constraints to be confronted with the formation and early evolution theories of O-stars.     

The low-mass initial mass function in the young cluster NGC 6611

NGC 6611 is the massive young cluster (2–3 Myr) that ionizes the Eagle Nebula. We present very deep photometric observations of the central region of NGC 6611 obtained with the Hubble Space Telescope and the following filters: ACS/WFC F775W and F850LP and NIC2 F110W and F160W, loosely equivalent to ground-based IZJH filters. This survey reaches down to   I ∼ 26 mag  . We construct the initial mass function (IMF) from  ∼1.5 M_⊙  well into the brown dwarf regime (down to  ∼0.02 M_⊙  ). We have detected 30–35 brown dwarf candidates in this sample. The low-mass IMF is combined with a higher-mass IMF constructed from the ground-based catalogue from Oliveira et al. We compare the final IMF with those of well-studied star-forming regions: we find that the IMF of NGC 6611 more closely resembles that of the low-mass star-forming region in Taurus than that of the more massive Orion Nebula Cluster. We conclude that there seems to be no severe environmental effect in the IMF due to the proximity of the massive stars in NGC 6611.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

The abundance of brown dwarfs

The amount of mass contained in low-mass objects is investigated anew. Instead of using a mass–luminosity relation to convert a luminosity function to a mass function, I predict the mass–luminosity relation from assumed mass functions and the luminosity functions of Jahreiss & Wielen and Gould, Bahcall & Flynn. Comparison of the resulting mass–luminosity relations with data for binary stars constrains the permissible mass functions. If the mass function is assumed to be a power law, the best-fitting slope lies either side of the critical slope, α =−2, below which the mass in low-mass objects is divergent, depending on the luminosity function adopted. If these power-law mass functions are truncated at 0.001 M_⊙, the contribution to the local density from stars lies between 0.013 and 0.10 M_⊙ pc⁻³ depending on the mass at which the mass function is normalized and the adopted value of α . Recent dynamical estimates of the local mass density rule out stellar mass densities above ∼0.05 M_⊙ pc⁻³. Hence, power laws steeper than α =−2 that extend down to 0.001 M_⊙ are allowed only if one adopts an implausible normalization of the mass function. If the mass function is generalized from a power law to a low-order polynomial in log( M ), the mass in stars with M <0.1 M_⊙ is either negligible or strongly divergent, depending on the order of the polynomial adopted.     

The formation of a bound star cluster: from the Orion nebula cluster to the Pleiades

Direct N -body calculations are presented of the formation of Galactic clusters using GasEx , which is a variant of the code Nbody6 . The calculations focus on the possible evolution of the Orion nebula cluster (ONC) by assuming that the embedded OB stars explosively drove out 2/3 of its mass in the form of gas about 0.4 Myr ago. A bound cluster forms readily and survives for 150 Myr despite additional mass loss from the large number of massive stars, and the Galactic tidal field. This is the very first time that cluster formation is obtained under such realistic conditions. The cluster contains about 1/3 of the initial 10⁴ stars, and resembles the Pleiades cluster to a remarkable degree, implying that an ONC-like cluster may have been a precursor of the Pleiades. This scenario predicts the present expansion velocity of the ONC, which will be measurable by upcoming astrometric space missions. These missions should also detect the original Pleiades members as an associated expanding young Galactic-field subpopulation. The results arrived at here suggest that Galactic clusters form as the nuclei of expanding OB associations.
The results have wide implications, also for the formation of globular clusters and the Galactic-field and halo stellar populations. In view of this, the distribution of binary orbital periods and the mass function within and outside the model ONC and Pleiades is quantified, finding consistency with observational constraints. Advanced mass segregation is evident in one of the ONC models. The calculations show that the primordial binary population of both clusters could have been much the same as is observed in the Taurus–Auriga star-forming region. The computations also demonstrate that the binary proportion of brown dwarfs is depleted significantly for all periods, whereas massive stars attain a high binary fraction.