Search

This paper presents the third data release of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 52 ultra-compact massive galaxies (UCMGs) observed with the X-Shooter spectrograph. We measure integrated stellar velocity dispersion, [Mg/Fe] abundances, ages, and metallicities for all the INSPIRE objects. We thus infer star formation histories and confirm the existence of a degree of relicness (DoR), defined in terms of the fraction of stellar mass formed by $z=2$, the time at which a galaxy has assembled 75\% of its mass, and the final assembly time. Objects with a high DoR assembled their stellar mass at early epochs, while low-DoR objects show a non-negligible fraction of later-formed populations and hence a spread in ages and metallicities. A higher DoR correlates with larger [Mg/Fe], super-solar metallicity, and larger velocity dispersion values. The 52 UMCGs span a large range of DoR from 0.83 to 0.06, with 38 of them having formed more than 75\% of their mass by $z=2$. Of these, nine are extreme relics (DoR$>0.7$), since they formed the totality ($>99\%$) of their stellar mass by redshift $z=2$. The remaining 14 UCMGs cannot be considered relics, as they are characterised by more extended star formation histories. With INSPIRE, we built the first sizeable sample of relics outside the local Universe, up to $z\sim0.4$, increasing the number of confirmed relics by a factor of $>10$, and opening up an important window to explain the mass assembly of massive galaxies in the high-z Universe., Comment: accepted for publication on MNRAS, 20 pages, 15 figures, 3 tables

In the last decade, growing evidence has emerged supporting a non-universal stellar Initial Mass Function (IMF) in massive galaxies, with a larger number of dwarf stars with respect to the Milky-Way (bottom-heavy IMF). However, a consensus about the mechanisms that cause IMF variations is yet to be reached. Recently, it has been suggested that stars formed early-on in cosmic time, via a star formation burst, could be characterised by a bottom-heavy IMF. A promising way to confirm this is to use relics, ultra-compact massive galaxies, almost entirely composed by these "pristine" stars. The INSPIRE Project aims at assembling a large sample of confirmed relics, that can serve as laboratory to investigate on the conditions of star formation in the first 1-3 Gyr of the Universe. In this third INSPIRE paper, we build a high signal-to-noise spectrum from five relics and one from five galaxies with similar sizes, masses, and kinematical properties, but characterised by a more extended star formation history (non-relics). Our detailed stellar population analysis suggests a systematically bottom-heavier IMF slope for relics than for non-relics, adding new observational evidence for the non-universality of the IMF at various redshifts and further supporting the above proposed physical scenario., Comment: Accepted for publication in MNRAS. 8 pages, 5 figures

This is the second data release (DR2) of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 21 new systems with observations completed before March 2022. For each system, we release four one-dimensional (1D) spectra to the ESO Science Archive, one spectrum for each arm of the X-Shooter spectrograph. In this paper, we focus on the line-of-sight velocity distribution, measuring integrated stellar velocity dispersions from the spectra, and assessing their robustness and the associated uncertainties. For each of the 21 new systems, we systematically investigated the effect of the parameters and set-ups of the full spectral fitting on the stellar velocity dispersion ($\sigma$) measurements. In particular, we tested how $\sigma$ changes when several parameters of the fit as well as the resolution and spectral coverage of the input spectra are varied. We found that the effect that causes the largest systematic uncertainties on $\sigma$ is the wavelength range used for the fit, especially for spectra with a lower signal-to-noise ratio (S/N $\leq$ 30). When using blue wavelengths (UVB arm) one generally underestimates the velocity dispersion (by $\sim$15 km/s). The values obtained from the near-IR (NIR) arm present a larger scatter because the quality of the spectra is lower. We finally compared our results with those in literature, finding a very good agreement overall. Joining results obtained in DR1 with those presented here, INSPIRE contains 40 ultra-compact massive galaxies, corresponding to 75% of the whole survey. By plotting these systems in a stellar mass-velocity dispersion diagram, we identify at least four highly reliable relic candidates among the new systems. Their velocity dispersion is larger than that of normal-sized galaxies of similar stellar mass., Comment: Accepted for publication, and in press on A&A

We apply a population-orbit superposition method to 16 galaxies in the Fornax cluster observed with MUSE/VLT in the context of the Fornax3D project. By fitting the luminosity distribution, stellar kinematics, and age and metallicity maps simultaneously, we obtained the internal stellar orbit distribution, as well as the age and metallicity distribution of stars on different orbits for each galaxy. Based on the model, we decompose each galaxy into a dynamically cold disk (orbital circularity $\lambda_z\ge0.8$) and a dynamically hot non-disk component (orbital circularity $\lambda_z<0.8$), and obtain the surface-brightness, age, and metallicity radial profiles of each component. The galaxy infall time into the cluster is strongly correlated with galaxy cold-disk age with older cold disks in ancient infallers. We quantify the infall time $t_{\rm infall}$ of each galaxy with its cold-disk age using a correlation calibrated with TNG50 cosmological simulations. For galaxies in the Fornax cluster, we found that the luminosity fraction of cold disk in galaxies with $t_{\rm infall}>8$ Gyr are a factor of $\sim 4$ lower than in more recent infallers while controlling for total stellar mass. Nine of the 16 galaxies have spatially extended cold disks, and most of them show positive or zero age gradients; stars in the inner disk are $\sim 2-5$ Gyr younger than that in the outer disk, in contrast to the expectation of inside-out growth. Our results indicate that the assembly of cold disks in galaxies is strongly affected by their infall into clusters, by either removal of gas in outer regions or even tidally stripping or heating part of the pre-existing disks. Star formation in outer disks can stop quickly after the galaxy falls into the cluster, while star formation in the inner disks can last for a few Gyrs more, building the positive age gradient measured in cold disks., Comment: 31 pages, 32 figures, Accepted to A&A

