Search

Over the last decades, evolutionary population synthesis models have powered an unmatched leap forward in our understanding of galaxies. From dating the age of the first galaxies in the Universe to detailed measurements of the chemical composition of nearby galaxies, the success of this approach built upon simple stellar population (SSP) spectro-photometric models is unquestionable. However, the internal constraints inherent to the construction of SSP models may hinder our ability to analyze the integrated spectra of galaxies in situations where the SSP assumption does not sufficiently hold. Thus, here we revisit the possibilities of fitting galaxy spectra as a linear combination of stellar templates without assuming any a priori knowledge on stellar evolution. We showcase the sensitivity of this alternative approach to changes in the stellar population properties, in particular the direct connection to variations in the stellar initial mass function, as well as its advantages when dealing with non-canonical integrated populations and semi-resolved observations. Furthermore, our analysis demonstrates that the absorption spectra of galaxies can be used to independently constrain tellar evolution theory beyond the limited conditions of the solar neighborhood., Comment: 5 pages, 4 figures. Accepted for publication in Astronomy & Astrophysics Letters

We studied the ionized gas in the inner region ($\sim680\times470$ pc$^2$) of the galaxy NGC 6868 using Gemini/GMOS integral field unit observations. Channel maps reveal complex kinematics and morphology, indicating multiple processes at work in NGC 6868. Through emission-line fitting, we identified two ubiquitous components in our data: a narrow ($\sigma\sim110$ km s$^{-1}$) tracing an ionized gas disc and a broad component ($\sigma\sim300$ km s$^{-1}$) mainly associated with inflowing/outflowing gas. The derived V-band reddening shows a spatial distribution consistent with that obtained from stellar population synthesis, although with generally higher values. For the first time, we measured the electron temperature in NGC 6868, finding values ranging from $\sim 14000$ K in the central region to $\gtrsim20000$ K with an outward increasing temperature gradient. The electron density map exhibits an inverse relationship, with central values reaching $N_e\sim4000$ cm$^{-3}$ for the broad component decreasing to $N_e\sim100$ cm$^{-3}$ towards the edges of the field of view. Using BPT diagrams, we found that all spaxels are consistent with both AGN and shock ionization. However, when this information is combined with our kinematic and temperature findings, and further supported by the WHAN diagram, we argue that an AGN is the dominant ionisation mechanism in the central region of NGC 6868, while the extended outer component is ionized by a combination of hot low-mass evolved stars and shocks. According to our findings, shocks play a significant role in the ionization balance of this galaxy., Comment: 16 pages, 16 figures, accepted for publication in MNRAS

The WHT Enhanced Area Velocity Explorer (WEAVE) is a new, massively multiplexing spectrograph. This new instrument will be exploited to obtain high S/N spectra of $\sim$25000 galaxies at intermediate redshifts for the WEAVE Stellar Population Survey (WEAVE-StePS). We test machine learning methods for retrieving the key physical parameters of galaxies from WEAVE-StePS-like spectra using both photometric and spectroscopic information at various S/Ns and redshifts. We simulated $\sim$105000 galaxy spectra assuming SFH with an exponentially declining star formation rate, covering a wide range of ages, stellar metallicities, sSFRs, and dust extinctions. We then evaluated the ability of the random forest and KNN algorithms to correctly predict such parameters assuming no measurement errors. We checked how much the predictive ability deteriorates for different S/Ns and redshifts, finding that both algorithms still accurately estimate the ages and metallicities with low bias. The dispersion varies from 0.08-0.16 dex for ages and 0.11-0.25 dex for metallicity, depending on the redshift and S/N. For dust attenuation, we find a similarly low bias and dispersion. For the sSFR, we find a very good constraining power for star-forming galaxies, log sSFR$\gtrsim$ -11, where the bias is $\sim$ 0.01 dex and the dispersion is $\sim$ 0.10 dex. For more quiescent galaxies, with log sSFR$\lesssim$ -11, we find a higher bias, 0.61-0.86 dex, and a higher dispersion, $\sim$ 0.4 dex, for different S/Ns and redshifts. Generally, we find that the RF outperforms the KNN. Finally, the retrieved sSFR was used to successfully classify galaxies as part of the blue cloud, green valley, or red sequence. We demonstrate that machine learning algorithms can accurately estimate the physical parameters of simulated galaxies even at relatively low S/N=10 per angstrom spectra with available ancillary photometric information., Comment: 19 pages, 10 + 2 figures, 4 tables, accepted in A&A

Near-infrared long-slit spectroscopy has been used to study the stellar population (SP) of the low luminosity active galactic nuclei (AGN) and matched analogues (LLAMA) sample. To perform the SP fits we have employed the X-shooter simple stellar population models together with the \st\ code. Our main conclusions are: The star formation history of the AGNs is very complex, presenting many episodes of star formation during their lifetimes. In general, AGN hosts have higher fractions of intermediate-age SP (light-weighted mean ages, $_L\lesssim$ 4.5 Gyr) when compared with their analogues ($_L\lesssim$ 8.0 Gyr). AGN are more affected by reddening and require significant fractions of featureless continuum and hot dust components. The ratio between the AGN radiated energy and the gravitational potential energy of the molecular gas ($E_{Rad}$/$E_{PG}$) for the AGN is compared with the \maL\ and a possible anti-correlation is observed. This suggests that the AGN is affecting the star formation in these galaxies, in the sense that more energetic AGN (log$(E_{Rad}$/$E_{PG}) \gtrsim 3$) tend to host nuclear younger SP ($ _L \lesssim$4Gyr). We found that the recent ($t<$2~Gyr) returned (recycled) stellar mass is higher in AGN than in the controls. We also provide evidence that the mass loss of stars would be enough to feed the AGN, thus providing observational constraints for models that predict that AGN feeding is partially due to the recycled gas from dying stars., Comment: Accepted for publication in MNRAS

Most of the observed galaxies cannot be resolved into individual stars and are studied through their integrated spectrum using simple stellar populations (SSPs) models, with stellar libraries being a key ingredient in building them. Spectroscopic observations are increasingly being directed towards the near-infrared (NIR), where much is yet to be explored. SSPs in the NIR are still limited, and there are inconsistencies between different sets of models. One of the ways to minimize this problem is to have reliable NIR stellar libraries. The main goal of this work is to present SMARTY (mileS Moderate resolution neAr-infRared sTellar librarY) a ~0.9-2.4$\mu$m stellar spectral library composed of 31 stars observed with the Gemini Near-IR Spectrograph (GNIRS) at the 8.1m Gemini North telescope and make it available to the community. The stars were chosen from the SMARTY library, for which the atmospheric parameters are reliable (and well tested), to populate different regions of the Hertzsprung-Russell (HR) diagram. Furthermore, five of these stars have NIR spectra available that we use to assess the quality of SMARTY. The remaining 26 stars are presented for the first time in the NIR. We compared the observed SMARTY spectra with synthetic and interpolated spectra, finding a mean difference of ~20% in the equivalent widths and ~1% in the overall continuum shape in both sets of comparisons. We computed the spectrophotometric broadband magnitudes and colours and compared them with the 2MASS ones, resulting in mean differences up to 0.07 and 0.10mag in magnitudes and colours, respectively. In general, a small difference was noted between the SMARTY spectra corrected using the continuum from the interpolated and the theoretical stars.

