As sites of some of the most efficient star formation in the Universe, globular clusters (GCs) have long been hypothesized to be the building blocks of young galaxies. Within the Milky Way, our best tracers of the contribution of GCs to the proto-Galaxy are stars with such anomalous overabundance in nitrogen and depletion in oxygen ("high-[N/O] stars") that they can be identified as having originated in a cluster long after they have escaped. We identify associations between these high-[N/O] field stars and GCs using integrals of motion and metallicities and compare to chemically typical halo stars to quantify any excess association, enabling a population-level exploration of the formation sites of the nitrogen-enhanced stars in the field. Relative to the halo as a whole, high-[N/O] stars show stronger associations with the most initially massive, inner Galaxy GCs, suggesting that many nitrogen-rich stars formed in these environments. However, when compared to a sample matched in orbital energy, the excess largely disappears: high-[N/O] stars are, on average, no more associated with surviving GCs than energy-matched halo stars, despite their [N/O] abundances indicating GC origins, consistent with a scenario in which a substantial fraction of low-energy inner-halo stars originate in GCs, so an energy-matched control dilutes any differential excess. We argue that associations between high-[N/O] stars and their parent GCs are further weakened because dynamical friction and the Galactic bar have altered integrals of motion, limiting the reliability of precise present-day associations and, especially, individual star-to-cluster tagging.

Recent James Webb Space Telescope (JWST) observations have revealed a growing population of galaxies at z > 4 with elevated nitrogen-to-oxygen ratios. These 'N/O-enhanced' galaxies (NOEGs) exhibit near to supersolar N/O at sub-solar O/H, clearly deviating from the well-established scaling relation between N/O and O/H observed in local galaxies. The origin of this abundance anomaly is unclear. Interestingly, local globular clusters also exhibit anomalous light-element abundances, whose origin remains debated. In this work, we compare the chemical abundance patterns of 22 known NOEGs at 0 less than or similar to z less than or similar to 12-primarily discovered with JWST-to those observed in local globular clusters. We find similarities in the abundances of C, N, O, Fe, and He between the two populations. The similar abundance patterns support the scenario in which globular cluster stars formed within proto-cluster environments-similar to those traced by NOEGs-that were self-enriched. Indeed, the enhancement in N/O in early galaxies appears to be only found in dense stellar environments with Sigma(*) >= 10(2.5 )M (R) pc(-2), as expected for the progenitors of globular clusters in the Milky Way, and similar to those of star clusters identified in strongly lensed high-redshift galaxies. Furthermore, we find a tentative positive correlation between N/O ratios and stellar mass among NOEGs. The apparent high occurrence rate of NOEGs at high redshift is consistent with the picture of cluster-dominated star formation during the early stages of galaxy evolution. Measuring chemical abundances across diverse stellar environments in high-redshift galaxies will be crucial for elucidating the connection between NOEGs and globular clusters.

The stellar halo of the Milky Way comprises an abundance of chemical signatures from accretion events and in-situ evolution, that form an interweaving tapestry in kinematic space. To untangle this, we consider the mixtures of chemical information, in a given region of integral of motion space, as a variant of the blind source separation problem using non-negative matrix factorisation (NMF). Specifically, we examine the variation in [Fe/H], [Mg/Fe], and [Al/Fe] distributions of APOGEE DR17 stars across the (E,L_z) plane of the halo. When 2 components are prescribed, the NMF algorithm splits stellar halo into low- and high-energy components in the (E,L_z) plane which approximately correspond to the accreted and in-situ halo respectively. We use these components to define a boundary between the in-situ and the accreted stellar halo, and calculate their fractional contribution to the stellar halo as a function of energy, galactocentric spherical radius (r), height (z), and galactocentric cylindrical radius (R). Using a stellar halo defined by kinematic cuts, we derive a boundary in (E,L_z) space where the halo transitions from in-situ dominated to accretion dominated. Spatially, we find that this transition happens at (r,z,R) ≈ (8.7, 3.0, 8.1) kpc. We find that between 34% to 53% of the stellar halo's content is of accreted origin. Upon prescribing more components to the NMF model, we find evidence for overlapping chemical evolution sequences. We examine features within these components that resemble known substructures in the halo, such as Eos and Aurora.