We study the evolution of kinematically-defined stellar discs in 10 Fornax-like clusters identified in the TNG50 run from the IllustrisTNG suite of cosmological simulations. We considered disc galaxies with present-day stellar mass $M_{\star}\geq 3 \times 10^{8} M_{\odot}$ and follow their evolution since first entering their host cluster. Very few stellar discs survive since falling in such dense environments, ranging from 40% surviving to all being disrupted. Such survival rates are consistent with what reported earlier for the two more massive, Virgo-like clusters in TNG50. In absolute terms, however, the low number of present-day disc galaxies in Fornax-like clusters could be at odds with the presence of three edge-on disc galaxies in the central regions of the actual Fornax cluster, as delineated by the Fornax3D survey. When looking at the Fornax analogues from random directions and with the same selection function of Fornax3D, the probability of finding three edge-on disc galaxies in any one Fornax-like cluster in TNG50 is rather low, albeit not impossible. We also compared the stellar-population properties near the equatorial plane derived from integral-field spectroscopy for the three edge-ons in Fornax to similar line-of-sight integrated values for present-day disc galaxies in TNG50. For one of these, the very old and metal-rich stellar population of its disc cannot be matched by any the disc galaxies in TNG50, including objects in the field. We discuss possible interpretations of these findings, while pointing to future studies on passive cluster spirals as a way to further test state-of-the-art cosmological simulations., Comment: 12 pages, 7 figures. Accepted on 21 of October, published on 25 of October

This work is based on high quality integral-field spectroscopic data obtained with the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT). The 21 brightest ($m_B\leq 15$ mag) early-type galaxies (ETGs) inside the virial radius of the Fornax cluster are observed out to distances of $\sim2-3\ R_{\rm e}$. Deep imaging from the VLT Survey Telescope (VST) is also available for the sample ETGs. We investigate the variation of the galaxy structural properties as a function of the total stellar mass and cluster environment. Moreover, we correlate the size scales of the luminous components derived from a multi-component decomposition of the VST surface-brightness radial profiles of the sample ETGs with the MUSE radial profiles of stellar kinematic and population properties. The results are compared with both theoretical predictions and previous observational studies and used to address the assembly history of the massive ETGs of the Fornax cluster. We find that galaxies in the core and north-south clump of the cluster, which have the highest accreted mass fraction, show milder metallicity gradients in their outskirts than the galaxies infalling into the cluster. We also find a segregation in both age and metallicity between the galaxies belonging to the core and north-south clump and the infalling galaxies. The new findings fit well within the general framework for the assembly history of the Fornax cluster., Comment: 23 pages, 7 figures, 3 tables. Accepted for publication in A&A

[arXiv Abridged] In this work, we explore new spatially-resolved measurements of the IMF for three edge-on lenticular galaxies in the Fornax cluster. Specifically, we utilise existing orbit-based dynamical models, which re-produce the measured stellar kinematics, in order to fit the new IMF maps within this orbital framework. We then investigate correlations between intrinsic orbital properties and the local IMF. We find that, within each galaxy, the high-angular-momentum, disk-like stars exhibit an IMF which is rich in dwarf stars. The centrally-concentrated pressure-supported orbits have IMF which are similarly rich in dwarf stars. Conversely, orbits at large radius which have intermediate angular momentum exhibit IMF which are markedly less dwarf-rich relative to the other regions of the same galaxy. Assuming that the stars which, in the present-day, reside on dynamically-hot orbits at large radii are dominated by accreted populations, we can interpret these findings as a correlation between the dwarf-richness of a population of stars, and the mass of the host in which it formed. Specifically, deeper gravitational potentials would produce more dwarf-rich populations, resulting in the relative deficiency of dwarf stars which originated in the lower-mass accreted satellites. Conversely, the central and high angular-momentum populations are likely dominated by in-situ stars, which were formed in the more massive host itself. There are also global differences between the three galaxies studied here, of up to $\sim 0.3\ \mathrm{dex}$ in the IMF parameter $\xi$. We find no local dynamical or chemical property which alone can fully account for the IMF variations., Comment: 11 pages. Accepted to MNRAS

The role played by environment in galaxy evolution is a current debate in astronomy. The degree to which environment can alter, re-shape, or drive galaxy evolution is a topic of discussion in both fronts, observations and simulations. This paper analyses the gas metallicity gradients for a sample of 10 Fornax cluster galaxies observed with MUSE as part of the Fornax3D project. Detailed maps of emission lines allowed a precise determination of gas metallicity and metallicity gradients. The integrated gas metallicity of our Fornax cluster galaxies show slightly higher metallicities (~0.045 dex) in comparison to a control sample. In addition, we find signs of a mass and metallicity segregation from the center to the outskirts of the cluster. By comparing our Fornax cluster metallicity gradients with a control sample we find a general median offset of ~0.04 dex/Re, with 8 of our galaxies showing flatter or more positive gradients. We find no systematic difference between the gradients of recent and intermediate infallers when considering the projected distance of each galaxy to the cluster center. To identify the origin of the observed offset in the metallicity gradients, we perform a similar analysis with data from the TNG50 simulation. We identify 12 subhalos in Fornax-like clusters and compare their metallicity gradients with a control sample of field subhalos. This exercise also shows a flattening in the metallicity gradients for galaxies in Fornax-like halos, with a median offset of ~0.05 dex/Re We also analyse the merger history, Mach numbers (M), and ram pressure stripping of our TNG50 sample. We conclude that the observed flattening in metallicity gradients is likely due to a combination of galaxies traveling at supersonic velocities (M>1) that are experiencing high ram pressure stripping and flybys., Comment: Accepted by A&A

