Assembly History Research Articles

Peeking beneath the precision floor − II. Probing the chemo-dynamical histories of the potential globular cluster siblings, NGC 288 and NGC 362

ABSTRACT The assembly history of the Milky Way (MW) is a rapidly evolving subject, with numerous small accretion events and at least one major merger proposed in the MW’s history. Accreted alongside these dwarf galaxies are globular clusters (GCs), which act as spatially coherent remnants of these past events. Using high precision differential abundance measurements from our recently published study, we investigate the likelihood that the MW clusters NGC 362 and NGC 288 are galactic siblings, accreted as part of the Gaia-Sausage-Enceladus (GSE) merger. To do this, we compare the two GCs at the 0.01 dex level for 20 + elements for the first time. Strong similarities are found, with the two showing chemical similarity on the same order as those seen between the three LMC GCs, NGC 1786, NGC 2210-, and NGC 2257. However, when comparing GC abundances directly to GSE stars, marked differences are observed. NGC 362 shows good agreement with GSE stars in the ratio of Eu to Mg and Si, as well as a clear dominance in the r- compared to the s-process, while NGC 288 exhibits only a slight r-process dominance. When fitting the two GC abundances with a GSE-like galactic chemical evolution model, NGC 362 shows agreement with both the model predictions and GSE abundance ratios (considering Si, Ni, Ba, and Eu) at the same metallicity. This is not the case for NGC 288. We propose that the two are either not galactic siblings, or GSE was chemically inhomogeneous enough to birth two similar, but not identical clusters with distinct chemistry relative to constituent stars.

The collapse of atomically cooled primordial haloes – I. High Lyman–Werner backgrounds

ABSTRACT Pristine, atomically cooled haloes may be the sites of primordial quasar formation because atomic cooling triggers rapid baryon collapse that can create 104–105 M⊙ black hole seeds. However, no numerical simulation has ever followed the collapse of these haloes for the times required to form supermassive stars and direct-collapse black holes (DCBHs). We have now modelled baryon collapse in atomically cooled haloes with a wide range of spin parameters and assembly histories for times that are sufficient for DCBH formation. Fragmentation of accretion discs after ∼500 kyr is nearly ubiquitous in these haloes and in most cases leads to the formation of binary or multiple supermassive stellar systems. They also confirm that rapid baryon collapse proceeds for the times required for these stars to form DCBHs. Our simulations suggest that binary or even multiple DCBH formation was the rule rather than the exception in the primordial Universe.

Late-formed haloes prefer to host quiescent central galaxies – I. Observational results

ABSTRACT The star formation and quenching of central galaxies are regulated by the assembly histories of their host haloes. In this work, we use the central stellar mass to halo mass ratio as a proxy of halo formation time, and we devise three different models, from the physical hydrodynamical simulation to the empirical statistical model, to demonstrate its robustness. With this proxy, we inferred the dependence of the central galaxy properties on the formation time of their host haloes using the SDSS main galaxy sample, where central galaxies are identified with the halo-based group finder. We found that central galaxies living in late-formed haloes have higher quiescent fractions and lower spiral fractions than their early-formed counterparts by $\lesssim 8~{{\ \rm per\ cent}}$ . Finally, we demonstrate that the group finding algorithm has a negligible impact on our results.

SDSS-IV MaNGA: The Effect of Stellar Mass and Halo Mass on the Assembly Histories of Satellite Galaxies

We combine an unprecedented MaNGA sample of over 3000 passive galaxies in the stellar mass range 109–1012 M ⊙ with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M *) and the mass of the host halo (M h ). After controlling for M * and M h , the stacked spectra of centrals and satellites beyond the effective radius (r e ) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf, a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on M h . We find that satellites in high-M h halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low-M h halos, especially for M * = 109.5–1010.5 M ⊙. These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation.

The mass distribution of quasars in optical time-domain surveys

ABSTRACTThe determination of supermassive black hole (SMBH) masses is the key to understanding the host galaxy build-up and the SMBH mass assembly histories. The SMBH masses of non-local quasars are frequently estimated via the single-epoch virial black-hole mass estimators, which may suffer from significant biases. Here we demonstrate a new approach to infer the mass distribution of SMBHs in quasars by modelling quasar ultraviolet /optical variability. Our inferred black hole masses are systematically smaller than the virial ones by 0.3–0.6 dex; the ∼0.3 dex offsets are roughly consistent with the expected biases of the virial black-hole mass estimators. In the upcoming time-domain astronomy era, our methodology can be used to constrain the cosmic evolution of quasar mass distributions.

