Assembly History Of Galaxies Research Articles

Half-mass Radii of Quiescent and Star-forming Galaxies Evolve Slowly from 0 ≲ z ≤ 2.5: Implications for Galaxy Assembly Histories*

Abstract We use high-resolution, multiband imaging of ∼16,500 galaxies in the CANDELS fields at 0 ≲ z ≤ 2.5 to study the evolution of color gradients and half-mass radii over cosmic time. We find that galaxy color gradients at fixed mass evolve rapidly between z ∼ 2.5 and z ∼ 1, but remain roughly constant below z ∼ 1. This result implies that the sizes of both star-forming and quiescent galaxies increase much more slowly than previous studies found using half-light radii. The half-mass radius evolution of quiescent galaxies is fully consistent with a model that uses observed minor merger rates to predict the increase in sizes due to the accretion of small galaxies. Progenitor bias may still contribute to the growth of quiescent galaxies, particularly if we assume a slower timescale for the minor merger growth model. The slower half-mass radius evolution of star-forming galaxies is in tension with cosmological simulations and semianalytic galaxy models. Further detailed, consistent comparisons with simulations are required to place these results in context.

The SLUGGS survey: measuring globular cluster ages using both photometry and spectroscopy

ABSTRACT Globular cluster ages provide both an important test of models of globular cluster formation and a powerful method to constrain the assembly history of galaxies. Unfortunately, measuring the ages of unresolved old stellar populations has proven challenging. Here, we present a novel technique that combines optical photometry with metallicity constraints from near-infrared spectroscopy in order to measure ages. After testing the method on globular clusters in the Milky Way and its satellite galaxies, we apply our technique to three massive early-type galaxies using data from the SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey. The three SLUGGS galaxies and the Milky Way show dramatically different globular cluster age and metallicity distributions, with NGC 1407 and the Milky Way showing mostly old globular clusters, while NGC 3115 and NGC 3377 show a range of globular ages. This diversity implies different galaxy formation histories and that the globular cluster optical colour–metallicity relation is not universal as is commonly assumed in globular cluster studies. We find a correlation between the median age of the metal-rich globular cluster populations and the age of the field star populations, in line with models where globular cluster formation is a natural outcome of high-intensity star formation.

A machine-vision method for automatic classification of stellar halo substructure

ABSTRACT Tidal debris structures formed from disrupted satellites contain important clues about the assembly histories of galaxies. To date, studies of these structures have been hampered by reliance on by-eye identification and morphological classification which leaves their interpretation significantly uncertain. In this work, we present a new machine-vision technique based on the Subspace-Constrained Mean Shift (SCMS) algorithm which can perform these tasks automatically. SCMS finds the location of the high-density ‘ridges’ that define substructure morphology. After identification, the coefficients of an orthogonal series density estimator are used to classify points on the ridges as part of a continuum between shell-like or stream-like debris, from which a global morphological classification can be determined. We dub this procedure Subspace-Constrained Unsupervised Detection of Structure (SCUDS). By applying this tool to controlled N-body simulations of minor mergers, we demonstrate that the extracted classifications correspond to the well-understood underlying physics of phase mixing. The application of SCUDS to resolved stellar population data from near-future surveys will inform our understanding of the buildup of galaxies’ stellar haloes.

The E-MOSAICS project: tracing galaxy formation and assembly with the age–metallicity distribution of globular clusters

