Initial Mass Function Research Articles

VISIONS: the VISTA Star Formation Atlas

VISIONS is an ESO public survey of five nearby (d < 500 pc) star-forming molecular cloud complexes that are canonically associated with the constellations of Chamaeleon, Corona Australis, Lupus, Ophiuchus, and Orion. The survey was carried out with the Visible and Infrared Survey Telescope for Astronomy (VISTA), using the VISTA Infrared Camera (VIRCAM), and collected data in the near-infrared passbands J (1.25 μm), H (1.65 μm), and KS (2.15 μm). With a total on-sky exposure time of 49.4h VISIONS covers an area of 650 deg2, it is designed to build an infrared legacy archive with a structure and content similar to the Two Micron All Sky Survey (2MASS) for the screened star-forming regions. Taking place between April 2017 and March 2022, the observations yielded approximately 1.15 million images, which comprise 19 TB of raw data. The observations undertaken within the survey are grouped into three different subsurveys. First, the wide subsurvey comprises shallow, large-scale observations and it has revisited the star-forming complexes six times over the course of its execution. Second, the deep subsurvey of dedicated high-sensitivity observations has collected data on areas with the largest amounts of dust extinction. Third, the control subsurvey includes observations of areas of low-to-negligible dust extinction. Using this strategy, the VISIONS observation program offers multi-epoch position measurements, with the ability to access deeply embedded objects, and it provides a baseline for statistical comparisons and sample completeness – all at the same time. In particular, VISIONS is designed to measure the proper motions of point sources, with a precision of 1 mas yr−1 or better, when complemented with data from the VISTA Hemisphere Survey (VHS). In this way, VISIONS can provide proper motions of complete ensembles of embedded and low-mass objects, including sources inaccessible to the optical ESA Gaia mission. VISIONS will enable the community to address a variety of research topics from a more informed perspective, including the 3D distribution and motion of embedded stars and the nearby interstellar medium, the identification and characterization of young stellar objects, the formation and evolution of embedded stellar clusters and their initial mass function, as well as the characteristics of interstellar dust and the reddening law.

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.

Exploring the impact of IMF and binary parameter stochasticity with a binary population synthesis code

ABSTRACT Low-mass star formation regions are unlikely to fully populate their initial mass functions (IMFs), leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 102 to 107 M⊙, using a binary stellar population synthesis code. For clusters with stellar masses less than 105 M⊙, observable quantities show substantial scatter and their mean properties reflect the expected deficit of massive stars. In common with previous work, we find that purely stochastic sampling of the IMF appears to underestimate the mass of the most massive star in known clusters. However, even with this constraint, the majority of clusters likely inject sufficient kinetic energy to clear their birth clusters of gas. For quantities that directly measure the impact of the most massive stars, such as Nion, ξion, and βUV, uncertainties due to stochastic sampling dominate over those from the IMF shape or distribution of binary parameters, while stochastic sampling has a negligible effect on the stellar continuum luminosity density.

Dynamics of intermediate mass black holes in globular clusters

Context. We recently introduced a new method for simulating collisional gravitational N-body systems with approximately linear time scaling with N. Our method is based on the multi-particle collision (MPC) scheme, previously applied in fluid dynamics and plasma physics. We were able to simulate globular clusters with a realistic number of stellar particles (at least up to several times 106) on a standard workstation. Aims. We simulated clusters hosting an intermediate mass black hole (IMBH), probing a broad range of BH-cluster and BH–average-star mass ratios, unrestricted by the computational constraints that affect direct N-body codes. Methods. We set up a grid of hybrid particle-in-cell-MPC N-body simulations using our implementation of the MPC method, MPCDSS. We used either single mass models or models with a Salpeter mass function (a single power law with an exponent of −2.35), with the IMBH initially sitting at the centre. The force exerted by and on the IMBH was evaluated with a direct sum scheme with or without softening. For all simulations we measured the evolution of the Lagrangian radii and core density and velocity dispersion over time. In addition, we also measured the evolution of the velocity anisotropy profiles. Results. We find that models with an IMBH undergo core collapse at earlier times, the larger the IMBH mass the shallower they are, with an approximately constant central density at core collapse. The presence of an IMBH tends to lower the central velocity dispersion. These results hold independently of the mass function of the model. For the models with Salpeter MF, we observed that equipartition of kinetic energies is never achieved, even long after core collapse. Orbital anisotropy at large radii appears to be driven by energetic escapers on radial orbits, triggered by strong collisions with the IMBH in the core. We measured the wander radius, that is the distance of the IMBH from the centre of mass of the parent system over time, finding that its distribution has positive kurtosis. Conclusions. Among the results we obtained, which mostly confirm or extend previously known trends that had been established over the range of parameters accessible to direct N-body simulations, we underline that the leptokurtic nature of the IMBH wander radius distribution might lead to IMBHs presenting as off-centre more frequently than expected, with implications on observational IMBH detection.