The Differential Assembly History of the Centers and Outskirts of Main-sequence Galaxies at z ∼ 2.3

We present a study of spatially resolved star formation histories (SFHs) for 60 z ∼ 2.3 main-sequence, star-forming galaxies selected from the MOSDEF spectroscopic survey in the GOODS-N field, with median stellar mass = 9.75 and spanning the range . Photometry is decomposed into a central and an outer spatial component using observed z F850LP − H F160W colors. The Prospector code is used to model spectral energy distributions for the center, outskirt, and integrated galaxy using Hubble Space Telescope/ACS and WFC3, Spitzer/IRAC, and ground-based photometry, with additional constraints on gas-phase metallicity and spectroscopic redshift from MOSDEF spectroscopy. For the low-resolution bands, spatially resolved photometry is determined with an iterative approach. The reconstructed SFHs indicate that the majority of galaxies with are observed while their central regions undergo relatively recent (<100 Myr) bursts of star formation, whereas the outskirts have a smooth, quasi-steady SFH that gently increases toward the redshift of observation. The enhanced star formation activity of the central parts is broadly consistent with the idea that it is produced by highly dissipative gas compaction and accretion. The wide range of central densities and sizes observed in the sample suggests that, for the selected galaxies, such a process has started but is still far from being completed. The implication would be that selecting star-forming galaxies at cosmic noon frequently includes systems in an “evolved” evolutionary phase where the centers have recently started a burst of star formation activity that will likely initiate inside-out quenching in the next several hundred million years.

Great balls of FIRE II: The evolution and destruction of star clusters across cosmic time in a Milky Way-mass galaxy

ABSTRACT The current generation of galaxy simulations can resolve individual giant molecular clouds, the progenitors of dense star clusters. But the evolutionary fate of these young massive clusters, and whether they can become the old globular clusters (GCs) observed in many galaxies, is determined by a complex interplay of internal dynamical processes and external galactic effects. We present the first star-by-star N-body models of massive (N ∼ 105–107) star clusters formed in a FIRE-2 MHD simulation of a Milky Way-mass galaxy, with the relevant initial conditions and tidal forces extracted from the cosmological simulation. We select 895 (∼30 per cent) of the YMCs with &amp;gt;6 × 104 M⊙ from Grudić et al. 2022 and integrate them to z = 0 using the cluster Monte Carlo code, CMC. This procedure predicts a MW-like system with 148 GCs, predominantly formed during the early, bursty mode of star formation. Our GCs are younger, less massive, and more core-collapsed than clusters in the Milky Way or M31. This results from the assembly history and age-metallicity relationship of the host galaxy: Younger clusters are preferentially born in stronger tidal fields and initially retain fewer stellar-mass black holes, causing them to lose mass faster and reach core collapse sooner than older GCs. Our results suggest that the masses and core/half-light radii of GCs are shaped not only by internal dynamical processes, but also by the specific evolutionary history of their host galaxies. These results emphasize that N-body studies with realistic stellar physics are crucial to understanding the evolution and present-day properties of GC systems.

The history of Shahrisabz city neighborhood assemblies during the years of independence.

The history of Shahrisabz city neighborhood assemblies during the years of independence.

Unravelling the mass spectrum of destroyed dwarf galaxies with the metallicity distribution function

ABSTRACTAccreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the ‘stellar haloes’ of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass–metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L &amp;gt; Lhost/100. We find that the MW stellar halo has N ∼ 1−3 massive progenitors (L ≳ 108L⊙) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of 2. Future data sets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass scale for both the MW and its satellites.

The SAMI survey: evidence for dynamical coupling of ionized gas and young stellar populations