We present 25 cosmological zoom-in simulations of Milky-Way-mass galaxies in the ‘MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE’ (E-MOSAICS) project. E-MOSAICS couples a detailed physical model for the formation, evolution, and disruption of star clusters to the EAGLE galaxy formation simulations. This enables following the co-formation and co-evolution of galaxies and their star cluster populations, thus realizing the long-standing promise of using globular clusters (GCs) as tracers of galaxy formation and assembly. The simulations show that the age–metallicity distributions of GC populations exhibit strong galaxy-to-galaxy variations, resulting from differences in their evolutionary histories. We develop a formalism for systematically constraining the assembly histories of galaxies using GC age–metallicity distributions. These distributions are characterized through 13 metrics that we correlate with 30 quantities describing galaxy formation and assembly (e.g. halo properties, formation/assembly redshifts, stellar mass assembly time-scales, galaxy merger statistics), resulting in 20 statistically (highly) significant correlations. The GC age–metallicity distribution is a sensitive probe of the mass growth, metal enrichment, and minor merger history of the host galaxy. No such relation is found between GCs and major mergers, which play a sub-dominant role in GC formation for Milky-Way-mass galaxies. Finally, we show how the GC age–metallicity distribution enables the reconstruction of the host galaxy’s merger tree, allowing us to identify all progenitors with masses |$M_*\gtrsim 10^8\,{\rm M}_\odot$| for redshifts 1 ≤ z ≤ 2.5. These results demonstrate that cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies successfully unlock the potential of GCs as quantitative tracers of galaxy formation and assembly.

The Auriga stellar haloes: connecting stellar population properties with accretion and merging history

We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies.

Dark matter response to galaxy assembly history

Aims. It is well known that the presence of baryons affects the dark matter host haloes. Exploring the galaxy assembly history together with the dark matter haloes properties through time can provide a way to measure these effects. Methods. We have studied the properties of four Milky Way mass dark matter haloes from the Aquarius project during their assembly history, between z = 0 − 4. In this work, we used a published SPH run and the dark matter only counterpart as case studies. To asses the robustness of our findings, we compared them with one of the haloes run using a moving-mesh technique and different sub-grid scheme. Results. Our results show that the cosmic evolution of the dark matter halo profiles depends on the assembly history of the baryons. We find that the dark matter profiles do not significantly change with time, hence they become stable, when the fraction of baryons accumulated in the central regions reaches 80 per cent of its present mass within the virial radius. Furthermore, the mass accretion history shows that the haloes that assembled earlier are those that contain a larger amount of baryonic mass aforetime, which in turn allows the dark matter halo profiles to reach a stable configuration earlier. For the SPH haloes, we find that the specific angular momentum of the dark matter particles within the five per cent of the virial radius at z = 0, remains approximately constant from the time at which 60 per cent of the stellar mass is gathered. We have explored different theoretical and empirical models for the contraction of the haloes through redshift. A model to better describe the contraction of the haloes through redshift evolution must depend on the stellar mass content in the inner regions.

Self-consistent population spectral synthesis with FADO