We use intergalactic medium (IGM) metallicity distributions from several state-of-the-art cosmological simulations of Milky Way analogs and a semi-analytic model of ultra-faint dwarf galaxy (UFD) formation to model the stellar metallicities of UFDs in MW-like environments. We study simulations with different treatments of star formation, stellar feedback, and Population III enrichment, and in all cases, we find that only a few percent of the IGM accretable by UFD progenitors is enriched to metallicities [Fe/H]≥-4. When the metallicity of accreted IGM in the semi-analytic galaxy formation model is set using these IGM metallicity distributions, the model underpredicts UFD metallicities and their scatter compared to the observed luminosity–metallicity relation. Our results indicate that IGM enrichment is not the dominant mechanism setting metallicities of UFD stars. Instead, UFD stellar metallicity is determined primarily by the interplay between internal enrichment and metal loss through feedback-driven outflows. We examine models with different values of the maximum outflow mass loading factor η_ max and show that the full range of average stellar metallicities of UFDs at M_V<-7 can be reproduced if the maximum mass loading factor varies in the range 200_ max≲ 2000. We also consider stellar metallicity distribution functions (MDFs) within individual model galaxies with different assumptions about IGM enrichment and η_ max. We find that all considered models are in reasonable agreement with observed UFD MDFs, with model differences less than the uncertainties of current metallicity measurements.

The particle mass of dark matter (DM) was previously constrained using kinematics of ultra-faint dwarf galaxies to m > 3 × 10^-19 eV. This constraint, which excludes the "fuzzy" range of ultra-light dark matter from comprising all of the DM, relies on an estimate of the heating rate from fuzzy dark matter (FDM) wave interference using linear perturbation theory. Here, we compare the results of this perturbative calculation to full Schrödinger-Poisson simulations of the evolution of star particles in FDM halos. This comparison confirms theoretical expectations that FDM heating is stronger in fully nonlinear simulations due to the formation of a dense central soliton whose fluctuations enhance gravitational perturbations, and that bounds on the DM particle mass using this perturbative method are indeed conservative. We also show that these bounds are not affected by possible tidal stripping, since for dwarf satellites like Segue 1, the tidal radius is much larger than the observed size of the galaxy. We further show that the constraints on the mass cannot be evaded by invoking DM self-interactions, due to constraints on the self-interaction from large-scale structure. Lastly, we show that if the recently discovered system Ursa Major III/UNIONS I is a galaxy, the observed properties of this object strengthen the lower bound on the DM mass by over an order of magnitude, to m > 8 × 10^-18 eV, at 95

We present model calculations of the reionization history of hydrogen using star formation histories, computed with a galaxy formation model which reproduces properties of local dwarf galaxies and UV luminosity functions of galaxies at $z=5-16$. We use the ionizing photon density functions predicted by the model along with different models for the escape fraction of ionizing photons, $f_{\rm esc}$, to study the effects of ionizing photons from faint galaxies and different assumptions about $f_{\rm esc}$ on the evolution of hydrogen ionized fraction with redshift, $Q_{\rm HII}(z)$. We show that accounting for the contribution of faint galaxies with UV luminosities $M_{1500}>-13$, and with a constant ionizing photon escape fraction of $f_{\rm esc}=0.1$ results in the hydrogen reionization history consistent with all current observational constraints. Comparing results of the $f_{\rm esc}=0.1$ model and two alternative models shows that the model with a strong luminosity dependence of $f_{\rm esc}$, which assigns high $f_{\rm esc}$ to faint galaxies, results in early reionization inconsistent with observational constraints. However, the model in which $f_{\rm esc}$ follows a universal redshift-independent correlation with the recent maximum specific star formation rate, motivated by the results of the SPHINX galaxy formation simulation, results in the reionization history in good agreement with existing observational constraints, even though this model produces a sizeable ionized hydrogen fraction of $Q_{\rm HII}\approx 0.15-0.2$ at redshifts $z=8-12$. Our results show that the relative contribution of faint dwarf galaxies to reionization depends sensitively on assumptions about the escape fraction for galaxies of different luminosities, and that this is the main source of uncertainty in modeling hydrogen reionization.