The stellar initial mass function (IMF) is central to our interpretation of astronomical observables and to our understanding of most baryonic processes within galaxies. The universality of the IMF, suggested by observations in our own Milky Way, has been thoroughly revisited due to the apparent excess of low-mass stars in the central regions of massive quiescent galaxies. As part of the efforts within the Fornax 3D project, we aim to characterize the two-dimensional IMF variations in a sample of 23 quiescent galaxies within the Fornax cluster. For each galaxy in the sample, we measured the mean age, metallicity, [Mg/Fe], and IMF slope maps from spatially resolved integrated spectra. The IMF maps show a variety of behaviors and internal substructures, roughly following metallicity variations. However, metallicity alone is not able to fully explain the complexity exhibited by the IMF maps. In particular, for relatively metal-poor stellar populations, the slope of the IMF seems to depend on the (specific) star formation rate at which stars were formed. Moreover, metallicity maps have systematically higher ellipticities than IMF slope ones. At the same time, both metallicity and IMF slope maps have at the same time higher ellipticities than the stellar light distribution in our sample of galaxies. In addition we find that, regardless of the stellar mass, every galaxy in our sample shows a positive radial [Mg/Fe] gradient. This results in a strong [Fe/H]-[Mg/Fe] relation, similar to what is observed in nearby, resolved galaxies. Since the formation history and chemical enrichment of galaxies are causally driven by changes in the IMF, our findings call for a physically motivated interpretation of stellar population measurements based on integrated spectra that take into account any possible time evolution of the stellar populations., Comment: 14 pages and 10 figures (plus Appendix including all the stellar population maps). Accepted for publication in Astronomy and Astrophysics

Extragalactic planetary nebulae (PNe) offer a way to determine the distance to their host galaxies thanks to the nearly universal shape of the planetary nebulae luminosity function (PNLF). Accurate PNe distance measurements rely on obtaining well-sampled PNLFs and the number of observed PNe scales with the encompassed stellar mass. This means either disposing of wide-field observations or focusing on the bright central regions of galaxies. In this work we take this second approach and conduct a census of the PNe population in the central regions of galaxies in the Fornax cluster, using VLT/MUSE data for the early-type galaxies observed over the course of the Fornax3D survey. Using such integral-field spectroscopic observations to carefully separate the nebular emission from the stellar continuum, we isolated [OIII] 5007 {\AA} sources of interest, filtered out unresolved impostor sources or kinematic outliers, and present a catalogue of 1350 unique PNe sources across 21 early-type galaxies, which includes their positions, [OIII] 5007 {\AA} line magnitudes, and line-of-sight velocities. Using the PNe catalogued within each galaxy, we present independently derived distance estimates based on the fit to the entire observed PNLF observed while carefully accounting for the PNe detection incompleteness. With these individual measurements, we arrive at an average distance to the Fornax cluster itself of 19.86 $\pm$ 0.32 Mpc ($\mu_{PNLF}$ = 31.49 $\pm$ 0.04 mag). Our PNLF distance measurements agree well with previous distances based on surface brightness fluctuations, finding no significant systematic offsets between the two methods as otherwise reported in previous studies., Comment: 26 pages, 16 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics Journal. Corrected Fig. 7 legend

Context. Extragalactic Planetary Nebulae (PNe) are useful distance indicators and are often used to trace the dark-matter content in external galaxies. At the same time, PNe can also be used as probes of their host galaxy stellar populations and to help understanding the later stages of stellar evolution. Previous works have indicated that specific number of PNe per stellar luminosity can vary across different galaxies and as a function of stellar-population properties, for instance increasing with decreasing stellar metallicity. Aims. In this study we further explore the importance of stellar metallicity in driving the properties of the PNe population in early-type galaxies, using three edge-on galaxies in the Fornax cluster offering a clear view into their predominantly metal-rich and metal-poor regions near the equatorial plane or both below and above it, respectively . Methods. Using VLT-MUSE integral-field observations and dedicated PNe detection procedures, we construct the PNe luminosity function and compute the luminosity-specific number of PNe alpha in both in- and off-plane regions of our edge-on systems. Results. Comparing these alpha values with metallicity measurements also based on the same MUSE data, we find no evidence for an increase in the specific abundance of PNe when transitioning between metal-rich and metal-poor regions. Conclusions. Our analysis highlights the importance of ensuring spatial consistency to avoid misleading results when investigating the link between PNe and their parent stellar populations and suggest that in passively-evolving systems variations in the specific number of PNe may pertain to rather extreme metallicity regimes found either in the innermost or outermost regions of galaxies.

