Slope Of Initial Mass Function Research Articles

Observable and ionizing properties of star-forming galaxies with very massive stars and different initial mass functions

Context. The presence of very massive stars (VMS) with masses > 100 M⊙ is now firmly established in the Local Group, nearby galaxies, and out to cosmological distances. If present, these stars could boost the UV luminosity and ionizing photon production of galaxies, helping to alleviate the overabundance of UV-bright galaxies found with JWST at high redshift. Aims. To examine these questions, we quantify the impact of VMS on properties of integrated stellar populations, exploring different stellar initial mass functions (IMFs) extending up to 400 M⊙, and with slopes between standard (Salpeter-like) and flatter, more top-heavy IMFs. Methods. Combing consistent stellar evolution and atmosphere models tailored to VMS at 1/2.5 Z⊙ metallicity with BPASS evolutionary synthesis models and including nebular emission, we computed integrated spectral energy distributions (SEDs) and derived quantities for a large set of models. Results. We find that VMS contribute significantly to the UV luminosity and Lyman continuum emission of young stellar populations, and they are characterized by strong stellar He II emission, with EW(He IIλ1640) up to 4–8 Å at young ages or ∼2.5 − 4 Å for a constant star formation rate (SFR) (for the IMFs considered here). For IMFs with a Salpeter slope, the boost of the UV luminosity is relatively modest (up to a factor of ∼1.6). However, small changes in the IMF slope (e.g., from α2 = −2.35 to −2) lead to large increases in LUV and the ionizing photon production, Q. The ionizing photon efficiency, ξion, is also increased with VMS, by typically 0.14–0.2 dex for a Salpeter slope, and by up to ∼0.4 dex when the IMF is extended up to 400 M⊙. Stronger H recombination lines are also predicted in the presence of VMS. Interestingly, SEDs including VMS show smaller Lyman breaks, and the shape of the ionizing spectra remain essentially unaltered up to ∼35 eV, but become softer at higher energies. We derive and discuss the maximum values that quantities such as LUV per stellar mass or unit SFR, ξion, Q, and others can reach when VMS are included, and we show that these values become essentially independent of the IMF. We propose observational methods to test for the presence of VMS and constrain the IMF in star-forming galaxies. Finally, using published JWST observations, we examine if high redshift (z ≳ 5 − 6) galaxies show some evidence of the presence of VMS and/or signs of non-standard IMFs. Very top-heavy IMFs can be excluded on average, but we find that the IMF could well extend into the regime of VMS and be flatter than Salpeter in the bulk of high-z galaxies. Conclusions. The predictions should improve our understanding of the stellar content and IMF of the most distant galaxies and allow us to establish if “extreme” or “unusual” populations extending well beyond 100 M⊙ existed in the early Universe.

A-sloth reveals the nature of the first stars

ABSTRACT The first generation of stars (Pop III) are too dim to be observed directly and probably too short-lived to have survived for local observations. Hence, we rely on simulations and indirect observations to constrain the nature of the first stars. In this study, we calibrate the semi-analytical model a-sloth (Ancient Stars and Local Observables by Tracing Haloes), designed for simulating star formation in the early Universe, using a likelihood function based on nine independent observables. These observables span Milky Way-specific and cosmologically representative variables, ensuring a comprehensive calibration process. This calibration methodology ensures that a-sloth provides a robust representation of the early Universe’s star formation processes, aligning simulated values with observed benchmarks across a diverse set of parameters. The outcome of this calibration process is best-fitting values and their uncertainties for 11 important parameters that describe star formation in the early Universe, such as the shape of the initial mass function (IMF) of Pop III stars or escape fractions of ionizing photons. Our best-fitting model has a Pop III IMF with a steeper slope, dN/d$M \propto M^{-1.77}$, than the log-flat models often proposed in the literature, and also relatively high minimum and maximum masses, $M_{\rm min} = 13.6~\, \mathrm{M}_\odot$ and $M_{\rm max} = 197~\, \mathrm{M}_\odot$. However, we emphasize that the IMF-generating parameters are poorly constrained and, e.g. the IMF slope could vary from log-flat to Salpeter. We also provide data products, such as delay time distribution, bubble size distributions for ionizing and metal-enriched bubbles at high redshift, and correlation plots between all 11 input parameters. Our study contributes to understanding the formation of early stars through a-sloth, providing valuable insights into the nature of Pop III stars and the intricate processes involved in the early Universe’s star formation.

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). VI. The High-mass Stellar Initial Mass Function of M33

