Initial Mass Function Research Articles

First Detection of CO Isotopologues in a High-redshift Main-sequence Galaxy: Evidence of a Top-heavy Stellar Initial Mass Function

Recent observations and theories have presented a strong challenge to the universality of the stellar initial mass function (IMF) in extreme environments. A notable example has been found for starburst conditions, where evidence favors a top-heavy IMF, i.e., there is a bias toward massive stars compared to the IMF that is responsible for the stellar mass function and elemental abundances observed in the Milky Way. Local starburst galaxies have star formation rates similar to those in high-redshift main-sequence galaxies, which appear to dominate the stellar mass budget at early epochs. However, the IMF of high-redshift main-sequence galaxies is yet to be probed. Since 13CO and C18O isotopologues are sensitive to the IMF, we have observed these lines toward four strongly lensed high-redshift main-sequence galaxies using the Atacama Large Millimeter/submillimeter Array. Of our four targets, SDSS J0901+1814, at z ≈ 2.26, is seen clearly in 13CO and C18O, the first detection of CO isotopologues in the high-redshift main-sequence galaxy population. The observed 13C/18O ratio, 2.4 ± 0.8, is significantly lower than that of local main-sequence galaxies. We estimate the isotope ratio, oxygen abundance, and stellar mass using a series of chemical evolution models with varying star formation histories and IMFs. All models favor an IMF that is more top-heavy than that of the Milky Way. Thus, as with starburst galaxies, main-sequence galaxies in the high-redshift Universe have a greater fraction of massive stars than a Milky Way IMF would imply.

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.

Herschel Gould Belt Survey in Taurus – II. A census of dense cores and filaments in the TMC1 region

ABSTRACT We present a catalogue of dense cores and filaments in a $3.8^\circ \times 2.4^\circ$ field around the TMC1 region of the Taurus molecular cloud. The catalogue was created using photometric data from the Herschel SPIRE and PACS instruments in the 70, 160, 250, 350, and 500 μm continuum bands. Extended structure in the region was reconstructed from a Herschel column density map. Power spectra and probability density functions (PDFs) of this structure are presented. The PDF splits into lognormal and power-law forms, with the high-density power-law component associated primarily with the central part of TMC1. The total mass in the mapped region is 2000 M$_\odot$, of which 34 per cent is above an extinction of $A_V\sim 3$ mag – a level that appears as a break in the PDF and as the minimum column density at which dense cores are found. A total of 35 dense filaments were extracted from the column density map. These have a typical full width at half-maximum (FWHM) width of 0.07 pc, but the TMC1 filament itself has a mean FWHM of $\sim 0.13$ pc. The thermally supercritical filaments in the region are aligned orthogonal to the prevailing magnetic field direction. Derived properties for the supercritical TMC1 filament support the scenario of it being relatively young. A catalogue of 44 robust and candidate prestellar cores is created and is assessed to be complete down to 0.1 M$_\odot$. The combined prestellar core mass function for the TMC1 and L1495 regions is well fit by a single lognormal distribution and is comparable to the standard initial mass function.

The Variation in the Galaxy-wide Initial Mass Function for Low-mass Stars: Modeling and Observational Insights

The stellar initial mass function (IMF) characterizes the mass distribution of newly formed stars in various cosmic environments, serving as a fundamental assumption in astrophysical research. Recent findings challenge the prevalent notion of a universal and static IMF, proposing instead that the IMF’s shape is contingent upon the star formation environment. In this study, we analyze the galaxy-wide variation in the IMF for low-mass stars in both dwarf and massive galaxies with diverse observational methods. Despite systematic discrepancies between different approaches, an IMF model with a metallicity-dependent slope for the low-mass stars aligns with the majority of observations, indicating a high degree of uniformity in the star formation processes across the Universe. We also emphasize the need for a more comprehensive understanding of the variation in the low-mass IMF, considering measurement biases and factors beyond metallicity.

Mg/Fe] and variable initial mass function: Revision of [α/Fe] for massive galaxies