Nuclear star clusters (NSCs) are the densest stellar systems in the Universe and are found in the centres of all types of galaxies. They are thought to form via mergers of star clusters such as ancient globular clusters (GCs) that spiral to the centre as a result of dynamical friction or through in-situ star formation directly at the galaxy centre. There is evidence that both paths occur, but the relative contribution of either channel and their correlation with galaxy properties are not yet constrained observationally. We aim to derive the dominant NSC formation channel for a sample of 25 nucleated galaxies, mostly in the Fornax galaxy cluster, with stellar masses between $M_\rm{gal} \sim 10^8$ and $10^{10.5} M_\odot$ and NSC masses between $M_\rm{NSC} \sim 10^5$ and $10^{8.5} M_\odot$. Using Multi-Unit Spectroscopic Explorer (MUSE) data from the Fornax 3D survey and the ESO archive, we derive star formation histories, mean ages and metallicities of NSCs, and compare them to the host galaxies. In many low-mass galaxies, the NSCs are significantly more metal-poor than the hosts with properties similar to GCs. In contrast, in the massive galaxies, we find diverse star formation histories and cases of ongoing or recent in-situ star formation. Massive NSCs ($> 10^7 M_\odot$) occupy a different region in the mass-metallicity diagram than lower mass NSCs and GCs, indicating a different enrichment history. We find a clear transition of the dominant NSC formation channel with both galaxy and NSC mass. We hypothesise that while GC-accretion forms the NSCs of the dwarf galaxies, central star formation is responsible for the efficient mass build up in the most massive NSCs in our sample. At intermediate masses, both channels can contribute. The transition between these formation channels seems to occur at galaxy masses $M_\rm{gal} \sim 10^9 M_\odot$ and NSC masses $M_\rm{NSC} \sim 10^7 M_\odot$., Comment: accepted for publication in Astronomy & Astrophysics, 18 pages, 10 figures (including appendix)