We present new H- and K-band spectroscopy for the bulge of M31, taken with the LUCI spectrograph at the Large Binocular Telescope (LBT). We studied radial trends of CO absorption features (namely, CO1.58, CO1.60, CO1.64, CO1.66, CO1.68, CO2.30, CO2.32, CO2.35) in the bulge of M31, out to a galactocentric distance of 100'' (380pc). We find that most COs do not exhibit a strong radial gradient, despite the strong metallicity gradient inferred from the optical spectral range, except for CO1.64, showing a steep increase in the center. We compared the observed line strengths to predictions of different state-of-the-art stellar population models, including an updated version of EMILES models, which also uses the extended IRTF spectral library. The observed COs are close to models' predictions, but in some models they turn out to be underestimated. We find that the lack of radial gradients is due to the combination of increasing CO strength with metallicity and C abundance, and decreasing CO strength with IMF slope and O abundance. We speculate that the steep gradient of CO1.64 might be due to Na overabundance. Remarkably, we were able to fit, at the same time, optical indices and all the NIR COs except for CO1.68, leaving abundance ratios (i.e., [C/Fe], [O/Fe], and [Mg/Fe]) as free-fitting parameters, imposing age and metallicity constraints from the optical, with no significant contribution from intermediate-age populations. For the majority of the bulge, we find [Mg/Fe]~0.15dex, [O/Fe] larger than [Mg/Fe] (by ~0.1dex), and C abundance consistent with that of Mg. In the central (few arcsec) region, we still find an enhancement of O and Mg, but significantly lower [C/Fe]. We find that the COs' line strengths of the bulge are significantly lower than those of massive galaxies, possibly because of a difference in carbon abundance, as well as, to some extent, total metallicity., Comment: 20 pages, 10 figures, accepted for publication on A&A

The goal of this work is to scrutinise the surface brightness fluctuation (SBF) calculation methodology. We analysed the SBF derivation procedure, measured the accuracy of the fitted SBF under controlled conditions, retrieved the uncertainty associated with the variability of a system that is inherently stochastic, and studied the SBF reliability under a wide range of conditions. Additionally, we address the possibility of an SBF gradient detection. We also examine the problems related with biased measurements of the SBF and low luminosity sources. All of this information allows us to put forward guidelines to ensure a valid SBF retrieval. To perform all these experiments, we carried out Monte Carlo simulations of mock galaxies as an ideal laboratory. Knowing its underlying properties, we attempted to retrieve SBFs under different conditions. We show how the SBF uncertainty can be obtained and we present a collection of qualitative recommendations for a safe SBF retrieval: it is important to model the instrumental noise, rather than fitting it; the target galaxies must be observed under appropriate observational conditions; in a traditional SBF derivation, one should avoid pixels with fluxes lower than ten times the SBF estimate to prevent biased results. We offer our computational implementation in the form of a simple code designed to estimate the uncertainty of the SBF measurement. This code can be used to predict the quality of future observations or to evaluate the reliability of those already conducted., Comment: 25 pages, 14 figures, 5 tables, accepted by A&A

In our previous work, we derive the CR-IGIMF: a new scenario for a variable stellar initial mass function (IMF), which combines numerical results on the role played by cosmic rays in setting the thermal state of star forming gas, with the analytical approach of the integrated galaxy-wide IMF. In this work, we study the implications of this scenario for the properties of local Early-Type galaxies (ETG), as inferred from dynamical, photometric and spectroscopic studies. We implement a library of CR-IGIMF shapes in the framework of the Galaxy Evolution and Assembly (GAEA) model. Our realization includes a derivation of synthetic spectral energy distribution for each model galaxy, allowing a direct derivation of the mass fraction in the mean IMF of low-mass stars (i.e. the dwarf-to-giant ratio - $\rm f_{dg}$) and a comparison with IMF sensitive spectral features. The predictions of the GAEA model implementing the CR-IGIMF confirm our previous findings: it correctly reproduces both the observed excess of z$\sim$0 dynamical mass (mass-to-light ratios) with respect to spectroscopic (photometric) estimates assuming a universal, MW-like, IMF, and the observed increase of [$\alpha$/Fe] ratios with stellar mass in spheroidal galaxies. Moreover, this realization reproduces the increasing trends of $\rm f_{dg}$, and IMF-sensitive line-strengths with velocity dispersion, although the predicted relations are significantly shallower than the observed ones. Our results show that the CR-IGIMF is a promising scenario that reproduces at the same time dynamical, photometric and spectroscopic indications of a varying IMF in local ETGs. The shallow relations found for spectral indices suggest that either a stronger variability as a function of galaxy properties or additional dependences (e.g. as a function of star forming gas metallicity) might be required to match the strength of the observed trends., Comment: 13 pages, 8 figures, 2 tables, A&A accepted

The upcoming generation of optical spectrographs on four meter-class telescopes, with their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage, will provide high-quality spectra for thousands of galaxies. These data will allow us to examine of the stellar population properties at intermediate redshift, an epoch that remains unexplored by large and deep surveys. We assess our capability to retrieve the mean stellar metallicity in galaxies at different redshifts and S/N, while simultaneously exploiting the UV and optical rest-frame wavelength coverage. The work is based on a comprehensive library of spectral templates of stellar populations, covering a wide range of age and metallicity values and built assuming various SFHs. We simulated realistic observations of a large sample of galaxies carried out with WEAVE at the WHT at different redshifts and S/N values. We measured all the reliable indices on the simulated spectra and on the comparison library. We then adopted a Bayesian approach to obtain the probability distribution of stellar metallicity. The analysis of the spectral indices has shown how some mid-UV indices can provide reliable constraints on stellar metallicity, along with optical indicators. The analysis of the mock observations has shown that even at S/N=10, the metallicity can be derived within 0.3 dex, in particular, for stellar populations older than 2 Gyr. Our results are in good agreement with other theoretical and observational works in the literature and show how the UV indicators can be advantageous in constraining metallicities. This is very promising for the upcoming surveys carried out with new, highly multiplexed, large-field spectrographs, such as StePS at the WEAVE and 4MOST, which will provide spectra of thousands of galaxies covering large spectral ranges at relatively high S/N., Comment: 16 pages, 14 figures. Submitted 31/03/2023, Accepted 20/07/2023

