The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0\sigma$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

We perform a comprehensive study of the stellar population properties of quiescent galaxies as a function of size and stellar mass to constrain the physical mechanism governing the stellar mass assembly and the likely evolutive scenarios that explain their growth in size. After selecting all the quiescent galaxies from the ALHAMBRA survey by the dust-corrected stellar mass$-$colour diagram, we built a shared sample of $\sim850$ quiescent galaxies with reliable sizes from the HST. The stellar population properties were retrieved using the SED-fitting code MUFFIT with various sets of composite stellar population models. Age, formation epoch, metallicity, and extinction were studied on the stellar mass$-$size plane as function of size through a Monte Carlo approach. This accounted for uncertainties and degeneracy effects amongst stellar population properties. The stellar population properties of quiescent galaxies and their stellar mass and size since $z\sim1$ are correlated. At fixed stellar mass, the more compact the quiescent galaxy, the older and richer in metals it is ($1$Gyr and $0.1$dex, respectively). In addition, more compact galaxies may present slight lower extinctions than their more extended counterparts at the same stellar mass ($<0.1$ mag). By means of studying constant regions of stellar population properties across the stellar mass$-$size plane, we obtained empirical relations to constrain the physical mechanism that governs the stellar mass assembly of the form $M_\star \propto r_\mathrm{c}^\alpha$, where $\alpha$ amounts to $0.50-0.55 \pm 0.09$. There are indications that support the idea that the velocity dispersion is tightly correlated with the stellar content of galaxies. The mechanisms driving the evolution of stellar populations can therefore be partly linked to the dynamical properties of galaxies, along with their gravitational potential.

The upcoming new generation of optical spectrographs on four-meter-class telescopes will provide invaluable information for reconstructing the history of star formation in individual galaxies up to redshifts of about 0.7. We aim at defining simple but robust and meaningful physical parameters that can be used to trace the coexistence of widely diverse stellar components: younger stellar populations superimposed on the bulk of older ones. We produce spectra of galaxies closely mimicking data from the forthcoming Stellar Populations at intermediate redshifts Survey (StePS), a survey that uses the WEAVE spectrograph on the William Herschel Telescope. First, we assess our ability to reliably measure both ultraviolet and optical spectral indices in galaxies of different spectral types for typically expected signal-to-noise levels. Then, we analyze such mock spectra with a Bayesian approach, deriving the probability density function of r- and u-band light-weighted ages as well as of their difference. We find that the ultraviolet indices significantly narrow the uncertainties in estimating the r- and u-band light-weighted ages and their difference in individual galaxies. These diagnostics, robustly retrievable for large galaxy samples even when observed at moderate signal-to-noise ratios, allow us to identify secondary episodes of star formation up to an age of ~0.1 Gyr for stellar populations older than ~1.5 Gyr, pushing up to an age of ~1 Gyr for stellar populations older than ~5 Gyr. The difference between r-band and u-band light-weighted ages is shown to be a powerful diagnostic to characterize and constrain extended star-formation histories and the presence of young stellar populations on top of older ones. This parameter can be used to explore the interplay between different galaxy star-formation histories and physical parameters such as galaxy mass, size, morphology, and environment.

We present the results of a 3-year long, medium-resolution spectroscopic campaign aimed at identifying very metal-poor stars from candidates selected with the CaHK, metallicity-sensitive Pristine survey. The catalogue consists of a total of 1007 stars, and includes 146 rediscoveries of metal-poor stars already presented in previous surveys, 707 new very metal-poor stars with [Fe/H] < -2.0, and 95 new extremely metal-poor stars with [Fe/H] < -3.0. We provide a spectroscopic [Fe/H] for every star in the catalogue, and [C/Fe] measurements for a subset of the stars (10% with [Fe/H] < -3 and 24% with -3 < [Fe/H] < -2) for which a carbon determination is possible, contingent mainly on the carbon abundance, effective temperature and S/N of the stellar spectra. We find an average carbon enhancement fraction ([C/Fe] >= +0.7) of 41 +- 4% for stars with -3 < [Fe/H] < -2 and 58 +- 14% for stars with [Fe/H] < -3, and report updated success rates for the Pristine survey of 56 % and 23 % to recover stars with [Fe/H] < -2.5 and [Fe/H] < -3, respectively. Finally, we discuss the current status of the survey and its preparation for providing targets to upcoming multi-object spectroscopic surveys such as WEAVE.

Context. There are more than 3000 known Galactic planetary nebulae (PNe), but only 492 central stars of Galactic planetary nebulae (CSPN) have known spectral types. It is vital to increase this number in order to have reliable statistics, which will lead to an increase of our understanding of these amazing objects. Aims. We aim to contribute to the knowledge of central stars of planetary nebulae and stellar evolution. Methods. This observational study is based on Gemini Multi-Object Spectrographs (GMOS) and with the Intermediate Dispersion Spectrograph (IDS) at the Isaac Newton Telescope (INT) spectra of 78 CSPN. The objects were selected because they did not have any previous classification, or the present classification is ambiguous. These new high quality spectra allowed us to identify the key stellar lines for determining spectral classification in the Morgan-Keenan (MK) system. Results. We have acquired optical spectra of a large sample of CSPN. From the observed targets, 50 are classified here for the first time while for 28 the existing classifications have been improved. In seven objects we have identified a P-Cygni profile at the He i lines. Six of these CSPN are late O-type. The vast majority of the stars in the sample exhibit an absorption-type spectrum, and in one case we have found wide emission lines typical of [WR] stars. We give a complementary, and preliminary, classification criterion to obtain the sub-type of the O(H)-type CSPN. Finally, we give a more realistic value of the proportion of CSPN that are rich or poor in hydrogen.