The most massive Population III stars

ABSTRACT Recent data from the JWST suggest that there are realistic prospects for detecting the earliest generation of stars at redshift ∼20. These metal-poor, gaseous Population III (Pop III) stars are likely in the mass range $10\!-\!10^3\, {\rm M}_{\odot }$ . We develop a framework for calculating the abundances of Pop III stars as well as the distribution of the most massive Pop III stars based on an application of extreme-value statistics. Our calculations use the star formation rate density from a recent simulation to calibrate the star formation efficiency from which the Pop III stellar abundances are derived. Our extreme-value modelling suggests that the most massive Pop III stars at redshifts 10 < z < 20 are likely to be $\gtrsim 10^3\!-\!10^4\, {\rm M}_\odot$ . Such extreme Pop III stars were sufficiently numerous to be the seeds of supermassive black holes at high redshifts and possibly source detectable gravitational waves. We conclude that the extreme-value formalism provides an effective way to constrain the stellar initial mass function.

The HH 24 Complex: Jets, Multiple Star Formation, and Orphaned Protostars

The HH 24 complex harbors five collimated jets emanating from a small protostellar multiple system. We have carried out a multiwavelength study of the jets, their driving sources, and the cloud core hosting the embedded stellar system, based on data from the Hubble Space Telescope, Gemini, Subaru, Apache Point Observatory 3.5 m, Karl G. Jansky Very Large Array, and Atacama Large Millimeter/submillimeter Array (ALMA) telescopes. The data show that the multiple system, SSV 63, contains at least 7 sources, ranging in mass from the hydrogen-burning limit to proto-Herbig Ae stars. The stars are in an unstable nonhierarchical configuration, and one member, a borderline brown dwarf, is moving away from the protostellar system with 25 km s−1, after being ejected ∼5800 yr ago as an orphaned protostar. Five of the embedded sources are surrounded by small, possibly truncated, disks resolved at 1.3 mm with ALMA. Proper motions and radial velocities imply jet speeds of 200–300 km s−1. The two main HH 24 jets, E and C, form a bipolar jet system that traces the innermost portions of parsec-scale chains of Herbig–Haro and H2 shocks with a total extent of at least 3 pc. H2CO and C18O observations show that the core has been churned and continuously fed by an infalling streamer. 13CO and 12CO trace compact, low-velocity, cavity walls carved by the jets and an ultracompact molecular outflow from the most embedded object. Chaotic N-body dynamics likely will eject several more of these objects. The ejection of stars from their feeding zones sets their masses. Dynamical decay of nonhierarchical systems can thus be a major contributor to establishing the initial mass function.

Formation of star clusters and enrichment by massive stars in simulations of low-metallicity galaxies with a fully sampled initial stellar mass function

ABSTRACT We present new griffin project hydrodynamical simulations that model the formation of galactic star cluster populations in low-metallicity (Z = 0.00021) dwarf galaxies, including radiation, supernova, and stellar wind feedback of individual massive stars. In the simulations, stars are sampled from the stellar initial mass function (IMF) down to the hydrogen-burning limit of 0.08 M⊙. Mass conservation is enforced within a radius of 1 pc for the formation of massive stars. We find that massive stars are preferentially found in star clusters and follow a correlation set at birth between the highest initial stellar mass and the star cluster mass that differs from pure stochastic IMF sampling. With a fully sampled IMF, star clusters lose mass in the galactic tidal field according to mass-loss rates observed in nearby galaxies. Of the released stellar feedback, 60 per cent of the supernova material and up to 35 per cent of the wind material reside either in the hot interstellar medium (ISM) or in gaseous, metal-enriched outflows. While stellar winds (instantaneously) and supernovae (delayed) start enriching the ISM right after the first massive stars form, the formation of supernova-enriched stars and star clusters is significantly delayed (by >50 Myr) compared to the formation of stars and star clusters enriched by stellar winds. Overall, supernova ejecta dominate the enrichment by mass, while the number of enriched stars is determined by continuous stellar winds. These results present a concept for the formation of chemically distinct populations of stars in bound star clusters, reminiscent of multiple populations in globular clusters.