We present spaxel-by-spaxel stellar population fits for the $\sim$10 thousand MaNGA datacubes. We provide multiple extension fits files, nominated as MEGACUBES, with maps of several properties as well as emission-line profiles that are provided for each spaxel. All the MEGACUBES are available through a web interface (\url{https://manga.linea.org.br/} or \url{http://www.if.ufrgs.br/~riffel/software/megacubes/}). We also defined a final Active Galactic Nuclei (AGN) sample, as well as a control sample matching the AGN host galaxy properties. We have analysed the stellar populations and spatially resolved emission-line diagnostic diagrams of these AGNs and compared them with the control galaxies sample. We find that the relative fractions of young ($t \leq $56 Myr) and intermediate-age (100 Myr $\leq t \leq$ 2 Gyr) show predominantly a positive gradient for both AGNs and controls. The relative fraction of intermediate-age stellar population is higher in AGN hosts when compared to the control sample, and this difference becomes larger for higher [O III] luminosity AGNs. We attribute this to the fact that extra gas is available in these more luminous sources and that it most likely originates from mass-loss from the intermediate-age stars. The spatially resolved diagnostic diagrams reveal that the AGN emission is concentrated in the inner 0.5 $R_e$ (effective radius) region of the galaxies, showing that the AGN classification is aperture dependent and that emission-line ratios have to be taken together with the H$\alpha$ equivalent width for proper activity classification. We present a composite ``BPT+WHAN" diagram that produces a more comprehensive mapping of the gas excitation., Comment: Accepted for publication in MNRAS; small typos corrected

We present a new library of semi-empirical stellar population models that are based on the empirical MILES and semi-empirical sMILES stellar libraries. The models span a large range of age and metallicity, in addition to an [$\alpha$/Fe] coverage from $-$0.2 to $+$0.6 dex, at MILES resolution (FWHM=2.5$ \mathring {\mathrm A}$) and wavelength coverage (3540.5-7409.6$ \mathring {\mathrm A}$). These models are aimed at exploring abundance ratios in the integrated light from stellar populations in star clusters and galaxies. Our approach is to build SSPs from semi-empirical stars at particular [$\alpha$/Fe] values, thus producing new SSPs at a range of [$\alpha$/Fe] values from sub-solar to super-solar. We compare these new SSPs with previously published and well-used models and find similar abundance pattern predictions, but with some differences in age indicators. We illustrate a potential application of our new SSPs, by fitting them to the high signal-to-noise data of stacked SDSS galaxy spectra. Age, metallicity and [$\alpha$/Fe] trends were measured for galaxy stacks with different stellar velocity dispersions and show systematic changes, in agreement with previous analyses of subsets of those data. These new SSPs are made publicly available., Comment: 21 pages, 18 figures plus Supplementary Material. Accepted for publication in MNRAS

The upcoming new generation of optical spectrographs on four-meter-class telescopes will provide valuable opportunities for forthcoming galaxy surveys through their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage. WEAVE is a new wide-field spectroscopic facility mounted on the 4.2 m William Herschel Telescope in La Palma. WEAVE-StePS is one of the five extragalactic surveys that will use WEAVE during its first five years of operations. It will observe galaxies using WEAVE MOS (~950 fibres across a field of view of ~3 deg2 on the sky) in low-resolution mode (R~5000, spanning the wavelength range 3660-9590 AA). WEAVE-StePS will obtain high-quality spectra (S/N ~ 10 per AA at R~5000) for a magnitude-limited (I_AB = 20.5) sample of ~25,000 galaxies, the majority selected at z>=0.3. The survey goal is to provide precise spectral measurements in the crucial interval that bridges the gap between LEGA-C and SDSS data. The wide area coverage of ~25 deg2 will enable us to observe galaxies in a variety of environments. The ancillary data available in each observed field (including X-ray coverage, multi-narrow-band photometry and spectroscopic redshift information) will provide an environmental characterisation for each observed galaxy. This paper presents the science case of WEAVE-StePS, the fields to be observed, the parent catalogues used to define the target sample, and the observing strategy chosen after a forecast of the expected performance of the instrument for our typical targets. WEAVE-StePS will go back further in cosmic time than SDSS, extending its reach to encompass more than ~6 Gyr, nearly half of the age of the Universe. The spectral and redshift range covered by WEAVE-StePS will open a new observational window by continuously tracing the evolutionary path of galaxies in the largely unexplored intermediate-redshift range., Comment: 15 pages, 9 figures, A&A in press

The X-shooter Spectral Library (XSL) is a large empirical stellar library used as a benchmark for the development of stellar population models. The inclusion of $\alpha$-elements abundances is crucial to disentangling the chemical evolution of any stellar system. The aim of this paper is to provide a catalogue of high-precision and accurate magnesium and calcium abundances from a wide variety of stars well distributed in the Hertzsprung-Russell (HR) diagram. We originally performed an analysis of the derived Mg and Ca abundances for medium-resolution spectra of 611 stars from the XSL Data Release 2. For this purpose, we used the GAUGUIN automated abundance estimation code to fit the ultraviolet-blue (UVB) and visible (VIS) spectra. We tested the consistency of the atmospheric parameters and chemical abundances with the Gaia DR3 and the AMBRE Project datasets. We have finally obtained precise [Mg/Fe] and [Ca/Fe] abundances for 192 and 217 stars respectively, from which 174 stars have measurements in both elements. The stars cover a broad range of effective temperature 4000 < T$_{\rm eff}$ < 6500K, surface gravity 0.3 < log(g) < 4.8 cm s$^{\rm -2}$, and metallicity -2.5 < [Fe/H] < +0.4 dex. We find an excellent agreement with the abundance estimates from the AMBRE:HARPS and the Gaia/RVS (Radial Velocity Spectrometer) analysis. Moreover, the resulting abundances reproduce a plateau in the metal-poor regime followed by a decreasing trend even at supersolar metallicities, as predicted by Galactic chemical evolution models. This catalogue is suitable for improving the modelling of evolutionary stellar population models with empirical $\alpha$-enhancements, which could significantly contribute to the analysis of external galaxies abundances in the near future., Comment: Accepted for publication in A&A

Near-ultraviolet (NUV) spectroscopic studies have suggested that passively evolving massive, early-type galaxies host sub-one percent fractions of young stars in their innermost regions. We shed light on the origin of these stars by analysing NGC 1277, a widely studied nearby prototypical massive compact relic galaxy. These are rare galaxies that have survived without experiencing significant size evolution via accretion and mergers since their formation at high redshift. We obtain a spectrum in the UV range within the central 1 kpc region of NGC 1277. We compare a carefully selected set of optical and NUV line-strengths to model predictions with star formation histories characteristic of massive galaxies. We find a 0.8% mass fraction of young stars in the centre of NGC 1277, similar to that found in massive early-type galaxies. Given the limited accretion history of NGC 1277, these results favour an intrinsic, in-situ, process triggering star formation at later epochs. Our results suggest a general constraint on the amount of young stars in the cores of massive early-type galaxies. This amount should be assumed as an upper limit for the young stellar contribution in massive galaxies, as there might be present other contributions from evolved stars., Comment: 6 pages, 3 figures

We perform a detailed study of the stellar populations in a sample of massive Fornax dwarf galaxies using a set of newly defined line indices. Using data from the Integral field spectroscopic data, we study abundance ratios of eight dEs with stellar mass ranging from 10$^8$ to 10$^{9.5}$ M$_\odot$ in the Fornax cluster. We present the definitions of a new set of high-resolution Lick-style indices to be used for stellar population studies of unresolved small stellar systems. We identify 23 absorption features and continuum regions, mainly dominated by 12 elements (Na, Ca, Sc, Ti, V, Cr, Mn, Fe, Ni, Y, Ba and Nd) in the wavelength range 4700 - 5400 \r{A} and characterise them as a function of age, metallicity and alpha element abundance ratios. We analyse eight dEs and interpret the line strengths, measured in our new high resolution system of indices, with the aid of stellar population models with high enough spectral resolution. We obtain abundance ratio proxies for a number of elements that have never been studied before for dwarf ellipticals outside the Local Group. These proxies represent relative deviations from predicted index-strengths of base stellar population models built-up following the abundance pattern of The Galaxy. The abundance proxy trend results are compared to abundance ratios from resolved stars in the Local Group, and indices from integrated light of larger early-type galaxies. We find that all our dwarfs show a pattern of abundance ratios consistent with the disk of the Milky Way, indicative of slow formation in comparison to their high mass counterparts., Comment: Accepted for publication in MNRAS, 31 pages, 30 figures

Massive Early-Type Galaxies (ETGs) show several strong CO absorption features in their H- and K-band spectra that cannot be explained by state-of-the-art stellar population models. For many years, the disagreement has been attributed to the presence of intermediate-age stellar components that are dominated by stars in the Asymptotic Giant Branch (AGB) phase. However, no robust evidence of this scenario has been provided so far. One way to test this claim is by comparison of CO indices for ETGs and for relic galaxies. Lacking the intermediate-age stellar populations, relic galaxies provide us with a unique opportunity to address the origin of strong CO absorptions in ETGs. Here, we utilize the prototype relic galaxy NGC 1277 and compare the CO absorption features of this galaxy with the ones of a representative sample of massive ETGs. We show that the CO lines in both systems have similar strengths, significantly stronger than the predictions of stellar population synthesis models. We conclude that intermediate-age stellar populations in massive ETGs are not the culprit of the strong CO absorptions., Comment: 6 pages, 3 figures, accepted for publication in MNRAS Letters

The leading method for the determination of relevant stellar population parameters of unresolved extragalactic Globular Clusters is through the study of their integrated spectroscopy, where Balmer line-strength indices are considered to be age sensitive. Previously, a splitting in the highly optimised spectral line-strength index H$\beta_o$ was observed in a sample of Galactic globular clusters at all metallicities resulting in an apparent "upper branch" and "lower branch" of globular clusters in the H$\beta_o$ - [MgFe] diagram. This was suggested to be caused by the presence of hot Blue straggler stars (BSSs), resulting in an underestimation of 'spectroscopic' ages in the upper branch. Over a decade on, we look to re-evaluate these findings. We make use of new, large Galactic Globular Cluster integrated spectroscopy datasets. To produce a large, homogeneously combined sample we have considered a number of factors including the radial dependence of Balmer and metal lines. Using this new sample, in disagreement with previous work, we find the splitting in H$\beta_o$ only occurs at intermediate to high metallicities ([M/H]$>-1$), and is not the result of an increased fraction of BSSs, but rather is due to an increased Helium abundance. We explore the possible impact of varying Helium on simple stellar population models to provide a theoretical basis for our hypothesis and then use the relationship between upper branch candidacy and enhanced Helium to predict the Helium content of three M31 clusters. We discuss what this can tell us about their mass and fraction of first generation stars., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

Massive Early-Type Galaxies (ETG) in the local Universe are believed to be the most mature stage of galaxy evolution. Their stellar population content reveals the evolutionary history of these galaxies. However, while state-of-the-art Stellar Population Synthesis (SPS) models provide an accurate description of observed galaxy spectra in the optical range, the modelling in the Near-Infrared (NIR) is still in its infancy. Here we focus on NIR CO absorption features to show, in a systematic and comprehensive manner, that for massive ETGs, all CO indices, from H through to K band, are significantly stronger than currently predicted by SPS models. We explore and discuss several possible explanations of this "CO mismatch", including the effect of intermediate-age, AGB-dominated, stellar populations, high metallicity populations, non-solar abundance ratios and the initial mass function. While none of these effects is able to reconcile models and observations, we show that ad-hoc "empirical" corrections, taking into account the effect of CO-strong giant stars in the low-temperature regime, provide model predictions that are closer to the observations. Our analysis points to the effect of carbon abundance as the most likely explanation of NIR CO line-strengths, indicating possible routes for improving the SPS models in the NIR., Comment: 24 pages, 12 figures, accepted for publication in MNRAS

We present simple stellar population models based on the empirical X-shooter Spectral Library (XSL) from NUV to NIR wavelengths. The unmatched characteristics of relatively high resolution and extended wavelength coverage ($350-2480$ nm, $R\sim10\,000$) of the XSL population models bring us closer to bridging optical and NIR studies of intermediate and old stellar populations. It is now common to find good agreement between observed and predicted NUV and optical properties of stellar clusters due to our good understanding of the main-sequence and early giant phases of stars. However, NIR spectra of intermediate-age and old stellar populations are sensitive to cool K and M giants. The asymptotic giant branch, especially the thermally pulsing asymptotic giant branch, shapes the NIR spectra of $0.5-2$ Gyr old stellar populations; the tip of the red giant branch defines the NIR spectra of populations with ages larger than that. We construct sequences of the average spectra of static giants, variable-rich giants, and C-rich giants to include in the models separately. The models span the metallicity range $-2.2<[Fe/H]<+0.2$ and ages above 50 Myr, a broader range in the NIR than in other models based on empirical spectral libraries. Our models can reproduce the integrated optical colours of the Coma cluster galaxies at the same level as other semi-empirical models found in the literature. In the NIR, there are notable differences between the colours of the models and Coma cluster galaxies. The XSL models expand the range of predicted values of NIR indices compared to other models based on empirical libraries. Our models make it possible to perform in-depth studies of colours and spectral features consistently throughout the optical and the NIR range to clarify the role of evolved cool stars in stellar populations., Comment: 30 pages, 26 figures, accepted to Astronomy & Astrophysics, models will be available on http://xsl.astro.unistra.fr/ upon publishing

We present the third data release (DR3) of the X-shooter Spectral Library (XSL). This moderate-to-high resolution, near-ultraviolet-to-near-infrared ($350-2480$ nm, R $\sim$ 10 000) spectral library is composed of 830 stellar spectra of 683 stars. DR3 improves upon the previous data release by providing the combined de-reddened spectra of the three X-shooter segments over the full $350-2480$ nm wavelength range. It also includes additional 20 M-dwarf spectra from the ESO archive. We provide detailed comparisons between this library and Gaia EDR3, MILES, NGSL, CaT library, and (E-)IRTF. The normalised rms deviation is better than $D=0.05$ or 5$\%$ for the majority of spectra in common between MILES (144 spectra of 180), NGSL (112$/$116), and (E-)IRTF (55$/$77) libraries. Comparing synthetic colours of those spectra reveals only negligible offsets and small rms scatter, such as the median offset(rms) 0.001$\pm$0.040 mag in the (box1-box2) colour of the UVB arm,-0.004$\pm$0.028 mag in (box3-box4) of the VIS arm, and -0.001$\pm$0.045 mag in (box2-box3) colour between the UVB and VIS arms, when comparing stars in common with MILES. We also find an excellent agreement between the Gaia published (BP-RP) colours and those measured from the XSL DR3 spectra, with a zero median offset and an rms scatter of 0.037 mag for 449 non-variable stars. The unmatched characteristics of this library, which combine a relatively high resolution, a large number of stars, and an extended wavelength coverage, will help us to bridge the gap between the optical and the near-IR studies of intermediate and old stellar populations, and to probe low-mass stellar systems., Comment: 26 pages, 25 figures, accepted to Astronomy & Astrophysics. The data are available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ or on the XSL web-page http://xsl.astro.unistra.fr

The aim of this work is to explore the potential of Surface Brightness Fluctuations (SBF) for studying composite stellar populations (CSP). To do so, we have computed the standard (mean) and SBF spectra with E-MILES stellar population synthesis code. We have created a set of models composed by different mass fractions of two single stellar populations (SSP), as a first approximation of a CSP scenario. With these models we present an ensemble of SBF colour-colour diagnostic diagrams that reveal different secondary populations depending on the bands used. For this work we focus on those colours capable of unveiling small fractions of metal-poor components in elliptical galaxies, which are dominated by old metal-rich stellar populations. We fit a set of synthetic models and a selection of nearby elliptical galaxies to our CSP models using both mean and SBF colours. We find that the results are highly improved and return small secondary components when mean and SBF values are applied simultaneously, instead of employing them separately or as a constraint. Finally, we explore the possibility of tracking chemical enrichment histories by including in the analysis a variety of SBF colours. For this purpose we present an example where, with two different SBF colour-colour diagrams, we untangle a small contribution of a young solar population and an old metal-poor component from an old solar principal population. The results we have found are promising, but limited by the available data. We highlight the urgent need for new, better and more consistent SBF observations., Comment: 25 pages, 19 figures, Published MNRAS

Using new, homogeneous, long-slit spectroscopy in the wavelength range from ~0.35 to ~1micron, we study radial gradients of optical and near-infrared (NIR) IMF-sensitive features along the major axis of the bulge of M31, out to a galacto-centric distance of ~200'' (~800pc). Based on state-of-the-art stellar population synthesis models with varying Na abundance ratio, we fit a number of spectral indices, from different chemical species (including TiO's, Ca, and Na indices), to constrain the low-mass (<0.5M_Sun) end slope (i.e. the fraction of low-mass stars) of the stellar IMF, as a function of galacto-centric distance. Outside a radial distance of ~10'', we infer an IMF similar to a Milky-Way-like distribution, while at small galacto-centric distances, an IMF radial gradient is detected, with a mildly bottom-heavy IMF in the few inner arcsec. We are able to fit Na features (both NaD and NaI8190), without requiring extremely high Na abundance ratios. [Na/Fe] is ~0.4dex for most of the bulge, rising up to ~0.6dex in the innermost radial bins. Our results imply an overall, luminosity-weighted, IMF and mass-to-light ratio for the M31 bulge, consistent with those for a Milky-Way-like distribution, in contrast to results obtained, in general, for most massive early-type galaxies., Comment: accepted for publication in MNRAS, 22 pages, 16 figures

We present a new library of semi-empirical stellar spectra that is based on the empirical MILES library. A new, high resolution library of theoretical stellar spectra is generated that is specifically designed for use in stellar population studies. We test these models across their full wavelength range against other model libraries and find reasonable agreement in their predictions of spectral changes due to atmospheric $\alpha$-element variations, known as differential corrections. We also test the models against the MILES and MaStar libraries of empirical stellar spectra and also find reasonable agreements, as expected from previous work. We then use the abundance pattern predictions of the new theoretical stellar spectra to differentially correct MILES spectra to create semi-empirical MILES (sMILES) star spectra with abundance patterns that differ from those present in the Milky Way. The final result is 5 families of 801 sMILES stars with [$\alpha$/Fe] abundances ranging from $-$0.20 to 0.60 dex at MILES resolution (FWHM=$2.5\,${\AA}) and wavelength coverage ($3540.5-7409.6\,${\AA}). We make the sMILES library publicly available., Comment: 26 pages, 18 figures plus Supplementary Material. Accepted for publication in MNRAS. Corrected a minor font size error and two figure caption placements from the previous version

Stellar population studies provide unique clues to constrain galaxy formation models. So far, detailed studies based on absorption line strengths have mainly focused on the optical spectral range although many diagnostic features are present in other spectral windows. In particular, the near-infrared (NIR) can provide a wealth of information about stars, such as evolved giants, that have less evident optical signatures. Due to significant advances in NIR instrumentation and extension of spectral libraries and stellar population synthesis (SPS) models to this domain, it is now possible to perform in-depth studies of spectral features in the NIR to a high level of precision. In the present work, taking advantage of state-of-the-art SPS models covering the NIR spectral range, we introduce a new set of NIR indices constructed to be maximally sensitive to the main stellar population parameters, namely age, metallicity and initial mass function (IMF). We fully characterize the new indices against these parameters as well as their sensitivity to individual elemental abundance variations, velocity dispersion broadening, wavelength shifts, signal-to-noise ratio and flux calibration. We also present, for the first time, a method to ensure that the analysis of spectral indices is not affected by sky contamination, which is a major challenge when dealing with NIR spectroscopy. Moreover, we discuss two main applications: (i) the ability of some NIR spectral indices to constrain the shape of the low-mass IMF and (ii) current issues in the analysis of NIR spectral indices for future developments of SPS modelling., Comment: 38 pages, 23 figures, published in MNRAS

The path towards robust near-infrared extensions of stellar population models involves the confrontation between empirical and synthetic stellar spectral libraries across the wavelength ranges of photospheric emission. [...] With its near-UV to near-IR coverage, the X-shooter Spectral Library (XSL) allows us to examine to what extent models succeed in reproducing stellar energy distributions (SEDs) and stellar absorption line spectra simultaneously. This study compares the stellar spectra of XSL with the PHOENIX spectra of the G\"ottingen Spectral Library. The comparison is carried out both separately in the three arms of the X-shooter spectrograph, and jointly across the whole spectrum. When adopting the stellar parameters published with data release DR2 of XSL, we find that the SEDs of the models are consistent with those of the data at Teff > 5000 K. Below 5000 K, there are significant discrepancies in the SEDs. When leaving the stellar parameters free to adjust, satisfactory representations of the SEDs are obtained down to about 4000 K. However, in particular below 5000 K and in the UVB spectral range, strong local residuals associated with intermediate resolution spectral features are then seen; the necessity of a compromise between reproducing the line spectra and reproducing the SEDs leads to dispersion between the parameters favored by various spectral ranges. We describe the main trends observed and we point out localized offsets between the parameters preferred in this global fit to the SEDs and the parameters in DR2. These depend in a complex way on position in the HR diagram (HRD). We estimate the effect of the offsets on bolometric corrections as a function of position in the HRD and use this for a brief discussion of their impact on the studies of stellar populations. [abridged], Comment: Accepted for publication in Astronomy and Astrophysics (24 Nov. 2020). 18 pages (10 figures) in the main text + 19 pages (22 figures) in the appendices (the appendices are spread out for better readability)

Central cluster galaxies are the largest and most massive galaxies in the Universe. Although they host very old stellar populations, several studies found the existence of blue cores in some BCGs indicating ongoing star formation. We analyse VLT/X-Shooter stacked spectra of 6 nearby massive central galaxies with high central velocity dispersions ($\sigma$>300 km/s) at different galactocentric distances. We quantify the young stellar population out to 4 kpc by fitting near-UV and optical absorption line indices with predictions of composite stellar populations modelled by an old and a young stellar component. We also use IMF-sensitive indices since these galaxies have been found to host a bottom-heavy IMF in their central regions. We derive negative young stellar populations gradients, with mass fractions of stars younger than 1 Gyr decreasing with galactocentric distance, from 0.70% within 0.8 kpc to zero beyond 2 kpc. We also measure the mass fraction in young stars for individual galaxies in the highest S/N central regions. All the galaxies have young components of less than one percent. Our results clearly suggest that the star formation in massive central cluster galaxies takes place in their galaxy cores (<2 kpc), which, with deeper gravitational potential wells, are capable of retaining more gas. Among the possible sources for the gas required to form these young stars, our results are consistent with an in-situ origin via stellar evolution, which is sufficient to produce the observed young stellar populations.

We report on the detection of a small contribution (around and below 1% in mass) from young stellar components with ages $\leq$20 Myr in low-mass galaxies purposely selected from the MaNGA survey to be already-quenched systems. Among the sample of 28 galaxies, eight of them show signatures of having suffered a very recent burst of star formation. The detection has been done through the analysis of line-strength indices in the red spectral range [5700,8800] \AA. The increasing contribution of red supergiants to this red regime is responsible for a deviation of the index measurements with respect to their position within the model grids in the standard spectral range [3600,5700] \AA. We demonstrate that a combination of red indices, as well as a qualitative assessment of the mean luminosity-weighted underlying stellar population, is required in order to distinguish between a true superyoung population and other possible causes of this deviation, such as abundance ratio variations. Our result implies that many presumably quenched low-mass galaxies actually contain gas that is triggering some level of star formation. They have, therefore, either accreted external gas, internally recycled enough gas from stellar evolution to trigger new star formation, or they kept a gas reservoir after the harassment or stripping process that quenched them in first place. Internal processes are favoured since we find no particular trends between our non-quenched galaxies and their environment, although more work is needed to fully discard an external influence., Comment: 12 pages, 6 figures, MNRAS in press

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

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., 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

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)

We present the second data release (DR2) of the X-Shooter Spectral Library (XSL), which contains all the spectra obtained over the six semesters of that program. This release supersedes our first data release from Chen et al. 2014, with a larger number of spectra (813 observations of 666 stars) and with a more extended wavelength coverage as the data from the near-infrared arm of the X-Shooter spectrograph are now included. The DR2 spectra then consist of three segments that were observed simultaneously and, if combined, cover the range between $\sim$300 nm and $\sim$2.45 $\mu$m at a spectral resolving power close to $R=10\,000$. The spectra were corrected for instrument transmission and telluric absorption, and they were also corrected for wavelength-dependent flux-losses in 85% of the cases. On average, synthesized broad-band colors agree with those of the MILES library and of the combined IRTF and Extended IRTF libraries to within $\sim\!1$%. The scatter in these comparisons indicates typical errors on individual colors in the XSL of 2$-$4 %. The comparison with 2MASS point source photometry shows systematics of up to 5% in some colors, which we attribute mostly to zero-point or transmission curve errors and a scatter that is consistent with the above uncertainty estimates. The final spectra were corrected for radial velocity and are provided in the rest-frame (with wavelengths in air). The spectra cover a large range of spectral types and chemical compositions (with an emphasis on the red giant branch), which makes this library an asset when creating stellar population synthesis models or for the validation of near-ultraviolet to near-infrared theoretical stellar spectra across the Hertzsprung-Russell diagram., Comment: Main paper: 16 pages, Full list of stars from pages 21 to 34, Typical XSL spectra on Figure 10 page 10

Early-type galaxies are considered to be the end-products of massive galaxy formation. Optical spectroscopic studies reveal that massive early-type galaxies formed the bulk of their stars over short timescales (12), followed by passive evolution to the present. However, their optical spectra are insensitive to constrain small episodes of recent star formation, since they are dominated by old stars. Fortunately, this problem can be tackled in the ultraviolet range. While recent studies that make use of ultraviolet absorption lines have suggested the presence of young stars in a few early-type galaxies, the age and mass fractions of young stars and their dependence on galaxy mass, is unknown. Here we report a detailed study of these young stellar populations, from high-quality stacked spectra of 28,663 galaxies from the BOSS survey, analysing optical and ultraviolet absorption lines simultaneously. We find that residual star formation is ubiquitous in massive early-type galaxies, measuring average mass fractions of ~0.5% in young stars in the last 2 Gyr of their evolution. This fraction shows a decreasing trend with galaxy stellar mass, consistent with a down-sizing scenario. We also find that synthetic galaxies from state-of-the-art cosmological numerical simulations significantly overproduce both intermediate and young stellar populations. Therefore, our results pose stringent constraints on numerical simulations of galaxy formation., Comment: Nature Astronomy, December 16, 2019. 21 pages, 8 figures, 2 tables. The published version is available at the following URL: https://www.nature.com/articles/s41550-019-0955-0

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., Comment: 18 pages, 20 figures, accepted for publication in A&A

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

Growing evidence in recent years suggests a systematic variation of the stellar initial mass function (IMF), being top-heavy for starburst galaxies and possibly bottom-heavy for massive ellipticals. Galaxy chemical evolution simulations adopting an invariant canonical IMF face difficulty in simultaneously reproducing the metallicity and $\alpha$-enhancement of the massive elliptical galaxies. Applying a variable IMF that changes with time is a promising solution, however, it is non-trivial to couple a variable IMF theory with the existing galaxy evolution codes. Here we present the first open source simulation code which recalculates the galaxy-wide IMF at each time step according to the Integrated-Galactic-IMF (IGIMF) theory where the galaxy-wide IMF depends on the galactic star formation rate and metallicity. The resulting galaxy-wide IMF and metal abundance evolve with time. With this pilot work, we explore the effect of the IGIMF theory on galaxy chemical evolution in comparison with an invariant IMF., Comment: 16 pages, 18 figures, accepted for publication in A&A

The X-shooter Spectral Library (XSL) is an empirical stellar library at medium spectral resolution covering the wavelength range from 3000 \AA to 24 800 \AA. This library aims to provide a benchmark for stellar population studies. In this work, we present a uniform set of stellar atmospheric parameters, effective temperatures, surface gravities, and iron abundances for 754 spectra of 616 XSL stars. We used the full-spectrum fitting package ULySS with the empirical MILES library as reference to fit the ultraviolet-blue (UVB) and visible (VIS) spectra. We tested the internal consistency and we compared our results with compilations from the literature. The stars cover a range of effective temperature 2900 < Teff < 38 000 K, surface gravity 0 < log g < 5.7, and iron abundance -2.5 < [Fe/H] < +1.0, with a couple of stars extending down to [Fe/H] = -3.9. The precisions of the measurements for the G- and K-type stars are 0.9%, 0.14, and 0.06 in Teff, log g, and [Fe/H], respectively. For the cool giants with log g < 1, the precisions are 2.1%, 0.21, and 0.22, and for the other cool stars these values are 1%, 0.14, and 0.10. For the hotter stars (Teff > 6500 K), these values are 2.6%, 0.20, and 0.10 for the three parameters., Comment: 15 pages, 12 figures, 3 tables, published in A&A

We construct Schwarzschild orbit-based models of NGC 7457, known as a peculiar low-mass lenticular galaxy. Our best-fitting model successfully retrieves most of the unusual kinematics behaviours of this galaxy, in which, the orbital distribution of stars is dominated by warm and hot orbits. The reconstructed surface brightness of the hot component matches fairly well the photometric bulge and the reconstructed LOSVD map of this component shows clear rotation around the major photometric axis of the galaxy. In the absence of a dominant cold component, the outer part of our model is dominated by warm orbits, representing an exponential thick disc. Our orbital analysis also confirm the existence of a counter-rotating orbital substructure in the very centre, reported in previous observational studies. By comparing our model with a variety of simulation studies, and considering the stellar kinematics and populations properties of this galaxy, we suggest that the thick disc is most likely a dynamically heated structure, formed through the interactions and accretion of satellite(s) with near-polar initial inclination. We also suggest a merger-driven process as the most plausible scenario to explain the observed and dynamically-modelled properties of the bulge of NGC 7457. We conclude that both the high level of cylindrical rotation and unusually low velocity dispersion reported for the NGC 7457 have most-likely external origins. Therefore, NGC 7457 could be considered as a candidate for merger-driven cylindrical rotation in the absence of a strong bar in disc galaxies., Comment: 15 pages, 12 figures, 2 tables. Accepted by MNRAS

Using samples drawn from the Sloan Digital Sky Survey, we study for the first time the relation between large-scale environments (Clusters, Groups and Voids) and the stellar Initial Mass Function (IMF). We perform an observational approach based on the comparison of IMF-sensitive indices of quiescent galaxies with similar mass in varying environments. These galaxies are selected within a narrow redshift interval ($ 0.020 < z < 0.055 $) and spanning a range in velocity dispersion from 100 to 200 kms$^{-1}$. The results of this paper are based upon analysis of composite spectra created by stacking the spectra of galaxies, binned by their velocity dispersion and redshift. The trends of spectral indices as measured from the stacked spectra, with respect to velocity dispersion, are compared in different environments. We find a lack of dependence of the IMF on the environment for intermediate-mass galaxy regime. We verify this finding by providing a more quantitative measurement of the IMF variations among galactic environments using MILES stellar population models with a precision of $\Delta\Gamma_{b}\sim0.2$., Comment: Accepted for publication in MNRAS

We analyze a set of optical-to-near-infrared long-slit nuclear spectra of 16 infrared-luminous spiral galaxies. All of the studied sources present H$_2$ emission, which reflects the star-forming nature of our sample, and they clearly display H I emission lines in the optical. Their continua contain many strong stellar absorption lines, with the most common features due to Ca I, Ca II, Fe I, Na I, Mg I, in addition to prominent absorption bands of TiO, VO, ZrO, CN and CO. We report a homogeneous set of equivalent width (EW) measurements for 45 indices, from optical to NIR species for the 16 star-forming galaxies as well as for 19 early type galaxies where we collected the data from the literature. This selected set of emission and absorption-feature measurements can be used to test predictions of the forthcoming generations of stellar population models. We find correlations among the different absorption features and propose here correlations between optical and NIR indices, as well as among different NIR indices, and compare them with model predictions. While for the optical absorption features the models consistently agree with the observations,the NIR indices are much harder to interpret. For early-type spirals the measurements agree roughly with the models, while for star-forming objects they fail to predict the strengths of these indices., Comment: accepted for publication in MNRAS

[KKS2000]04 (NGC1052-DF2) has become a controversial and well-studied galaxy after the claims suggesting a lack of dark matter and the presence of an anomalously bright globular cluster (GC) system around it. A precise determination of its overall star formation history (SFH) as well as a better characterisation of its GC or planetary nebulae (PN) systems are crucial aspects to: i) understand its real nature, in particular placing it within the family of ultra diffuse galaxies; ii) shed light on its possible formation, evolution, and survival in the absence of dark matter. With this purpose we expand on the knowledge of [KKS2000]04 from the analysis of OSIRIS@GTC spectroscopic data. On the one hand, we claim the possible detection of two new PNe and confirm membership of 5 GCs. On the other hand, we find that the stars shaping [KKS2000]04 are intermediate-age to old (90\% of its stellar mass older than 5 Gyr, average age of 8.7 $\pm$ 0.7 Gyr) and metal-poor ([M/H] $\sim$ -1.18 $\pm$ 0.05), in general agreement with previous results. We do not find any clear hints of significant changes in its stellar content with radius. In addition, the possibility of [KKS2000]04 being a tidal dwarf galaxy with no dark matter is highly disfavoured., Comment: 10 pages, 4 figures, accepted for publication in MNRAS

