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.

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.

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.

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.

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.

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.

Nuclear star clusters (NSCs) are found in at least 70% of all galaxies, but their formation path is still unclear. In the most common scenarios, NSCs form in-situ from the galaxy's central gas reservoir, through merging of globular clusters (GCs), or through a combination of the two. As the scenarios pose different expectations for angular momentum and stellar population properties of the NSC in comparison to the host galaxy and the GC system, it is necessary to characterise the stellar light, NSC and GCs simultaneously. The large NSC (r$_\rm{eff} = 66$ pc) and rich GC system of the early-type Fornax cluster galaxy FCC47 (NGC1336) render this galaxy an ideal laboratory to constrain NSC formation. Using MUSE science verification data assisted by adaptive optics, we obtained maps for the stellar kinematics and for stellar-population properties of FCC47. We extracted the spectra of the central NSC and determined line-of-sight velocities of 24 GCs and metallicities of five. FCC47 shows two decoupled components (KDCs): a rotating disk and the NSC. Our orbit-based dynamical Schwarzschild model revealed that the NSC is a distinct kinematic feature and it constitutes the peak of metallicity and old ages in the galaxy. The main body consists of two counter-rotating populations and is dominated by a more metal-poor population. The GC system is bimodal with a dominant metal-poor population and the total GC system mass is $\sim 17\%$ of the NSC mass ($\sim$ 7 $\times$ $10^8 M_\odot$). The rotation, high metallicity and high mass of the NSC cannot be uniquely explained by GC-inspiral and most likely requires additional, but quickly quenched, in-situ formation. The presence of two KDCs most probably are evidence of a major merger that has altered the structure of FCC47 significantly, indicating the important role of galaxy mergers in forming the complex kinematics in the galaxy-NSC system.

The formation and assembly process of massive galaxies is a combination of two phases: an initial in-situ-dominated one followed by an ex-situ-dominated evolution. Separating these two contributions is therefore crucial to understand the baryonic cycle within massive halos. A recently discovered population of so-called naked red nuggets, galaxies that shortcut the ex-situ stage preserving their pristine properties, presents a unique opportunity to study in detail star formation in massive galaxies without the confounding effect of later accretion. We investigate the spatially resolved star formation histories of a sample of 12 naked red nuggets. We measure how their radial light distributions, star formation rates and central densities evolved in time. We find that, while forming stars, red nuggets become gradually more concentrated, reaching a maximum concentration at quenching. After being quenched, they kept forming stars in a more disky-like configuration. Our measurements suggest that super-massive black holes and host galaxies grow their mass in a self-regulated way until a characteristic Mbh/Mhalo is reached. Once black holes are massive enough, red nuggets get quenched and depart from the star formation main sequence. While in the main sequence, red nuggets evolve at roughly constant star formation rate. This can explain up to ~0.3 dex of the scatter of the star formation main sequence, as well as its higher normalization observed in the early Universe. Hence, our results suggest that the main sequence is composed of populations of galaxies at different evolutionary stages, and that the scatter is therefore due to secular processes.

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.

The $M_{\bullet} - \sigma$ relation establishes a connection between central black holes (BHs) and their host spheroids. Supported by observations at $M_{\bullet} \gtrsim 10^5 \, M_{\odot}$, there is limited data on its validity at lower masses. Employing a semi-analytical model to simulate the combined evolution of BHs and their host galaxies, we predict the observational consequences of assuming a bimodality in the accretion efficiency of BHs, with low-mass BHs ($M_{\bullet} \lesssim 10^5 \, M_{\odot}$) accreting inefficiently. We predict a departure from the $M_{\bullet} - \sigma$ relation at a transitional BH mass $\sim 10^5 \, M_{\odot}$, with lower-mass BHs unable to reach the mass dictated by the relation and becoming disconnected from the evolution of the host galaxy. This prediction is an alternative to previous works suggesting a flattening of the relation at $\sim 10^5-10^6 \, M_{\odot}$. Furthermore, we predict a deficit of BHs shining at bolometric luminosities $\sim 10^{42} \, \mathrm{erg \, s^{-1}}$. Joined with a detection bias, this could partly explain the scarce number of intermediate-mass BHs detected. Conversely, we predict an increase in source density at lower bolometric luminosities, $<10^{42} \, \mathrm{erg \, s^{-1}}$. Because our predictions assume a bimodal population of high-redshift BH seeds, future observations of fainter BHs will be fundamental for constraining the nature of these seeds.

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]

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.

We use Keck/DEIMOS spectroscopy to confirm the cluster membership of 16 ultra-diffuse galaxies (UDGs) in the Coma cluster, bringing the total number of spectroscopically con- firmed UDGs to 24. We also identify a new cluster background UDG. In this pilot study of Coma UDGs in velocity phase-space, we find evidence that most present-day Coma UDGs have a recent infall epoch while a few may be ancient infalls. These recent infall UDGs have higher absolute relative line-of-sight velocities, bluer optical colors, and are smaller in size, unlike the ancient infalls. The kinematics of the spectroscopically confirmed Coma UDG sample is similar to that of the cluster late-type galaxy population. Our velocity phase-space analysis suggests that present-day cluster UDGs have a predominantly accretion origin from the field, acquire velocities corresponding to the mass of the cluster at accretion as they are accelerated towards the cluster center, and become redder and bigger as they experience the various physical processes at work within the cluster.

Dedicated photometric and spectroscopic surveys have provided unambiguous evidence for a strong stellar mass-size evolution of galaxies within the last 10 Gyr. The likely progenitors of today's most massive galaxies are remarkably small, disky, passive and have already assembled much of their stellar mass at redshift z=2. An in-depth analysis of these objects, however, is currently not feasible due to the lack of high-quality, spatially-resolved photometric and spectroscopic data. In this paper, we present a sample of nearby compact elliptical galaxies (CEGs), which bear resemblance to the massive and quiescent galaxy population at earlier times. Hubble Space Telescope (HST) and wide-field integral field unit (IFU) data have been obtained, and are used to constrain orbit-based dynamical models and stellar population synthesis (SPS) fits, to unravel their structural and dynamical properties. We first show that our galaxies are outliers in the present-day stellar mass-size relation. They are, however, consistent with the mass-size relation of compact, massive and quiescent galaxies at redshift z=2. The compact sizes of our nearby galaxies imply high central stellar mass surface densities, which are also in agreement with the massive galaxy population at higher redshift, hinting at strong dissipational processes during their formation. Corroborating evidence for a largely passive evolution within the last 10 Gyr is provided by their orbital distribution as well as their stellar populations, which are difficult to reconcile with a very active (major) merging history. This all supports that we can use nearby CEGs as local analogues of the high-redshift, massive and quiescent galaxy population, thus providing additional constraints for models of galaxy formation and evolution.

We confirm two new local massive relic galaxies, i.e. untouched survivors of the early universe massive population: Mrk1216 and PGC032873. Both show early and peaked formation events within very short timescales (<1 Gyr) and thus old mean mass-weighted ages (~13 Gyr). Their star formation histories remain virtually unchanged out to several effective radii, even when considering the steeper IMF values inferred out to ~3 effective radii. Their morphologies, kinematics and density profiles are like those found in the z>2 massive population, setting them apart of the typical z~0 massive early-type galaxies. We find that there seems to exist a "degree of relic" that is related on how far into the path to become one of these typical z~0 massive galaxies the compact relic has undergone. This path is partly dictated by the environment the galaxy lives in. For galaxies in rich environments, such as the previously reported relic galaxy NGC1277, the most extreme properties (e.g. sizes, short formation timescales, larger super-massive black holes) are expected, while lower density environments will have galaxies with delayed and/or extended star formations, slightly larger sizes and not that extreme black hole masses. The confirmation of 3 relic galaxies up to a distance of 106Mpc implies a lower limit in the number density of these red nuggets in the local universe of 6x10^{-7} Mpc^{3}, which is within the theoretical expectations.

Coevolution between supermassive black holes (BHs) and their host galaxies is universally adopted in models for galaxy formation. In the absence of feedback from active galactic nuclei, simulated massive galaxies keep forming stars in the local Universe. From an observational point of view, however, such coevolution remains unclear. We present a stellar population analysis of galaxies with direct BH mass measurements and the BH mass-{\sigma} relation as a working framework. We find that over-massive BH galaxies, i.e., galaxies lying above the best-fitting BH mass-{\sigma} line, tend to be older and more {\alpha}-element enhanced than under-massive BH galaxies. The scatter in the BH mass-{\sigma}-[{\alpha}/Fe] plane is significantly lower than in the standard BH mass-{\sigma} relation. We interpret this trend as an imprint of active galactic nucleus feedback on the star formation histories of massive galaxies.

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.

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.

We investigate the compact, early-type galaxy NGC 1281 with integral field unit observations to map the stellar LOSVD out to 5 effective radii and construct orbit-based dynamical models to constrain its dark and luminous matter content. Under the assumption of mass-follows-light, the H-band stellar mass-to-light ratio (M/L) is {\Upsilon} = 2.7(+-0.1) {\Upsilon}_{sun}, higher than expected from our stellar population synthesis fits with either a canonical Kroupa ({\Upsilon} = 1.3 {\Upsilon}_{sun}) or Salpeter ({\Upsilon} = 1.7 {\Upsilon}_{sun}) stellar initial mass function. Such models also cannot reproduce the details of the LOSVD. Models with a dark halo recover the kinematics well and indicate that NGC 1281 is dark matter dominated, making up ~ 90 per cent of the total enclosed mass within the kinematic bounds. Parameterised as a spherical NFW profile, the dark halo mass is 11.5 < log(M_{DM}/M_{sun}) < 11.8 and the stellar M/L is 0.6 < {\Upsilon} < 1.1. However, this stellar M/L is lower than predicted by its old stellar population. Moreover, the halo mass within the kinematic extent is ten times larger than expected based on {\Lambda}CDM predictions, and an extrapolation yields cluster sized dark halo masses. Adopting {\Upsilon} = 1.7 {\Upsilon}_{sun} yields more moderate dark halo virial masses, but these models fit the kinematics worse. A non-NFW model might solve the discrepancy between the unphysical consequences of the best-fitting dynamical models and models based on more reasonable assumptions for the dark halo and stellar mass-to-light ratio, which are disfavoured according to our parameter estimation.

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.

Variations in the stellar initial mass function (IMF) have been invoked to explain the spectroscopic and dynamical properties of early-type galaxies. However, no observations have yet been able to disentangle the physical driver. We analyse here a sample of 24 early-type galaxies drawn from the CALIFA survey, deriving in a homogeneous way their stellar population and kinematic properties. We find that the local IMF is tightly related to the local metallicity, becoming more bottom-heavy towards metal-rich populations. Our result, combined with the galaxy mass-metallicity relation, naturally explains previous claims of a galaxy mass-IMF relation, derived from non-IFU spectra. If we assume that - within the star formation environment of early-type galaxies - metallicity is the main driver of IMF variations, a significant revision of the interpretation of galaxy evolution observables is necessary.

Massive relic galaxies formed the bulk of their stellar component before z~2 and have remained unaltered since then. Therefore, they represent a unique opportunity to study in great detail the frozen stellar population properties of those galaxies that populated the primitive Universe. We have combined optical to near-infrared line-strength indices in order to infer, out to 1.5 Reff, the IMF of the nearby relic massive galaxy NGC 1277. The IMF of this galaxy is bottom-heavy at all radii, with the fraction of low-mass stars being at least a factor of two larger than that found in the Milky Way. The excess of low-mass stars is present throughout the galaxy, while the velocity dispersion profile shows a strong decrease with radius. This behaviour suggests that local velocity dispersion is not the only driver of the observed IMF variations seen among nearby early-type galaxies. In addition, the excess of low-mass stars shown in NGC 1277 could reflect the effect on the IMF of dramatically different and intense star formation processes at z~2, compared to the less extreme conditions observed in the local Universe.

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.

The hypothesis of a universal initial mass function (IMF) -- motivated by observations in nearby stellar systems -- has been recently challenged by the discovery of a systematic variation of the IMF with the central velocity dispersion, {\sigma}, of early-type galaxies (ETGs), towards an excess of low-mass stars in high-{\sigma} galaxies. This trend has been derived so far from integrated spectra, and remains unexplained at present. To test whether such trend depends on the local properties within a galaxy, we have obtained new, extremely deep, spectroscopic data, for three nearby ETGs, two galaxies with high {\sigma} (~300 km/s), and one lower mass system, with {\sigma} ~ 100 km/s. From the analysis of IMF-sensitive spectral features, we find that the IMF depends significantly on galactocentric distance in the massive ETGs, with the enhanced fraction of low-mass stars f mostly confined to their central regions. In contrast, the low-{\sigma} galaxy does not show any significant radial gradient in the IMF, well described by a shallower distribution, relative to the innermost regions of massive galaxies, at all radii. Such a result indicates that the IMF should be regarded as a local (rather than global) property, and suggests a significant difference between the formation process of the core and the outer regions of massive ETGs.

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.

The mechanism causing breaks in the radial surface brightness distribution of spiral galaxies is not yet well known. Despite theoretical efforts, there is not a unique explanation for these features and the observational results are not conclusive. In an attempt to address this problem, we have selected a sample of 34 highly inclined spiral galaxies present both in the Sloan Digital Sky Survey and in the Spitzer Survey of Stellar Structure in Galaxies. We have measured the surface brightness profiles in the five Sloan optical bands and in the 3.6$\mu m$ Spitzer band. We have also calculated the color and stellar surface mass density profiles using the available photometric information, finding two differentiated features: an innermost break radius at distances of $\sim 8 \pm 1$ kpc [$0.77 \pm 0.06$ $R_{25}$] and a second characteristic radius, or truncation radius, close to the outermost optical extent ($\sim 14 \pm 2$ kpc [$1.09 \pm 0.05$ $R_{25}$]) of the galaxy. We propose in this work that the breaks might be a phenomena related to a threshold in the star formation, while truncations are more likely a real drop in the stellar mass density of the disk associated with the maximum angular momentum of the stars.