Although the standard A & thorn; cold dark matter (ACDM) model is well tested on large scales, the primordial power spectrum may deviate from the ACDM spectrum on small scales due to specific dark matter properties or alternative inflationary models. These deviations affect the formation of dark matter structure, which subsequently leads to different observable properties of galaxies. In this work, we study the impact of a blue and red tilted power spectrum on the central density of dwarf galaxies. To do this, we model densities of dwarf galaxies using a combination of high-resolution numerical simulations and galaxy formation model. The model galaxies in ACDM are consistent with observations of 41 faint dwarf satellite galaxies of the Milky Way. The deviations from the ACDM power spectrum are constrained by the central matter densities of dwarf galaxies, which set stringent constraints on the possible small-scale tilt of the primordial power spectrum, improving on the current limits. Moreover, similar analysis can be applied to test any feature in the power spectrum at small scales between k similar to 10-100 Mpc-1.

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which include radiative transfer, a non-equilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H$_2$ abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the HI-H$_2$ transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the interstellar medium drops drastically at $Z \lesssim 0.2 \, Z_\odot$ as the formation time of H$_2$ becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected $f_\mathrm{H_2}$ measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy the overall molecular gas mass is reproduced to better than a factor of $\lesssim 1.5$ across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous $f_\mathrm{H_2}$ models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

The innermost portions of the Milky Way's stellar halo have avoided scrutiny until recently. The lack of wide-area survey data, made it difficult to reconstruct an uninterrupted view of the density distribution of the metal-poor stars inside the Solar radius. In this study, we utilize red giant branch (RGB) stars from Gaia, with metallicities estimated using spectro-photometry from Gaia Data Release 3. Accounting for Gaia's selection function, we examine the spatial distribution of metal-poor ([M/H]<-1.3) RGB stars, from the Galactic centre (r 1 kpc) out to beyond the Solar radius (r 18 kpc). Our best-fitting single-component cored power-law model shows a vertical flattening of 0.5 and a slope -3.4, consistent with previous studies. Motivated by the mounting evidence for two distinct stellar populations in the inner halo, we additionally test a range of two-component models. One of the components models the tidal debris from the Gaia Sausage/Enceladus merger, while the other captures the Aurora population – stars that predate the Galactic disk formation. Our best-fit two-component model suggests that both populations contribute equally around the Solar radius, but Aurora dominates the inner halo with a steeper power-law index of -4.5, in agreement with the nitrogen-rich star distribution measured by Horta et al. (2021).

We use a well-motivated galaxy formation framework to predict stellar masses, star formation rates (SFR), and ultraviolet (UV) luminosities of galaxy populations at redshifts z∈ 5-16, taking into account stochasticity of SFR in a controlled manner. We demonstrate that the model can match observational estimates of UV luminosity functions (LFs) at 51 indicating that SFR stochasticity cannot be higher. We discuss several testable consequences of the increased SFR stochasticity at z>10. The increase of SFR stochasticity with increasing z, for example, prevents steepening of UV LF and even results in some flattening of UV LF at z≳ 13. The median stellar ages of model galaxies at z≈ 11-16 are predicted to decrease from ≈ 20-30 Myr for M_ UV≳ -21 galaxies to ≈ 5-10 Myr for brighter ones. Likewise, the scatter in median stellar age is predicted to decrease with increasing luminosity. The scatter in the ratio of star formation rates averaged over 10 and 100 Myr should increase with redshift. Fluctuations of ionizing flux should increase at z>10 resulting in the increasing scatter in the line fluxes and their ratios for the lines sensitive to ionization parameter.

