Initial Mass Function Research Articles

Exponential amplification of the magnetic field in the primordial star-forming cloud

In the study of the initial mass function associated with the first generation of stars, known as Population III (Pop III) stars, a fundamental yet unresolved question pertains to the ultimate destiny of the secondary protostars emerging within the accretion disk – specifically, their likelihood of either merging or persisting as distinct entities. Our research concentrates on the magnetic influences affecting the genesis of these first stars under the conditions set by the cosmological initial magnetic field strength. We employ ideal magnetohydrodynamic simulations, utilizing a stiff equation-of-state (EOS) model, to accurately depict the magnetic field structure interconnecting these protostars. We observe that the magnetic field experiences rapid intensification due to the gas near the protostar completing multiple tens of orbital rotations in the initial decade following the formation of the protostar. Concurrently, as mass accretion continues, the region influenced by the significant magnetic field expands outward. This process of magnetic braking effectively curtails the disk fragmentation that would typically occur without a magnetic field. The resulting exponential augmentation of the magnetic field is posited to facilitate the formation of supermassive first stars.

Cold Dark Matter and Self-interacting Dark Matter Interpretations of the Strong Gravitational Lensing Object JWST-ER1

van Dokkum et al. reported the discovery of JWST-ER1, a strong lensing object at redshift z ≈ 2, using data from the James Webb Space Telescope. The lens mass within the Einstein ring is 5.9 times higher than the expected stellar mass from a Chabrier initial mass function, indicating a high dark matter density. In this work, we show that a cold dark matter halo, influenced by gas-driven adiabatic contraction, can account for the observed lens mass. We interpret the measurement of JWST-ER1 in the self-interacting dark matter scenario and show that the cross section per particle mass σ/m ≈ 0.1 cm2 g−1 is generally favored. Intriguingly, σ/m ≈ 0.1 cm2 g−1 can also be consistent with the strong lensing observations of early-type galaxies at redshift z ≈ 0.2, where adiabatic contraction is not observed overall.

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.

Star clusters forming in a low-metallicity starburst – rapid self-enrichment by (very) massive stars

ABSTRACT Stellar winds of massive ($\gtrsim 9\, \mathrm{M_\odot }$) and very massive ($\gtrsim 100\, \mathrm{M_\odot }$) stars may play an important role in the metal-enrichment during the formation of star clusters. With novel high-resolution hydrodynamical griffin-project simulations, we investigate the rapid recycling of stellar wind-material during the formation of massive star clusters up to $M_\mathrm{cluster}\sim 2\times 10^5\, \mathrm{M_\odot }$ in a low-metallicity dwarf galaxy starburst. The simulation realizes new stars from a stellar initial mass function (IMF) between $0.08$ and $\sim 400\, \mathrm{M_\odot }$ and follows stellar winds, radiation and supernova-feedback of single massive stars with evolution tracks. Star clusters form on time-scales less than ∼5 Myr, and their supernova-material is very inefficiently recycled. Stellar wind-material, however, is trapped in massive clusters resulting in the formation of stars self-enriched in Na, Al, and N within only a few Myr. Wind-enriched (second population, 2P) stars can be centrally concentrated in the most massive clusters ($\gtrsim 10^4\, \mathrm{M_\odot }$) and the locked wind-material increases approximately as $M_\mathrm{cluster}^{2}$. These trends resemble the characteristics of observed 2P stars in globular clusters (GCs). We fit scaling relations to the lognormal distributed wind-mass fractions and extrapolate to possible GC progenitors of $M_\mathrm{cluster}=10^7\, \mathrm{M_\odot }$ to investigate whether a dominant 2P could form. This can only happen if the IMF is well-sampled, single massive stars produce at least a factor of a few more enriched winds, for example, through a top-heavy IMF, and a significant fraction of the first population (unenriched) stars is lost during cluster evolution.

Simulating emission line galaxies for the next generation of large-scale structure surveys