We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging the Hubble Space Telescope’s high resolving power, we fully model the IMF probabilistically. We first model the optical color–magnitude diagram of each cluster to constrain its power-law slope Γ, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of mass function (MF) slopes for a highly strict cluster sample of nine clusters more massive than log(Mass/M ⊙) = 3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find that the ensemble IMF is best described by a mean high-mass slope of Γ¯=1.49±0.18 , with an intrinsic scatter of σΓ2=0.020.00+0.16 , consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate (SFR) or galactocentric radius within M33, which serves as a proxy for metallicity. This Γ¯ measurement is consistent with similar measurements in M31, despite M33 having a much higher SFR intensity. While this measurement is formally consistent with the canonical Kroupa (Γ = 1.30) IMF, as well as the Salpeter (Γ = 1.35) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical subfields.

INSPIRE: INvestigating Stellar Population In RElics – VI. The low-mass end slope of the stellar initial mass function and chemical composition

ABSTRACT The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.1 < z < 0.4 and obtained their star formation histories (SFHs). Due to their preserved very old stellar populations, relics are the perfect systems to constrain the earliest epochs of mass assembly in the Universe and the formation of massive early-type galaxies. The goal of this work is to investigate whether a correlation exists between the degree of relicness (DoR), quantifying the fraction of stellar mass formed at z > 2, and the other stellar population parameters. We use the Full-Index-Fitting method to fit the INSPIRE spectra to single stellar population (SSP) models. This allows us to measure, for the first time, the slope of the IMF, as well as stellar metallicity [M/H], [Mg/Fe], [Ti/Fe], and [Na/Fe] ratios, and study correlations between them and the DoR. Similarly to normal-sized galaxies, UCMGs with larger stellar masses have overall higher metallicities. We found a correlation between the IMF slope and the DoR, that, however, breaks down for systems with a more extended SFH. An even stronger dependency is found between the IMF and the fraction of mass formed at high-z. At equal velocity dispersion and metallicity, galaxies with a higher DoR have a larger dwarf-to-giant ratio, i.e. a bottom heavy IMF, than that of low-DoR counterparts. This might indicate that the cosmic epoch and therefore different formation scenarios influence the fragmentation of the star formation cloud and hence might be the explanation for IMF variations detected in massive ETGs.

The universal variability of the stellar initial mass function probed by the TIMER survey

The debate about the universality of the stellar initial mass function (IMF) revolves around two competing lines of evidence. While measurements in the Milky Way, an archetypal spiral galaxy, seem to support an invariant IMF, the observed properties of massive early-type galaxies (ETGs) favor an IMF somehow sensitive to the local star-formation conditions. However, the fundamental methodological and physical differences between the two approaches have hampered a comprehensive understanding of IMF variations. Here, we describe an improved modeling scheme that, for the first time, allows consistent IMF measurements across stellar populations with different ages and complex star-formation histories (SFHs). Making use of the exquisite MUSE optical data from the TIMER survey and powered by the MILES stellar population models, we show the age, metallicity, [Mg/Fe], and IMF slope maps of the inner regions of NGC 3351, a spiral galaxy with a mass similar to that of the Milky Way. The measured IMF values in NGC 3351 follow the expectations from a Milky Way-like IMF, although they simultaneously show systematic and spatially coherent variations, particularly for low-mass stars. In addition, our stellar population analysis reveals the presence of metal-poor and Mg-enhanced star-forming regions that appear to be predominantly enriched by the stellar ejecta of core-collapse supernovae. Our findings therefore showcase the potential of detailed studies of young stellar populations to provide the means to better understand the early stages of galaxy evolution and, in particular, the origin of the observed IMF variations beyond and within the Milky Way.