Context.Spectral population synthesis (PS) is a fundamental tool in extragalactic research that aims to decipher the assembly history of galaxies from their spectral energy distribution (SED). Whereas this technique has led to key insights into galaxy evolution in recent decades, star formation histories (SFHs) inferred therefrom have been plagued by considerable uncertainties stemming from inherent degeneracies and the fact that until recently allPScodes were restricted to purely stellar fits, neglecting the essential contribution of nebular emission (ne). With the advent ofFADO(Fitting Analysis using Differential evolution Optimisation), the now possible self-consistent modelling of stellar andneopens new routes to the exploration of galaxy SFHs.Aims.The main goal of this study is to quantitatively explore the accuracy to whichFADOcan recover physical and evolutionary properties of galaxies and compare its output with that from purely stellarPScodes.Methods.FADOandSTARLIGHTwere applied to synthetic SEDs that track the spectral evolution of stars and gas in extinction-free mock galaxies of solar metallicity that form their stellar mass (M⋆) according to different parametric SFHs. Spectral fits were computed for two different set-ups that approximate the spectral range of SDSS and CALIFA (V500) data, using up to seven libraries of simple stellar population spectra in the 0.005–2.5Z⊙metallicity range.Results.Our analysis indicates thatFADOcan recover the key physical and evolutionary properties of galaxies, such asM⋆and mass- and light-weighted mean age and metallicity, with an accuracy better than 0.2 dex. This is the case even in phases of strongly elevated specific star formation rate (sSFR) and thus with considerablenecontamination (EW(Hα) > 103Å). Likewise, population vectors fromFADOadequately recover the mass fraction of stars younger than 10 Myr and older than 1 Gyr (M⋆<10Myr/M⋆totalandM⋆>1Gyr/M⋆total, respectively) and reproduce with a high fidelity the observed Hαluminosity. As forSTARLIGHT, our analysis documents a moderately good agreement with theoretical values only for evolutionary phases for whichnedrops to low levels (EW(Hα) ≤ 60 Å) which, depending on the assumed SFH, correspond to an age between ∼0.1 Gyr and 2–4 Gyr. However, fits withSTARLIGHTduring phases of high sSFR severely overestimate bothM⋆and the mass-weighted stellar age, whereas strongly underestimate the light-weighted age and metallicity. Furthermore, our analysis suggests a subtle tendency ofSTARLIGHTto favour a bi-modal SFH, as well a slightly overestimatedM⋆<10Myr/M⋆total, regardless of galaxy age. Whereas the amplitude of these biases can be reduced, depending on the specifics of the fitting procedure (e.g. accuracy and completeness of flagging emission lines, omission of the Balmer and Paschen jump from the fit), they persist even in the idealised case of a line-free SED comprising only stellar and nebular continuum emission.Conclusions.The insights from this study suggest that the neglect of nebular continuum emission inSTARLIGHTand similar purely stellarPScodes could systematically impactM⋆and SFH estimates for star-forming galaxies. We argue that these biases can be relevant in the study of a range of topics in extragalactic research, including the redshift-dependent slope of the star formation (SF) main sequence, the SF frosting hypothesis, and the regulatory role of supermassive black holes on the global SFH of galaxies.

KMTNet Nearby Galaxy Survey. I. Optimal Strategy for Low Surface Brightness Imaging with KMTNet

Abstract In hierarchical galaxy formation models, galaxies evolve through mergers and accretions. Tidally disrupted debris from these processes can remain as diffuse, faint structures, which can provide useful insight into the assembly history of galaxies. To investigate the properties of the faint structures in outskirts of nearby galaxies, we conduct deep and wide-field imaging survey with Korea Microlensing Telescope Network (KMTNet). We present our observing strategy and optimal data reduction process to recover the faint extended features in the imaging data of NGC 1291 taken with KMTNet. Through the dark sky flat-fielding and optimal sky subtraction, we can effectively remove inhomogeneous patterns. In the combined images, the peak-to-peak global sky gradients were reduced to less than ∼0.5% and ∼0.3% of the original B- and R-band sky levels, respectively. However, we find local spatial fluctuations in the background sky that can affect the precise measurement of the sky value. Consequently, we can reach the surface brightness of and mag arcsec−2 in azimuthally averaged one-dimensional surface brightness profiles, which is mainly limited by the uncertainty in the sky determination. These results suggest that the deep imaging data produced by KMTNet are suitable to study the faint features of nearby galaxies such as outer disks and dwarf companions, but unideal (not impossible) to detect stellar halos. The one-dimensional profile revealed that NGC 1291 appeared to have a Type I disk out to R ∼30 kpc with no obvious color gradient, and excess light due to a stellar halo was undetected.

Testing the uniqueness of gravitational lens mass models