Observations of nearby massive galaxies have revealed that they are older and richer in metals and magnesium than their low-mass counterparts. In particular, the overabundance of magnesium compared to iron, [Mg/Fe], is interpreted to reflect the short star formation history that the current massive galaxies underwent early in the Universe. We present a systematic revision of the [Mg/Fe] – velocity dispersion (σ) relation based on stacked spectra of early-type galaxies with a high signal-to-noise ratio from the Sloan Digital Sky Survey. Using the penalized pixel-fitting (pPXF) method and the MILES single stellar population models, we fit a wide optical wavelength range to measure the net α-abundance. The combination of pPXF and α-enhanced MILES models incorrectly leads to an apparently decreasing trend of [α/Fe] with velocity dispersion. We interpret this result as a consequence of variations in the individual abundances of the different α-elements. This warrants caution for a naive use of full spectral fitting algorithms paired with stellar population models that do not take individual elemental abundance variations into account, especially when deriving averaged quantities such as the mean [α/Fe] of a stellar population. In addition, and based on line-strength measurements, we quantify the impact of a non-universal initial mass function on the recovered abundance pattern of galaxies. In particular, we find that a simultaneous fit of the slope of the initial mass function and the [Mg/Fe] results in a shallower [Mg/Fe]–σ relation. Therefore, our results suggest that star formation in massive galaxies lasted longer than what has been reported previously, although it still occurred significantly faster than in the solar neighbourhood.

Exploring the Evolution of a Dwarf Spheroidal Galaxy with Smoothed-particle Hydrodynamics Simulations. I. Stellar Feedback

A fundamental question regarding the evolution of dwarf spheroidal galaxies is the identification of the key physical mechanisms responsible for gas depletion. Here, we focus on the study of stellar feedback in isolated dwarf spheroidal galaxies by performing numerical simulations using a modified version of the smoothed-particle hydrodynamics code GADGET-3. The Milky Way satellite Leo II (PGC 34176) in the Local Group was considered as our default model dwarf galaxy. The parameter space for the stellar feedback models was explored to match observational constraints of Leo II, such as residual gas mass, total mass within the tidal radius, star formation history, final stellar mass, stellar ages, and metallicity. Additionally, we examined the impact of the binary fraction of stars, initial mass function, dark matter halo mass, and initial gas reservoir. Many simulations revealed recent star formation quenching due to stellar feedback. In general, the gas depletion, expected star formation history, total mass of stars, and total mass within the tidal radius were adequately reproduced in the simulations when compared to observational estimates. However, there were discrepancies in the distribution of stellar ages and metallicities, which suggested that the cosmic gas infall would play a more complex role in our dwarf spheroidal galaxy than captured by a monolithic infall scenario. Our results suggest that currently quenched dwarf galaxies may not necessarily need to evolve within clusters or groups and that stellar feedback alone could be a sufficient factor in shaping at least some of these galaxies as we observe them today.

ALMA-IMF

Context. A crucial aspect in addressing the challenge of measuring the core mass function (CMF), that is pivotal for comprehending the origin of the initial mass function (IMF), lies in constraining the temperatures of the cores. Aims. We aim to measure the luminosity, mass, column density and dust temperature of star-forming regions imaged by the ALMA-IMF large program. These fields were chosen to encompass early evolutionary stages of massive protoclusters. High angular resolution mapping is required to capture the properties of protostellar and pre-stellar cores within these regions, and to effectively separate them from larger features, such as dusty filaments. Methods. We employed the point process mapping (PPMAP) technique, enabling us to perform spectral energy distribution fitting of far-infrared and submillimeter observations across the 15 ALMA-IMF fields, at an unmatched 2.5″ angular resolution. By combining the modified blackbody model with near-infrared data, we derived bolometric luminosity maps. We estimated the errors impacting values of each pixel in the temperature, column density, and luminosity maps. Subsequently, we employed the extraction algorithm getsf on the luminosity maps in order to detect luminosity peaks and measure their associated masses. Results. We obtained high-resolution constraints on the luminosity, dust temperature, and mass of protoclusters, that are in agreement with previously reported measurements made at a coarser angular resolution. We find that the luminosity-to-mass ratio correlates with the evolutionary stage of the studied regions, albeit with intra-region variability. We compiled a PPMAP source catalog of 313 luminosity peaks using getsf on the derived bolometric luminosity maps. The PPMAP source catalog provides constraints on the mass and luminosity of protostars and cores, although one source may encompass several objects. Finally, we compare the estimated luminosity-to-mass ratio of PPMAP sources with evolutionary tracks and discuss the limitations imposed by the 2.5″ beam.