Recovery of the low- and high-mass end slopes of the IMF in massive early-type galaxies using detailed elemental abundances

ABSTRACT Star formation in the early Universe has left its imprint on the chemistry of observable stars in galaxies. We derive elemental abundances and the slope of the low-mass end of the initial mass function (IMF) for a sample of 25 very massive galaxies, separated into brightest cluster galaxies (BCGs) and their massive satellites. The elemental abundances of BGCs and their satellites are similar, but for some elements, satellite galaxies show a correlation with the global velocity dispersion. Using a subset of derived elemental abundances, we model the star formation histories of these galaxies with chemical evolution models, and predict the high-mass end slope of the IMF and star formation time-scales. The high-mass end IMF slope of the satellite galaxies correlates with the global velocity dispersion. The low- and the high-mass end IMF slopes are weakly correlated in a general sense that top heavy IMFs are paired with bottom heavy IMFs. Our results do not necessarily imply that the IMF was simultaneously bottom and top heavy. Instead, our findings can be considered consistent with a temporal variation in the IMF, where, for massive galaxies, the high-mass end IMF slope is representative of the very early age and the low-mass end slope of the later star formation. The small but noticeable differences between the BCGs and the satellites in terms of their elemental abundances and IMF slopes, together with their stellar kinematical properties, suggest somewhat different formation pathways, where BCGs experience more major, gas-free mergers.

Constraining the Initial Mass Function in the Epoch of Reionization from Astrophysical and Cosmological Data

We aim to constrain the stellar initial mass function (IMF) during the epoch of reionization. To this purpose, we build up a semi-empirical model for the reionization history of the Universe based on various ingredients: the latest determination of the UV galaxy luminosity function from JWST out to redshift z≲12; data-inferred and simulation-driven assumptions on the redshift-dependent escape fraction of ionizing photons from primordial galaxies; a simple yet flexible parameterization of the IMF ϕ(m⋆)∼m⋆ξe−m⋆,c/m⋆ in terms of a high-mass end slope ξ<0 and a characteristic mass m⋆,c, below which a flattening or a bending sets in (allowing description of a variety of IMF shapes from the classic Salpeter to top-heavy ones); the PARSEC stellar evolution code to compute the UV and ionizing emission from different stars’ masses as a function of age and metallicity; and a few physical constraints related to stellar and galaxy formation in faint galaxies at the reionization redshifts. We then compare our model outcomes with the reionization observables from different astrophysical and cosmological probes and perform Bayesian inference on the IMF parameters via a standard MCMC technique. We find that the IMF slope ξ is within the range from −2.8 to −2.3, consistent with direct determination from star counts in the Milky Way, while appreciably flatter slopes are excluded at great significance. However, the bestfit value of the IMF characteristic mass m⋆,c∼a few M⊙ implies a suppression in the formation of small stellar masses at variance with the IMF in the local Universe. This may be induced by the thermal background of ∼20–30 K provided by CMB photons at the reionization redshifts. We check that our results are robust against different parameterizations for the redshift evolution of the escape fraction. Finally, we investigate the implications of our reconstructed IMF for the recent JWST detections of massive galaxies at and beyond the reionization epoch, showing that any putative tension with the standard cosmological framework is substantially alleviated.

A statistical and multiwavelength photometric analysis of a young embedded open star cluster: IC 1590