Star cluster formation and survival in the first galaxies

ABSTRACT Using radiation-hydrodynamic cosmological simulations, we present a detailed (0.1 pc resolution), physically motivated portrait of a typical-mass dwarf galaxy before the epoch of reionization, resolving the formation, and evolution of star clusters into individual 10 M⊙ star particles. In the rest-frame ultraviolet, the galaxy has an irregular morphology with no bulge or disc, dominated by light emitted from numerous, compact, and gravitationally-bound star clusters. This is especially interesting in light of recent James Webb Space Telescope observations that − aided by the magnifying power of gravitational lenses – have imaged, at parsec-scale resolution, individual young star clusters forming in similar galaxies at z> 6. Because of their low metallicities and high temperatures, star-forming gas clouds in this galaxy have densities ∼100 times higher than typical giant molecular clouds; hence, their expected star formation efficiencies (SFEs) are high enough (around 10 − 70 per cent) to produce a sizeable population of potential globular cluster progenitors, but typically smaller (a few 100 − 2 ×104M⊙, half-mass radii of up to 3 pc) and of lower metallicities (10−3.5– 10−2.5 Z⊙). The initial mass function of the star-forming clouds is log-normal, whereas the bound star cluster mass function is a power-law with a slope that depends mainly on SFE but also on the temporal proximity to a major starburst. We find slopes between −0.5 and −2.5 depending on the assumed sub-grid SFE. Star formation is self-regulated on galactic scales; however, the multimodal metallicity distribution of the star clusters and the fraction of stars locked into surviving bound star clusters depends on SFE.

Abundances of CNO elements in z ∼ 0.3–0.4 Lyman continuum leaking galaxies

ABSTRACT We present observations with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope (HST) of 11 Lyman continuum (LyC) leaking galaxies at redshifts, z, in the range 0.29–0.43, with oxygen abundances 12+log(O/H) = 7.64–8.16, stellar masses M⋆ ∼ 107.8–109.8 M⊙, and O32 = [O iii] λ5007/[O ii] λ3727 of ∼ 5–20 aiming to detect the C iii] λ1908 emission line. We combine these observations with the optical Sloan Digital Sky Survey (SDSS) spectra for the determination of carbon, nitrogen, and oxygen abundances. Our sample was supplemented by 31 galaxies from the literature, for which carbon, nitrogen, and oxygen abundances can be derived from the HST and SDSS spectra. These additional galaxies, however, do not have LyC observations. We find that log(C/O) for the entire sample at 12+log(O/H) < 8.1 does not depend on metallicity, with a small dispersion of ∼0.13 dex around the average value of ∼−0.75 dex. On the other hand, the log(N/O) in galaxies at z > 0.1, including LyC leakers, is systematically higher compared to the rest of the sample with lower metallicity. We find that log(C/O) slightly decreases with increasing M⋆ from ∼ −0.65 at M⋆ = 106 M⊙ to ∼ −0.80 at M⋆ = 109–1010 M⊙, whereas log(N/O) is considerably enhanced at M⋆ > 108 M⊙. The origin of these trends remains basically unknown. A possible solution would be to assume that the upper mass limit of the stellar initial mass function in more massive galaxies is higher. This would result in a higher production of oxygen and a larger fraction of massive stars with stellar wind polluting the interstellar medium with nitrogen.

Local stellar formation history from the 40 pc white dwarf sample

ABSTRACT We derive the local stellar formation history from the Gaia-defined 40 pc white dwarf sample. This is currently the largest volume-complete sample of white dwarfs for which spectroscopy is available, allowing for classification of the chemical abundances at the photosphere, and subsequently accurate determination of the atmospheric parameters. We create a population synthesis model and show that a uniform stellar formation history for the last ≈10.5 Gyr provides a satisfactory fit to the observed distribution of absolute Gaia G magnitudes. To test the robustness of our derivation, we vary various assumptions in the population synthesis model, including the initial mass function, initial-to-final mass relation, kinematic evolution, binary fraction, and white dwarf cooling time-scales. From these tests, we conclude that the assumptions in our model have an insignificant effect on the derived relative stellar formation rate as a function of look-back time. However, the onset of stellar formation (age of Galactic disc) is sensitive to a variety of input parameters including the white dwarf cooling models. Our derived stellar formation history gives a much better fit to the absolute Gaia G magnitudes than most previous studies.

Global survey of star clusters in the Milky Way