ABSTRACT We explore local and global dynamical differences between the kinematics of ionized gas and stars in a sample of galaxies from Data Release 3 of the SAMI Galaxy Survey. We find better agreement between local (i.e. comparing on a spaxel-to-spaxel basis) velocities and dispersion of gas and stars in younger systems as with previous work on the asymmetric drift in galaxies, suggesting that the dynamics of stars and ionized gas are initially coupled. The intrinsic scatter around the velocity and dispersion relations increases with increasing stellar age and mass, suggesting that subsequent mechanisms, such as internal processes, divergent star formation, and assembly histories, also play a role in setting and altering the dynamics of galaxies. The global (flux-weighted) dynamical support of older galaxies is hotter than in younger systems. We find that the ionized gas in galaxies is almost always dynamically colder than the stars with a steeper velocity gradient. In absolute terms, the local difference in velocity dispersion is more pronounced than the local difference in velocity, possibly reflecting inherent differences in the impact of turbulence, inflow and/or feedback on gas compared to stars. We suggest how these findings may be taken into account when comparing high and low redshift galaxy samples to infer dynamical evolution.

The fragility of thin discs in galaxies – I. Building tailored N-body galaxy models

ABSTRACT Thin stellar discs on both galactic and nuclear, sub-kpc scales are believed to be fragile structures that would be easily destroyed in major mergers. In turn, this makes the age dating of their stellar populations a useful diagnostics for the assembly history of galaxies. We aim at carefully exploring the fragility of such stellar discs in intermediate- and low-mass encounters, using high-resolution N-body simulations of galaxy models with structural and kinematic properties tailored to actually observed galaxies. As a first but challenging step, we create a dynamical model of FCC 170, a nearly edge-on galaxy in the Fornax cluster with multiple galactic components and including both galactic-scale and nuclear stellar discs (NSDs), using detailed kinematic data from the Multi Unit Spectroscopic Explorer and a novel method for constructing distribution function-based self-consistent galaxy models. We then create N-body realizations of this model and demonstrate that it remains in equilibrium and preserves its properties over many Gyr, when evolved with a sufficiently high particle number. However, the NSD is more prone to numerical heating, which gradually increases its thickness by up to 22 per cent in 10 Gyr even in our highest resolution runs. Nevertheless, these N-body models can serve as realistic representations of actual galaxies in merger simulations.

On the edge: the relation between stellar and dark matter haloes of Milky Way-mass galaxies

ABSTRACTWe investigate the build-up of the accreted stellar and dark matter haloes of Milky Way-like galaxies in APOSTLE suite of cosmological hydrodynamics simulations. We show that the stellar halo is made up primarily of stars stripped from a small number of massive dwarfs, most of which are disrupted by the present day. The dark matter halo, on the other hand, is made up primarily of small unresolved subhaloes (≲106 M⊙) and a ‘smooth’ component consisting of particles which were never bound to a subhalo. Despite these differences, the massive dwarfs that make up the majority of the stellar halo also contribute a significant fraction of the dark matter. The stars and dark matter stripped from these dwarfs are related through their kinematics and this leaves imprints in the phase-space structure of the haloes. We examine the relation between the location of features, such as caustics, in the phase space of the stars and dark halo properties. We show that the ‘edge’ of the stellar halo is a probe of dark matter halo mass and assembly history. The edges of Milky Way-mass galaxies should be visible at a surface brightness of 31–36 mag arcsec−2.

Relating galaxies across different redshift to study galaxy evolution

ABSTRACT We propose a general framework leveraging the galaxy–halo connection to link galaxies observed at different redshift in a statistical way, and use the link to infer the redshift evolution of the galaxy population. Our tests based on hydrodynamic simulations show that our method can accurately recover the stellar mass assembly histories up to z ∼ 3 for present star-forming and quiescent galaxies down to 1010 h−1 M⊙. Applying the method to observational data shows that the stellar mass evolution of the main progenitors of galaxies depends strongly on the properties of descendants, such as stellar mass, halo mass, and star formation states. Galaxies hosted by low-mass groups/haloes at the present time have since z ∼ 1.8 grown their stellar mass ∼2.5 times as fast as those hosted by massive clusters. This dependence on host halo mass becomes much weaker for descendant galaxies with similar star formation states. Star-forming galaxies grow about 2–4 times faster than their quiescent counterparts since z ∼ 1.8. Both TNG and EAGLE simulations overpredict the progenitor stellar mass at z &amp;gt; 1, particularly for low-mass descendants.

The accretion history of the Milky Way – I. How it shapes globular clusters and dwarf galaxies