ABSTRACTWe present a statistical and multiwavelength photometric studies of young open cluster IC 1590. We identified 91 cluster members using Gaia DR3 astrometry data using ensemble-based unsupervised machine learning techniques. From Gaia EDR3 data, we estimate the best-fitting parameters for IC 1590 using the Automated Stellar Cluster Analysis package (asteca) yielding the distance d ∼ 2.87 ± 0.02 kpc, age ∼ 3.54 ± 0.05 Myr, metallicity z ∼ 0.0212 ± 0.003, binarity value of ∼ 0.558, and extinction Av ∼ 1.252 ± 0.4 mag for an Rv value of ∼ 3.322 ± 0.23. We estimate the initial mass function slope of the cluster to be α = 1.081 ± 0.112 for single stars and α = 1.490 ± 0.051 for a binary fraction of ∼ 0.558 in the mass range 1 M⊙ ≤ m (M⊙) ≤ 100 M⊙. The G-band luminosity function slope is estimated to be ∼ 0.33 ± 0.09. We use (J − H) versus (H − Ks) colour–colour diagram to identify young stellar objects (YSOs). We found that all the identified YSOs have ages ≤ 2 Myr and masses ∼ 0.35 – 5.5 M⊙. We also fit the radial surface density profile. Using the galpy, we performed orbit analysis of the cluster. The extinction map for the cluster region has been generated using the PNICER technique, and it is almost similar to the dust structure obtained from the 500 μm dust continuum emissions map of Herschel SPIRE. We finally at the end discussed the star formation scenario in the cluster region.

A study of extreme C iii]1908 & [O iii]88/[C ii]157 emission in Pox 186: Implications for JWST+ALMA (FUV+FIR) studies of distant galaxies

Abstract Carbon spectral features are ubiquitous in the ultraviolet (UV) and far-infrared (FIR) spectra of the reionization-era galaxies. We probe the ionized carbon content of a dwarf galaxy Pox 186 using the UV, optical, mid-infrared and FIR data taken with Hubble, Gemini, Spitzer and Herschel, respectively. This local (z∼0.0040705) galaxy is likely an analogue of reionization-era galaxies, as revealed by its extreme FIR emission line ratio, [O iii] 88μm/[C ii] 157μm (>10). The UV spectra reveal extreme C iii] λλ 1907, 1909 emission with the strongest equivalent width (EW) = 35.85 ± 0.73 Å detected so far in the local (z∼0) Universe, a relatively strong C iv λλ 1548, 1550 emission with EW = 7.95 ±0.45 Å, but no He ii λ 1640 detection. Several scenarios are explored to explain the high EW of carbon lines, including high effective temperature, high carbon-to-oxygen ratio, slope and upper mass of top-heavy initial mass function, hard ionizing radiation and in-homogeneous dust distribution. Both C iii] and C iv line profiles are broadened with respect to the O iii] λ 1660 emission line. Each emission line of C iv λλ 1548, 1550 shows the most distinct double-peak structure ever detected, which we model via two scenarios, firstly a double-peaked profile that might emerge from resonant scattering and secondly, a single nebular emission line along with a weaker interstellar absorption. The study demonstrates that galaxies with extreme FIR properties may also show extreme UV properties, hence paving a promising avenue of using FIR+UV in the local (via Hubble+Herschel/SOFIA) and distant (via JWST+ALMA) Universe for unveiling the mysteries of the reionization-era.

Variation of the High-mass Slope of the Stellar Initial Mass Function: Theory Meets Observations

We present observational evidence of a correlation between the high-mass slope of the stellar initial mass function (IMF) in young star clusters and their stellar surface density, σ *. When the high-mass end of the IMF is described by a power law of the form , the value of Γ is seen to decrease weakly with increasing σ *, following a relation. We also present a model that can explain these observations. The model is based on the idea that the coalescence of protostellar cores in a protocluster-forming clump is more efficient in high-density environments where cores are more closely packed. The efficiency of the coalescence process is calculated as a function of the parental clump properties, in particular the relation between its mass and radius as well as its core formation efficiency. The main result of this model is that the increased efficiency of the coalescence process leads to shallower slopes of the IMF, in agreement with the observations of young clusters, and the observations are best reproduced with compact protocluster-forming clumps. These results have significant implications for the shape of the IMF in different Galactic and extragalactic environments and have very important consequences for galactic evolution.