The X-Shooter Spectral Library (XSL) contains more than 800 spectra of stars across the color-magnitude diagram, that extend from near-UV to near-IR wavelengths (320-2450 nm). We summarize properties of the spectra of O-rich Long Period Variables in XSL, such as phase-related features, and we confront the data with synthetic spectra based on static and dynamical stellar atmosphere models. We discuss successes and remaining discrepancies, keeping in mind the applications to population synthesis modeling that XSL is designed for., Comment: (4 pages, 2 figures, in press in the proceedings of IAU Symp. 343, Why Galaxies Care about AGB Stars (Aug. 2018). Ed. F. Kerschbaum, M. Groenewegen, H. Olofsson.)

The stellar initial mass function (IMF) is commonly assumed to be an invariant probability density distribution function of initial stellar masses being represented by the canonical IMF. As a consequence the galaxy-wide IMF (gwIMF), defined as the sum of the IMFs of all star forming regions, should also be invariant. Recent observational and theoretical results challenge the hypothesis that the gwIMF is invariant. In order to study the possible reasons for this variation we use the IMF determined in resolved star clusters and apply the IGIMF-theory to calculate a grid of gwIMF models for metallicities, -3<[Fe/H]<1, and galaxy-wide star formation rates, $10^{-5}$ 1\,M_\odot$/yr, which is a common condition in the early Universe, we find that the gwIMF is top-heavy (more massive stars), when compared to the canonical IMF. For a SFR $< 1\,\mathrm{M_{\odot}/yr}$ the gwIMF becomes top-light regardless of the metallicity. For metallicities $\mathrm{[Fe/H]} > 0$ the gwIMF can become bottom-heavy regardless of the SFR. The IGIMF models predict that massive elliptical galaxies should have formed with a gwIMF that is top-heavy within the first few hundred Myr of the galaxy's life and that it evolves into a bottom-heavy gwIMF in the metal-enriched galactic center. We study the SFR$-$H$\alpha$ relation, its dependency on metallicity and the SFR, the correction factors to the Kennicutt SFR$_{\rm K}-$H$\alpha$ relation, and provide new fitting functions Late-type dwarf galaxies show significantly higher SFRs with respect to Kennicutt SFRs, while star forming massive galaxies have significantly lower SFRs than hitherto thought. This has implications for the gas-consumption time scales and for the main sequence of galaxies. The Leo P and ultra-faint dwarf galaxies are discussed explicitly. [abridged], Comment: Astronomy and Astrophysics (A&A) in press. 15 pages, 8 figures