ABSTRACT We investigate emission line galaxies across cosmic time by combining the modified L-Galaxies semi-analytical galaxy formation model with the JiuTian cosmological simulation. We improve the tidal disruption model of satellite galaxies in L-Galaxies to address the time dependence problem. We utilize the public code cloudy to compute emission line ratios for a grid of H ii region models. The emission line models assume the same initial mass function as that used to generate the spectral energy distribution of semi-analytical galaxies, ensuring a coherent treatment for modelling the full galaxy spectrum. By incorporating these emission line ratios with galaxy properties, we reproduce observed luminosity functions for H α, H β, [O ii], and [O iii] in the local Universe and at high redshifts. We also find good agreement between model predictions and observations for autocorrelation and cross-correlation functions of [O ii]-selected galaxies, as well as their luminosity dependence. The bias of emission line galaxies depends on both luminosity and redshift. At lower redshifts, it remains constant with increasing luminosity up to around $\sim 10^{42.5}\, {\rm erg\, s^{-1}}$ and then rises steeply for higher luminosities. The transition luminosity increases with redshift and becomes insignificant above z = 1.5. Generally, galaxy bias shows an increasing trend with redshift. However, for luminous galaxies, the bias is higher at low redshifts, as the strong luminosity dependence observed at low redshifts diminishes at higher redshifts. We provide a fitting formula for the bias of emission line galaxies as a function of luminosity and redshift, which can be utilized for large-scale structure studies with future galaxy surveys.

The minimum mass for star formation by dynamical fragmentation: dependence on epoch, dust abundance, and environment

ABSTRACT We estimate the minimum mass of a star formed by dynamical collapse and fragmentation, as a function of epoch, dust abundance, and environment. Epoch is parametrized by redshift, zred, through the variation in the temperature of the cosmic microwave background. The dust abundance is parametrized by the mass-fraction in dust, ZD, with the additional simplifying assumption that the intrinsic properties of dust do not change with ZD, only the amount of dust. Environment is parametrized by the energy-density of the ambient suprathermal radiation fields through a dilution factor ω⋆ (applied to a blackbody radiation field at $T_{\star }=10^4\, {\rm K}$). The critical condition is that a spherical proto-fragment should be able to cool, and therefore contract, fast enough to detach from neighbouring proto-fragments. The minimum mass increases with increasing redshift, increasing dust abundance, and increasing suprathermal background. Values in the range from $M_{{\rm MIN}}\sim 0.002\, {\rm M}_{\odot }$ to $M_{{\rm MIN}}\sim 0.2\, {\rm M}_{\odot }$ are obtained at the extremes of the parameter ranges we have considered (0 ≤ zred ≤ 8, 0.00016 < ZD < 0.04, and 10−15 ≤ ω⋆ ≤ 10−8). Our results agree quite well with the predictions of detailed numerical simulations invoking similar redshifts and dust abundances, but our estimates are somewhat lower; we attribute this difference to resolution issues and the small-number statistics from the simulations. The increased minimum masses predicted at high redshift and/or high suprathermal background result in significantly bottom-light initial mass functions, and therefore low mass-to-light ratios, provided that the dust abundance is not too low. The changes due to high suprathermal background may be particularly important for star formation in galactic nuclei and at high redshift.

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.

Detailed chemical composition of the globular cluster Sextans A GC-1 on the outskirts of the Local Group

The chemical composition of globular clusters (GCs) across the Local Group provides information on chemical abundance trends. Studying GCs in isolated systems in particular provides us with important initial conditions plausibly unperturbed by mergers and tidal forces from the large Local Group spirals. We present a detailed chemical abundance analysis of Sextans A GC-1. The host galaxy, Sextans A, is a low-surface-brightness dwarf irregular galaxy located on the edge of the Local Group. We derive the dynamical mass of the GC together with the mass-to-light ratio and the abundances of the alpha , Fe-peak, and heavy elements. Abundance ratios were determined from the analysis of an optical integrated-light spectrum of Sextans A GC-1, obtained with UVES on the VLT. We apply non-local thermodynamic equilibrium (NLTE) corrections to Mg, Ca, Ti, Fe, and Ni. The GC appears to be younger and more metal-poor than the majority of the GCs of the Milky Way, with an age of 8.6pm 2.7 Gyr and $ Fe/H dex. The calculated dynamical mass is dyn odot $, which results in an atypically high value of the mass-to-light ratio, 4.35pm 1.40 M$_ odot $/L$_ V odot $. Sextans A GC-1 has varying alpha elements -- the Mg abundance is extremely low, Ca and Ti are solar-scaled or mildly enhanced, and Si is enhanced. The measured values are $ Mg/Fe Ca/Fe Ti/Fe and $ Si/Fe which makes the mean alpha abundance (excluding Mg) to be enhanced $ <Si,Ca,Ti>/Fe NLTE The Fe-peak elements are consistent with scaled-solar or slightly enhanced abundances: $ Cr/Fe Mn/Fe Sc/Fe and $ Ni/Fe The heavy elements measured are Ba, Cu, Zn, and Eu. Ba and Cu have sub-solar abundance ratios ($ Ba/Fe and $ Cu/Fe <-0.343$), while Zn and Eu are consistent with their upper limits being solar-scaled and enhanced, $ Zn/Fe <+0.171$ and $ Eu/Fe The composition of Sextans A GC-1 resembles the overall pattern and behaviour of GCs in the Local Group. The anomalous values are the mass-to-light ratio and the depleted abundance of Mg. There is no definite explanation for such an extreme abundance value. Variations in the initial mass function or the presence of an intermediate-mass black hole might explain the high mass-to-light ratio value.