The merger of hard binaries in globular clusters as the primary channel for the formation of second-generation stars

ABSTRACT We have recently presented observational evidence which suggests that the origin of the second-generation (G2) stars in globular clusters (GCs) is due to the binary-mediated collision of primordial (G1) low-mass main-sequence (MS) stars. This mechanism avoids both the mass budget problem and the need of external gas for dilution. Here, we report on another piece of evidence supporting this scenario: (1) the fraction of MS binaries is proportional to the fraction of G1 stars in GCs and, at the same time, (2) the smaller the fraction of G1 stars is, the more deficient binaries of higher mass ratio (q>0.7) are. They are, on average, harder than their smaller mass-ratio counterparts due to higher binding energy at a given primary mass. Then (2) implies that (1) is due to the merging/collisions of hard binaries rather than to their disruption. These new results complemented by the present-day data on binaries lead to the following conclusions: (i) the mass-ratio distribution of binaries, particularly short-period ones, with low-mass primaries, MP < 1.5 M⊙, is strongly peaked close to q=1.0, whereas (ii) dynamical processes at high stellar density tend to destroy softer binaries and make hard (nearly) twin binaries to become even harder and favour their mergers and collisions. G2 stars formed this way gain mass that virtually doubles the primary one, 2MP, at which the number of G1 stars is approximately five times smaller than at MP according to the slope of a Milky Way-like initial mass function at MMS < 1.0 M⊙.

Initial mass function variability from the integrated light of diverse stellar systems

ABSTRACT We present a uniform analysis of the stellar initial mass function (IMF) from integrated light spectroscopy of 15 compact stellar systems (11 globular clusters in M31 and 4 ultra compact dwarfs in the Virgo cluster, UCDs) and two brightest Coma cluster galaxies (BCGs), covering a wide range of metallicities (−1.7 < [Fe/H] < 0.01) and velocity dispersions (7.4 km s−1 <σ < 275 km s−1). The S/N ∼100 Å−1 Keck LRIS spectra are fitted over the range 4000 < λ/Å < 10 000 with flexible full-spectrum stellar population synthesis models. We use the models to fit simultaneously for ages, metallicities, and individual elemental abundances of the population, allowing us to decouple abundance variations from variations in IMF slope. We show that compact stellar systems do not follow the same trends with physical parameters that have been found for early-type galaxies. Most globular clusters in our sample have an IMF consistent with that of the Milky Way, over a wide range of [Fe/H] and [Mg/Fe]. There is more diversity among the UCDs, with some showing evidence for a bottom-heavy IMF, but with no clear correlation with metallicity, abundance, or velocity dispersion. The two Coma BCGs have similar velocity dispersion and metallicity, but we find the IMF of NGC 4874 is consistent with that of the Milky Way while NGC 4889 presents evidence for a significantly bottom-heavy IMF. For this sample, the IMF appears to vary between objects in a way that is not explained by a single metallicity-dependent prescription.

Why should models of dwarf galaxy evolution care about the initial mass function at low star-formation rates?