Abridged for arXiv: In this work, we apply a powerful new technique in order to observationally derive accurate assembly histories through a self-consistent combined stellar dynamical and population galaxy model. We present this approach for three edge-on lenticular galaxies from the Fornax3D project -- FCC 153, FCC 170, and FCC 177 -- in order to infer their mass assembly histories individually and in the context of the Fornax cluster. The method was tested on mock data from simulations to quantify its reliability. We find that the galaxies studied here have all been able to form dynamically-cold (intrinsic vertical velocity dispersion $\sigma_z \lesssim 50\ {\rm km}\ {\rm s}^{-1}$) stellar disks after cluster infall. Moreover, the pre-existing (old) high angular momentum components have retained their angular momentum (orbital circularity $\lambda_z > 0.8$) through to the present day. Comparing the derived assembly histories with a comparable galaxy in a low-density environment -- NGC 3115 -- we find evidence for cluster-driven suppression of stellar accretion and merging. We measured the intrinsic stellar age--velocity-dispersion relation and find that the shape of the relation is consistent with galaxies in the literature across redshift. There is tentative evidence for enhancement in the luminosity-weighted intrinsic vertical velocity dispersion due to the cluster environment. But importantly, there is an indication that metallicity may be a key driver of this relation. We finally speculate that the cluster environment is responsible for the S0 morphology of these galaxies via the gradual external perturbations, or `harassment', generated within the cluster., Comment: Accepted for publication in A&A

Studying the stellar kinematics of galaxies is a key tool in the reconstruction of their evolution. However, the current measurements of the stellar kinematics are complicated by several factors, including dust extinction and the presence of multiple stellar populations. We use integral field spectroscopic data of four galaxies from the TIMER survey to explore and compare the kinematics measured in different spectral regions that are sensitive to distinct stellar populations. We derive the line-of-sight velocity and velocity dispersion of both a young (<2 Gyr) and an old stellar population from the spectral regions around the H$\beta$ line and the Ca II Triplet. In addition we obtain colour excess, mean age, and metallicity. We report a correlation of the colour excess with the difference in the kinematic parameters of the H$\beta$ line and the Ca II Triplet range, which are dominated by young and old stellar populations, respectively. Young stellar populations, located primarily in nuclear rings, have higher velocity dispersions than old ones. These differences in the rings are typically 10 km/s in velocity dispersion, but up to a mean value of 24 km/s in the most extreme case. Trends with age exist in the nuclear rings but are less significant than those with dust extinction. We report different degrees of correlation of these trends among the galaxies in the sample, which are related to the size of the Voronoi bins in their rings. No clear trends for the difference of line-of-sight velocity are observed. The absence of these trends can be explained as a consequence of the masking process of the H$\beta$ line during the kinematic extraction, as confirmed by dedicated simulations. Our study demonstrates that kinematic differences caused by different stellar populations can be identified in the central regions of nearby galaxies even from intermediate resolution spectroscopy., Comment: 22 pages, 24 figures

Based on very deep photometry, Surface Brightness Fluctuations (SBF) have been traditionally used to determine galaxy distances. We have recently computed SBF spectra of stellar populations at moderately high resolution, which are fully based on empirical stars. We show that the SBF spectra provide an unprecedented potential for stellar population studies that, so far, have been tackled on the basis of the mean fluxes. We find that the SBFs are able to unveil metal-poor stellar components at the one percent level, which are not possible to disentangle with the standard analysis. As these metal-poor components correspond to the first stages of the chemical enrichment, the SBF analysis provides stringent constrains on the quenching epoch., Comment: Proceedings of the IAU Symposium No. 359: Galaxy evolution and feedback across different environments. Editors: T. Storchi-Bergmann, R. Overzier, W. Forman & R. Riffel

Current instruments and spectral analysis programs are now able to decompose the integrated spectrum of a stellar system into distributions of ages and metallicities. The reliability of these methods have rarely been tested on nearby systems with resolved stellar ages and metallicities. Here we derive the age-metallicity distribution of M54, the nucleus of the Sagittarius dwarf spheroidal galaxy, from its integrated MUSE spectrum. We find a dominant old (8-14 Gyr), metal-poor (-1.5 dex) and a young (1 Gyr), metal-rich (+0.25 dex) component - consistent with the complex stellar populations measured from individual stars in the same MUSE data set. There is excellent agreement between the (mass-weighted) average age and metallicity of the resolved and integrated analyses. Differences are only 3% in age and 0.2 dex metallicitiy. By co-adding individual stars to create M54's integrated spectrum, we show that the recovered age-metallicity distribution is insensitive to the magnitude limit of the stars or the contribution of blue horizontal branch stars - even when including additional blue wavelength coverage from the WAGGS survey. However, we find that the brightest stars can induce the spurious recovery of an old ($>8$ Gyr), metal-rich (+0.25 dex) stellar population, which is otherwise not expected from our understanding of chemical enrichment in M54. The overall derived stellar mass-to-light ratio of M54 is M/L$_{\mathrm{V}}=1.46$ with a scatter of 0.22 across the field-of-view, which we attribute to the stochastic contribution of a young, metal-rich component. These findings provide strong evidence that complex stellar population distributions can be reliably recovered from integrated spectra of extragalactic systems., Comment: 28 pages, 15 figures incl. appendix, accepted for publication in ApJ

Globular clusters (GCs) are found ubiquitously in massive galaxies and due to their old ages, they are regarded as fossil records of galaxy evolution. Spectroscopic studies of GC systems are often limited to the outskirts of galaxies, where GCs stand out against the galaxy background and serve as bright tracers of galaxy assembly. In this work, we use the capabilities of the Multi Unit Explorer Spectrograph (MUSE) to extract a spectroscopic sample of 722 GCs in the inner regions ($\lesssim 3 R_\text{eff}$) of 32 galaxies in the Fornax cluster. These galaxies were observed as part of the Fornax 3D project, a MUSE survey that targets early and late-type galaxies within the virial radius of Fornax. After accounting for the galaxy background in the GC spectra, we extracted line-of-sight velocities and determined metallicities of a sub-sample of 238 GCs. We found signatures of rotation within GC systems, and comparing the GC kinematics and that of the stellar body shows that the GCs trace the spheroid of the galaxies. While the red GCs prove to closely follow the metallicity profile of the host galaxy, the blue GCs show a large spread of metallicities but they are generally more metal-poor than the host., Comment: 13 pages + appendix, accepted for publication in A&A, updated to published version

Globular cluster (GC) systems of massive galaxies often show a bimodal colour distribution. This has been interpreted as a metallicity bimodality, created by a two-stage galaxy formation where the red, metal-rich GCs were formed in the parent halo and the blue metal-poor GCs were accreted. This interpretation, however, crucially depends on the assumption that GCs are exclusively old stellar systems with a linear colour-metallicity relation (CZR). The shape of the CZR and range of GC ages are currently under debate, because their study requires high quality spectra to derive reliable stellar population properties. We determined metallicities with full spectral fitting from a sample of 187 GCs with high spectral signal-to-noise ratio in 23 galaxies of the Fornax cluster that were observed as part of the Fornax 3D project. The derived CZR from this sample is non-linear and can be described by a piecewise linear function with a break point at ($g - z$) $\sim$ 1.1 mag. The less massive galaxies in our sample ($M_\ast < 10^{10} M_\odot$) appear to have slightly younger GCs, but the shape of the CZR is insensitive to the GC ages. Although the least massive galaxies lack red, metal-rich GCs, a non-linear CZR is found irrespective of the galaxy mass, even in the most massive galaxies ($M_\ast \geq 10^{11} M_\odot$). Our CZR predicts narrow unimodal GC metallicity distributions for low mass and broad unimodal distributions for very massive galaxies, dominated by a metal-poor and metal-rich peak, respectively, and bimodal distributions for galaxies with intermediate masses (10$^{10}$ $\leq$ $M_\ast < 10^{11} M_\odot$) as a consequence of the relative fraction of red and blue GCs. The diverse metallicity distributions challenge the simple differentiation of GC populations solely based on their colour., Comment: 12 pages + appendix, updated to published version

Extragalactic planetary nebulae (PNe) are detectable through relatively strong nebulous [OIII] emission and act as direct probes into the local stellar population. Because they have an apparently universal invariant magnitude cut-off, PNe are also considered to be a remarkable standard candle for distance estimation. Through detecting PNe within the galaxies, we aim to connect the relative abundances of PNe to the properties of their host galaxy stellar population. By removing the stellar background components from FCC 167 and FCC 219, we aim to produce PN luminosity functions (PNLF) of these galaxies, and thereby also estimate the distance modulus to these two systems. Finally, we test the reliability and robustness of our novel detection and analysis method. It detects unresolved point sources by their [OIII] 5007{\AA} emission within regions that have previously been unexplored. We model the [OIII] emissions in the spatial and spectral dimensions together, as afforded to us by the Multi Unit Spectroscopic Explorer (MUSE), and we draw on data gathered as part of the Fornax3D survey. For each source, we inspect the properties of the nebular emission lines to remove other sources that might hinder the safe construction of the PNLF, such as supernova remnants and HII regions. As a further step, we characterise any potential limitations and draw conclusions about the reliability of our modelling approach through a set of simulations. By applying this novel detection and modelling approach to integral field unit observations, we report for the distance estimates and luminosity-specific PNe frequency values for the two galaxies. Furthermore, we include an overview of source contamination, galaxy differences, and possible effects on the PNe populations in the dense stellar environments., Comment: 16 pages, 17 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics Journal. Replaced with published version

We present a new set of Surface Brightness Fluctuation spectra computed with the E-MILES stellar population synthesis models. The model SBF spectra cover the range 1680-50000 Angstrom at moderately high resolution, all based on extensive empirical stellar libraries. The models span the metallicity range -2.3<[M/H<+0.26 for a suite of IMF types with varying slopes. These predictions can complement and aid fluctuation magnitudes studies, permitting a first order approximation by applying filter responses to the SBF spectra to obtain spectroscopic SBF magnitudes. We provide a recipe for obtaining the latter and discuss their uncertainties and limitations. We compare our spectroscopic SBF magnitudes to photometric data of a sample of early-type galaxies. We also show that the SBF spectra can be very useful for constraining relevant stellar population parameters. We find small (<5%) mass-fractions of extremely metal-poor components ([M/H]<-1) on the top of the dominant, old and metal-rich stellar population. These results put stringent constraints on the early stages of galaxy formation in massive elliptical galaxies. This is remarkable given the high degree of degeneracy of the standard spectral analysis to such metal-poor stellar populations in the visible and in the near-IR. The new SBF models show great potential for exploiting ongoing surveys, particularly those based on narrow-band filters., Comment: MNRAS, accepted. 25 pages, 13 figures. Models can be downloaded from http://miles.iac.es/ (OTHER PREDICTIONS/DATA)

Using new long-slit spectroscopy obtained with X-Shooter at ESO-VLT, we study, for the first time, radial gradients of optical and Near-Infrared IMF-sensitive features in a representative sample of galaxies at the very high-mass end of the galaxy population. The sample consists of seven early-type galaxies (ETGs) at $z\sim0.05$, with central velocity dispersion in the range $300<\sigma<350$km/s. Using state-of-art stellar population synthesis models, we fit a number of spectral indices, from different chemical species (including TiO's and Na indices), to constrain the IMF slope (i.e. the fraction of low-mass stars), as a function of galactocentric distance, over a radial range out to $\sim4$kpc. ETGs in our sample show a significant correlation of IMF slope and surface mass density. The bottom-heavy population (i.e. an excess of low-mass stars in the IMF) is confined to central galaxy regions with surface mass density above $\sim 10^{10} M_\odot kpc^{-2}$, or, alternatively, within a characteristic radius of $\sim2$~kpc. Radial distance, in physical units, and surface mass density, are the best correlators to IMF variations, with respect to other dynamical (e.g. velocity dispersion) and stellar population (e.g. metallicity) properties. Our results for the most massive galaxies suggest that there is no single parameter} that fully explains variations in the stellar IMF, but IMF radial profiles at z$\sim$0 rather result from the complex formation and mass accretion history of galaxy inner and outer regions., Comment: 23 pages, 13 figures, accepted for publication in MNRAS

The 31 brightest galaxies (m_B < 15 mag) inside the virial radius of the Fornax cluster were observed from the centres to the outskirts with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. These observations provide detailed high-resolution maps of the line-of-sight kinematics and line strengths of the stars and ionised gas reaching 2-3 Re for 21 early-type galaxies and 1-2 Re for 10 late-type galaxies. The majority of the galaxies are regular rotators, with eight hosting a kinematically distinct core. Only two galaxies are slow rotators. The mean age, total metallicity, and [Mg/Fe] abundance ratio in the bright central region inside 0.5 Re and in the galaxy outskirts are presented. Extended emission-line gas is detected in 13 galaxies, most of them are late-type objects with wide-spread star formation. The measured structural properties are analysed in relation to the galaxies' position in the projected phase space of the cluster. This shows that the Fornax cluster appears to consist of three main groups of galaxies inside the virial radius: the old core; a clump of galaxies, which is aligned with the local large-scale structure and was accreted soon after the formation of the core; and a group of galaxies that fell in more recently., Comment: Accepted for publication in Astronomy and Astrophysics Journal. The high-resolution pdf file, including all maps in the Appendix, is available at the following link: https://www.dropbox.com/s/oukfd56dogf6v9t/F3D_sample_paper.pdf?dl=0

We used deep MUSE observations to perform a stellar-kinematic and population analysis of FCC 153 and FCC 177, two edge-on S0 galaxies in the Fornax cluster. The geometrical definition of the different structural components of these two galaxies allows us to describe the nature of their thick disks. These are both old, relatively metal poor and [Mg/Fe]-enhanced, and their star formation history (SFH) reveals a minor younger component whose chemical properties suggest its later accretion. Moreover, the outer regions of these geometrically defined thick disks show higher values of metallicity and lower values of [Mg/Fe]. These stars probably formed in the thin-disk region and they were dynamically heated to form the flares present in these two galaxies. We propose different formation scenarios for the three populations of these thick disks: in-situ formation, accretion and disk heating. A clear distinction in age is found between the metal poor and [Mg/Fe]-enhanced thick disks (old, $\sim 12-13$ Gyr), and the metal rich and less [Mg/Fe]-enhanced thin disks (young, $\sim 4-5$ Gyr). These two galaxies show signs of relatively recent star formation in their thin disks and nuclear regions. While the thin disks show more continuous SFHs, the nuclei display a rather bursty SFH. These two galaxies are located outside of the densest region of the Fornax cluster where FCC 170 resides. This other edge-on S0 galaxy was studied by \citet{Pinna2019}. We compare and discuss our results with this previous study. The differences between these three galaxies, at different distances from the cluster center, suggest that the environment can have a strong effect on the galaxy evolutionary path., Comment: 24 pages, 18 figures

The stellar initial mass function (IMF) regulates the baryonic cycle within galaxies, and is a key ingredient to translate observations into physical quantities. Although for decades it was assumed to be universal, there is now growing observational evidence showing that the center of massive early-type galaxies host an enhanced population of low-mass stars compared to the expectations from the Milky Way. Moreover, these variations in the IMF have been found to be related to the radial metallicity variations in massive galaxies. We present here a two-dimensional stellar population analysis of the massive lenticular galaxy FCC 167 (NGC 1380) as part of the Fornax3D project. Using a newly developed stellar population fitting scheme, we derive a full two-dimensional IMF map of an early-type galaxy. This two-dimensional analysis allows us go further than a radial analysis, showing how the metallicity changes along a disc-like structure while the IMF follows a distinct, less disky distribution. Thus, our findings indicate that metallicity cannot be the sole driver of the observed radial IMF variations. In addition, a comparison with the orbital decomposition shows suggestive evidence of a coupling between stellar population properties and the internal dynamical structure of FCC 167, where metallicity and IMF maps seem to track the distribution of cold and warm orbits, respectively., Comment: 12 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

Context. Most of our knowledge on the stellar component of galaxies is based on the analysis of distant systems and comes from integrated light data. It is important to test whether the results of the star formation histories (SFH) obtained with standard full-spectrum fitting methods are in agreement with those obtained through colour-magnitude diagram (CMD) fitting (usually considered the most reliable approach). Aims. We compare SFHs recovered from both techniques in Leo~A, a Local Group dwarf galaxy whose majority of stars formed during the last 8 Gyrs. This complements our previous findings in a field in the Large Magellanic Cloud bar, where star formation has been on-going since early epochs though at varying rates. Methods. We have used GTC/OSIRIS in long-slit mode to obtain a high-quality integrated light spectrum by scanning a selected region within Leo~A, for which a CMD reaching the old main mequence turn-off (oMSTO) is available from HST. We compared the SFH obtained from the two datasets, using state-of-art methods of integrated light ({\tt STECKMAP}) and resolved stellar population analysis. In the case of the CMD, we computed the SFH both from a deep CMD (observed with HST/ACS), and from a shallower one (archival data from HST/WFPC2). Results. The agreement between the SFHs recovered from the oMSTO CMD and from full spectrum fitting is remarkable, particularly regarding the time evolution of the star formation rate. The overall extremely low metallicity of Leo~A is recovered up to the last 2 Gyrs, when some discrepancies appear. A relatively high metallicity found for the youngest stars from the integrated data is a recurring feature that might indicate that the current models or synthesis codes should be revised, but that can be significantly mitigated using a more restrictive metallicity range... [Abridged], Comment: 15 pages, 12 figures, accepted for publication in A&A

Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star formation histories (SFH) and mean population parameters, such as ages and metallicities, for a sample of five UDG candidates in the Coma cluster. From the radial velocities, we confirm the Coma membership of these galaxies. We find that their rotation properties, if detected at all, are compatible with dwarf-like galaxies. The SFHs of the UDG are dominated by old (~ 7 Gyr), metal-poor ([M/H] ~ -1.1) and alpha-enhanced ([Mg/Fe] ~ 0.4) populations followed by a smooth or episodic decline which halted ~ 2 Gyr ago, possibly a sign of cluster-induced quenching. We find no obvious correlation between individual SFH shapes and any UDG morphological properties. The recovered stellar properties for UDGs are similar to those found for DDO44, a local UDG analogue resolved into stars. We conclude that the UDGs in our sample are extended dwarfs whose properties are likely the outcome of both internal processes, such as bursty SFHs and/or high-spin haloes, as well as environmental effects within the Coma cluster., Comment: Accepted for publication in MNRAS

In this letter we describe how we use stellar dynamics information to constrain the shape of the stellar IMF in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, $\Upsilon_\mathrm{dyn}$ and $\Upsilon_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $\Upsilon_\mathrm{dyn} \ge \Upsilon_{\ast}$ to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power law IMF shape for 75% of the galaxies in our sample. Conversely, we find that a double power law IMF shape with a varying high-mass end slope is compatible (within 1$\sigma$) with 95% of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies., Comment: 6 pages, 3 figures, accepted for publication in MNRAS Letters

In the local Universe, massive early-type galaxies exhibit enhanced [Mg/Fe] ratios, which has been traditionally interpreted as the result of a rapid ($\tau \lesssim 1$ Gyr) collapse. However, recent claims of a non-universal, steep initial mass function call for a revision of this standard interpretation. In the present work we show how the simultaneous consideration of a high [Mg/Fe] and a steep IMF slope would imply unreasonably short ($\tau \sim 7$ Myr) and intense (SFR $\sim 10^{5}$ Msun yr$^{-1}$) formation events for massive early-type galaxies. We discuss possible caveats and explanations to this apparent inconsistency, and we suggest that further IMF determinations, both in the local Universe and at high redshift, are necessary to better understand the problem., Comment: 5 pages, MNRAS Letters accepted

At present, the main challenge to the interpretation of variations in gravity-sensitive line strengths as driven by a non-universal initial mass function (IMF), lies in understanding the effect of other parameters describing unresolved stellar populations, such as elemental abundance ratios. We combine various TiO-based, IMF-sensitive indicators in the optical and NIR spectral windows, along with the FeH-based Wing-Ford band to break this degeneracy. We obtain a significant radial trend of the IMF slope in XSG1, a massive early-type galaxy (ETG), with velocity dispersion sigma~300km/s, observed with the VLT/X-SHOOTER instrument. In addition, we constrain both the shape and normalization of the IMF based only on a stellar population analysis. We robustly rule out a single power-law to describe the IMF, whereas a power law tapered off to a constant value at low masses (defined as a bimodal IMF) is consistent with all the observational spectroscopic data and with the stellar M/L constraints based on the Jeans Anisotropic Modelling method. The IMF in XSG1 is bottom-heavy in the central regions (corresponding to a bimodal IMF slope Gb~3, or a mass normalization mismatch parameter alpha~2), changing towards a standard Milky-Way like IMF (Gb~1.3; alpha~1) around half of the effective radius. This result, combined with previous observations of local IMF variations in massive ETGs, reflects the varying processes underlying the formation of this type of galaxies, between the central core and the outer regions., Comment: 24pages, 16 figures, accepted for publication in MNRAS

Several authors have reported that the dynamical masses of massive compact galaxies ($M_\star \gtrsim 10^{11} \ \mathrm{M_\odot}$, $r_\mathrm{e} \sim 1 \ \mathrm{kpc}$), computed as $M_\mathrm{dyn} = 5.0 \ \sigma_\mathrm{e}^2 r_\mathrm{e} / G$, are lower than their stellar masses $M_\star$. In a previous study from our group, the discrepancy is interpreted as a breakdown of the assumption of homology that underlie the $M_\mathrm{dyn}$ determinations. Here, we present new spectroscopy of six redshift $z \approx 1.0$ massive compact ellipticals from the Extended Groth Strip, obtained with the 10.4 m Gran Telescopio Canarias. We obtain velocity dispersions in the range $161-340 \ \mathrm{km \ s^{-1}}$. As found by previous studies of massive compact galaxies, our velocity dispersions are lower than the virial expectation, and all of our galaxies show $M_\mathrm{dyn} < M_\star$ (assuming a Salpeter initial mass function). Adding data from the literature, we build a sample covering a range of stellar masses and compactness in a narrow redshift range $\mathit{z \approx 1.0}$. This allows us to exclude systematic effects on the data and evolutionary effects on the galaxy population, which could have affected previous studies. We confirm that mass discrepancy scales with galaxy compactness. We use the stellar mass plane ($M_\star$, $\sigma_\mathrm{e}$, $r_\mathrm{e}$) populated by our sample to constrain a generic evolution mechanism. We find that the simulations of the growth of massive ellipticals due to mergers agree with our constraints and discard the assumption of homology., Comment: 17 pages, 13 figures, 3 tables. Updated to match the published version in MNRAS

We explore the stellar initial mass function (IMF) of a sample of 49 massive quiescent galaxies (MQGs) at 0.9$<$z$<$1.5. We base our analysis on intermediate resolution spectro-photometric data in the GOODS-N field taken in the near-infrared and optical with the HST/WFC3 G141 grism and the Survey for High-z Absorption Red and Dead Sources (SHARDS). To constrain the slope of the IMF, we have measured the TiO$_2$ spectral feature, whose strength depends strongly on the content of low-mass stars, as well as on stellar age. Using ultraviolet to near-infrared individual and stacked spectral energy distributions, we have independently estimated the stellar ages of our galaxies. Knowing the age of the stellar population, we interpret the strong differences in the TiO$_2$ feature as an IMF variation. In particular, for the heaviest z$\sim$1 MQGs (M$>$10$^{11}$Msun) we find an average age of 1.7$\pm$0.3 Gyr and a bottom-heavy IMF ($\Gamma_b$=3.2$\pm$0.2). Lighter MQGs (2$\times$10$^{10}$$<$M$<$10$^{11}$ Msun) at the same redshift are younger on average (1.0$\pm$0.2 Gyr) and present a shallower IMF slope ($\Gamma_b=2.7^{+0.3}_{-0.4}$). Our results are in good agreement with the findings about the IMF slope in early-type galaxies of similar mass in the present-day Universe. This suggests that the IMF, a key characteristic of the stellar populations in galaxies, is bottom-heavier for more massive galaxies and has remained unchanged in the last $\sim$8 Gyr., Comment: 6 pages, 5 figures, accepted for publication in ApJL

The absence of stellar disc truncations in low-inclined spiral galaxies has been a matter of debate in the last decade. Disc truncations are often observed in highly inclined galaxies but no obvious detection of this feature has so far been made in face-on spirals. Here we show, using a simple exponential disc plus stellar halo model based on current observational constraints, that truncations in face-on projections occur at surface brightness levels comparable to the brightness of stellar haloes at the same radial distance. In this sense, stellar haloes outshine the galaxy disc at the expected position of the truncations, forcing their studies only in highly inclined (edge-on) orientations., Comment: 6 pages, 2 figures. Submitted to MNRAS. Comments are welcome