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.     

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

Evolutionary synthesis of stellar populations: a modular tool

A new tool for the evolutionary synthesis of stellar populations is presented, which is based on three independent matrices, giving respectively (1) the fuel consumption during each evolutionary phase as a function of stellar mass, (2) the typical temperatures and gravities during such phases, and (3) the colours and bolometric corrections as functions of gravity and temperature. The modular structure of the code allows one easily to assess the impact on the synthetic spectral energy distribution of the various assumptions and model ingredients, such as, for example, uncertainties in stellar evolutionary models, the mixing length, the temperature distribution of horizontal branch stars, asymptotic giant branch mass loss, and colour–temperature transformations. The so-called 'AGB phase transition' in Magellanic Cloud clusters is used to calibrate the contribution of the thermally pulsing asymptotic giant branch phase to the synthetic integrated luminosity. As an illustrative example, solar-metallicity ( Y  = 0.27, Z  = 0.02) models, with ages ranging between 30 Myr and 15 Gyr and various choices for the slope of the initial mass function, are presented. Synthetic broad-band colours and the luminosity contributions of the various evolutionary stages are compared with Large Magellanic Cloud and Galactic globular cluster data. In all these cases, a good agreement is found. Finally, the evolution is presented of stellar mass-to-light ratios in the bolometric and U B V R K passbands, in which the contribution of stellar remnants is accounted for.     

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.     

Galaxy bimodality versus stellar mass and environment

We analyse a   z < 0.1  galaxy sample from the Sloan Digital Sky Survey focusing on the variation in the galaxy colour bimodality with stellar mass     and projected neighbour density Σ, and on measurements of the galaxy stellar mass functions. The characteristic mass increases with environmental density from about  10^10.6  to     (Kroupa initial mass function,   H ₀= 70  ) for Σ in the range  0.1–10 Mpc⁻²  . The galaxy population naturally divides into a red and blue sequence with the locus of the sequences in colour–mass and colour–concentration indices not varying strongly with environment. The fraction of galaxies on the red sequence is determined in bins of 0.2 in  log Σ  and     bins). The red fraction   f _r   generally increases continuously in both Σ and     such that there is a unified relation:     . Two simple functions are proposed which provide good fits to the data. These data are compared with analogous quantities in semi-analytical models based on the Millennium N -body simulation: the Bower et al. and Croton et al. models that incorporate active galactic nucleus feedback. Both models predict a strong dependence of the red fraction on stellar mass and environment that is qualitatively similar to the observations. However, a quantitative comparison shows that the Bower et al. model is a significantly better match; this appears to be due to the different treatment of feedback in central galaxies.     

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.