The positions of images produced by the gravitational lensing of background sources provide unique insight in to galaxy-lens mass distribution. However, even quad images of extended sources are not able to fully characterize the central regions of the host galaxy. Most previous work has focused either on the radial density profile of the lenses or localized substructure clumps. Here, we concentrate on the azimuthal mass asymmetries near the image circle. The motivation for considering such mass inhomogeneities is that the transition between the central stellar dominated region and the outer dark matter dominated region, though well represented by a power law density profile, is unlikely to be featureless, and encodes information about the dynamical state and assembly history of galaxies. It also happens to roughly coincide with the Einstein radius. We ask if galaxies that have mass asymmetries beyond ellipticity can be modeled with simpler lenses, i.e., can complex mass distributions masquerade as simple elliptical+shear lenses? Our preliminary study indicates that for galaxies with elliptical stellar and dark matter distributions, but with no mass asymmetry, and an extended source filling the diamond caustic, an elliptical+shear lens model can reproduce the images well, thereby hiding the potential complexity of the actual mass distribution. For galaxies with non-zero mass asymmetry, the answer depends on the size and brightness distribution of the source, and its location within the diamond caustic. In roughly half of the cases we considered the mass asymmetries can easily evade detection.

Metallicity gradients in the globular cluster systems of early-type galaxies: in situ and accreted components?

Massive early-type galaxies typically have two subpopulations of globular clusters (GCs) which often reveal radial colour (metallicity) gradients. Collating gradients from the literature, we show that the gradients in the metal-rich and metal-poor GC subpopulations are the same, within measurement uncertainties, in a given galaxy. Furthermore, these GC gradients are similar in strength to the {\it stellar} metallicity gradient of the host galaxy. At the very largest radii (e.g. greater than 8 galaxy effective radii) there is some evidence that the GC gradients become flat with near constant mean metallicity. Using stellar metallicity gradients as a proxy, we probe the assembly histories of massive early-type galaxies with hydrodynamical simulations from the Magneticum suite of models. In particular, we measure the stellar metallicity gradient for the in-situ and accreted components over a similar radial range as those observed for GC subpopulations. We find that the in-situ and accreted stellar metallicity gradients are similar but have a larger scatter than the metal-rich and metal-poor GC subpopulations gradients in a given galaxy. We conclude that although metal-rich GCs are predominately formed during the in-situ phase and metal-poor GCs during the accretion phase of massive galaxy formation, they do not have a strict one-to-one connection.

On the connectivity of the cosmic web: theory and implications for cosmology and galaxy formation

ABSTRACT Cosmic connectivity and multiplicity, i.e. the number of filaments globally or locally connected to a given cluster is a natural probe of the growth of structure and in particular of the nature of dark energy. It is also a critical ingredient driving the assembly history of galaxies as it controls mass and angular momentum accretion. The connectivity of the cosmic web is investigated here via the persistent skeleton. This tool identifies topologically the ridges of the cosmic landscape which allows us to investigate how the nodes of the cosmic web are connected together. When applied to Gaussian random fields corresponding to the high-redshift universe, it is found that on average the nodes are connected to exactly κ = 4 neighbours in two dimensions and ∼6.1 in three dimensions. Investigating spatial dimensions up to d = 6, typical departures from a cubic lattice κ = 2d are shown to scale like the power 7/4 of the dimension. These numbers strongly depend on the height of the peaks: the higher the peak the larger the connectivity. Predictions from first principles based on peak theory are shown to reproduce well the connectivity and multiplicity of Gaussian random fields and cosmological simulations. As an illustration, connectivity is quantified in galaxy lensing convergence maps and large dark haloes catalogues. As a function of redshift and scale the mean connectivity decreases in a cosmology-dependent way. As a function of halo mass, it scales like 10/3 times the log of the mass. Implications on galactic scales are discussed.

Field spheroid-dominated galaxies in a Λ-CDM Universe