The finding that massive galaxies grow with cosmic time fired the starting gun for the search of objects which could have survived up to the present day without suffering substantial changes (neither in their structures, neither in their stellar populations). Nevertheless, and despite the community efforts, up to now only one firm candidate to be considered one of these relics is known: NGC 1277. Curiously, this galaxy is located at the centre of one of the most rich near galaxy clusters: Perseus. Is its location a matter of chance? Should relic hunters focus their search on galaxy clusters? In order to reply this question, we have performed a simultaneous and analogous analysis using simulations (Millennium I-WMAP7) and observations (New York University Value-Added Galaxy Catalogue). Our results in both frameworks agree: it is more probable to find relics in high density environments.

We study the preferred environments of $z \sim 0$ massive relic galaxies ($M_\star \gtrsim 10^{10}~\mathrm{M_\odot}$ galaxies with little or no growth from star formation or mergers since $z \sim 2$). Significantly, we carry out our analysis on both a large cosmological simulation and an observed galaxy catalogue. Working on the Millennium I-WMAP7 simulation we show that the fraction of today massive objects which have grown less than 10 per cent in mass since $z \sim 2$ is ~0.04 per cent for the whole massive galaxy population with $M_\star > 10^{10}~\mathrm{M_\odot}$. This fraction rises to ~0.18 per cent in galaxy clusters, confirming that clusters help massive galaxies remain unaltered. Simulations also show that massive relic galaxies tend to be closer to cluster centres than other massive galaxies. Using the New York University Value-Added Galaxy Catalogue, and defining relics as $M_\star \gtrsim 10^{10}~\mathrm{M_\odot}$ early-type galaxies with colours compatible with single-stellar population ages older than 10 Gyr, and which occupy the bottom 5-percentile in the stellar mass-size distribution, we find $1.11 \pm 0.05$ per cent of relics among massive galaxies. This fraction rises to $2.4 \pm 0.4$ per cent in high-density environments. Our findings point in the same direction as the works by Poggianti et al. and Stringer et al. Our results may reflect the fact that the cores of the clusters are created very early on, hence the centres host the first cluster members. Near the centres, high-velocity dispersions and harassment help cluster core members avoid the growth of an accreted stellar envelope via mergers, while a hot intracluster medium prevents cold gas from reaching the galaxies, inhibiting star formation.

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.

We address the problem that dynamical masses of high-redshift massive galaxies, derived using virial scaling, often come out lower than stellar masses inferred from population fitting to multi-band photometry. We compare dynamical and stellar masses for various samples spanning ranges of mass, compactness and redshift, including the SDSS. The discrepancy between dynamical and stellar masses occurs both at low and high redshifts, and systematically increases with galaxy compactness. Because it is unlikely that stellar masses show systematic errors with galaxy compactness, the correlation of mass discrepancy with compactness points to errors in the dynamical mass estimates which assume homology with massive, nearby ellipticals. We quantify the deviations from homology and propose specific non-virial scaling of dynamical mass with effective radius and velocity dispersion.

For many massive compact galaxies, their dynamical masses ($M_\mathrm{dyn} \propto \sigma^2 r_\mathrm{e}$) are lower than their stellar masses ($M_\star$). We analyse the unphysical mass discrepancy $M_\star / M_\mathrm{dyn} > 1$ on a stellar-mass-selected sample of early-type galaxies ($M_\star \gtrsim 10^{11} \ \mathrm{M_\odot}$) at redshifts $z \sim 0.2$ to $z \sim 1.1$. We build stacked spectra for bins of redshift, size and stellar mass, obtain velocity dispersions, and infer dynamical masses using the virial relation $M_\mathrm{dyn} \equiv K \ \sigma_\mathrm{e}^2 r_\mathrm{e} / G$ with $K = 5.0$; this assumes homology between our galaxies and nearby massive ellipticals. Our sample is completed using literature data, including individual objects up to $z \sim 2.5$ and a large local reference sample from the Sloan Digital Sky Survey (SDSS). We find that, at all redshifts, the discrepancy between $M_\star$ and $M_\mathrm{dyn}$ grows as galaxies depart from the present-day relation between stellar mass and size: the more compact a galaxy, the larger its $M_\star / M_\mathrm{dyn}$. Current uncertainties in stellar masses cannot account for values of $M_\star / M_\mathrm{dyn}$ above 1. Our results suggest that the homology hypothesis contained in the $M_\mathrm{dyn}$ formula above breaks down for compact galaxies. We provide an approximation to the virial coefficient $K \sim 6.0 \left[ r_\mathrm{e} / (3.185 \ \mathrm{kpc}) \right]^{-0.81} \left[ M_\star / (10^{11} \ \mathrm{M_\odot}) \right]^{0.45}$, which solves the mass discrepancy problem. A rough approximation to the dynamical mass is given by $M_\mathrm{dyn} \sim \left[ \sigma_\mathrm{e} / (200 \ \mathrm{km \ s^{-1}}) \right]^{3.6} \left[ r_\mathrm{e} / (3 \ \mathrm{kpc}) \right]^{0.35} 2.1 \times 10^{11} \ \mathrm{M_\odot}$.