The COSMOS-Web ring: In-depth characterization of an Einstein ring lensing system at z ∼ 2

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z ≈ 2 that we serendipitously discovered during the data reduction of the COSMOS-Web survey and that could be the most distant lens discovered to date. Methods. We extracted the visible and near-infrared photometry of the source and the lens from more than 25 bands. We combined these observations with far-infrared detections to study the dusty nature of the source and we derived the photometric redshifts and physical properties of both the lens and the source with three different spectral energy distribution (SED) fitting codes. Using JWST/NIRCam images, we also produced two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. We find the lens to be a very massive elliptical galaxy at z = 2.02 ± 0.02 with a total mass within the Einstein radius of Mtot(<θEin = (3.66 ± 0.36) × 1011 M⊙ and a total stellar mass of M⋆ = 1.37−0.11+0.14 × 1011 M⊙. We also estimate it to be compact and quiescent with a specific star formation rate below 10−13 yr. Compared to stellar-to-halo mass relations from the literature, we find that the total mass of the lens within the Einstein radius is consistent with the presence of a dark matter (DM) halo of total mass Mh = 1.09−0.57+1.46 × 1013 M⊙. In addition, the background source is a M⋆ = (1.26 ± 0.17) × 1010 M⊙ star-forming galaxy (SFR ≈ (78 ± 15) M⊙ yr) at z = 5.48 ± 0.06. The morphology reconstructed in the source plane shows two clear components with different colors. Dust attenuation values from SED fitting and nearby detections in the far infrared also suggest that the background source could be at least partially dust-obscured. Conclusions. We find the lens at z ≈ 2. Its total, stellar, and DM halo masses are consistent within the Einstein ring, so we do not need any unexpected changes in our description of the lens such as changing its initial mass function or including a non-negligible gas contribution. The most likely solution for the lensed source is at z ≈ 5.5. Its reconstructed morphology is complex and highly wavelength dependent, possibly because it is a merger or a main sequence galaxy with a heterogeneous dust distribution.

JADES

Finding the first generation of stars formed out of pristine gas in the early Universe, known as Population III (PopIII) stars, is one of the most important goals of modern astrophysics. Recent models have suggested that PopIII stars may form in pockets of pristine gas in the halo of more evolved galaxies. We present NIRSpec integral field spectroscopy and micro-shutter array spectroscopic observations of the region around GN-z11, an exceptionally luminous galaxy at z = 10.6, that reveal a greater than 5σ detection of a feature consistent with being HeIIλ1640 emission at the redshift of GN-z11. The very high equivalent width of the putative HeII emission in this clump (log(EWrest(HeII)/Å) = 1.79−0.25+0.15) and a lack of metal lines can be explained in terms of photoionisation by PopIII stars, while photoionisation by PopII stars is inconsistent with the data. The high equivalent width would also indicate that the putative PopIII stars likely have an initial mass function with an upper cutoff reaching at least 500 M⊙. The PopIII bolometric luminosity inferred from the HeII line would be ∼7 × 109 L⊙, which would imply a total stellar mass formed in the burst of ∼2 × 105 M⊙. We find that photoionisation by the active galactic nucleus (AGN) in GN-z11 cannot account for the HeII luminosity observed in the clump but can potentially be responsible for an additional HeII emission observed closer to GN-z11. We also consider the possibility of in situ photoionisation by an accreting direct collapse black hole hosted by the HeII clump. We find that this scenario is less favoured, but it remains a possible alternative interpretation. We also report the detection of a Lyα halo stemming out of GN-z11 and extending out to ∼2 kpc as well as resolved funnel-shaped CIII emission likely tracing the ionisation cone of the AGN.

Intermediate-mass Black Hole Progenitors from Stellar Collisions in Dense Star Clusters

Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes, which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies have investigated how different cluster initial conditions affect the formation of a VMS, including mass segregation, stellar collisions, and binaries, among others. In this study, we investigate the growth of VMSs with a new grid of Cluster Monte Carlo star cluster simulations—the most expansive to date. The simulations span a wide range of initial conditions, varying the number of stars, cluster density, stellar initial mass function (IMF), and primordial binary fraction. We find a gradual shift in the mass of the most massive collision product across the parameter space; in particular, denser clusters born with top-heavy IMFs provide strong collisional regimes that form VMSs with masses easily exceeding 1000 M ⊙. Our results are used to derive a fitting formula that can predict the typical mass of a VMS formed as a function of the star cluster properties. Additionally, we study the stochasticity of this process and derive a statistical distribution for the mass of the VMS formed in one of our models, recomputing the model 50 times with different initial random seeds.

Surveying the Onset and Evolution of Supermassive Black Holes at High-z with AXIS

The nature and origin of supermassive black holes (SMBHs) remain an open matter of debate within the scientific community. While various theoretical scenarios have been proposed, each with specific observational signatures, the lack of sufficiently sensitive X-ray observations hinders the progress of observational tests. In this white paper, we present how AXIS will contribute to solving this issue. With an angular resolution of 1.5″ on-axis and minimal off-axis degradation, we designed a deep survey capable of reaching flux limits in the [0.5–2] keV range of approximately 2 × 10−18 erg s−1 cm−2 over an area of 0.13 deg2 in approximately 7 million seconds (7 Ms). Furthermore, we planned an intermediate depth survey covering approximately 2 deg2 and reaching flux limits of about 2 × 10−17 erg s−1 cm−2 in order to detect a significant number of SMBHs with X-ray luminosities (LX) of approximately 1042 erg s−1 up to z∼10. These observations will enable AXIS to detect SMBHs with masses smaller than 105 M⊙, assuming Eddington-limited accretion and a typical bolometric correction for Type II AGN. AXIS will provide valuable information on the seeding and population synthesis models of SMBHs, allowing for more accurate constraints on their initial mass function (IMF) and accretion history from z∼0–10. To accomplish this, AXIS will leverage the unique synergy of survey telescopes such as the JWST, Roman, Euclid, Vera Rubin Telescope, and the new generation of 30 m class telescopes. These instruments will provide optical identification and redshift measurements, while AXIS will discover the smoking gun of nuclear activity, particularly in the case of highly obscured AGN or peculiar UV spectra as predicted and recently observed by the JWST in the early Universe.

Stellar Population near NGC 2021: Procession of Star Formation in the South Rim of Supergiant Shell LMC 4

Supergiant shells (SGSs) are the largest interstellar structures where heated and enriched gas flows into the host galaxy’s halo. The SGSs in the Large Magellanic Cloud (LMC) are so close that their stars can be resolved with ground-based telescopes to allow studies of star-formation history. Aiming to study the star formation history and energy budget of LMC 4, we have conducted a pilot study of the cluster NGC 2021 and the OB associations in its vicinity near the south rim of LMC 4. We use the Magellanic Cloud Photometric Survey data of the LMC to establish a methodology to examine the stellar population and assess the massive star formation history. We find a radial procession of massive star formation from the northwest part of the OB association LH79 through NGC 2021 to the OB association LH78 in the south. Using the stellar content of NGC 2021 and the assumption of Salpeter’s initial mass function, we estimate that ∼4 supernovae have occurred in NGC 2021, injecting at least 4 × 1051 erg of kinetic energy into the interior of LMC 4.

No Top-heavy Stellar Initial Mass Function Needed: The Ionizing Radiation of GS9422 Can Be Powered by a Mixture of an Active Galactic Nucleus and Stars

JWST is producing high-quality rest-frame optical and UV spectra of faint galaxies at z > 4 for the first time, challenging models of galaxy and stellar populations. One galaxy recently observed at z = 5.943, GS9422, has nebular line and UV continuum emission that appears to require a high ionizing photon production efficiency. This has been explained with an exotic stellar initial mass function (IMF; Cameron et al. 2023a), 10–30x more top-heavy than a Salpeter IMF. Here we suggest an alternate explanation to this exotic IMF. We use a new flexible neural net emulator for CLOUDY, Cue, to infer the shape of the ionizing spectrum directly from the observed emission line fluxes. By describing the ionizing spectrum with a piecewise power law, Cue is agnostic to the source of the ionizing photons. Cue finds that the ionizing radiation from GS9422 can be approximated by a double power law characterized by QHeIIQH=−1.5 , which can be interpreted as a combination of young, metal-poor stars and a low-luminosity active galactic nucleus with F ν ∝ λ 2 in a 65%/35% ratio. This suggests a significantly lower nebular continuum contribution to the observed UV flux (24%) than a top-heavy IMF (≳80%), and hence, necessitates a damped Lyα absorber to explain the continuum turnover bluewards of ∼1400 Å. While current data cannot rule out either scenario, given the immense impact the proposed top-heavy IMF would have on models of galaxy formation, it is important to propose viable alternative explanations and to further investigate the nature of peculiar high-z nebular emitters.

Nebular dominated galaxies: insights into the stellar initial mass function at high redshift

ABSTRACT We identify a low-metallicity ($12+\log ({\rm O}/{\rm H})=7.59$) Ly $\alpha$-emitting galaxy at $z=5.943$ with evidence of a strong Balmer jump, arising from nebular continuum. While Balmer jumps are sometimes observed in low-redshift star-forming galaxies, this galaxy also exhibits a steep turnover in the UV continuum. Such turnovers are typically attributed to absorption by a damped Ly $\alpha$ system (DLA); however, the shape of the turnover and the high observed Ly $\alpha$ escape fraction ($f_{\rm esc,Ly\alpha }~\sim 27~{{\ \rm per\ cent}}$) is also consistent with strong nebular two-photon continuum emission. Modelling the UV turnover with a DLA requires extreme column densities ($N_{\rm HI}\,\,\gt\,\, 10^{23}$ cm$^{-2}$), and simultaneously explaining the high $f_{\rm esc,Ly\alpha }$ requires a fine-tuned geometry. In contrast, modelling the spectrum as primarily nebular provides a good fit to both the continuum and emission lines, motivating scenarios in which (a) we are observing only nebular emission or (b) the ionizing source is powering extreme nebular emission that outshines the stellar emission. The nebular-only scenario could arise if the ionizing source has ‘turned off’ more recently than the recombination time-scale ($\sim$1000 yr), hence we may be catching the object at a very specific time. Alternatively, hot stars with $T_{\rm eff}\gtrsim 10^5$ K (e.g. Wolf–Rayet or low-metallicity massive stars) produce enough ionizing photons such that the two-photon emission becomes visible. While several stellar SEDs from the literature fit the observed spectrum well, the hot-star scenario requires that the number of $\gtrsim 50~{\rm M}_\odot$ stars relative to $\sim 5\!-\!50~{\rm M}_\odot$ stars is significantly higher than predicted by typical stellar initial mass functions (IMFs). The identification of more galaxies with similar spectra may provide evidence for a top-heavy IMF at high redshift.

The Complete CEERS Early Universe Galaxy Sample: A Surprisingly Slow Evolution of the Space Density of Bright Galaxies at z ∼ 8.5–14.5

We present a sample of 88 candidate z ∼ 8.5–14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science survey. These data cover ∼90 arcmin2 (10 NIRCam pointings) in six broadband imaging filters and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z ∼ 9, 11, and 14 and show that the space density of bright (M UV = −20) galaxies changes only modestly from z ∼ 14 to z ∼ 9, compared to a steeper increase from z ∼ 8 to z ∼ 4. While our candidates are photometrically selected, spectroscopic follow-up has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z > 10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models that naturally have enhanced star formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic follow-up of a large sample of early galaxies can distinguish between these competing scenarios.

ASTRAEUS

Context. Observations with the James Webb Space Telescope (JWST) have revealed an abundance of bright z > 10 galaxy candidates, challenging the predictions of most theoretical models at high redshifts. Aims. Since massive stars dominate the observable ultraviolet (UV) emission, we explore whether a stellar initial mass function (IMF) that becomes increasingly top-heavy towards higher redshifts and lower gas-phase metallicities results in a higher abundance of bright objects in the early universe and how it influences the evolution of galaxy properties compared to a constant Salpeter IMF. Methods. We parameterised the IMF based on the findings from hydrodynamical simulations that track the formation of stars in differently metal-enriched gas clouds in the presence of the cosmic microwave background (CMB) at different redshifts. We incorporated this evolving IMF into the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionisation in N-body dArk mattEr simUlationS) framework, which couples galaxy evolution and reionisation in the first billion years. Our implementation accounts for the IMF dependence of supernova (SN) feedback, metal enrichment, and ionising and UV radiation emission. We conducted two simulations: one with a Salpeter IMF and the other with the evolving IMF. In both, we adjusted the free model parameters to reproduce key observables. Results. Compared to a constant Salpeter IMF, we find that (i) the higher abundance of massive stars in the evolving IMF results in more light per unit stellar mass, resulting in a slower build-up of the stellar mass and lower stellar-to-halo mass ratio; (ii) due to the self-similar growth of the underlying dark matter (DM) halos, the evolving IMF’s star formation main sequence scarcely deviates from that of the Salpeter IMF; (iii) the evolving IMF’s stellar mass to gas-phase metallicity relation shifts to higher metallicities, while its halo mass to gas-phase metallicity relation remains unchanged; (iv) the evolving IMF’s median dust-to-metal mass ratio is lower due to its stronger SN feedback; and (v) the evolving IMF requires lower values of the escape fraction of ionising photons and exhibits a flatter median relation and smaller scatter between the ionising photons emerging from galaxies and the halo mass. However, the ionising emissivities of the galaxies mainly driving reionisation (Mh ∼ 1010 M⊙) are comparable to those of a Salpeter IMF, resulting in minimal changes to the topology of the ionised regions. Conclusions. These results suggest that a top-heavier IMF alone is unlikely to explain the higher abundance of bright z > 10 sources, since the lower mass-to-light ratio driven by the greater abundance of massive stars is counteracted by stronger stellar feedback.

