The initial mass function of the rich young cluster NGC 1818 in the Large Magellanic Cloud

We use deep Hubble Space Telescope photometry of the rich, young (∼20- to 45-Myr old) star cluster NGC 1818 in the Large Magellanic Cloud to derive its stellar mass function (MF) down to  ∼0.15 M_⊙  . This represents the deepest robust MF thus far obtained for a stellar system in an extragalactic, low-metallicity  ([Fe/H]≃−0.4 dex)  environment. Combining our results with the published MF for masses above  1.0 M_⊙  , we obtain a complete present-day MF. This is a good representation of the cluster's initial MF (IMF), particularly at low masses, because our observations are centred on the cluster's uncrowded half-mass radius. Therefore, stellar and dynamical evolution of the cluster will not have affected the low-mass stars significantly. The NGC 1818 IMF is well described by both a lognormal and a broken power-law distribution with slopes of  Γ= 0.46 ± 0.10  and  Γ≃−1.35  (Salpeter-like) for masses in the range from 0.15 to  0.8 M_⊙  and greater than  0.8 M_⊙  , respectively. Within the uncertainties, the NGC 1818 IMF is fully consistent with both the Kroupa solar neighbourhood and the Chabrier lognormal mass distributions.     

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⁻¹  .     

Modelling the components of binaries in the Hyades: the dependence of the mixing-length parameter on stellar mass

We present our findings based on a detailed analysis of the binaries of the Hyades, in which the masses of the components are well known. We fit the models of the components of a binary system to observations so as to give the observed total V and B − V of that system and the observed slope of the main sequence in the corresponding parts. According to our findings, there is a very definite relationship between the mixing-length parameter and the stellar mass. The fitting formula for this relationship can be given as  α= 9.19( M /M_⊙− 0.74)^0.053− 6.65  , which is valid for stellar masses greater than  0.77 M_⊙  . While no strict information is gathered for the chemical composition of the cluster, as a result of degeneracy in the colour–magnitude diagram, by adopting   Z = 0.033  and using models for the components of 70 Tau and θ² Tau we find the hydrogen abundance to be   X = 0.676  and the age to be 670 Myr. If we assume that   Z = 0.024  , then   X = 0.718  and the age is 720 Myr. Our findings concerning the mixing-length parameter are valid for both sets of the solution. For both components of the active binary system V818 Tau, the differences between radii of the models with   Z = 0.024  and the observed radii are only about 4 per cent. More generally, the effective temperatures of the models of low-mass stars in the binary systems studied are in good agreement with those determined by spectroscopic methods.     

Monte Carlo simulations of star clusters – III. A million-body star cluster

A revision of Stodółkiewicz's Monte Carlo code is used to simulate the evolution of million-body star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. All models consist of 1 000 000 stars. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods. During the evolution, the initial mass function (IMF) changes significantly. The local mass function around the half-mass radius closely resembles the actual global mass function. At the late stages of evolution, the mass of the evolved stars inside the core can be as high as 97 per cent of the total mass in this region. For the whole system, the evolved stars can compose up to 75 per cent of the total mass. The evolution of cluster anisotropy strongly depends on initial cluster concentration, IMF and the strength of the tidal field. The results presented are the first step in the direction of simulating the evolution of real globular clusters by means of the Monte Carlo method.     

On the binarity of Herbig Ae/Be stars

We present high-resolution spectro-astrometry of a sample of 28 Herbig Ae/Be and three F-type pre-main-sequence stars. The spectro-astrometry, which is essentially the study of unresolved features in long-slit spectra, is shown from both empirical and simulated data to be capable of detecting binary companions that are fainter by up to 6 mag at separations larger than ∼0.1 arcsec. The nine targets that were previously known to be binary are all detected. In addition, we report the discovery of six new binaries and present five further possible binaries. The resulting binary fraction is 68 ± 11 per cent. This overall binary fraction is the largest reported for any observed sample of Herbig Ae/Be stars, presumably because of the exquisite sensitivity of spectro-astrometry for detecting binary systems. The data hint that the binary frequency of the Herbig Be stars is larger than that of the Herbig Ae stars. The Appendix presents model simulations to assess the capabilities of spectro-astrometry and reinforces the empirical findings. Most spectro-astrometric signatures in this sample of Herbig Ae/Be stars can be explained by the presence of a binary system. Two objects, HD 87643 and Z CMa, display evidence for asymmetric outflows. Finally, the position angles of the binary systems have been compared with available orientations of the circumprimary disc and these appear to be coplanar. The alignment between the circumprimary discs and the binary systems strongly suggests that the formation of binaries with intermediate-mass primaries is due to fragmentation as the alternative, stellar capture, does not naturally predict aligned discs. The alignment extends to the most massive B-type stars in our sample. This leads us to conclude that formation mechanisms that do result in massive stars, but predict random angles between the binaries and the circumprimary discs, such as stellar collisions, are also ruled out for the same reason.     

The influence of gas expulsion and initial mass segregation on the stellar mass function of globular star clusters

