Top-heavy Initial Mass Function Research Articles

Osaka Feedback Model. III. Cosmological Simulation CROCODILE

We introduce our new cosmological simulation data set CROCODILE, executed using the GADGET4-Osaka smoothed particle hydrodynamics code. This simulation incorporates an updated supernova (SN) feedback model of Y. Oku et al. and an active galactic nuclei (AGN) feedback model. A key innovation in our SN feedback model is the integration of a metallicity- and redshift-dependent, top-heavy initial mass function. Our SN model introduces a new consideration that results in an order of magnitude difference in the energy injection rate per unit stellar mass formed at high redshift. The CROCODILE data set is comprehensive, encompassing a variety of runs with diverse feedback parameters. This allows for an in-depth exploration of the relative impacts of different feedback processes in galactic evolution. Our initial comparisons with observational data, spanning the galaxy stellar mass function, the star formation main sequence, and the mass–metallicity relation, show promising agreement, especially for the Fiducial run. These results establish a solid foundation for our future work. We find that SN feedback is a key driver in the chemical enrichment of the intergalactic medium (IGM). Additionally, the AGN feedback creates metal-rich, bipolar outflows that extend and enrich the circumgalactic medium and IGM over a few Mpc scales.

True unicorns and false positives: Simulated probabilities of dark massive companions to bright stars

Abstract Many compact objects (black holes and neutron stars) exist in binaries. These binaries are normally discovered through their interactions, either from accretion as an X-ray binary or collisions as a gravitational wave source. However, the majority of compact objects in binaries should be non-interacting. Recently proposed discoveries have used radial velocities of a bright star (main sequence or evolved) that are indicative of a massive but dark companion, which is inferred to be a compact object. Unfortunately, this burgeoning new field has been hindered by false positives, including the “Unicorn” (V723 Mon) which was initially believed to be a red giant/black hole binary before being refuted. In this work, we investigate the evolution of stellar binary populations over time, using the binary evolution code COSMIC to simulate binary populations and determine the probability of a candidate object being either a “true Unicorn” (actual compact objects in binaries) or a false positive. We find that main sequence stars have a higher true Unicorn probability than red giants or naked helium stars (an exposed core of an evolved star), particularly if the companion is more massive and is ≥3 times less luminous than the MS star. We also find that a top-heavy initial mass function raises the true Unicorn probability further, that super-solar metallicity reduces the probability, and that most true Unicorns are found at periods ≤100 days. Finally, we find that a significant fraction of true Unicorns do not evolve into x-ray binaries during the age of the universe.

Photometric classification of stars around the Milky Way’s central black hole

Context. The presence of young massive stars in the Galactic Centre (GC) raises questions about how such stars could form near the massive black hole Sagittarius A* (Sgr A*). Furthermore, the shape of the initial mass function (IMF) in this region seems to differ from its standard Salpeter/Kroupa law. Due to observational challenges such as extreme extinction and crowding, our understanding of the stellar population in this region remains limited, with spectroscopic data available only for selected small and comparably bright sources. Aims. We aim to improve our knowledge about the distribution and the IMF of young, massive, stars in the vicinity of Sgr A*. Methods. We used intermediate band (IB) photometry to identify candidates for massive young stars. To ensure robust classification, we applied three different, but complementary methods: Bayesian inference, a basic neural network, and a fast gradient-boosted trees algorithm. Results. We obtain spectral energy distributions for 6590 stars, 1181 of which have been previously classified spectroscopically. We identify 351 stars that are classified as early types by all three classification methods, with 155 of them being newly identified candidates. The radial density profiles for late and early-type stars fit well with broken power laws, revealing a break radius of 9.2 ± 0.6″ for early-type stars. The late-type stars show a core-like distribution around Sgr A* while the density of the early-type stars increases steeply towards the black hole, consistent with previous work. We infer a top-heavy IMF of the young stars near Sgr A* (R < 9″), with a power-law of 1.6 ± 0.1. At greater distances from Sgr A* a standard Salpeter/Kroupa IMF can explain the data. Additionally, we demonstrate that IB photometry can also constrain the metallicities of late-type stars, estimating metallicities for over 600 late-type stars. Conclusions. The variation of the IMF with radial distance from Sgr A* suggests that different mechanisms of star formation may have been at work in this region. The top-heavy IMF in the innermost region is consistent with star formation in a disc around Sgr A*.

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.

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.

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 Interplay between the Initial Mass Function and Star Formation Efficiency through Radiative Feedback at High Stellar Surface Densities