In this paper, we present a new derivation of the shape and evolution of the integrated galaxy-wide initial mass function (IGIMF), incorporating explicitly the effects of cosmic rays (CRs) as regulators of the chemical and thermal state of the gas in the dense cores of molecular clouds. We predict the shape of the IGIMF as a function of star formation rate (SFR) and CR density, and show that it can be significantly different with respect to local estimates. In particular, we focus on the physical conditions corresponding to IGIMF shapes that are simultaneously shallower at high-mass end and steeper at the low-mass end than a Kroupa IMF. These solutions can explain both the levels of $\alpha$-enrichment and the excess of low-mass stars as a function of stellar mass, observed for local spheroidal galaxies. As a preliminary test of our scenario, we use idealized star formation histories to estimate the mean IMF shape for galaxies of different $z=0$ stellar mass. We show that the fraction of low-mass stars as a function of galaxy stellar mass predicted by these mean IMFs agrees with the values derived from high-resolution spectroscopic surveys., Comment: 8 pages, 5 figures, MNRAS accepted

The claimed detection of a diffuse galaxy lacking dark matter represents a possible challenge to our understanding of the properties of these galaxies and galaxy formation in general. The galaxy, already identified in photographic plates taken in the summer of 1976 at the UK 48-in Schmidt telescope, presents normal distance-independent properties (e.g. colour, velocity dispersion of its globular clusters). However, distance-dependent quantities are at odds with those of other similar galaxies, namely the luminosity function and sizes of its globular clusters, mass-to-light ratio and dark matter content. Here we carry out a careful analysis of all extant data and show that they consistently indicate a much shorter distance (13 Mpc) than previously indicated (20 Mpc). With this revised distance, the galaxy appears to be a rather ordinary low surface brightness galaxy (R_e=1.4+-0.1 kpc; M*=6.0+-3.6x10^7 Msun) with plenty of room for dark matter (the fraction of dark matter inside the half mass radius is >75% and M_halo/M*>20) corresponding to a minimum halo mass >10^9 Msun. At 13 Mpc, the luminosity and structural properties of the globular clusters around the object are the same as those found in other galaxies., Comment: MNRAS, in press. This version includes an extra new independent distance indicator confirming the 13 Mpc distance. New tests also reject the possibility that blends are causing a misidentification of the location of the TRGB

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