Exploring the nature of UV-bright z ≳ 10 galaxies detected by JWST: star formation, black hole accretion, or a non-universal IMF?

ABSTRACT We use the Cosmic Archaeology Tool (CAT) semi-analytical model to explore the contribution of Population (Pop) III/II stars and active galactic nuclei (AGNs) to the galaxy ultraviolet (UV) luminosity function (LF) evolution at 4 ≤ z ≤ 20. We compare in particular with recent JWST data in order to explore the apparent tension between observations and theoretical models in the number density of bright galaxies at z ≳ 10. The model predicts a star formation history dominated by UV faint (MUV > −18) galaxies, with a Pop III contribution of $\lesssim 10~{{\ \rm per\ cent}}$ ($\lesssim 0.5~{{\ \rm per\ cent}}$) at z ≃ 20 (z ≃ 10). Stars are the primary sources of cosmic reionization, with $5~{{\ \rm per\ cent}}\!-\!10~{{\ \rm per\ cent}}$ of ionizing photons escaping into the intergalatic medium at 5 ≤ z ≤ 10, while the contribution of unobscured AGNs becomes dominant only at z ≲ 5. The predicted stellar and AGN UV LFs reproduce the observational data at 5 ≲ z ≲ 9–10. At higher redshift, CAT predicts a steeper evolution in the faint-end slope (MUV > –18), and a number density of bright galaxies (MUV ≃ −20) consistent with data at z ∼ 10–11, but smaller by 0.8 dex at z ∼ 12–13, and 1.2 dex at z ∼ 14–16, when compared to the values estimated by recent studies. Including the AGN emission does not affect the above findings, as AGNs contribute at most to $\lesssim 10~{{\ \rm per\ cent}}$ of the total UV luminosity at MUV < –19 and z ≳ 10. Interestingly, considering a gradual transition in the stellar initial mass function, modulated by metallicity and redshift as suggested by recent simulations, the model agrees with JWST data at z ∼ 12–13, and the disagreement at z ∼ 14–16 is reduced to 0.5 dex.

How Population III Supernovae Determined the Properties of the First Galaxies

Massive Population III stars can die as energetic supernovae that enrich the early Universe with metals and determine the properties of the first galaxies. With masses of about 109 M ⊙ at z ≳ 10, these galaxies are believed to be the ancestors of the Milky Way. This paper investigates the impact of Population III supernova remnants (SNRs) from both Salpeter-like and top-heavy initial mass functions (IMFs) on the formation of first galaxies with high-resolution radiation-hydrodynamical simulations with the ENZO code. Our findings indicate that SNRs from a top-heavy Population III IMF produce more metals, leading to more efficient gas cooling and earlier Population II star formation in the first galaxies. From a few hundred to a few thousand Population II stars can form in the central regions of these galaxies. These stars have metallicities of 10−3–10−2, Z ⊙, greater than those of extremely metal-poor (EMP) stars. Their mass function follows a power-law distribution with dN(M*)/dM*∝M*α , where M * is stellar mass, and α = 2.66–5.83 and is steeper for a top-heavy IMF. We thus find that EMP stars were not typical of most primitive galaxies.

Quantifying the Effects of Known Unknowns on Inferred High-redshift Galaxy Properties: Burstiness, IMF, and Nebular Physics