Variation of the stellar initial mass function in semi-analytical models

Context. In our previous work, we derived the CR-IGIMF, which is a new scenario for a variable stellar initial mass function (IMF) that combines numerical results on the role played by cosmic rays (CR) in setting the thermal state of star-forming gas with the analytical approach of the integrated galaxy-wide IMF (IGIMF). Aims. In this work, we study the implications of this scenario for the properties of local early-type galaxies (ETG) as inferred from dynamical, photometric, and spectroscopic studies. Methods. We implemented a library of CR-IGIMF shapes in the framework of the galaxy evolution and assembly (GAEA) model. GAEA provides predictions for the physical and photometric properties of model galaxies and for their chemical composition. Our realization includes a self-consistent derivation of the synthetic spectral energy distribution for each model galaxy, which allows a direct derivation of the mass fraction in the mean IMF of low-mass stars (i.e., the dwarf-to-giant ratio, fdg) and a comparison with IMF-sensitive spectral features. Results. The predictions of the GAEA model implementing the CR-IGIMF confirm our previous findings: It correctly reproduces both the observed excess of z ∼ 0 dynamical mass (mass-to-light ratio) with respect to spectroscopic (photometric) estimates assuming a universal MW-like IMF, and the observed increase in [α/Fe] ratios with stellar mass in spheroidal galaxies. Moreover, this realization reproduces the increasing trends of fdg and IMF-sensitive line strengths with velocity dispersion, although the predicted relations are significantly shallower than the observed ones. Conclusions. Our results show that the CR-IGIMF is a promising scenario that reproduces at the same time dynamical, photometric, and spectroscopic indications of a varying IMF in local ETGs. The shallow relations found for spectral indices suggest that either a stronger variability as a function of galaxy properties or additional dependences (e.g., as a function of star forming gas metallicity) might be required to match the strength of the observed trends.

ALMA-IMF

Context. One of the central questions in astrophysics is the origin of the initial mass function (IMF). It is intrinsically linked to the processes from which it originates, and hence its connection with the core mass function (CMF) must be elucidated. Aims. We aim to measure the CMF in the evolved W33-Main star-forming protocluster to compare it with CMF recently obtained in other Galactic star-forming regions, including the ones that are part of the ALMA-IMF program. Methods. We used observations from the ALMA-IMF large programme: ~2′ × 2′ maps of emission from the continuum and selected lines at 1.3 mm and 3 mm observed by the ALMA 12m only antennas. Our angular resolution was typically 1″, that is, ~2400 au at a distance of 2.4 kpc. The lines we analysed are CO (2–1), SiO (5–4), N2H+ (1–0), H41α as well as He41α blended with C41α. We built a census of dense cores in the region, and we measured the associated CMF based on a core-dependent temperature value. Results. We confirmed the ‘evolved’ status of W33-Main by identifiying three H II regions within the field, and to a lesser extent based on the number and extension of N2H+ filaments. We produced a filtered core catalogue of 94 candidates that we refined to take into account the contamination of the continuum by free-free and line emission, obtaining 80 cores with masses that range from 0.03 to 13.2 M⊙. We fitted the resulting high-mass end of the CMF with a single power law of the form N(log(M)) ∝ Mα, obtaining α = −1.44−0.22+0.16, which is slightly steeper but consistent with the Salpeter index. We categorised our cores as prestellar and protostellar, mostly based on outflow activity and hot core nature. We found the prestellar CMF to be steeper than a Salpeter-like distribution, and the protostellar CMF to be slightly top heavy. We found a higher proportion of cores within the H II regions and their surroundings than in the rest of the field. We also found that the cores’ masses were rather low (maximum mass of ~13 M⊙). Conclusions. We find that star formation in W33-Main could be compatible with a ‘clump-fed’ scenario of star formation in an evolved cloud characterised by stellar feedback in the form of H II regions, and under the influence of massive stars outside the field. Our results differ from those found in less evolved young star-forming regions in the ALMA-IMF program. Further investigations are needed to elucidate the evolution of late CMFs towards the IMF over statistically significant samples.

Evidence for very massive stars in extremely UV-bright star-forming galaxies at z ∼ 2.2–3.6

We present a comprehensive analysis of the presence of very massive stars (VMS > 100 M⊙) in the integrated spectra of 13 UV-bright star-forming galaxies at 2.2 ≲ z ≲ 3.6 taken with the Gran Telescopio Canarias (GTC). These galaxies have very high UV absolute magnitudes (MUV ≃ −24), intense star formation (star formation rate ≃100 − 1000 M⊙ yr−1), and metallicities in the range of 12 + log(O/H) ≃ 8.10 − 8.50 inferred from strong rest-optical lines. The GTC rest-UV spectra reveal spectral features indicative of very young stellar populations with VMS, such as strong P-Cygni line profiles in the wind lines N Vλ1240 and C IVλ1550 along with intense and broad He IIλ1640 emission with equivalent width (EW0) ≃ 1.40 − 4.60 Å, and full width half maximum (FWHM) ≃1150 − 3170 km s−1. A Comparison with known VMS-dominated sources and typical galaxies without VMS reveals that some UV-bright galaxies closely resemble VMS-dominated clusters (e.g., R136 cluster). The presence of VMS is further supported by a quantitative comparison of the observed strength of the He II emission with population synthesis models with and without VMS, where models with VMS are clearly preferred. Employing an empirical threshold for EW0 (He II) ≥ 3.0 Å, along with the detection of other VMS-related spectral profiles (N IVλ1486, 1719), we classify nine out of 13 UV-bright galaxies as VMS-dominated sources. This high incidence of VMS-dominated sources in the UV-bright galaxy population (≈70%) contrasts significantly with the negligible presence of VMS in typical LUV∗ LBGs at similar redshifts (< 1%). Our results thus indicate that VMS are common in UV-bright galaxies, suggesting a different initial mass function (IMF) with upper mass limits between 175 M⊙ and 475 M⊙.

The Low-mass Stellar Initial Mass Function in Nearby Ultrafaint Dwarf Galaxies

The stellar initial mass function (IMF) describes the distribution of stellar masses that form in a given star formation event. The long main-sequence lifetimes of low-mass stars mean that the IMF in this regime (below ∼ 1 M ⊙) can be investigated through star counts. Ultrafaint dwarf galaxies are low-luminosity systems with ancient, metal-poor stellar populations. We investigate the low-mass IMF in four such systems (Reticulum II, Ursa Major II, Triangulum II, and Segue 1), using Hubble Space Telescope imaging data that reaches to ≲ 0.2 M ⊙ in each galaxy. The analysis techniques that we adopt depend on the number of low-mass stars in each sample. We use Kolmogorov–Smirnov tests for all four galaxies to determine whether their observed apparent magnitude distributions can reject a given combination of IMF parameters and binary fraction for the underlying population. We forward model 1000 synthetic populations for each combination of parameters, and reject those parameters only if each of the 1000 realizations reject the null hypothesis. We find that all four galaxies reject a variety of IMFs, and the IMFs that they cannot reject include those that are identical, or similar, to that of the stellar populations of the Milky Way. We determine the best-fit parameter values for the IMF in Reticulum II and Ursa Major II and find that the IMF in Reticulum II is generally consistent with that of the Milky Way, while the IMF in Ursa Major II is more bottom heavy. The interpretation of the results for Ursa Major II is complicated by possible contamination from two known background galaxy clusters.