It is generally recognized that massive galaxies form through a combination of in-situ collapse and ex-situ accretion. The in-situ component forms early, where gas collapse and compaction leads to the formation of massive compact systems (blue and red "nuggets") seen at $z>1$. The subsequent accretion of satellites brings in ex-situ material, growing these nuggets in size and mass to appear as the massive early-type galaxies (ETGs) we see locally. Due to stochasticity in the accretion process, in a few rare cases a red nugget will evolve to the present day having undergone little ex-situ mass accretion. The resulting massive, compact and ancient objects have been termed "relic galaxies". Detailed stellar population and kinematic analyses are required to characterise these systems. However, an additional crucial aspect lies in determining the fraction of ex-situ mass they have accreted since their formation. Globular cluster systems can be used to constrain this fraction, since the oldest and most metal-poor globular clusters in massive galaxies are primarily an accreted, ex-situ population. Models for the formation of relic galaxies and their globular cluster systems suggest that, due to their early compaction and limited accretion of dark-matter dominated satellites, relic galaxies should have characteristically low dark-matter mass fractions compared to ETGs of the same stellar mass.

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.

Observationally, constraining the baryonic cycle within massive galaxies has proven to be quite difficult. In particular, the role of black hole feedback in regulating star formation, a key process in our theoretical understanding of galaxy formation, remains highly debated. We present here observational evidence showing that, at fixed stellar velocity dispersion, the temperature of the hot gas is higher for those galaxies hosting more massive black holes in their centers. Analyzed in the context of well-established scaling relations, particularly the mass-size plane, the relation between the mass of the black hole and the temperature of the hot gas around massive galaxies provides further observational support to the idea that baryonic processes within massive galaxies are regulated by the combined effects of the galaxy halo virial temperature and black hole feedback, in agreement with the expectations from the EAGLE cosmological numerical simulation.

Super-massive black holes are fundamental ingredients in our theoretical understanding of galaxy formation. They are likely the only sources energetic enough to regulate star formation within massive dark matter halos, but observational evidence of this process remains elusive. The effect of black hole feedback is expected to be a strong function of halo mass, and galaxy groups and clusters are among the most massive structures in the Universe. At fixed halo mass, we find an enhanced fraction of quiescent satellite galaxies and a hotter X-ray intragroup and intracluster medium in those groups and clusters hosting more massive black holes in their centers. These results indicate that black hole feedback makes quenching processes more efficient through a cumulative heating of the gaseous intragroup and intracluster medium.

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.

Ultra diffuse galaxies are a population of extended galaxies but with relatively low luminosities. The origin of these objects remains unclear, largely due to the observational challenges of the low surface brightness Universe. We present here a detailed stellar population analysis of a relatively isolated UDG, DGSAT I, based on spectroscopic data from the Keck Cosmic Web Imager integral field unit. The star formation history of DGSAT I seems to be extended, with a mean luminosity-weighted age of ~3 Gyr, in agreement with previous photometric studies. However, we find a very high [Mg/Fe] abundance ratio, which is extreme even in the context of the highly alpha-enhanced massive ellipticals and ultra-faint dwarfs. The [Mg/Fe]-enhancement of DGSAT I appears to be 10 times higher than the most magnesium-enhanced stellar systems discovered to date, and suggests that the chemical enrichment of this object was dominated by core-collapse supernovae. Intriguingly, this breaks the canonical relation between [Mg/Fe] and star formation time-scale. With a measured velocity dispersion of 56 +/- 10 km/s, DGSAT I also shows a high mass-to-light ratio, which indicates that it is highly dark matter-dominated. The metal-poor conditions of DGSAT I may have enhanced the formation of massive stars, while at the same time, additional mechanisms are needed to prevent iron-rich yields from being recycled into stars. These results suggest that some ultra-diffuse galaxies could have experienced chemical enrichment episodes similar to the first building blocks of galaxies.

NGC7252, which is one of the nearest major-merger galaxy remnants, is an ideal laboratory to study the processes inherent to the transformation of disc galaxies to ellipticals as observed about ~1Gyr after the collision. We obtained wide-field IFU spectroscopy with the VLT-VIMOS integral-field spectrograph covering the central 50"x50" of NGC7252 to map the stellar and ionised gas kinematics, and the distribution and conditions of the ionised gas, revealing the extent of ongoing star formation and recent star formation history. We find that the inner gas disc is not counter-rotating with respect to the stars and that the stellar kinematics appear complex with a clear indication of a prolate-like rotation component suggesting a polar merger configuration. The ongoing star formation rate is 2.2+-0.6 M_sun/yr and implies a typical depletion time of 2Gyr given the molecular gas content. Furthermore, the spatially-resolved star formation history suggests a slight radial dependence, moving outwards at later times. We confirm a large AGN-ionised gas cloud previously discovered 5kpc south of the nucleus, and find higher ionisation also at the galaxy centre relative to the surrounding gas disc. Although the higher ionisation towards the centre is potentially degenerate within the central star forming ring, it may be associated with a low-luminosity AGN. Although NGC7252 has been classified as post-starburst galaxy at the centre, the elliptical-like major-merger remnant still appears very active. A central kpc-scale gas disc has presumably re-formed quickly within the last 100Myr after final coalescence. The disc features ongoing star formation, implying Gyr long timescale to reach the red sequence through gas consumption alone. While NGC7252 is useful to probe the transformation from discs to ellipticals, it is not well-suited to study the transformation from blue to red at this point.