The era of the James Webb Space Telescope ushers stellar population models into uncharted territories, particularly at the high-redshift frontier. In a companion paper, we apply the Prospector Bayesian framework to jointly infer galaxy redshifts and stellar population properties from broadband photometry as part of the UNCOVER survey. Here we present a comprehensive error budget in spectral energy distribution (SED) modeling. Using a sample selected to have photometric redshifts higher than 9, we quantify the systematic shifts stemming from various model choices in inferred stellar mass, star formation rate (SFR), and age. These choices encompass different timescales for changes in the star formation history (SFH), nonuniversal stellar initial mass functions (IMF), and the inclusion of variable nebular abundances, gas density, and ionizing photon budget. We find that the IMF exerts the strongest influence on the inferred properties: the systematic uncertainties can be as much as 1 dex, 2–5 times larger than the formal reported uncertainties in mass and SFR, and importantly, exceed the scatter seen when using different SED fitting codes. Although the assumptions on the lower end of the IMF induce degeneracy, our findings suggest that a common practice in the literature of assessing uncertainties in SED-fitting processes by comparing multiple codes is substantively underestimating the true systematic uncertainty. Highly stochastic SFHs change the inferred SFH by much larger than the formal uncertainties, and introduce ∼0.8 dex systematics in SFR averaged over a short timescale and ∼0.3 dex systematics in average age. Finally, employing a flexible nebular emission model causes ∼0.2 dex systematic increase in mass and SFR, comparable to the formal uncertainty. This paper constitutes an initial step toward a complete uncertainty estimate in SED modeling.

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Impact of Galaxies on the Circumgalactic Medium Metal Enrichment at z > 6 Using the JWST and VLT

We characterize the multiphase circumgalactic medium (CGM) and galaxy properties at z = 6.0–6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305–3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (O i, Mg ii, Fe ii, and C ii). They are combined with the published absorbing systems in Davies et al. at the same redshift range to form a sample of nine metal absorbers at z = 6.03–6.49. We identify eight galaxies within 1000 km s−1 and 350 kpc around the absorbing gas from the ASPIRE spectroscopic data, with their redshifts secured by [O iii] (λ λ4959, 5007) doublets and Hβ emission lines. Our spectral energy distribution fitting indicates that the absorbing galaxies have stellar masses ranging from 107.2 to 108.8 M ⊙ and metallicity between 0.02 and 0.4 solar. Notably, the z = 6.2713 system in the J0305–3150 field resides in a galaxy overdensity region, which contains two (tentatively) merging galaxies within 350 kpc and seven galaxies within 1 Mpc. We measure the relative abundances of α elements to iron ([α/Fe]) and find that the CGM gas in the most overdense region exhibits a lower [α/Fe] ratio. Our modeling of the galaxy’s chemical abundance favors a top-heavy stellar initial mass function and hints that we may be witnessing the contribution of the first generation of Population III stars to the CGM at the end of the reionization epoch.

The stellar ‘Snake’ – II. The mass function

ABSTRACT We present a comprehensive investigation on the mass function (MF) of a snake-like stellar structure in the solar neighbourhood, building on our previous discovery. To ensure the reliability of the data, we reselect the member stars of the Stellar ‘Snake’ in the latest Gaia Data Release 3 using the same approach as the initial series of articles. We also precisely measure the physical parameters of the clusters within the Stellar Snake. In light of the high completeness of the member stars in the cluster regions, we develop a simulated model colour–magnitude diagram-based inference method to derive the mass function, binary fraction, and mass-ratio distribution of the clusters in the Stellar Snake. Notably, despite their similar ages and metallicity, we discover systematic variations in the MFs along the elongation direction of the Snake in the mass range of 0.5 to 2.0 M⊙. The ‘head’ of the Snake conforms to a canonical initial mass function with a power-law slope of α ∼ −2.3. Extending towards the ‘tail’, the MF becomes more top-light, indicating a deficiency of massive stars within these clusters. This result provides evidence for the delayed formation of massive stars in the clusters. Such clues give support to the hypothesis that the Stellar Snake constitutes as a hierarchically primordial structure.

Stellar populations of massive early-type galaxies observed by MUSE