Recently, De Marchi, Paresce & Pulone studied a sample of 20 globular clusters and found that all clusters with high concentrations have steep stellar mass functions while clusters with low concentration have comparatively shallow mass functions. No globular clusters were found with a flat mass function and high concentration. This seems curious since more concentrated star clusters are believed to be dynamically more evolved and should have lost more low-mass stars via evaporation, which would result in a shallower mass function in the low-mass part.
We show that this effect can be explained by residual-gas expulsion from initially mass segregated star clusters, and is enhanced further through unresolved binaries. If gas expulsion is the correct mechanism to produce the observed trend in the   c –α  -plane, then observation of these parameters would allow to constrain cluster starting conditions such as star formation efficiency and the time-scale of gas expulsion.     

On the infant weight loss of low- to intermediate-mass star clusters

Star clusters are born in a highly compact configuration, typically with radii of less than about 1 pc roughly independently of mass. Since the star formation efficiency is less than 50 per cent by observation and because the residual gas is removed from the embedded cluster, the cluster must expand. In the process of doing so it only retains a fraction f _st of its stars. To date there are no observational constraints for f _st, although N -body calculations by Kroupa, Aarseth & Hurley suggest it to be about 20–30 per cent for Orion-type clusters. Here we use the data compiled by Testi et al., Testi, Palla & Natta and Testi, Palla & Natta for clusters around young Ae/Be stars and by de Wit et al. and de Wit et al. around young O stars and the study of de Zeeuw et al. of OB associations and combine these measurements with the expected number of stars in clusters with primary Ae/Be and O stars, respectively, using the empirical correlation between maximal stellar mass and star cluster mass of Weidner & Kroupa. We find that   f _st < 50  per cent with a decrease to higher cluster masses/more massive primaries. The interpretation would be that cluster formation is very disruptive. It appears that clusters with a birth stellar mass in the range  10–10³ M_⊙  keep at most 50 per cent of their stars.     

Accretion in stellar clusters and the initial mass function

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses M ∗≳0.5 M_⊙ . The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a ρ  ∝  R ⁻² distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of γ =−3/2 (where Salpeter is γ =−2.35) . Once the stars dominate the potential and are virialized, which occurs first in the cluster core, the accretion radius is the Bondi–Hoyle radius. The resultant mass spectrum has an asymptotic limit of γ =−2 with slightly steeper slopes ( γ ≈−2.5) if the stars are already mass-segregated. Simulations of accretion on to clusters containing 1000 stars show that, as expected, the low-mass stars accumulate the majority of their masses during the gas-dominated phase whereas the high-mass stars accumulate the majority of their masses during the stellar-dominated phase. This results in a mass spectrum with a relatively shallow γ ≈3/2 power law for low-mass stars and a steeper power law for high-mass stars −2.5≲ γ ≤−2 . This competitive accretion model also results in a mass-segregated cluster.     

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.     

On the formation of massive stars

We present a model for the formation of massive ( M ≳10 M⊙) stars through accretion-induced collisions in the cores of embedded dense stellar clusters. This model circumvents the problem of accreting on to a star whose luminosity is sufficient to reverse the infall of gas. Instead, the central core of the cluster accretes from the surrounding gas, thereby decreasing its radius until collisions between individual components become sufficient. These components are, in general, intermediate-mass stars that have formed through accretion on to low-mass protostars. Once a sufficiently massive star has formed to expel the remaining gas, the cluster expands in accordance with this loss of mass, halting further collisions. This process implies a critical stellar density for the formation of massive stars, and a high rate of binaries formed by tidal capture.     

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.     

Stellar contents and star formation in the young star cluster Be 59

We present   UBV  I _c   CCD photometry of the young open cluster Be 59 with the aim to study the star formation scenario in the cluster. The radial extent of the cluster is found to be ∼10 arcmin (2.9 pc). The interstellar extinction in the cluster region varies between   E ( B − V ) ≃ 1.4  to 1.8 mag. The ratio of total-to-selective extinction in the cluster region is estimated as  3.7 ± 0.3  . The distance of the cluster is found to be  1.00 ± 0.05 kpc  . Using near-infrared (NIR) colours and slitless spectroscopy, we have identified young stellar objects (YSOs) in the open cluster Be 59 region. The ages of these YSOs range between <1 and ∼2 Myr, whereas the mean age of the massive stars in the cluster region is found to be ∼2 Myr. There is evidence for second-generation star formation outside the boundary of the cluster, which may be triggered by massive stars in the cluster. The slope of the initial mass function, Γ, in the mass range  2.5 < M /M_⊙≤ 28  is found to be  −1.01 ± 0.11  which is shallower than the Salpeter value (−1.35), whereas in the mass range  1.5 < M /M_⊙≤ 2.5  the slope is almost flat. The slope of the K -band luminosity function is estimated as  0.27 ± 0.02  , which is smaller than the average value (∼0.4) reported for young embedded clusters. Approximately 32 per cent of Hα emission stars of Be 59 exhibit NIR excess indicating that inner discs of the T Tauri star (TTS) population have not dissipated. The Midcourse Space Experiment (MSX) and IRAS-HIRES images around the cluster region are also used to study the emission from unidentified infrared bands and to estimate the spatial distribution of optical depth of warm and cold interstellar dust.     