ABSTRACT We present a new scheme for the classification of the in-situ and accreted globular clusters (GCs). The scheme uses total energy E and z-component of the orbital angular momentum and is calibrated using the [Al/Fe] abundance ratio. We demonstrate that this classification results in two GC populations with distinct spatial, kinematic, and chemical abundance distributions. The in-situ GCs are distributed within the central 10 kpc of the Galaxy in a flattened configuration aligned with the Milky Way (MW) disc, while the accreted GCs have a wide distribution of distances and a spatial distribution close to spherical. In-situ and accreted GCs have different $\rm [Fe/H]$ distributions with the well-known bimodality present only in the metallicity distribution of the in-situ GCs. Furthermore, the accreted and in-situ GCs are well separated in the plane of $\rm [Al/Fe]-[Mg/Fe]$ abundance ratios and follow distinct sequences in the age–$\rm [Fe/H]$ plane. The in-situ GCs in our classification show a clear disc spin-up signature – the increase of median Vϕ at metallicities −1.3 < [Fe/H] < −1 similar to the spin-up in the in-situ field stars. This signature signals the MW’s disc formation, which occurred ≈11.7−12.7 Gyr ago (or at z ≈ 3.1−5.3) according to in-situ GC ages. In-situ GCs with metallicities of $\rm [Fe/H]\gtrsim -1.3$ were thus born in the MW disc, while lower metallicity in-situ GCs were born during early, turbulent, pre-disc stages of the evolution of the Galaxy and are part of its Aurora stellar component.

We analyze high-resolution hydrodynamics simulations of an isolated disk dwarf galaxy with an explicit model for unresolved turbulence and turbulence-based star formation prescription. We examine the characteristic values of the star formation efficiency per free-fall time, $\epsilon_\mathrm{ff}$, and its variations with local environment properties, such as metallicity, UV flux, and surface density. We show that the star formation efficiency per free-fall time in $\approx 10$ pc star-forming regions of the simulated disks has values in the range $\epsilon_\mathrm{ff}\approx 0.01-0.1$, similar to observational estimates, with no trend with metallicity and only a weak trend with the UV flux. Likewise, $\epsilon_{\rm ff}$ estimated using projected patches of 500 pc size does not vary with metallicity and shows only a weak trend with average UV flux and gas surface density. The characteristic values of $\epsilon_\mathrm{ff}\approx 0.01-0.1$ arise naturally in the simulations via the combined effect of dynamical gas compression and ensuing stellar feedback that injects thermal and turbulent energy. The compression and feedback regulate the virial parameter, $\alpha_\mathrm{vir}$, in star-forming regions, limiting it to $\alpha_\mathrm{vir}\approx 3-10$. Turbulence plays an important role in the universality of $\epsilon_\mathrm{ff}$ because turbulent energy and its dissipation are not sensitive to metallicity and UV flux that affect thermal energy. Our results indicate that the universality of observational estimates of $\epsilon_\mathrm{ff}$ can be plausibly explained by the turbulence-driven and feedback-regulated properties of star-forming regions.

We present estimates of the ultraviolet (UV) and Lyman continuum flux density contributed by galaxies of luminosities from $M_{\rm UV}\approx -25$ to $M_{\rm UV}=-4$ at redshifts $5\leq z\leq 10$ using a galaxy formation model that reproduces properties of local dwarf galaxies down to the luminosities of the ultra-faint satellites. We characterize the UV luminosity function (LF) of galaxies and their abundance as a function of the ionizing photon emission rate predicted by our model and present accurate fitting functions describing them. Although the slope of the LF becomes gradually shallower with decreasing luminosity due to feedback-driven outflows, the UV LF predicted by the model remains quite steep at the luminosities $M_{\rm UV}\lesssim -14$. After reionization, the UV LF flattens at $M_{\rm UV}\gtrsim -12$ due to UV heating of intergalactic gas. However, before reionization, the slope of the LF remains steep and approximately constant from $M_{\rm UV}\approx -14$ to $M_{\rm UV}=-4$. We show that for a constant ionizing photon escape fraction the contribution of faint galaxies with $M_{\rm UV}>-14$ to the UV flux and ionizing photon budget is $\approx 40-60\%$ at $z>7$ and decreases to $\approx 20\%$ at $z=6$. Before reionization, even ultra-faint galaxies of $M_{\rm UV}>-10$ contribute $\approx 10-25\%$ of ionizing photons. If the escape fraction increases strongly for fainter galaxies, the contribution of $M_{\rm UV}>-14$ galaxies before reionization increases to $\approx 60-75\%$. Our results imply that dwarf galaxies fainter than $M_{\rm UV}=-14$, beyond the James Webb Space Telescope limit, contribute significantly to the UV flux density and ionizing photon budget before reionization alleviating requirements on the escape fraction of Lyman continuum photons.