Aims. We built Galactic open star cluster mass functions (CMFs) for different age sub-samples and spatial locations in the wider solar neighbourhood. Here, we present a simple cluster formation and evolution model to reproduce the main features of the CMFs. Methods. We used an unbiased working sample of 2227 clusters of the Milky Way Star Cluster (MWSC) catalogue, which occupy the heliocentric cylinders with magnitude-dependent completeness radii of 1–5 kpc. The MWSC survey provides an extended set of open star cluster parameters, including tidal radii, distances, and ages. From an analytic three-component Galaxy model, we derived tidal masses of clusters with a typical accuracy of about 70%. Our simple model includes a two-section cluster initial mass function, constant cluster formation rate, supervirial phase after a sudden expulsion of the remaining gas, and cluster mass loss due to stellar evolution and the clusters’ gradual destruction in the Galactic tidal field. The dynamical evolution model is based on previous N-body simulations. Results. The obtained tidal masses have been added to the MWSC catalogue. A general CMF (GCMF), built for all cluster ages around the Sun, has a bell-like shape and extends over four decades in mass. The high-mass slope found for tidal mass log mt/M⊙ ≥ 2.3 is equal to 1.14 ± 0.07. The CMFs for different age groups show the same high-mass slopes, while the low-mass slope is nearly flat for the youngest sub-sample (clusters younger than 20 Myr) and about −0.7 for the others. The inner and outer sub-samples covering Galactocentric radii R = 4.2–8.1 kpc and 8.9–13.5 kpc, respectively, are consistent with the GCMF, once the exponential decline of the Galactic disc density is taken into account. The model suggests star formation with low efficiency of 15–20%, where only 10% of stars remain bound in a cluster after gas expulsion and subsequent violent relaxation. The cluster formation rate required to reproduce the observed distributions in age and mass is about 0.4 M⊙ pc−2 Gyr−1. Conclusions. The obtained high-mass slope of the GCMF for the wide neighbourhood of the Sun is similar to slopes determined earlier in nearby galaxies for more luminous clusters with log m/M⊙ > 3.8. The MWSC catalogue supports models with a low star-formation efficiency, where 90% of stars are lost quickly after gas expulsion. The obtained cluster formation rate corresponds to open clusters’ contribution to the stellar content of the thin disc at the level of 30%.

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

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

An Analytical Theory for the Growth from Planetesimals to Planets by Polydisperse Pebble Accretion

Pebble accretion is recognized as a significant accelerator of planet formation. Yet only formulae for single-sized (monodisperse) distribution have been derived in the literature. These can lead to significant underestimates for Bondi accretion, for which the best accreted pebble size may not be the one that dominates the mass distribution. We derive in this paper the polydisperse theory of pebble accretion. We consider a power-law distribution in pebble radius, and we find the resulting surface and volume number density distribution functions. We derive also the exact monodisperse analytical pebble accretion rate for which 3D accretion and 2D accretion are limits. In addition, we find analytical solutions to the polydisperse 2D Hill and 3D Bondi limits. We integrate the polydisperse pebble accretion numerically for the MRN distribution, finding a slight decrease (by an exact factor 3/7) in the Hill regime compared to the monodisperse case. In contrast, in the Bondi regime, we find accretion rates 1–2 orders of magnitude higher compared to monodisperse, also extending the onset of pebble accretion to 1–2 orders of magnitude lower in mass. We find megayear timescales, within the disk lifetime, for Bondi accretion on top of planetary seeds of masses 10−6 to 10−4 M ⊕, over a significant range of the parameter space. This mass range overlaps with the high-mass end of the planetesimal initial mass function, and thus pebble accretion is possible directly following formation by streaming instability. This alleviates the need for mutual planetesimal collisions as a major contribution to planetary growth.

GMC collisions as triggers of star formation – VIII. The core mass function

ABSTRACTCompression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger the formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetized GMCs and examine the properties of pre-stellar cores, selected from projected mass surface density maps, including after synthetic ALMA observations. We then examine core properties, including mass, size, density, velocity, velocity dispersion, temperature, and magnetic field strength. After 4 Myr, ∼1000 cores have formed in the GMC collision, and the high-mass end of the core mass function (CMF) can be fit by a power-law dN/dlogM ∝ M−α with α ≃ 0.7, i.e. relatively top heavy compared to a Salpeter mass function. Depending on how cores are identified, a break in the power law can appear around a few $\times 10\, \mathrm{M}_\odot$. The non-colliding GMCs form fewer cores with a CMF with α ≃ 0.8–1.2, i.e. closer to the Salpeter index. We compare the properties of these CMFs to those of several observed samples of cores. Considering other properties, cores formed from colliding clouds are typically warmer, have more disturbed internal kinematics, and are more likely to be gravitational unbound, than cores formed from non-colliding GMCs. The dynamical state of the protocluster of cores formed in the GMC–GMC collision is intrinsically subvirial but can appear to be supervirial if the total mass measurement is affected by observations that miss mass on large scales or at low densities.