In this second paper of the series, we study with new Keck/DEIMOS spectra the stellar populations of seven spectroscopically confirmed ultra-diffuse galaxies (UDGs) in the Coma cluster. We find intermediate to old ages (~ 7Gyr), low metallicities ([Z/H] ~ -0.7dex) and mostly super-solar abundance patterns ([Mg/Fe] ~ 0.13dex). These properties are similar to those of low-luminosity (dwarf) galaxies inhabiting the same area in the cluster and are most consistent with being the continuity of the stellar mass scaling relations of more massive galaxies. These UDGs' star formation histories imply a relatively recent infall into the Coma cluster, consistent with the theoretical predictions for a dwarf-like origin. However, considering the scatter in the resulting properties and including other UDGs in Coma, together with the results from the velocity phase-space study of the Paper I in this series, a mixed-bag of origins is needed to explain the nature of all UDGs. Our results thus reinforce a scenario in which many UDGs are field dwarfs that become quenched through their later infall onto cluster environments, whereas some UDGs could be genuine primordial galaxies that failed to develop due to an early quenching phase. The unknown proportion of dwarf-like to primordial-like UDGs leaves the enigma of the nature of UDGs still open.

We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the prospector fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). We find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. In the future, our approach can be applied to a large and diverse sample of UDGs down to faint surface brightness limits, with the goal of constraining their stellar ages, stellar metallicities, and circumstellar and diffuse interstellar dust content.

Energy feedback, either from active galactic nuclei (AGN) or from supernovae, is required to understand galaxy formation within a $\Lambda$-Cold Dark Matter cosmology. We study a sample of 127 low-mass galaxies, comparing their stellar populations properties to the mass of the central supermassive black hole, in order to investigate the effect of AGN feedback. We find a loose coupling between star formation history and black hole mass, which seems to suggest that AGN activity does not dominate baryonic cooling in low-mass galaxies. We also find that a break in the $M_\bullet$-$\sigma$ relation marks a transitional stellar mass, M$_\mathrm{trans}=3.4\pm2.1 \times 10^{10}$ M$_{\odot}$, remarkably similar to M$_\star$. Our results are in agreement with a bi-modal star formation process where the AGN-dominated feedback of high-mass galaxies transitions towards a supernovae-driven regime in low-mass systems, as suggested by numerical simulations.

To investigate star formation and assembly processes of massive galaxies, we present here a spatially-resolved stellar populations analysis of a sample of 45 elliptical galaxies (Es) selected from the CALIFA survey. We find rather flat age and [Mg/Fe] radial gradients, weakly dependent on the effective velocity dispersion of the galaxy within half-light radius. However, our analysis shows that metallicity gradients become steeper with increasing galaxy velocity dispersion. In addition, we have homogeneously compared the stellar populations gradients of our sample of Es to a sample of nearby relic galaxies, i.e., local remnants of the high-z population of red nuggets. This comparison indicates that, first, the cores of present-day massive galaxies were likely formed in gas-rich, rapid star formation events at high redshift (z>2). This led to radial metallicity variations steeper than observed in the local Universe, and positive [Mg/Fe] gradients. Second, our analysis also suggests that a later sequence of minor dry mergers, populating the outskirts of early-type galaxies (ETGs), flattened the pristine [Mg/Fe] and metallicity gradients. Finally, we find a tight age-[Mg/Fe] relation, supporting that the duration of the star formation is the main driver of the [Mg/Fe] enhancement in massive ETGs. However, the star formation time-scale alone is not able to fully explain our [Mg/Fe] measurements. Interestingly, our results match the expected effect that a variable stellar initial mass function would have on the [Mg/Fe] ratio.

Super-massive black holes, with masses larger than a million times that of the Sun, appear to inhabit the centers of all massive galaxies. Cosmologically-motivated theories of galaxy formation need feedback from these super-massive black holes to regulate star formation. In the absence of such feedback, state-of-the-art numerical simulations dramatically fail to reproduce the number density and properties of massive galaxies in the local Universe. However, there is no observational evidence of this strongly coupled co-evolution between super-massive black holes and star formation, impeding our understanding of baryonic processes within galaxies. Here we show that the star formation histories (SFHs) of nearby massive galaxies, as measured from their integrated optical spectra, depend on the mass of the central super-massive black hole. Our results suggest that black hole mass growth scales with gas cooling rate in the early Universe. The subsequent quenching of star formation takes place earlier and more efficiently in galaxies hosting more massive central black holes. The observed relation between black hole mass and star formation efficiency applies to all generations of stars formed throughout a galaxy's life, revealing a continuous interplay between black hole activity and baryon cooling.

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.