ABSTRACT When star clusters are formed at low star-formation rates (SFRs), their stellar initial mass function (IMF) can hardly be filled continuously with stars at each mass. This lack holds for massive stars and is verified observationally by the correlation between star-cluster mass and its most massive cluster star. Since galaxy evolution is strongly affected by massive stars, numerical models should account for this lack. Because a filled IMF is mostly applied even when only fractions of massive stars form, here we investigate, by means of 3D chemo-dynamical simulations of isolated dwarf galaxies, how deviations from a standard IMF in star clusters affect the evolution. We compare two different IMF recipes, a filled IMF and one truncated at a maximum mass at which a single complete star forms. Attention is given to energetic and chemical feedback by massive stars. Since their energy release is mass-dependent but steeper than the negative IMF slope, the energetic feedback retains a positive mass dependence, so that a filled IMF regulates star formation (SF) more strongly than truncated IMFs, though only stellar number fractions exist. The higher SFR of the truncated IMF in the simulation leads to more Type II supernovae (SNeII), driving galactic winds. Whether this results from the model-inherent larger SFR is questioned and therefore explored analytically. This shows the expected result for the Lyman continuum, but that the total SNII energy release is equal for both IMF modes, while the power is smaller for the truncated IMF. Reasonably, the different IMFs leave fingerprints in the abundance ratios of massive to intermediate-mass star elements.

The formation of globular clusters with top-heavy initial mass functions

ABSTRACT We study the formation of globular clusters (GCs) in massive compact clouds with the low metallicity of Z = 10−3 Z⊙ by performing three-dimensional radiative-hydrodynamic simulations. Considering the uncertainty of the initial mass function (IMF) of stars formed in low-metallicity and high-density clouds, we investigate the impacts of the IMF on the cloud condition for the GC formation with the range of the power-law index of IMF as γ = 1−2.35. We find that the threshold surface density (Σthr) for the GC formation increases from 800 M⊙ pc−2 at γ = 2.35 to 1600 M⊙ pc−2 at γ = 1.5 in the cases of clouds with Mcl = 106 M⊙ because the emissivity of ionizing photons per stellar mass increases as γ decreases. For γ < 1.5, Σthr saturates with ∼2000 M⊙ pc−2 that is quite rare and observed only in local starburst galaxies due to e.g. merger processes. Thus, we suggest that formation sites of low-metallicity GCs could be limited only in the very high-surface density regions. We also find that Σthr can be modelled by a power-law function with the cloud mass (Mcl) and the emissivity of ionizing photons (s*) as $\propto M_{\rm cl}^{-1/5} s_{*}^{2/5}$. Based on the relation between the power-law slope of IMF and Σthr, future observations with e.g. the JWST can allow us to constrain the IMF of GCs.

ALMA-IMF

Context.Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF.Aims.The ALMA-IMF Large Program aims to investigate the variations in the core distributions (CMF and mass segregation) with cloud characteristics, such as the density and kinematic of the gas, as diagnostic observables of the formation process and evolution of clouds. The present study focuses on the W43-MM2&MM3 mini-starburst, whose CMF has recently been found to be top-heavy with respect to the Salpeter slope of the canonical IMF.Methods.W43-MM2&MM3 is a useful test case for environmental studies because it harbors a rich cluster that contains a statistically significant number of cores (specifically, 205 cores), which was previously characterized in Paper III. We applied a multi-scale decomposition technique to the ALMA 1.3 mm and 3 mm continuum images of W43-MM2&MM3 to define six subregions, each 0.5–1 pc in size. For each subregion we characterized the probability distribution function of the high column density gas,η-PDF, using the 1.3 mm images. Using the core catalog, we investigate correlations between the CMF and cloud and core properties, such as theη-PDF and the core mass segregation.Results.We classify the W43-MM2&MM3 subregions into different stages of evolution, from quiescent to burst to post-burst, based on the surface number density of cores, number of outflows, and ultra-compact HII presence. The high-mass end (>1M⊙) of the subregion CMFs varies from close to the Salpeter slope (quiescent) to top-heavy (burst and post-burst). Moreover, the second tail of theη-PDF varies from steep (quiescent) to flat (burst and post-burst), as observed for high-mass star-forming clouds. We find that subregions with flat secondη-PDF tails display top-heavy CMFs.Conclusions.In dynamical environments such as W43-MM2&MM3, the high-mass end of the CMF appears to be rooted in the cloud structure, which is at high column density and surrounds cores. This connection stems from the fact that cores and their immediate surroundings are both determined and shaped by the cloud formation process, the current evolutionary state of the cloud, and, more broadly, the star formation history. The CMF may evolve from Salpeter to top-heavy throughout the star formation process from the quiescent to the burst phase. This scenario raises the question of if the CMF might revert again to Salpeter as the cloud approaches the end of its star formation stage, a hypothesis that remains to be tested.

