Competitive accretion in embedded stellar clusters

We investigate the physics of gas accretion in young stellar clusters. Accretion in clusters is a dynamic phenomenon as both the stars and the gas respond to the same gravitational potential. Accretion rates are highly non-uniform with stars nearer the centre of the cluster, where gas densities are higher, accreting more than others. This competitive accretion naturally results in both initial mass segregation and a spectrum of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal-lobe radius instead of the commonly used Bondi–Hoyle accretion radius. This works as both the stellar and gas velocities are under the influence of the same gravitational potential and are thus comparable. The low relative velocity which results means that R _tidal< R _BH in these systems. In contrast, when the stars dominate the potential and are virialized, R _BH< R _tidal and Bondi–Hoyle accretion is a better fit to the accretion rates.     

Mass segregation in young stellar clusters

We investigate the evolutionary effect of dynamical mass segregation in young stellar clusters. Dynamical mass segregation acts on a time-scale of order the relaxation time of a cluster. Although some degree of mass segregation occurs earlier, the position of massive stars in rich young clusters generally reflects the cluster's initial conditions. In particular, the positions of the massive stars in the Trapezium cluster in Orion cannot be due to dynamical mass segregation, but indicate that they formed in, or near, the centre of the cluster. Implications of this for cluster formation and for the formation of high-mass stars are discussed.     

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.     

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.     

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

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.     

Binary formation in stellar clusters

We consider how the tidal potential of a stellar cluster or a dense molecular cloud affects the fragmentation of gravitationally unstable molecular cloud cores. We find that molecular cloud cores which would collapse to form a single star in the absence of tidal shear, can be forced to fragment if they are subjected to tides. This may enhance the frequency of binaries in star-forming regions such as Ophiuchus and the frequency of binaries with separations ≲100 au in the Orion Trapezium Cluster. We also find that clouds which collapse to form binary systems in the absence of a tidal potential will form bound binary systems if exposed to weak tidal shear. However, if the tidal shear is sufficiently strong, even though the cloud still collapses to form two fragments, the fragments are pulled apart while they are forming by the tidal shear and two single stars are formed. This sets an upper limit for the separation of binaries that form near dense molecular clouds or in stellar clusters.     

The initial–final mass relationship of white dwarfs revisited: effect on the luminosity function and mass distribution

The initial–final mass relationship connects the mass of a white dwarf with the mass of its progenitor in the main sequence. Although this function is of fundamental importance to several fields in modern astrophysics, it is not well constrained either from the theoretical or from the observational points of view. In this work, we revise the present semi-empirical initial–final mass relationship by re-evaluating the available data. The distribution obtained from grouping all our results presents a considerable dispersion, which is larger than the uncertainties. We have carried out a weighted least-squares linear fit of these data and a careful analysis to give some clues on the dependence of this relationship on some parameters such as metallicity or rotation. The semi-empirical initial–final mass relationship arising from our study covers the range of initial masses from 1.0 to  6.5 M_⊙  , including in this way the low-mass domain, poorly studied until recently. Finally, we have also performed a test of the initial–final mass relationship by studying its effect on the luminosity function and on the mass distribution of white dwarfs. This was done by using different initial–final mass relationships from the literature, including the expression derived in this work, and comparing the results obtained with the observational data from the Palomar Green Survey and the Sloan Digital Sky Survey. We find that the semi-empirical initial–final mass relationship derived here gives results in good agreement with the observational data, especially in the case of the white dwarf mass distribution.