ABSTRACT Stellar population studies of massive early-type galaxies (ETGs) suggest that the stellar initial mass function may not be universal. In particular, the centres of ETGs seem to contain an excess of low-mass dwarf stars compared to our own Galaxy. Through high-resolution data from MUSE (Multi Unit Spectroscopic Explorer), we carry out a detailed study of the stellar populations of eight massive ETGs. We use full spectrum fitting to determine ages, element abundances, and initial mass function (IMF) slopes for spatially binned spectra. We measure flat gradients in age and [Mg/Fe] ratio, as well as negative gradients in metallicity and [Na/Fe]. We detect IMF gradients in some galaxies, with the centres hosting bottom-heavy IMFs and mass excess factors between 1.5 and 2.5 compared to a Kroupa IMF. The IMF slope below 0.5 M$_\odot$varies for our galaxy sample between 1 and 2.8, with negative radial gradients, while the IMF slope between 0.5 and 1M$_\odot$ has a steep value of $\sim$3 with mildly positive gradients for most galaxies. For M87, we find excellent agreement with the dynamical mass-to-light ratio (M/L) as a function of radius. For the other galaxies, we find systematically higher M/L from stellar populations compared to orbit-based dynamical analysis of the same data. This discrepancy increases with NaI strength, suggesting a combination of calibration issues of this line and correlated uncertainties.

Effects of physical conditions on the stellar initial mass function: The low-metallicity star-forming region Sh 2-209

Recent work suggested that the variation of the initial mass function (IMF) of stars depends on the physical conditions, notably, the metallicity and gas density. We investigated the properties of two clusters, namely the main cluster (MC) and the subcluster (SC), in the low-metallicity HII region Sh 2-209 (S209) based on recently derived IMFs. We tested three previously published correlations using previous observations: the top-heaviness of the IMF in dependence on metallicity, the half-mass radius, and the most massive star in dependence on the stellar mass of the embedded clusters. For this region, two different galactocentric distances, namely $10.5\ kpc $ and $18\ kpc $, were considered, where an age-distance-degeneracy was found for the previously determined IMF to be consistent with other formulated metallicity and density dependent IMFs. The determined half-mass radius $r_ h pc $ and the embedded cluster density $ ecl 0.1)\,10^6 M_ pc $ for the MC with an age of 0.5 Myr in S209 assuming a galactocentric distance of $18\ kpc $ support the assumption that a low-metallicity environment results in a denser cluster, which leads to a top-heavy IMF. Thus, all three tests are consistent with the previously published correlations. The results for S209 are placed in the context with the IMF determination within the metal-poor cluster in the star-forming region NGC 346 in the Small Magellanic Cloud.

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

ABSTRACT We 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.

Predicting the heaviest black holes below the pair instability gap

ABSTRACT Traditionally, the pair instability (PI) mass gap is located between 50 and 130 M⊙, with stellar mass black holes (BHs) expected to ‘pile up’ towards the lower PI edge. However, this lower PI boundary is based on the assumption that the star has already lost its hydrogen (H) envelope. With the announcement of an ‘impossibly’ heavy BH of 85 M⊙ as part of GW 190521 located inside the traditional PI gap, we realized that blue supergiant (BSG) progenitors with small cores but large hydrogen envelopes at low metallicity (Z) could directly collapse to heavier BHs than had hitherto been assumed. The question of whether a single star can produce such a heavy BH is important, independent of gravitational wave events. Here, we systematically investigate the masses of stars inside the traditional PI gap by way of a grid of 336 detailed mesa stellar evolution models calculated across a wide parameter space, varying stellar mass, overshooting, rotation, semiconvection, and Z. We evolve low Z stars in the range 10−3 < Z/Z⊙ < ZSMC, making no prior assumption regarding the mass of an envelope, but instead employing a wind mass-loss recipe to calculate it. We compute critical carbon–oxygen and helium core masses to determine our lower limit to PI physics, and we provide two equations for Mcore and Mfinal that can also be of use for binary population synthesis. Assuming the H envelope falls into the BH, we confirm the maximum BH mass below PI is MBH ≃ 93.3 M⊙. Our grid allows us to populate the traditional PI gap, and we conclude that the distribution of BHs above the traditional boundary is not solely due to the shape of the initial mass function, but also to the same stellar interior physics (i.e. mixing) that which sets the BH maximum.

Cosmic evolution early release science survey (CEERS): multiclassing galactic dwarf stars in the deep JWST/NIRCam