Context. Understanding the formation and evolution of early-type, spheroid-dominated galaxies is an open question within the context of the hierarchical clustering scenario, particularly in low-density environments. Aims. Our goal is to study the main structural, dynamical, and stellar population properties and assembly histories of field spheroid-dominated galaxies formed in a Λ-cold dark matter (Λ-CDM) scenario to assess to what extent they are consistent with observations. Methods. We selected spheroid-dominated systems from a Λ-CDM simulation that includes star formation (SF), chemical evolution, and supernova feedback. The sample is made up of 18 field systems with MStar ≲ 6 × 1010M⊙ that are dominated by the spheroid component. For this sample we estimated the fundamental relations of ellipticals and compared them with current observations. Results. The simulated spheroid galaxies have sizes that are in good agreement with observations. The bulges follow a Sersic law with Sersic indexes that correlate with the bulge-to-total mass ratios. The structural-dynamical properties of the simulated galaxies are consistent with observed Faber–Jackson, fundamental plane, and Tully–Fisher relations. However, the simulated galaxies are bluer and with higher star formation rates (SFRs) than the observed isolated early-type galaxies. The archaeological mass growth histories show a slightly delayed formation and more prominent inside-out growth mode than observational inferences based on the fossil record method. Conclusions. The main structural and dynamical properties of the simulated spheroid-dominated galaxies are consistent with observations. This is remarkable since our simulation has not been calibrated to match them. However, the simulated galaxies are blue and star-forming, and with later stellar mass growth histories compared to observational inferences. This is mainly due to the persistence of extended discs in the simulations. The need for more efficient quenching mechanisms able to avoid further disc growth and SF is required in order to reproduce current observational trends.

SDSS-IV MaNGA: stellar angular momentum of about 2300 galaxies: unveiling the bimodality of massive galaxy properties

We measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

First Results on the Cluster Galaxy Population from the Subaru Hyper Suprime-Cam Survey. III. Brightest Cluster Galaxies, Stellar Mass Distribution, and Active Galaxies

Abstract The unprecedented depth and area surveyed by the Subaru Strategic Program with the Hyper Suprime-Cam (HSC-SSP) have enabled us to construct and publish the largest distant cluster sample out to to date. In this exploratory study of cluster galaxy evolution from z = 1 to z = 0.3, we investigate the stellar mass assembly history of brightest cluster galaxies (BCGs), the evolution of stellar mass and luminosity distributions, the stellar mass surface density profile, as well as the population of radio galaxies. Our analysis is the first high-redshift application of the top N richest cluster selection, which is shown to allow us to trace the cluster galaxy evolution faithfully. Over the 230 deg2 area of the current HSC-SSP footprint, selecting the top 100 clusters in each of the four redshift bins allows us to observe the buildup of galaxy population in descendants of clusters whose mass is about . Our stellar mass is derived from a machine-learning algorithm, which is found to be unbiased and accurate with respect to the COSMOS data. We find very mild stellar mass growth in BCGs (about 35% between z = 1 and 0.3), and no evidence for evolution in both the total stellar mass–cluster mass correlation and the shape of the stellar mass surface density profile. We also present the first measurement of the radio luminosity distribution in clusters out to , and show hints of changes in the dominant accretion mode powering the cluster radio galaxies at .

Individual stellar haloes of massive galaxies measured to 100 kpc at 0.3 &lt; z &lt; 0.5 using Hyper Suprime-Cam

Massive galaxies display extended light profiles that can reach several hundreds of kilo parsecs. These stellar halos provide a fossil record of galaxy assembly histories. Using data that is both wide (~100 square degree) and deep (i>28.5 mag/arcsec^2 in i-band), we present a systematic study of the stellar halos of a sample of more than 3000 galaxies at 0.3 < z < 0.5 with $\log M_{\star}/M_{\odot} > 11.4$. Our study is based on high-quality (0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP), which enables us to individually estimate surface mass density profiles to 100 kpc without stacking. As in previous work, we find that more massive galaxies exhibit more extended outer profiles. When this extended light is not properly accounted for as a result of shallow imaging or inadequate profile modeling, the derived stellar mass function can be significantly underestimated at the highest masses. Across our sample, the ellipticity of outer light profiles increases substantially as we probe larger radii. We show for the first time that these ellipticity gradients steepen dramatically as a function of galaxy mass, but we detect no mass-dependence in outer color gradients. Our results support the two-phase formation scenario for massive galaxies in which outer envelopes are built up at late times from a series of merging events. We provide surface mass surface mass density profiles in a convenient tabulated format to facilitate comparisons with predictions from numerical simulations of galaxy formation.

