No Arabic abstract
Some ultra-compact dwarf galaxies have large dynamical mass to light (M/L) ratios and also appear to contain an overabundance of LMXB sources, and some Milky Way globular clusters have a low concentration and appear to have a deficit of low-mass stars. These observations can be explained if the stellar IMF becomes increasingly top-heavy with decreasing metallicity and increasing gas density of the forming object. The thus constrained stellar IMF then accounts for the observed trend of metallicity and M/L ratio found amongst M31 globular star clusters. It also accounts for the overall shift of the observationally deduced galaxy-wide IMF from top-light to top-heavy with increasing star formation rate amongst galaxies. If the IMF varies similarly to deduced here, then extremely young very massive star-burst clusters observed at a high redshift would appear quasar-like (Jerabkova et al. 2017) .
We address the turbulent fragmentation scenario for the origin of the stellar initial mass function (IMF), using a large set of numerical simulations of randomly driven supersonic MHD turbulence. The turbulent fragmentation model successfully predicts the main features of the observed stellar IMF assuming an isothermal equation of state without any stellar feedback. As a test of the model, we focus on the case of a magnetized isothermal gas, neglecting stellar feedback, while pursuing a large dynamic range in both space and timescales covering the full spectrum of stellar masses from brown dwarfs to massive stars. Our simulations represent a generic 4 pc region within a typical Galactic molecular cloud, with a mass of 3000 Msun and an rms velocity 10 times the isothermal sound speed and 5 times the average Alfven velocity, in agreement with observations. We achieve a maximum resolution of 50 au and a maximum duration of star formation of 4.0 Myr, forming up to a thousand sink particles whose mass distribution closely matches the observed stellar IMF. A large set of medium-size simulations is used to test the sink particle algorithm, while larger simulations are used to test the numerical convergence of the IMF and the dependence of the IMF turnover on physical parameters predicted by the turbulent fragmentation model. We find a clear trend toward numerical convergence and strong support for the model predictions, including the initial time evolution of the IMF. We conclude that the physics of isothermal MHD turbulence is sufficient to explain the origin of the IMF.
Galaxy formation simulations frequently use Initial Mass Function (IMF) averaged feedback prescriptions, where star particles are assumed to represent single stellar populations that fully sample the IMF. This approximation breaks down at high mass resolution, where stochastic variations in stellar populations become important. We discuss various schemes to populate star particles with stellar masses explicitly sampled from the IMF. We use Monte Carlo numerical experiments to examine the ability of the schemes to reproduce an input IMF in an unbiased manner while conserving mass. We present our preferred scheme which can easily be added to pre-existing star formation prescriptions. We then carry out a series of high resolution isolated simulations of dwarf galaxies with supernovae, photoionization and photoelectric heating to compare the differences between using IMF averaged feedback and explicitly sampling the IMF. We find that if supernovae are the only form of feedback, triggering individual supernovae from IMF averaged rates gives identical results to IMF sampling. However, we find that photoionization is more effective at regulating star formation when IMF averaged rates are used, creating more, smaller H II regions than the rare, bright sources produced by IMF sampling. We note that the increased efficiency of the IMF averaged feedback versus IMF sampling is not necessarily a general trend and may be reversed depending on feedback channel, resolution and other details. However, IMF sampling is always the more physically motivated approach. We conservatively suggest that it should be used for star particles less massive than $sim500,mathrm{M_odot}$.
(abridged) We develop a model which describes the coevolution of the mass function of dense cores and of the IMF in a protocluster clump. In the model, cores injected in the clump evolve under the effect of gas accretion. Accretion onto the cores follows a time-dependent accretion rate that describes accretion in a turbulent medium. Once the accretion timescales of cores exceed their contraction timescales, they are turned into stars. We include the effect of feedback by the newly formed massive stars through their stellar winds. A fraction of the winds energy is assumed to counter gravity and disperse the gas from the protocluster and as a consequence, quench further star formation. The latter effect sets the final IMF of the cluster. We apply our model to a clump that is expected to resemble the progenitor clump of the Orion Nebula Cluster (ONC). Our model is able to reproduce both the shape and normalization of the ONCs IMF and the mass function of dense cores in Orion. The complex features of the ONCs IMF,i.e., a shallow slope in the mass range ~0.3-2.5 Msol,a steeper slope in the mass range ~2.5-12 Msol, and a nearly flat tail at the high mass end are reproduced. The model predicts a rapid star formation process with an age spread for the stars of 2.3 10^5 yr which is consistent with the fact that 80% of the ONCs stars have ages of <=0.3 Myr. The model predicts a primordial mass segregation with the most massive stars being born in the region between 2-4 times the core radius of the cluster. In parallel, the model also reproduces, simultaneously, the mass function of dense cores in Orion. We study the effects of varying the model parameters on the resulting IMF and show that the IMF of stellar clusters is expected to show significant variations, provided variations in the clumps and cores properties exist.
The stellar initial mass function (IMF), which is often assumed to be universal across unresolved stellar populations, has recently been suggested to be bottom-heavy for massive ellipticals. In these galaxies, the prevalence of gravity-sensitive absorption lines (e.g. Na I and Ca II) in their near-IR spectra implies an excess of low-mass ($m <= 0.5$ $M_odot$) stars over that expected from a canonical IMF observed in low-mass ellipticals. A direct extrapolation of such a bottom-heavy IMF to high stellar masses ($m >= 8$ $M_odot$) would lead to a corresponding deficit of neutron stars and black holes, and therefore of low-mass X-ray binaries (LMXBs), per unit near-IR luminosity in these galaxies. Peacock et al. (2014) searched for evidence of this trend and found that the observed number of LMXBs per unit $K$-band luminosity ($N/L_K$) was nearly constant. We extend this work using new and archival Chandra X-ray Observatory (Chandra) and Hubble Space Telescope (HST) observations of seven low-mass ellipticals where $N/L_K$ is expected to be the largest and compare these data with a variety of IMF models to test which are consistent with the observed $N/L_K$. We reproduce the result of Peacock et al. (2014), strengthening the constraint that the slope of the IMF at $m >= 8$ $M_odot$ must be consistent with a Kroupa-like IMF. We construct an IMF model that is a linear combination of a Milky Way-like IMF and a broken power-law IMF, with a steep slope ($alpha_1=$ $3.84$) for stars < 0.5 $M_odot$ (as suggested by near-IR indices), and that flattens out ($alpha_2=$ $2.14$) for stars > 0.5 $M_odot$, and discuss its wider ramifications and limitations.
In this contributed talk I present recent results on the connection between stellar population properties and the normalisation of the stellar initial mass function (IMF) measured using stellar dynamics, based on a large sample of 260 early-type galaxies observed as part of the Atlas3D project. This measure of the IMF normalisation is found to vary non-uniformly with age- and metallicity-sensitive absorption line strengths. Applying single stellar population models, there are weak but measurable trends of the IMF with age and abundance ratio. Accounting for the dependence of stellar population parameters on velocity dispersion effectively removes these trends, but subsequently introduces a trend with metallicity, such that `heavy IMFs favour lower metallicities. The correlations are weaker than those found from previous studies directly detecting low-mass stars, suggesting some degree of tension between the different approaches of measuring the IMF. Resolving these discrepancies will be the focus of future work.