ABSTRACT Halo inhabitants are individual stars, stellar streams, star, and globular clusters, and dwarf galaxies. Here we compare the two last categories that include objects of similar stellar mass, which are often studied as self-dynamical equilibrium systems. We discover that the half-light radius of globular clusters depends on their orbital pericentre and total energy, and that Milky Way (MW) tides may explain the observed correlation. We also suggest that the accretion epoch of stellar systems in the MW halo can be calibrated by the total orbital energy, and that such a relation is due to both the mass growth of the MW and dynamical friction affecting mostly satellites with numerous orbits. This calibration starts from the bulge to Kraken, Gaia Sausage Enceladus, Sagittarius stellar systems, and finally to the new coming dwarfs, either or not linked to the vast-polar structure. The most eccentric globular clusters, and dwarfs have their half-light radius scaling as the inverse of their binding energy, and this over more than two decades. This means that earlier arriving satellites are smaller due to the tidal effects of the MW. Therefore, most halo inhabitants appear to have their structural parameters shaped by MW tides, and also by ram-pressure for the most recent arrivals, the dwarf galaxies. The correlations found in this study can be used as tools to further investigate the origin of globular clusters and dwarfs, as well as the assembly history of our Galaxy.

Applying unsupervised learning to resolve evolutionary histories and explore the galaxy--halo connection in IllustrisTNG

ABSTRACT We examine the effectiveness of identifying distinct evolutionary histories in IllustrisTNG-100 galaxies using unsupervised machine learning with Gaussian mixture models. We focus on how clustering compressed metallicity histories and star formation histories produces sub-population of galaxies with distinct evolutionary properties (for both halo mass assembly and merger histories). By contrast, clustering with photometric colours fails to resolve such histories. We identify several populations of interest that reflect a variety of evolutionary scenarios supported by the literature. Notably, we identify a population of galaxies inhabiting the upper red sequence, M* &amp;gt; 1010 M⊙, that has a significantly higher ex-situ merger mass fraction present at fixed masses and a star formation history that has yet to fully quench, in contrast to an overlapping, satellite-dominated population along the red sequence, which is distinctly quiescent. Extending the clustering to study four clusters instead of three further divides quiescent galaxies, whereas star-forming ones are mostly contained in a single cluster, demonstrating a variety of supported pathways to quenching. In addition to these populations, we identify a handful of populations from our other clusters that are readily applicable to observational surveys, including a population related to post-starburst galaxies, allowing for possible extensions of this work in an observational context, and to corroborate results within the IllustrisTNG ecosystem.

On the choice of the most suitable indicator for the assembly state of dark matter haloes through cosmic time

ABSTRACT The dynamical state and morphological features of galaxies and galaxy clusters, and their high-redshift precursors, are tightly connected with their assembly history, encoding crucial information about the formation and evolution of such cosmic structures. As a first step towards finding an optimal indicator of the assembly state of observed structures, we use a cosmological simulation of a moderate volume to critically examine the best definition of an indicator that is able to discriminate dark matter haloes undergoing mergers and/or strong accretion from haloes experimenting a relaxed evolution. Using a combination of centre offset, virial ratio, mean radial velocity, sparsity, and ellipticity of the dark matter halo, we study how the thresholds on these parameters, as well as their relative weights, should evolve with redshift to provide the best classification possible. This allows us to split a sample of haloes in a totally relaxed, a marginally relaxed and an unrelaxed subsamples. The resulting classification strongly correlates with the merging activity obtained from the analysis of complete merger trees extracted from whole simulation data. The results on how the different indicators depend on redshift and halo mass, and their optimal combination to better match the true assembly history of haloes, could constitute relevant hints to find a suitable set of indicators applicable to observational data.

DS+: A method for the identification of cluster substructures

Context. The study of cluster substructures is important for the determination of the cluster dynamical status, assembly history, and the evolution of cluster galaxies, and it allows us to set constraints on the nature of dark matter and cosmological parameters. Aims. We present and test DS+, a new method for the identification and characterization of group-sized substructures in clusters. Methods. Our new method is based on the projected positions and line-of-sight (l.o.s. hereafter) velocities of cluster galaxies, and it is an improvement and extension of the traditional method of Dressler &amp; Shectman (1988, AJ, 95, 985). We tested it on cluster-size cosmological halos extracted from the IllustrisTNG simulations, with virial masses 14 ≲ log(M200/M⊙) ≲ 14.6 that contain ~190 galaxies on average. We also present an application of our method on a real data set, the Bullet cluster. Results. DS+ is able to identify ~80% of real group galaxies as members of substructures, and at least 60% of the galaxies assigned to substructures belong to real groups. The physical properties of the real groups are significantly correlated with those of the corresponding detected substructures, but with significant scatter, and they are overestimated on average. Application of the DS+ method to the Bullet cluster confirms the presence and main properties of the high-speed collision and identifies other substructures along the main cluster axis. Conclusions. DS+ proves to be a reliable method for the identification of substructures in clusters. The method is made freely available to the community as a Python code.

