No Arabic abstract
We examine the cosmic evolution of a stellar initial mass function (IMF) in galaxies that varies with the Jeans mass in the interstellar medium, paying particular attention to the K-band stellar mass to light ratio (M/L_K) of present-epoch massive galaxies. We calculate the typical Jeans mass using high-resolution hydrodynamic simulations coupled with a fully radiative model for the ISM, which yields a parameterisation of the IMF characteristic mass as a function of galaxy star formation rate (SFR). We then calculate the star formation histories of galaxies utilising an equilibrium galaxy growth model coupled with constraints on the star formation histories set by abundance matching models. We find that at early times, energetic coupling between dust and gas drive warm conditions in the ISM, yielding bottom-light/top- heavy IMFs associated with large ISM Jeans masses for massive star-forming galaxies. Owing to the remnants of massive stars that formed during the top-heavy phases at early times, the resultant M/L_K(sigma) in massive galaxies at the present epoch is increased relative to the non- varying IMF case. At late times, lower cosmic ray fluxes allow for cooler ISM temperatures in massive galaxies, and hence newly formed clusters will exhibit bottom-heavy IMFs, further increasing M/L_K(sigma). Our central result is hence that a given massive galaxy may go through both top-heavy and bottom-heavy IMF phases during its lifetime, though the bulk of the stars form during a top-heavy phase. Qualitatively, the variations in M/L_K(sigma) with galaxy mass are in agreement with observations, however, our model may not be able to account for bottom-heavy mass functions as indicated by stellar absorption features.
We tested the implementation of different IMFs in our model for the chemical evolution of ellipticals, with the aim of reproducing the observed relations of [Fe/H] and [Mg/Fe] abundances with galaxy mass in a sample of early-type galaxies selected from the SPIDER-SDSS catalog. Abundances in the catalog were derived from averaged spectra, obtained by stacking individual spectra according to central velocity dispersion, as a proxy of galaxy mass. We tested initial mass functions already used in a previous work, as well as two new models, based on low-mass tapered (bimodal) IMFs, where the IMF becomes either (1) bottom-heavy in more massive galaxies, or (2) is time-dependent, switching from top-heavy to bottom-heavy in the course of galactic evolution. We found that observations could only be reproduced by models assuming either a constant, Salpeter IMF, or a time-dependent distribution, as other IMFs failed. We further tested the models by calculating their M/L ratios. We conclude that a constant, time-independent bottom-heavy IMF does not reproduce the data, especially the increase of the $[alpha/Fe]$ ratio with galactic stellar mass, whereas a variable IMF, switching from top to bottom-heavy, can match observations. For the latter models, the IMF switch always occurs at the earliest possible considered time, i.e. $t_{text{switch}}= 0.1$ Gyr.
We use the relations between aperture stellar velocity dispersion (sigma_ap), stellar mass (M_sps), and galaxy size (R_e) for a sample of sim 150,000 early-type galaxies from SDSS/DR7 to place constraints on the stellar initial mass function (IMF) and dark halo response to galaxy formation. We build LCDM based mass models that reproduce, by construction, the relations between galaxy size, light concentration and stellar mass, and use the spherical Jeans equations to predict sigma_ap. Given our model assumptions (including those in the stellar population synthesis models), we find that reproducing the median sigma_ap vs M_sps relation is not possible with {it both} a universal IMF and a universal dark halo response. Significant departures from a universal IMF and/or dark halo response are required, but there is a degeneracy between these two solutions. We show that this degeneracy can be broken using the strength of the correlation between residuals of the velocity-mass (Delta log sigma_ap) and size-mass (Delta log R_e) relations. The slope of this correlation, d_vr equiv Delta log sigma_ap/Delta log R_e, varies systematically with galaxy mass from d_vr simeq -0.45 at M_sps sim 10^{10}M_sun, to d_vr simeq -0.15 at M_sps sim 10^{11.6} M_sun. The virial fundamental plane (FP) has d_vr=-1/2, and thus we find the tilt of the observed FP is mass dependent. Reproducing this tilt requires {it both} a non-universal IMF and a non-universal halo response. Our best model has mass-follows-light at low masses (Msps < 10^{11.2}M_sun) and unmodified NFW haloes at M_sps sim 10^{11.5} M_sun. The stellar masses imply a mass dependent IMF which is lighter than Salpeter at low masses and heavier than Salpeter at high masses.
Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called small-scale turbulent dynamo and the large-scale mean-field dynamo. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.
The characteristic mass that sets the peak of the stellar initial mass function (IMF) is closely linked to the thermodynamic behaviour of interstellar gas, which controls how gas fragments as it collapses under gravity. As the Universe has grown in metal abundance over cosmic time, this thermodynamic behaviour has evolved from a primordial regime dominated by the competition between compressional heating and molecular hydrogen cooling to a modern regime where the dominant process in dense gas is protostellar radiation feedback, transmitted to the gas by dust-gas collisions. In this paper we map out the primordial-to-modern transition by constructing a model for the thermodynamics of collapsing, dusty gas clouds at a wide range of metallicities. We show the transition from the primordial regime to the modern regime begins at metallicity $Zsim 10^{-4} rm{Z_odot}$, passes through an intermediate stage where metal line cooling is dominant at $Z sim 10^{-3},rm{Z_{odot}}$, and then transitions to the modern dust- and feedback-dominated regime at $Zsim 10^{-2} rm{Z_odot}$. In low pressure environments like the Milky Way, this transition is accompanied by a dramatic change in the characteristic stellar mass, from $sim 50,rm{M_odot}$ at $Z sim 10^{-6},rm{Z_{odot}}$ to $sim 0.3,rm{M_odot}$ once radiation feedback begins to dominate, which marks the appearance of the modern bottom-heavy Milky Way IMF. In the high pressure environments typical of massive elliptical galaxies, the characteristic mass for the modern, dust-dominated regime falls to $sim 0.1,rm{M_{odot}}$, thus providing an explanation for the brown dwarf rich population observed in these galaxies. We conclude that metallicity is a key driver of variations in the characteristic stellar mass, and by extension, the IMF.
The next generation of large aperture ground based telescopes will offer the opportunity to perform accurate stellar photometry in very crowded fields. This future capability will allow one to study in detail the stellar population in distant galaxies. In this paper we explore the effect of photometric errors on the stellar metallicity distribution derived from the color distribution of the Red Giant Branch stars in the central regions of galaxies at the distance of the Virgo cluster. We focus on the analysis of the Color-Magnitude Diagrams at different radii in a typical giant Elliptical galaxy obtained from synthetic data constructed to exemplify observations of the European Extremely Large Telescope. The simulations adopt the specifications of the first light high resolution imager MICADO and the expected performance of the Multi-Conjugate Adaptive Optics Module MAORY. We find that the foreseen photometric accuracy allows us to recover the shape of the metallicity distribution with a resolution $lesssim 0.4$ dex in the inner regions ($mu_{rm B}$ = 20.5 mag arcsec$^{-2}$) and $simeq 0.2$ dex in regions with $mu_{rm B}$ = 21.6 mag arcsec$^{-2}$, that corresponds to approximately half of the effective radius for a typical giant elliptical in Virgo. At the effective radius ($mu_{rm B} simeq 23$ mag arcsec$^{-2}$), the metallicity distribution is recovered with a resolution of $simeq 0.1$ dex. It will thus be possible to study in detail the metallicity gradient of the stellar population over (almost) the whole extension of galaxies in Virgo. We also evaluate the impact of moderate degradations of the Point Spread Function from the assumed optimal conditions and find similar results, showing that this science case is robust.