Type-III S0 galaxies present tight scaling relations between their surface brightness photometric and structural parameters. Several evolutionary models have been proposed for the formation of Type-III S0 galaxies but the observations of are usually limited to the local Universe. We study the evolution of the photometric and structural scaling relations found between the parameters of the surface brightness profiles in the rest-frame R-band of Type-III S0 galaxies with z and the possible differences between the rest-frame (B-R) colours of the inner and outer disc profiles. We make use of a sample of 14 Type-III E/S0--S0 galaxies at 0.2

In this second paper of T$h$e role of $E$nvironment in shaping $L$ow-mass $E$arly-type $N$earby g$a$laxies (hELENa) series we study [Mg/Fe] abundance distribution trends of early-type galaxies observed with the SAURON integral field unit, spanning a wide range in mass and local environment densities: 20 low-mass early-types (dEs) of Sybilska et al. (2017) and 258 massive early types (ETGs) of the $ATLAS^{3D}$ project, all homogeneously reduced and analyzed. We show that the [Mg/Fe] ratios scale with velocity dispersion ($\sigma$) at fixed [Fe/H] and that they evolve with [Fe/H] along similar paths for all early-types, grouped in bins of increasing local and global {$\sigma$}, as well as the second velocity moment $V_{rms}$, indicating a common inside-out formation pattern. We then place our dEs on the [Mg/Fe] $vs.$ [Fe/H] diagram of Local Group galaxies and show that dEs occupy the same region and show a similar trend line slope in the diagram as the high-metallicity stars of the Milky Way and the Large Magellanic Cloud. This finding extends the similar trend found for dwarf spheroidal $vs.$ dwarf irregular galaxies and supports the notion that dEs have evolved from late-type galaxies that have lost their gas at a point of their evolution, which likely coincided with them entering denser environments., Comment: Accepted for publication in MNRAS