ABSTRACT Low-mass (sub)stellar objects represent the low end of the initial mass function, the transition to free-floating planets and a prominent interloper population in the search for high-redshift galaxies. To what accuracy can photometry only classify these? JWST/NIRCam has several advantages over Hubble Space Telescope (HST)/WFC3 near-infrared (NIR): more filters, a greater wavelength range, and greater spatial resolution. Here, we present a catalogue of (sub)stellar dwarfs identified in the cosmic evolution early release science survey (CEERS). We identify 518 stellar objects down to mF200W 28 using half-light radius, a full three magnitudes deeper than typical HST/WFC3 images. A k-means nearest neighbour (kNN) algorithm identifies and types these sources, using four HST/WFC3 and four NIRCam filters, trained on IRTF/spex spectra of nearby brown dwarfs. The kNN with four neighbours classifies well within two subtypes: e.g. M2±2, achieving ∼95 per cent precision and recall. In CEERS, we find 9 M8±2, 2 L6±2, 1 T4±2, and 15 T8±2. We compare the observed long wavelength NIRCam colours – not used in the kNN – to those expected for brown dwarf atmospheric models. The NIRCam F356W−F444W and F410M−F444W colours are redder by a magnitude for the type assigned by the kNN, hinting at a wider variety of atmospheres for these objects. We find a 300–350 pc scale height for M6±2 dwarfs plus a second structural component and a 150–200 pc scale height for T6±2 type dwarfs, consistent with literature values. A single M8±2 dwarf is spectroscopically confirmed at 4 kpc distance.

True Pair-instability Supernova Descendant: Implications for the First Stars’ Mass Distribution

The initial mass function (IMF) of the first Population III (Pop III) stars remains a persistent mystery. Their predicted massive nature implies the existence of stars exploding as pair-instability supernovae (PISNe), but no observational evidence had been found. Now, the LAMOST survey claims to have discovered a pure PISN descendant, J1010+2358, at [Fe/H] = − 2.4. Here we confirm that a massive 250–260 M ⊙ PISN is needed to reproduce the abundance pattern of J1010+2358. However, the PISN contribution can be as low as 10%, since key elements are missing to discriminate between scenarios. We investigate the implications of this discovery for the Pop III IMF, by statistical comparison with the predictions of our cosmological galaxy formation model, NEFERTITI. First, we show that the nondetection of mono-enriched PISN descendants at [Fe/H] < − 2.5 allows us to exclude (i) a flat IMF at a 90% confidence level; and (ii) a Larson-type IMF with characteristic mass m ch/M ⊙ > 191.16x − 132.44, where x is the slope, at a 75% confidence level. Second, we show that if J1010+2358 has only inherited <70% of its metals from a massive PISN, no further constraints can be put on the Pop III IMF. If, instead, J1010+2358 will be confirmed to be a nearly pure (>90%) PISN descendant, it will offer strong and complementary constraints on the Pop III IMF, excluding the steepest and bottom-heaviest IMFs: m ch/M ⊙ < 143.21x − 225.94. Our work shows that even a single detection of a pure PISN descendant can be crucial to our understanding of the mass distribution of the first stars.

Kindling the First Stars. I. Dependence of Detectability of the First Stars with JWST on the Population III Stellar Masses

The first Population III (Pop III) stars formed out of primordial, metal-free gas, in minihalos at z > 20, and kickstarted the cosmic processes of reionization and enrichment. While these stars are likely more massive than their enriched counterparts, the current unknowns of their astrophysics include when the first Pop III stars ignited, how massive they were, and when and how the era of the first stars ended. Investigating these questions requires an exploration of a multidimensional parameter space, including the slope of the Pop III stellar initial mass function (IMF) and the strength of the nonionizing UV background. In this work, we present a novel model which treats both the slope and maximum mass of Pop III stars as truly free parameters while including the physics of the fragmentation of primordial gas. Our results also hint at a nonuniversal Pop III IMF which is dependent on the efficiency of primordial gas fragmentation. Our relatively simple model reproduces the results from hydrodynamic simulations, but with a computational efficiency which allows us to investigate the observable differences between a wide range of potential Pop III IMFs. In addition, the evolution of the number density of Pop III stars may provide insight into the evolution of the H2 dissociating background. While the slope of the Pop III IMF does not significantly affect the predicted number density of the first stars, more top-heavy IMFs produce Pop III star clusters which are 2–3 magnitudes brighter than their more bottom-heavy counterparts. While the Pop III star clusters are too dim for direct detection by JWST, we find they are within the reach of gravitational lensing.