The Initial Mass Function and Other Stellar Properties Across the Core of the Hydra I Cluster* * This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory, Chile.

The Hydra I cluster offers an excellent opportunity to study and compare the relic old stellar populations in the core of its two brightest galaxies. In addition, the differing kinematics of the two galaxies allows a test of the local validity of general scaling relations. In this work, we present a direct comparison employing full spectral fitting of new high-quality long-slit optical and near-infrared spectroscopic data. We retrieve age, metallicity, and 19 elemental abundances out to ∼12 kpc within each galaxy, as well as the Initial Mass Function (IMF) in their central regions. Our results suggest that the inner ∼5 kpc regions of both galaxies, despite their different masses, formed at the same time and evolved with a similar star formation timescale and chemical enrichment, confirming their early formation in the cluster buildup. Only the overall metallicity and IMF radial profiles show differences connected with their different velocity dispersion profiles. The radial trends of the IMF positively correlate with both [Z/H] and σ. While the trends of the IMF slope values with metallicity agree with a global trend for both galaxies, the trends with the velocity dispersion exhibit differences. The outer regions show signs of mixed stellar populations with large differences in chemical content compared to the centers, but with similarly old ages.

Multimass modelling of Milky Way globular clusters – I. Implications on their stellar initial mass function above 1 M⊙

ABSTRACT The distribution of stars and stellar remnants (white dwarfs, neutron stars, and black holes) within globular clusters holds clues about their formation and long-term evolution, with important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from various data sets, including proper motions from Gaia EDR3 and HST, line-of-sight velocities from ground-based spectroscopy and deep stellar mass functions from HST. We use metallicity-dependent stellar evolution recipes to obtain present-day mass functions of stars and remnants from the IMF. By dynamically probing the present-day mass function of all objects in a cluster, including the mass distribution of remnants, these models allow us to explore in detail the stellar (initial) mass functions of a large sample of Milky Way GCs. We show that, while the low-mass mass function slopes are strongly dependent on the dynamical age of the clusters, the high-mass slope (α3; m > 1 M⊙) is not, indicating that the mass function in this regime has generally been less affected by dynamical mass loss. Examination of this high-mass mass function slope suggests an IMF in this mass regime consistent with a Salpeter IMF is required to reproduce the observations. This high-mass IMF is incompatible with a top-heavy IMF, as has been proposed recently. Finally, based on multimass model fits to our sample of Milky Way GCs, no significant correlation is found between the high-mass IMF slope and cluster metallicity.

Stellar Initial Mass Function (IMF) Probed with Supernova Rates and Neutrino Background: Cosmic-average IMF Slope Is ≃2–3 Similar to the Salpeter IMF

The stellar initial mass function (IMF) is expressed by ϕ(m) ∝ m − α with the slope α, and known as a poorly constrained but very important function in studies of star and galaxy formation. There are no sensible observational constraints on the IMF slopes beyond the Milky Way and nearby galaxies. Here we combine two sets of observational results, (1) cosmic densities of core-collapse supernova (CCSN) explosion rates and (2) cosmic far-UV radiation (and infrared reradiation) densities, which are sensitive to massive (≃8–50 M ⊙) and moderately massive (≃2.5–7 M ⊙) stars, respectively, and constrain the IMF slope at m > 1 M ⊙ with a freedom of redshift evolution. Although no redshift evolution is identified beyond the uncertainties, we find that the cosmic-average IMF slope at z = 0 is α = 1.8–3.2 at the 95% confidence level that is comparable with the Salpeter IMF, α = 2.35, which marks the first constraint on the cosmic-average IMF. We show a forecast for the Nancy Grace Roman Space Telescope supernova survey that will provide significantly strong constraints on the IMF slope with δ α ≃ 0.5 over z = 0–2. Moreover, as for an independent IMF probe instead of (1), we suggest to use diffuse supernovae neutrino background (DSNB), relic neutrinos from CCSNe. We expect that the Hyper-Kamiokande neutrino observations over 20 yr will improve the constraints on the IMF slope and the redshift evolution significantly better than those obtained today, if the systematic uncertainties of DSNB production physics are reduced in the future numerical simulations.