The observed rest-UV luminosity function at cosmic dawn (z ∼ 8–14) measured by JWST revealed an excess of UV-luminous galaxies relative to many prelaunch theoretical predictions. A high star formation efficiency (SFE) and a top-heavy initial mass function (IMF) are among the mechanisms proposed for explaining this excess. Although a top-heavy IMF has been proposed for its ability to increase the light-to-mass ratio (ΨUV), the resulting enhanced radiative pressure from young stars could decrease the SFE, potentially driving galaxy luminosities back down. In this Letter, we use idealized radiation hydrodynamic simulations of star cluster formation to explore the effects of a top-heavy IMF on the SFE of clouds typical of the high-pressure conditions found at these redshifts. We find that the SFE in star clusters with solar-neighborhood-like dust abundance decreases with increasingly top-heavy IMFs—by ∼20% for an increase of a factor of 4 in ΨUV and by 50% for a factor of ∼10 in ΨUV. However, we find that an expected decrease in the dust-to-gas ratio (∼0.01 × solar) at these redshifts can completely compensate for the enhanced light output. This leads to a (cloud-scale; ∼10 pc) SFE that is ≳70% even for a factor of 10 increase in ΨUV, implying that highly efficient star formation is unavoidable for high surface density and low-metallicity conditions. Our results suggest that a top-heavy IMF, if present, likely coexists with efficient star formation in these galaxies.

Contribution of very massive stars to the sulfur abundance in star-forming galaxies: Role of pair-instability supernovae

Context. Recent work presented increasing evidence of high non-constant S/O abundance ratios observed in star-forming metal-poor galaxies that deviated from the constant canonical S/O across a wide range of O/H abundances. Similar peculiarly high Fe/O ratios have also recently been detected. Aims. We investigate whether these high S/O ratios at low metallicities could be explained when the process of pair-instability supernovae (PISN) in chemical modelling is included, through which a similar behaviour of the Fe/O ratios was reproduced successfully. Methods. We used chemical evolution models that considered the stages of PISN in the previously published yields and adopted a suitable initial mass function (IMF) to characterize this evolutionary stage appropriately. Results. The peculiarly high values and the behaviour of the observed S/O versus O/H relation can be reproduced when the ejecta of very massive stars that go through the process of PISN are taken into account. Additionally, a bimodal top-heavy IMF and an initial strong burst of star formation are required to reach the reported high S/O values. Conclusions. We show that the role of very massive stars going through the process of PISN should be taken into account to explain the chemical enrichment of sulphur and oxygen in metal-poor star-forming regions.

Heavy black hole seed formation in high-z atomic cooling halos

Context. Halos with masses in excess of the atomic limit are believed to be ideal environments in which to form heavy black hole seeds with masses above 103 M⊙. In cases where the H2 fraction is suppressed, this is expected to lead to reduced fragmentation of the gas and the generation of a top-heavy initial mass function. In extreme cases this can result in the formation of massive black hole seeds. Resolving the initial fragmentation scale and the resulting protostellar masses has, until now, not been robustly tested. Aims. We run zoom-in simulations of atomically cooled halos in which the formation of H2 is suppressed to assess whether they can truly resist fragmentation at high densities and tilt the initial mass function towards a more top-heavy form and the formation of massive black hole seeds. Methods. Cosmological simulations were performed with the moving mesh code AREPO, using a primordial chemistry network until z ∼ 11. Three haloes with masses in excess of the atomic cooling mass were then selected for detailed examination via zoom-ins. A series of zoom-in simulations, with varying levels of maximum spatial resolution, captured the resulting fragmentation and formation of metal-free stars using the sink particle technique. The highest resolution simulations resolved densities up to 10−6 g cm−3 (1018 cm−3) and captured a further 100 yr of fragmentation behaviour at the centre of the halo. Lower resolution simulations were then used to model the future accretion behaviour of the sinks over longer timescales. Results. Our simulations show intense fragmentation in the central region of the halos, leading to a large number of near-solar mass protostars. Even in the presence of a super-critical Lyman-Werner radiation field (JLW > 105J21), H2 continues to form within the inner ∼2000 au of the halo. Despite the increased fragmentation, the halos produce a protostellar mass spectrum that peaks at higher masses relative to standard Population III star-forming halos. The most massive protostars have accretion rates of 10−3–10−1 M⊙ yr−1 after the first 100 years of evolution, while the total mass of the central region grows at 1 M⊙ yr−1. Lower resolution zoom-ins show that the total mass of the system continues to accrete at ∼1 M⊙ yr−1 for at least 104 yr, although how this mass is distributed amongst the rapidly growing number of protostars is unclear. However, assuming that a fraction of stars can continue to accrete rapidly, the formation of a sub-population of stars with masses in excess of 103 M⊙ is likely in these halos. In the most optimistic case, we predict the formation of heavy black hole seeds with masses in excess of 104 M⊙, assuming an accretion behaviour in line with expectations from super-competitive accretion and/or frequent mergers with secondary protostars.

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.

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.

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.

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.

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.

Pure Spectroscopic Constraints on UV Luminosity Functions and Cosmic Star Formation History from 25 Galaxies at z spec = 8.61–13.20 Confirmed with JWST/NIRSpec

We present pure spectroscopic constraints on the UV luminosity functions and cosmic star formation rate (SFR) densities from 25 galaxies at z spec = 8.61–13.20. By reducing the JWST/NIRSpec spectra taken in multiple programs of Early Release Observation, Early Release Science, General Observer, and Director’s Discretionary Time observations with our analysis technique, we independently confirm 16 galaxies at z spec = 8.61–11.40, including new redshift determinations, and a bright interloper at z spec = 4.91 that was claimed as a photometric candidate at z ∼ 16. In conjunction with nine galaxies at redshifts up to z spec = 13.20 in the literature, we make a sample of 25 spectroscopically confirmed galaxies in total and carefully derive the best estimates and lower limits of the UV luminosity functions. These UV luminosity function constraints are consistent with the previous photometric estimates within the uncertainties and indicate mild redshift evolution toward z ∼ 12, showing tensions with some theoretical models of rapid evolution. With these spectroscopic constraints, we obtain firm lower limits of the cosmic SFR densities and spectroscopically confirm a high SFR density at z ∼ 12 beyond the constant star formation efficiency models, which supports earlier claims from the photometric studies. While there are no spectroscopically confirmed galaxies with very large stellar masses violating the ΛCDM model due to the removal of the bright interloper, we confirm star-forming galaxies at z spec = 11–13 with stellar masses much higher than model predictions. Our results indicate possibilities of high star formation efficiency (>5%), a hidden active galactic nucleus, a top-heavy initial mass function (possibly with Population III), and large scatter/variance. Having these successful and unsuccessful spectroscopy results, we suggest observational strategies for efficiently removing low-redshift interlopers for future JWST programs.

Modeling the JWST High-redshift Galaxies with a General Formation Scenario and the Consistency with the ΛCDM Model