Cosmic chronometers may be used to measure the age difference between passively evolving galaxy populations to calculate the Hubble parameter H(z). The age estimator emerges from the relationship between the amplitude of the rest frame Balmer break at 4000 angstroms and the age of a galaxy, assuming that there is one single stellar population within each galaxy. However, recent literature has shown possible contamination (up to 2.4% of the stellar mass in a high redshift sample) of a young component embedded within the predominantly old population of the quiescent galaxy. We compared the data with the predictions of each model, using a new approach of distinguishing between systematic and statistical errors (in previous works, these had incorrectly been added in quadrature) and evaluating the effects of contamination by a young stellar component. The ages inferred using cosmic chronometers represent a galaxy-wide average rather than a characteristic of the oldest population alone. The average contribution from the young component to the rest luminosity at 4000 angstroms may constitute a third of the luminosity in some samples, which means that this is far from negligible. This ratio is significantly dependent on stellar mass, proportional to M^{-0.7}. Consequently, the measurements of the absolute value of the age or the differential age between different redshifts are incorrect and make the previous calculations of H(z) very inaccurate. Some cosmological models, such as the Einstein-de Sitter model or quasi-steady state cosmology, which are rejected under the assumption of a purely old population, can be made compatible with the predicted ages of the Universe as a function of redshift if we take this contamination into account. However, the static Universe models are rejected by these H(z) measurements, even when this contamination is taken into account., Comment: 12 pages, accepted to be published in A&A