Identifying the progenitors of present-day early-type galaxies in observational surveys: correcting ‘progenitor bias’ using the Horizon-AGN simulation

As endpoints of the hierarchical mass-assembly process, the stellar populations of local early-type galaxies encode the assembly history of galaxies over cosmic time. We use Horizon-AGN, a cosmological hydrodynamical simulation, to study the merger histories of local early-type galaxies and track how the morphological mix of their progenitors evolves over time. We provide a framework for alleviating `progenitor bias' -- the bias that occurs if one uses only early-type galaxies to study the progenitor population. Early-types attain their final morphology at relatively early epochs -- by $z\sim1$, around 60 per cent of today's early-types have had their last significant merger. At all redshifts, the majority of mergers have one late-type progenitor, with late-late mergers dominating at $z>1.5$ and early-early mergers becoming significant only at $z<0.5$. Progenitor bias is severe at all but the lowest redshifts -- e.g. at $z\sim0.6$, less than 50 per cent of the stellar mass in today's early-types is actually in progenitors with early-type morphology, while, at $z\sim2$, studying only early-types misses almost all (80 per cent) of the stellar mass that eventually ends up in local early-type systems. At high redshift, almost all massive late-type galaxies, regardless of their local environment or star-formation rate, are progenitors of local early-type galaxies, as are lower-mass (M$_\star$ $<$ 10$^{10.5}$ M$_{\odot}$) late-types as long as they reside in high density environments. In this new era of large observational surveys (e.g. LSST, JWST), this study provides a framework for studying how today's early-type galaxies have been built up over cosmic time.

The E-MOSAICS project: simulating the formation and co-evolution of galaxies and their star cluster populations

We introduce the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) project. E-MOSAICS incorporates models describing the formation, evolution, and disruption of star clusters into the EAGLE galaxy formation simulations, enabling the examination of the co-evolution of star clusters and their host galaxies in a fully cosmological context. A fraction of the star formation rate of dense gas is assumed to yield a cluster population; this fraction and the population's initial properties are governed by the physical properties of the natal gas. The subsequent evolution and disruption of the entire cluster population are followed accounting for two-body relaxation, stellar evolution, and gravitational shocks induced by the local tidal field. This introductory paper presents a detailed description of the model and initial results from a suite of 10 simulations of ∼L⋆ galaxies with disc-like morphologies at z = 0. The simulations broadly reproduce key observed characteristics of young star clusters and globular clusters (GCs), without invoking separate formation mechanisms for each population. The simulated GCs are the surviving population of massive clusters formed at early epochs (z ≳ 1–2), when the characteristic pressures and surface densities of star-forming gas were significantly higher than observed in local galaxies. We examine the influence of the star formation and assembly histories of galaxies on their cluster populations, finding that (at similar present-day mass) earlier-forming galaxies foster a more massive and disruption-resilient cluster population, while galaxies with late mergers are capable of forming massive clusters even at late cosmic epochs. We find that the phenomenological treatment of interstellar gas in EAGLE precludes the accurate modelling of cluster disruption in low-density environments, but infer that simulations incorporating an explicitly modelled cold interstellar gas phase will overcome this shortcoming.

COSMOS2015 photometric redshifts probe the impact of filaments on galaxy properties

The variations of galaxy stellar masses and colour-types with the distance to projected cosmic filaments are quantified using the precise photometric redshifts of the COSMOS2015 catalogue extracted from COSMOS field (2 deg$^{2}$). Realistic mock catalogues are also extracted from the lightcone of the cosmological hydrodynamical simulation Horizon-AGN. They show that the photometric redshift accuracy of the observed catalogue ($\sigma_z<0.015$ at $M_*>10^{10}{\rm M}_{\odot}$ and $z<0.9$) is sufficient to provide 2D filaments that closely match their projected 3D counterparts. Transverse stellar mass gradients are measured in projected slices of thickness 75 Mpc between $0.5< z <0.9$, showing that the most massive galaxies are statistically closer to their neighbouring filament. At fixed stellar mass, passive galaxies are also found closer to their filament while active star-forming galaxies statistically lie further away. The contributions of nodes and local density are removed from these gradients to highlight the specific role played by the geometry of the filaments. We find that the measured signal does persist after this removal, clearly demonstrating that proximity to a filament is not equivalent to proximity to an over-density. These findings are in agreement with gradients measured both in 2D or 3D in the Horizon-AGN simulation and those observed in the spectroscopic VIPERS survey (which rely on the identification of 3D filaments). They are consistent with a picture in which the influence of the geometry of the large-scale environment drives anisotropic tides which impact the assembly history of galaxies, and hence their observed properties.

Galaxies with Shells in the Illustris Simulation: Metallicity Signatures

Stellar shells are low surface brightness arcs of overdense stellar regions, extending to large galactocentric distances. In a companion study, we identified 39 shell galaxies in a sample of 220 massive ellipticals ( M 200 crit &gt; 6 × 10 12 M ⊙ ) from the Illustris cosmological simulation. We used stellar history catalogs to trace the history of each individual star particle inside the shell substructures, and we found that shells in high-mass galaxies form through mergers with massive satellites (stellar mass ratios μ stars ≳ 1 : 10 ). Using the same sample of shell galaxies, the current study extends the stellar history catalogs in order to investigate the metallicity of stellar shells around massive galaxies. Our results indicate that outer shells are often times more metal-rich than the surrounding stellar material in a galaxy’s halo. For a galaxy with two different satellites forming z = 0 shells, we find a significant difference in the metallicity of the shells produced by each progenitor. We also find that shell galaxies have higher mass-weighted logarithmic metallicities ([Z/H]) at 2– 4 R eff compared to galaxies without shells. Our results indicate that observations comparing the metallicities of stars in tidal features, such as shells, to the average metallicities in the stellar halo can provide information about the assembly histories of galaxies.

Fitting Analysis using Differential evolution Optimization (FADO):

The goal of population spectral synthesis (PSS) is to decipher from the spectrum of a galaxy the mass, age and metallicity of its constituent stellar populations. This technique has been established as a fundamental tool in extragalactic research. It has been extensively applied to large spectroscopic data sets, notably the SDSS, leading to important insights into the galaxy assembly history. However, despite significant improvements over the past decade, all current PSS codes suffer from two major deficiencies that inhibit us from gaining sharp insights into the star-formation history (SFH) of galaxies and potentially introduce substantial biases in studies of their physical properties (e.g., stellar mass, mass-weighted stellar age and specific star formation rate). These are i) the neglect of nebular emission in spectral fits, consequently, ii) the lack of a mechanism that ensures consistency between the best-fitting SFH and the observed nebular emission characteristics of a star-forming (SF) galaxy. In this article, we present FADO (Fitting Analysis using Differential evolution Optimization): a conceptually novel, publicly available PSS tool with the distinctive capability of permitting identification of the SFH that reproduces the observed nebular characteristics of a SF galaxy. This so-far unique self-consistency concept allows us to significantly alleviate degeneracies in current spectral synthesis. The innovative character of FADO is further augmented by its mathematical foundation: FADO is the first PSS code employing genetic differential evolution optimization. This, in conjunction with other unique elements in its mathematical concept (e.g., optimization of the spectral library using artificial intelligence, convergence test, quasi-parallelization) results in key improvements with respect to computational efficiency and uniqueness of the best-fitting SFHs.