I compare the dark matter content within stellar half-mass radius expected in a $\Lambda$CDM-based galaxy formation model with existing observational estimates for the observed dwarf satellites of the Milky Way and ultra-diffuse galaxies (UDGs). The model reproduces the main properties and scaling relations of dwarf galaxies, in particular their stellar mass-size relation. I show that the model also reproduces the relation between the dark matter mass within the half-mass radius, $M_{\rm dm}(

We present a study illustrating the effects of the passage of a Large Magellanic Cloud (LMC) mass satellite on the distance and velocity distributions of satellites in Λ+Cold Dark Matter simulations of Milky Way (MW) sized halos. In agreement with previous studies, we find that during such a passage the velocity distribution develops a high-velocity tail, which can bias velocity-based virial halo mass estimates. When the velocity distribution of MW satellites is corrected for effects of the LMC passage, it is consistent with the distributions in halos of masses as low as and as high as 1.5×1012M⊙. We present a new halo mass estimator , where c is the coefficient calibrated using satellite systems in the simulated MW-sized halos, is the variance of 3D velocities taken with the sign of the radial velocity of each satellite, and is the median halocentric distance of the satellites. We show that the estimator has only s=8% scatter around the median relation of the estimated and true halo masses and deviates by <2s from the median during the pericentric passage of an LMC-like subhalo. This is because and deviate in the opposite directions during such passages. We apply the estimator to the MW satellite system and estimate the virial mass of the Milky Way of , in good agreement with several recent estimates using other methods.

Anomalously high nitrogen-to-oxygen abundance ratios [N/O] are observed in globular clusters (GCs), among the field stars of the Milky Way (MW), and even in the gas in a $z\approx 11$ galaxy. Using data from the APOGEE Data Release 17 and the Gaia Data Release 3, we present several independent lines of evidence that most of the MW's high-[N/O] stars were born in situ in massive bound clusters during the early, pre-disk evolution of the Galaxy. Specifically, we show that distributions of metallicity [Fe/H], energy, the angular momentum $L_z$, and distance of the low-metallicity high-[N/O] stars match the corresponding distributions of stars of the Aurora population and of the in-situ GCs. We also show that the fraction of in-situ field high-[N/O] stars, $f_{\rm N/O}$, increases rapidly with decreasing metallicity. During epochs when metallicity evolves from $\rm [Fe/H]=-1.5$ to $\rm [Fe/H]=-0.9$, the Galaxy spins up and transitions from a turbulent Aurora state to a coherently rotating disk. This transformation is accompanied by many qualitative changes. In particular, we show that high N/O abundances similar to those observed in GN-z11 were common before the spin-up ($\rm [Fe/H]\lesssim -1.5$) when up to $\approx 50\%-70\%$ of the in-situ stars formed in massive bound clusters. The dramatic drop of $f_{\rm N/O}$ at $\rm [Fe/H]\gtrsim -0.9$ indicates that after the disk emerges the fraction of stars forming in massive bound clusters decreases by two orders of magnitude.

ABSTRACT We examine the spatial distribution and orbital pole correlations of satellites in a suite of zoom-in high-resolution dissipationless simulations of Milky Way (MW)-sized haloes. We use the measured distribution to estimate the incidence of satellite configurations as flattened and as correlated in their orbital pole distribution as the satellite system of the MW. We confirm that this incidence is sensitive to the radial distribution of subhaloes and thereby to the processes that affect it, such as artificial disruption due to numerical effects and disruption due to the central disc. Controlling for the resolution effects and bracketing the effects of the disc, we find that the MW satellite system is somewhat unusual (at the ≈2–3σ level) but is statistically consistent with the Lambda cold dark matter model, in general agreement with results and conclusions of other recent studies.

We investigate the impact of bursty star formation on several galaxy scaling relations of dwarf galaxies using the $\texttt{GRUMPY}$ galaxy formation model. While this model reproduces the star formation rate (SFR)-stellar mass, stellar mass-gas mass, and stellar mass-metallicity relations, the scatter of these relations in the original model is smaller than observed. We explore the effects of additional stochasticity of SFR on the scaling relations using a model that reproduces the level of SFR burstiness in high-resolution zoom-in simulations. The additional SFR stochasticity increases the scatter in the SFR-stellar mass relation to a level similar to that exhibited by most nearby dwarf galaxies. The most extreme observed starbursting dwarfs, however, require higher levels of SFR stochasticity. We find that bursty star formation increases the scatter in the colour-magnitude distribution (CMD) for brighter dwarf galaxies $(M_V < -12)$ to the observed level, but not for fainter ones for which scatter remains significantly smaller than observed. This is due to the predominant old stellar populations in these faint model galaxies and their generally declining SFR over the past 10 Gyrs, rather than quenching caused by reionization. We examine the possibility that the colour scatter is due to scatter in metallicity, but show that the level of scatter required leads to an overestimation of scatter in the metallicity-mass relation. This illustrates that the scatter of observed scaling relations in the dwarf galaxy regime represents a powerful constraint on the properties of their star formation.

ABSTRACT We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 \, \mathrm{M}_\odot \lesssim M_\star \lesssim 10^{11} \, \mathrm{M}_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulations, which was originally designed to model Milky Way-like host galaxies, resemble the populations in the Exploration of Local VolumE Satellites (ELVES) Survey and the Satellites Around Galactic Analogs (SAGA) survey in their luminosity function in the luminosity range −12 ≲ MV ≲ −15 and resemble ELVES in their quenched fraction and colour–magnitude distribution in the luminosity range −12 ≲ Mg ≲ −15. We find that satellites transition from blue colours to red colours at the luminosity range −15 ≲ Mg ≲ −12 in both the simulations and observations and we show that this shift is driven by environmental effects in the simulations. We demonstrate also that the colour distribution in both simulations and observations can be decomposed into two statistically distinct populations based on their morphological type or star-forming status that are statistically distinct. In the simulations, these two populations also have statistically distinct infall time distributions. The comparison presented here seems to indicate that this tension is resolved by the improved target selection of ELVES, but there are still tensions in understanding the colours of faint galaxies, of which ELVES appears to have a significant population of faint blue satellites not recovered in Auriga.

ABSTRACT We use the GRUMPY galaxy formation model based on a suite of zoom-in, high-resolution, dissipationless Λ Cold Dark Matter (ΛCDM) simulations of the Milky Way (MW) sized haloes to examine total matter density within the half-mass radius of stellar distribution, ρtot(< r1/2), of satellite dwarf galaxies around the MW hosts and their mass assembly histories. We compare model results to ρtot(< r1/2) estimates for observed dwarf satellites of the Milky Way spanning their entire luminosity range. We show that observed MW dwarf satellites exhibit a trend of decreasing total matter density within a half-mass radius, ρtot(< r1/2), with increasing stellar mass. This trend is in general agreement with the trend predicted by the model. None of the observed satellites are overly dense compared to the results of our ΛCDM-based model. We also show that although the halo mass of many satellite galaxies is comparable to the halo mass of the MW progenitor at z ≳ 10, at these early epochs halos that survive as satellites to z = 0 are located many virial radii away from the MW progenitors and thus do not have a chance to merge with it. Our results show that neither the densities estimated in observed Milky Way satellites nor their mass assembly histories pose a challenge to the ΛCDM model. In fact, the broad agreement between density trends with the stellar mass of the observed and model galaxies can be considered as yet another success of the model.