On Fixed-Accuracy Confidence Intervals for the Parameters of Lindley Distribution and Its Extensions

The purpose of the present paper is to deal with sequential estimation of the parameter θ in a Lindley distribution. A fixed-accuracy confidence interval for θ with a preassigned confidence coefficient is developed. It is established that, no fixed sample size procedure can solve the estimation problem and hence a purely sequential methodology is proposed to deal with the situation. The first-order asymptotic efficiency and consistency properties associated with our purely sequential strategy are derived. Similar estimation strategies are also outlined for a few other extensions of the Lindley distribution. Extensive simulation analysis is carried out to validate the theoretical findings. We also provide a real data example, where we estimate the parameter related to the “initial mass function" for a particular cluster of stars.

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.

Evidence for a bottom-light initial mass function in massive star clusters

ABSTRACT We have determined stellar mass functions of 120 Milky Way globular clusters and massive Large Magellanic Cloud/Small Magellanic Cloud star clusters based on a comparison of archival Hubble Space Telescope photometry with a large grid of direct N-body simulations. We find a strong correlation of the global mass function slopes of star clusters with both their internal relaxation times and their lifetimes. Once dynamical effects are being accounted for, the mass functions of most star clusters are compatible with an initial mass function described by a broken power-law distribution N(m) ∼ mα with break masses at 0.4 and 1.0 M⊙ and mass function slopes of αLow = −0.3 for stars with masses m < 0.4 M⊙, αHigh = −2.30 for stars with m > 1.0 M⊙, and αMed = −1.65 for intermediate-mass stars. Alternatively, a lognormal mass function with a characteristic mass log MC = −0.36 and width σC = 0.28 for low-mass stars and a power-law mass function for stars with m > 1 M⊙ also fit our data. We do not find a significant environmental dependence of the initial mass function on cluster mass, density, global velocity dispersion, or metallicity. Our results lead to a larger fraction of high-mass stars in globular clusters compared to canonical Kroupa/Chabrier mass functions, increasing the efficiency of self-enrichment in clusters and helping to alleviate the mass budget problem of multiple stellar populations in globular clusters. By comparing our results with direct N-body simulations, we finally find that only simulations in which most black holes are ejected by natal birth kicks correctly reproduce the observed correlations.

Dust-buried Compact Sources in the Dwarf Galaxy NGC 4449

Multiwavelength images from the Hubble Space Telescope covering the wavelength range 0.27–1.6 μm show that the central area of the nearby dwarf galaxy NGC 4449 contains several tens of compact sources that are emitting in the hydrogen recombination line Paβ (1.2818 μm) but are only marginally detected in Hα (0.6563 μm) and undetected at wavelengths λ ≤ 0.55 μm. An analysis of the spectral energy distributions (SEDs) of these sources indicates that they are likely relatively young stellar clusters heavily attenuated by dust. The selection function used to identify the sources prevents meaningful statistical analyses of their age, mass, and dust extinction distributions. However, these cluster candidates have ages ∼5–6 Myr and A V > 6 mag, according to their SED fits, and are extremely compact, with typical deconvolved radii of 1 pc. The dusty clusters are located at the periphery of the dark clouds within the galaxy and appear to be partially embedded. Density and pressure considerations indicate that the H ii regions surrounding these clusters may be stalled, and that pre-supernova (pre-SN) feedback has not been able to clear the clusters of their natal cocoons. These findings are in potential tension with existing models that regulate star formation with pre-SN feedback, since pre-SN feedback acts on short timescales, ≲4 Myr, for a standard stellar initial mass function. The existence of a population of dusty stellar clusters with ages >4 Myr, if confirmed by future observations, paints a more complex picture for the role of stellar feedback in controlling star formation.

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.

Early evolution and three-dimensional structure of embedded star clusters

ABSTRACT We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the Amuse framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in Flash. We simulate a suite of 104 M⊙ clouds at 0.0683 pc resolution for ∼2 Myr after the onset of star formation, with virial parameters α = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial α having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over ∼0.01 Myr time-scales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses ≲1000 M⊙ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation.