Top-heavy Initial Mass Function Research Articles

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.

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.

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 (≳ 9 M⊙) and very massive (≳ 100 M⊙) 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 Mcluster ∼ 2 × 105 M⊙ in a low-metallicity dwarf galaxy starburst. The simulation realises new stars from a stellar initial mass function (IMF) between 0.08 M⊙ and ∼400 M⊙ and follows stellar winds, radiation and supernova-feedback of single massive stars with evolution tracks. Star clusters form on timescales 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) stars can be centrally concentrated in the most massive clusters (≳ 104 M⊙) and the locked wind-material increases approximately as $M_\mathrm{cluster}^{2}$. These trends resemble the characteristics of observed second population stars in globular clusters. We fit scaling relations to the log-normal distributed wind-mass fractions and extrapolate to possible globular cluster progenitors of Mcluster = 107 M⊙ to investigate whether a dominant second population 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 e.g. 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.

CEERS Key Paper. I. An Early Look into the First 500 Myr of Galaxy Formation with JWST

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z ∼ 0.5 to >10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 arcmin2, to search for candidate galaxies at z > 9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multiband photometry across seven NIRCam broad- and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z ∼ 9–16. These objects are compact with a median half-light radius of ∼0.5 kpc. We present an early estimate of the z ∼ 11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M UV ∼ −20 appears to evolve very little from z ∼ 9 to 11. We also find that the abundance (surface density [arcmin−2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z > 10, star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultrahigh-redshift galaxies with which to further explore these conclusions.

A Comprehensive Study of Galaxies at z ∼ 9–16 Found in the Early JWST Data: Ultraviolet Luminosity Functions and Cosmic Star Formation History at the Pre-reionization Epoch

We conduct a comprehensive study on dropout galaxy candidates at z ∼ 9–16 using the first 90 arcmin2 James Webb Space Telescope (JWST) Near Infrared Camera images taken by the early release observations (ERO) and early release science programs. With the JWST simulation images, we find that a number of foreground interlopers are selected with a weak photo-z determination (Δχ 2 > 4). We thus carefully apply a secure photo-z selection criterion (Δχ 2 > 9) and conventional color criteria with confirmations of the ERO Near Infrared Spectrograph spectroscopic redshifts, and obtain a total of 23 dropout galaxies at z ∼ 9–16, including two candidates at and . We perform thorough comparisons of dropout galaxies found in our work with recent JWST studies, and conclude that our galaxy sample is reliable enough for statistical analyses. We derive the UV luminosity functions at z ∼ 9–16, and confirm that our UV luminosity functions at z ∼ 9 and 12 agree with those determined by other Hubble Space Telescope and JWST studies. The cosmic star formation rate (SFR) density decreases from z ∼ 9 to 12, and perhaps to 16, but the densities at z ∼ 12–16 are higher than the constant star formation efficiency model. Interestingly, there are six bright galaxy candidates at z ∼ 10–16 with M UV < −19.5 mag and M * ∼ 108−9 M ⊙. Because a majority (∼80%) of these galaxies show no signatures of active galactic nuclei in their morphologies, the high cosmic SFR densities and the existence of these UV-luminous galaxies are explained by the lack of suppression of star formation by the UV background radiation at the pre-reionization epoch and/or an efficient UV radiation production by a top-heavy initial mass function with Population III–like star formation.

A closer look at supernovae as seeds for galactic magnetization

Context. Explaining the currently observed magnetic fields in galaxies requires relatively strong seeding in the early Universe. One of the current theories proposes that magnetic seeds on the order of μG were expelled by supernova (SN) explosions after primordial fields of nG strength or weaker were amplified in stellar interiors. Aims. In this work, we take a closer look at this theory and calculate the maximum magnetic energy that can be injected in the interstellar medium by a stellar cluster of mass Mcl based on what is currently known about stellar magnetism. Methods. We consider early-type stars and adopt either a Salpeter or a top-heavy initial mass function. For their magnetic fields, we adopt either a Gaussian or a bimodal distribution. The Gaussian model assumes that all massive stars are magnetized with 103 &lt; ⟨B*⟩&lt; 104 G, while the bimodal, consistent with observations of Milky Way stars, assumes only 5 − 10% of OB stars have 103 &lt; ⟨B*⟩&lt; 104 G, while the rest have 10 &lt; ⟨B*⟩&lt; 102 G. We ignore the effect of magnetic diffusion and assume no losses of magnetic energy. Results. We find that the maximum magnetic energy that can be injected by a stellar population is between 10−10 and 10−7 times the total SN energy. The highest end of these estimates is about five orders of magnitude lower than what is usually employed in cosmological simulations, where about 10−2 of the SN energy is injected as magnetic. Conclusions. Pure advection of the stellar magnetic field by SN explosions is a good candidate for seeding a dynamo, but not enough to magnetize galaxies. Assuming SNe as the main mechanism for galactic magnetization, the magnetic field cannot exceed an intensity of 10−7 G in the best-case scenario for a population of 105 solar masses in a superbubble of 300 pc radius, while more typical values are between 10−10 and 10−9 G. Therefore, other scenarios for galactic magnetization at high redshift need to be explored.