Early results from the James Webb Space Telescope (JWST) observations have hinted at two traces beyond the standard cosmological framework. One is the extraordinarily high stellar masses and their density at z = 7.5 ∼ 9.1; another is the unexpected abundance of ultraviolet (UV) bright galaxies at z ≥ 10. Nevertheless, both pieces of evidence are not statistically robust yet. In this work, we construct rest-frame UV luminosity functions (LFs) based on a general formation model for these high-redshift galaxy candidates, since UV LFs always carry the information of stellar formation efficiency (SFE), initial mass function (IMF), dust attenuation, and other crucial elements for galaxy evolution. By updating the massive galaxies candidates with spectroscopic observations and exploring the parameter space of SFE, we are able to reasonably explain the cumulative stellar mass density within the redshift range of 7.5–9.1, with only one galaxy exhibiting unusual characteristics. We also reveal a potential nonmonotonic trend of SFE with the increasing redshift. At higher redshift (z ∼ 13), bright UV LFs can be well fitted with non–dust attenuation or top-heavy IMF for Population III stars. The Population III star scenario can also naturally account for the possible dip of the peak SFE evolution curve at z ∼ 9.

Star Formation in Self-gravitating Disks in Active Galactic Nuclei. III. Efficient Production of Iron and Infrared Spectral Energy Distributions

Strong iron lines are a common feature of the optical spectra of active galactic nuclei (AGNs) and quasars from z ∼ 6−7 to the local universe, and [Fe/Mg] ratios do not show cosmic evolution. During active episodes, accretion disks surrounding supermassive black holes (SMBHs) inevitably form stars in the self-gravitating part, and these stars accrete with high accretion rates. In this paper, we investigate the population evolution of accretion-modified stars (AMSs) to produce iron and magnesium in AGNs. The AMSs, as a new type of star, are allowed to have any metallicity but without significant loss from stellar winds, since the winds are choked by the dense medium of the disks and return to the core stars. Mass functions of the AMS population show a pile-up or cutoff pile-up shape in top-heavy or top-dominant forms if the stellar winds are strong, consistent with the narrow range of supernovae (SNe) explosions driven by the known pair-instability. This provides an efficient way to produce metals. Meanwhile, SN explosions support an inflated disk as a dusty torus. Furthermore, the evolving top-heavy initial mass functions lead to bright luminosity in infrared bands in dusty regions. This contributes a new component in infrared bands, which is independent of the emissions from the central part of accretion disks, appearing as a long-term trending of the NIR continuum compared to optical variations. Moreover, the model can be further tested through reverberation mapping of emission lines, including LIGO/LISA detections of gravitational waves and signatures from spatially resolved observations of GRAVITY+/VLTI.

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.

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.

Population III X-ray binaries and their impact on the early universe

ABSTRACT The first population of X-ray binaries (XRBs) is expected to affect the thermal and ionization states of the gas in the early Universe. Although these X-ray sources are predicted to have important implications for high-redshift observable signals, such as the hydrogen 21-cm signal from cosmic dawn and the cosmic X-ray background, their properties are poorly explored, leaving theoretical models largely uninformed. In this paper we model a population of X-ray binaries arising from zero metallicity stars. We explore how their properties depend on the adopted initial mass function (IMF) of primordial stars, finding a strong effect on their number and X-ray production efficiency. We also present scaling relations between XRBs and their X-ray emission with the local star formation rate, which can be used in sub-grid models in numerical simulations to improve the X-ray feedback prescriptions. Specifically, we find that the uniformity and strength of the X-ray feedback in the intergalactic medium is strongly dependant on the IMF. Bottom-heavy IMFs result in a smoother distribution of XRBs, but have a luminosity orders of magnitude lower than more top-heavy IMFs. Top-heavy IMFs lead to more spatially uneven, albeit strong, X-ray emission. An intermediate IMF has a strong X-ray feedback while sustaining an even emission across the intergalactic medium. These differences in X-ray feedback could be probed in the future with measurements of the cosmic dawn 21-cm line of neutral hydrogen, which offers us a new way of constraining population III IMF.