Digging into the Galactic Bulge: Stellar Population and Structure of the Poorly Studied Cluster NGC 6316

High-resolution Hubble Space Telescope optical observations have been used to analyze the stellar population and the structure of the poorly investigated bulge globular cluster NGC 6316. We constructed the first high-resolution reddening map in the cluster direction, which allowed us to correct the evolutionary sequences in the color–magnitude diagram (CMD) for the effects of differential reddening. A comparison between the CMDs of NGC 6316 and 47 Tucanae revealed strikingly similar stellar populations, with the two systems basically sharing the same turnoff, subgiant branch, and horizontal branch morphologies, indicating comparable ages. The red giant branch in NGC 6316 appears slightly bluer than in 47 Tucanae, suggesting a lower metal content. This has been confirmed by the isochrone fitting of the observed CMD, which provided us with updated values of the cluster age, distance, average color excess, and metallicity. We estimated an absolute age of 13.1 ± 0.5 Gyr, consistent with the age of 47 Tucanae, an average color excess E(B − V) = 0.64 ± 0.01, and a true distance modulus (m − M)0 = 15.27 ± 0.03 that sets the cluster distance at 11.3 kpc from the Sun. In addition, the photometric estimate of the cluster metallicity suggests [Fe/H] ≈ −0.9, which is ∼0.2 dex smaller than that of 47 Tucanae. We also determined the gravitational center and the density profile of the system from resolved stars. The latter is well reproduced by a King model. Our results confirm that NGC 6316 is another extremely old relic of the assembly history of the Galaxy.

Community assembly history alters relationships between biodiversity and ecosystem functions during restoration.

Relationships between biodiversity and ecosystem functioning depend on the processes structuring community assembly. However, predicting biodiversity-ecosystem functioning (BEF) relationships based on community assembly remains challenging because assembly outcomes are often contingent on history and the consequences of history for ecosystem functions are poorly understood. In a grassland restoration experiment, we isolated the role of history for the relationships between plant biodiversity and multiple ecosystem functions by initiating assembly in three different years, while controlling for all other aspects of community assembly. We found that two aspects of assembly history-establishment year and succession-altered species and trait community trajectories, which in turn altered net primary productivity, decomposition rates, and floral resources. Moreover, history altered BEF relationships (which ranged from positive to negative), both within and across functions, by modifying the causal pathways linking species identity, traits, diversity, and ecosystem functions. Our results show that the interplay of deterministic succession and environmental stochasticity during establishment mediate historical contingencies that cause variation in biodiversity and ecosystem functions, even under otherwise identical assembly conditions. An explicit attention to history is needed to understand why biodiversity-ecosystem function relationships vary in natural ecosystems: a critical question at the intersection of fundamental theory and applications to environmental change biology and ecosystem restoration.

Dust contribution to the panchromatic galaxy emission

ABSTRACT We have developed a pipeline called mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the Dusty SAGEsemi-analytic galaxy formation model. Dusty SAGEincorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically motivated attenuation model from the computed dust properties in Dusty SAGE, so each galaxy has a self-consistent set of attenuation parameters based on the complicated dust physics that occurred across the galaxy’s assembly history. Then, we explore several dust emission templates to produce infrared spectra. Our results show that a physically-motivated attenuation model is better for obtaining a consistent multiwavelength description of galaxy formation and evolution, compared to using a constant attenuation. We compare our predictions with a compilation of observations and find that the fiducial model is in reasonable agreement with: (i) the observed z= 0 luminosity functions from the far-ultraviolet to far-infrared simultaneously, and hence (ii) the local cosmic SED in the same range, (iii) the rest-frame K-band luminosity function across 0 &amp;lt; z&amp;lt; 3, and (iv) the rest-frame far-ultraviolet luminosity function across 0 &amp;lt; z&amp;lt; 1. Our model underproduces the far-ultraviolet emission at z= 2 and z= 3, which can be improved by altering the AGN feedback and dust processes in Dusty SAGE, However, this combination thus worsens the agreement at z= 0, which suggests that more detailed treatment of such processes is required.