In this paper we investigate whether the stellar initial mass function of early-type galaxies depends on their host environment. To this purpose, we have selected a sample of early-type galaxies from the SPIDER catalogue, characterized their environment through the group catalogue of Wang et al. and used their optical SDSS spectra to constrain the IMF slope, through the analysis of IMF-sensitive spectral indices. To reach a high enough signal-to-noise ratio, we have stacked spectra in velocity dispersion ($\sigma_0$) bins, on top of separating the sample by galaxy hierarchy and host halo mass, as proxies for galaxy environment. In order to constrain the IMF, we have compared observed line strengths to predictions of MIUSCAT/EMILES synthetic stellar population models, with varying age, metallicity, and "bimodal" (low-mass tapered) IMF slope ($\rm \Gamma_b$). Consistent with previous studies, we find that $\rm \Gamma_b$ increases with $\sigma_0$, becoming bottom-heavy (i.e. an excess of low-mass stars with respect to the Milky-Way-like IMF) at high $\sigma_0$. We find that this result is robust against the set of isochrones used in the stellar population models, as well as the way the effect of elemental abundance ratios is taken into account. We thus conclude that it is possible to use currently state-of-the-art stellar population models and intermediate resolution spectra to consistently probe IMF variations. For the first time, we show that there is no dependence of $\Gamma_b$ on environment or galaxy hierarchy, as measured within the $3"$ SDSS fibre, thus leaving the IMF as an intrinsic galaxy property, possibly set already at high redshift.

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., Comment: Accepted for publication in MNRAS