INSPIRE: INvestigating Stellar Population In RElics – IV. The initial mass function slope in relics

ABSTRACTIn the last decade, growing evidence has emerged supporting a non-universal stellar initial mass function (IMF) in massive galaxies, with a larger number of dwarf stars with respect to the Milky Way (bottom-heavy IMF). However, a consensus about the mechanisms that cause IMF variations is yet to be reached. Recently, it has been suggested that stars formed early-on in cosmic time, via a star formation burst, could be characterized by a bottom-heavy IMF. A promising way to confirm this is to use relics, ultra-compact massive galaxies, almost entirely composed by these ‘pristine’ stars. The INvestigating Stellar Population In RElics (INSPIRE) Project aims at assembling a large sample of confirmed relics, that can serve as laboratory to investigate on the conditions of star formation in the first 1–3 Gyr of the Universe. In this third INSPIRE paper, we build a high signal-to-noise spectrum from five relics, and one from five galaxies with similar sizes, masses, and kinematical properties, but characterized by a more extended star formation history (non-relics). Our detailed stellar population analysis suggests a systematically bottom-heavier IMF slope for relics than for non-relics, adding new observational evidence for the non-universality of the IMF at various redshifts and further supporting the above proposed physical scenario.

How much metal did the first stars provide to the ultra-faint dwarfs?

Numerical simulations of dwarf galaxies have so far failed to reproduce the observed metallicity-luminosity relation, down to the regime of ultra-faint dwarfs (UFDs). We address this issue by exploring how the first generations of metal-free stars (Pop III) could help increase the mean metallicity ([Fe/H]) of those small and faint galaxies. We ran zoom-in chemo-dynamical simulations of 19 halos extracted from a Λ Cold Dark Matter (CDM) cosmological box and followed their evolution down to redshiftz = 0. Models were validated not only on the basis of galaxy global properties, but also on the detailed investigation of the stellar abundance ratios ([α/Fe]). We identified the necessary conditions for the formation of the first stars in mini-halos and derived constraints on the metal ejection schemes. The impact of Pop III stars on the final metallicity of UFDs was evaluated by considering different stellar mass ranges for their initial mass function (IMF), the influence of pair-instability supernovae (PISNe), and their energetic feedback, as well as the metallicity threshold that marks the transition from the first massive stars to the formation of low-mass long-lived stars. The inclusion of Pop III stars with masses below 140 M⊙, and a standard IMF slope of −1.3 does increase the global metallicity of UFDs, although these are insufficient to resolve the tension with observations. The PISNe with progenitor masses above 140 M⊙do allow the metal content of UFDs to further increase. However, as PISNe are very rare and sometimes absent in the faintest UFDs, they have a limited impact on the global faint end of the metallicity-luminosity relation. Despite a limited number of spectroscopically confirmed members in UFDs, which make the stellar metallicity distribution of some UFDs uncertain, our analysis reveals that this is essentially the metal-rich tail that is missing in the models. The remaining challenges are thus both observational and numerical: (i) to extend high-resolution spectroscopy data samples and confirm the mean metallicity of the faintest UFDs; and (ii) to explain the presence of chemically enriched stars in galaxies with very short star formation histories.