Magnetic activity and evolution of Algol-type stars – II

We examine the possibility of probing dynamo action in mass-losing stars, components of Algol-type binaries. Our analysis is based on the calculation of non-conservative evolution of these systems. We model the systems U Sge and β Per where the more massive companion fills its Roche lobe at the main sequence (case AB) and where it has a small helium core (early case B) respectively. We show that to maintain evolution of these systems at the late stages which are presumably driven by stellar 'magnetic braking', an efficient mechanism for producing large-scale surface magnetic fields in the donor star is needed. We discuss the relevance of dynamo operation in the donor star to the accelerated mass transfer during the late stages of evolution of Algol-type binaries. We suggest that the observed X-ray activity in Algol-type systems may be a good indicator of their evolutionary status and internal structure of the mass-losing stellar components.     

Models of α Centauri A and B with and without seismic constraints: time dependence of the mixing-length parameter

The α Centauri (α Cen) binary system is a well-known stellar system with very accurate observational constraints on the structure of its component stars. In addition to the classical non-seismic constraints, there are also seismic constraints for the interior models of α Cen A and B. These two types of constraint give very different values for the age of the system. While we obtain 8.9 Gyr for the age of the system from the non-seismic constraints, the seismic constraints imply that the age is about 5.6–5.9 Gyr. There may be observational or theoretical reasons for this discrepancy, which can be found by careful consideration of similar stars. The α Cen binary system, with its solar-type components, is also suitable for testing the stellar mass dependence of the mixing-length parameter for convection derived from the binaries of Hyades. The values of the mixing-length parameter for α Cen A and B are 2.10 and 1.90 for the non-seismic constraints. If we prioritize the seismic constraints, we obtain 1.64 and 1.91 for α Cen A and B, respectively. By taking into account these two contrasting cases for stellar mass dependence of the mixing-length parameter, we derive two expressions for its time dependence, which are also compatible with the mass dependence of the mixing-length parameter derived from the Hyades stars. For assessment, these expressions should be tested in other stellar systems and clusters.     

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.     

Binary stars and the fundamental initial mass function

We have carried out Monte Carlo simulations in which we generate a random pairing of objects drawn from a pre-assumed single-star power-law initial mass function (IMF), which we call the fundamental IMF. We show how the mass functions of primary stars and secondary stars and the mass function of the total mass of systems (if we could resolve them) differ from the underlying fundamental IMF for different slopes of this IMF. We also compare our results with the observed IMF, the binary frequency and the binary mass-ratio distributions for field stars and conclude that the fundamental IMF of subsolar mass stars could be steeper than is currently believed. In other words, the low-mass turn-over of the observed ('apparent') IMF could be spurious, if the main-sequence binary fraction of field stars is close to 100 per cent (perhaps owing to invisible companions).     

The distribution of ejected brown dwarfs in clusters

We examine the spatial distribution of brown dwarfs produced by the decay of small‐N stellar systems as expected from the embryo ejection scenario. We model a cluster of several hundred stars grouped into ‘cores’ of a few stars/brown dwarfs. These cores decay, preferentially ejecting their lowest‐mass members. Brown dwarfs are found to have a wider spatial distribution than stars, however once the effects of limited survey areas and unresolved binaries are taken into account it can be difficult to distinguish between clusters with many or no ejections. A large difference between the distributions probably indicates that ejections have occurred, however similar distributions sometimes arise even with ejections. Thus the spatial distribution of brown dwarfs is not necessarily a good discriminator between ejection and non‐ejection scenarios. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the galaxy stellar mass function, the mass–metallicity relation and the implied baryonic mass function

A comparison between published field galaxy stellar mass functions (GSMFs) shows that the cosmic stellar mass density is in the range 4–8 per cent of the baryon density (assuming  Ω_b= 0.045  ). There remain significant sources of uncertainty for the dust correction and underlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass function. We determine the   z < 0.05  GSMF using the New York University Value-Added Galaxy Catalog sample of 49 968 galaxies derived from the Sloan Digital Sky Survey and various estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is fitted with a double Schechter function that has  α₂≃−1.6  . At masses below  ∼10^8.5 M_⊙  , the GSMF may be significantly incomplete because of missing low-surface-brightness galaxies. One interpretation of the stellar mass–metallicity relation is that it is primarily caused by a lower fraction of available baryons converted to stars in low-mass galaxies. Using this principle, we determine a simple relationship between baryonic mass and stellar mass and present an 'implied baryonic mass function'. This function has a faint-end slope,  α₂≃−1.9  . Thus, we find evidence that the slope of the low-mass end of the galaxy mass function could plausibly be as steep as the halo mass function. We illustrate the relationship between halo baryonic mass function → galaxy baryonic mass function → GSMF. This demonstrates the requirement for peak galaxy formation efficiency at baryonic masses  ∼10¹¹ M_⊙  corresponding to a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower masses.