Showing 1–50 of 77 results for author: Martín-Navarro, I

Deriving the star formation histories of galaxies from spectra with simulation-based inference

Authors: Patricia Iglesias-Navarro, Marc Huertas-Company, Ignacio Martín-Navarro, Johan H. Knapen, Emilie Pernet

Abstract: High-resolution galaxy spectra encode information about the stellar populations within galaxies. The properties of the stars, such as their ages, masses, and metallicities, provide insights into the underlying physical processes that drive the growth and transformation of galaxies over cosmic time. We explore a simulation-based inference (SBI) workflow to infer from optical absorption spectra th… ▽ More High-resolution galaxy spectra encode information about the stellar populations within galaxies. The properties of the stars, such as their ages, masses, and metallicities, provide insights into the underlying physical processes that drive the growth and transformation of galaxies over cosmic time. We explore a simulation-based inference (SBI) workflow to infer from optical absorption spectra the posterior distributions of metallicities and the star formation histories (SFHs) of galaxies (i.e. the star formation rate as a function of time). We generated a dataset of synthetic spectra to train and test our model using the spectroscopic predictions of the MILES stellar population library and non-parametric SFHs. We reliably estimate the mass assembly of an integrated stellar population with well-calibrated uncertainties. Specifically, we reach a score of $0.97\,R^2$ for the time at which a given galaxy from the test set formed $50\%$ of its stellar mass, obtaining samples of the posteriors in only $10^{-4}$\,s. We then applied the pipeline to real observations of massive elliptical galaxies, recovering the well-known relationship between the age and the velocity dispersion, and show that the most massive galaxies ($σ\sim300$ km/s) built up to 90\% of their total stellar masses within $1$\,Gyr of the Big Bang. The inferred properties also agree with the state-of-the-art inversion codes, but the inference is performed up to five orders of magnitude faster. This SBI approach coupled with machine learning and applied to full spectral fitting makes it possible to address large numbers of galaxies while performing a thick sampling of the posteriors. It will allow both the deterministic trends and the inherent uncertainties of the highly degenerated inversion problem to be estimated for large and complex upcoming spectroscopic surveys, such as DESI, WEAVE, or 4MOST. △ Less

Submitted 28 June, 2024; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 12 pages, 9 figures, 1 table. Accepted for publication on A&A

Unveiling the (in)consistencies among the galaxy stellar mass function, star formation histories, satellite abundances and intracluster light from a semi-empirical perspective

Authors: Hao Fu, Francesco Shankar, Mohammadreza Ayromlou, Ioanna Koutsouridou, Andrea Cattaneo, Caroline Bertemes, Sabine Bellstedt, Ignacio Martín-Navarro, Joel Leja, Viola Allevato, Mariangela Bernardi, Lumen Boco, Paola Dimauro, Carlotta Gruppioni, Andrea Lapi, Nicola Menci, Iván Muñoz Rodríguez, Annagrazia Puglisi, Alba V. Alonso-Tetilla

Abstract: In a hierarchical, dark matter-dominated Universe, stellar mass functions (SMFs), galaxy merger rates, star formation histories (SFHs), satellite abundances, and intracluster light, should all be intimately connected observables. However, the systematics affecting observations still prevent universal and uniform measurements of, for example, the SMF and the SFHs, inevitably preventing theoretical… ▽ More In a hierarchical, dark matter-dominated Universe, stellar mass functions (SMFs), galaxy merger rates, star formation histories (SFHs), satellite abundances, and intracluster light, should all be intimately connected observables. However, the systematics affecting observations still prevent universal and uniform measurements of, for example, the SMF and the SFHs, inevitably preventing theoretical models to compare with multiple data sets robustly and simultaneously. We here present our holistic semi-empirical model DECODE (Discrete statistical sEmi-empiriCal mODEl) that converts via abundance matching dark matter merger trees into galaxy assembly histories, using different SMFs in input and predicting all other observables in output in a fully data-driven and self-consistent fashion with minimal assumptions. We find that: 1) weakly evolving or nearly constant SMFs below the knee ($M_\star \lesssim 10^{11} \, M_\odot$) are the best suited to generate star formation histories aligned with those inferred from MaNGA, SDSS, GAMA, and, more recently, JWST; 2) the evolution of satellites after infall only affects the satellite abundances and star formation histories of massive central galaxies but not their merger histories; 3) the resulting SFR-$M_\star$ relation is lower in normalization by a factor of $\sim 2$ with respect to observations, with a flattening at high masses more pronounced in the presence of mergers; 4) the latest data on intracluster light can be reproduced if mass loss from mergers is included in the models. Our findings are pivotal in acting as pathfinder to test the self-consistency of the high-quality data from, e.g., JWST and Euclid. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 21 pages, 12 figures, MNRAS accepted

[Mg/Fe] and variable initial mass function: Revision of [$α$/Fe] for massive galaxies

Authors: Emilie Pernet, Alina Böcker, Ignacio Martín-Navarro

Abstract: Observations of nearby massive galaxies have revealed that they are older and richer in metals and magnesium than their low-mass counterparts. In particular, the overabundance of magnesium compared to iron, [Mg/Fe], is interpreted to reflect the short star formation history that the current massive galaxies underwent early in the Universe. We present a systematic revision of the [Mg/Fe] - velocity… ▽ More Observations of nearby massive galaxies have revealed that they are older and richer in metals and magnesium than their low-mass counterparts. In particular, the overabundance of magnesium compared to iron, [Mg/Fe], is interpreted to reflect the short star formation history that the current massive galaxies underwent early in the Universe. We present a systematic revision of the [Mg/Fe] - velocity dispersion ($σ$) relation based on stacked spectra of early-type galaxies with a high signal-to-noise ratio from the Sloan Digital Sky Survey (SDSS). Using the penalized pixel-fitting (pPXF) method of Cappellari & Emsellem 2004 and the Vazdekis et al. 2015 MILES single stellar population (SSP) models, we fit a wide optical wavelength range to measure the net $α$-abundance. The combination of pPXF and $α$-enhanced MILES models incorrectly leads to an apparently decreasing trend of [$α$/Fe] with velocity dispersion. We interpret this result as a consequence of variations in the individual abundances of the different $α$-elements. This warrants caution for a naïve use of full spectral fitting algorithms paired with stellar population models that do not take individual elemental abundance variations into account, especially when deriving averaged quantities such as the mean [$α$/Fe] of a stellar population. In addition, and based on line-strength measurements, we quantify the impact of a non-universal initial mass function on the recovered abundance pattern of galaxies. In particular, we find that a simultaneous fit of the slope of the initial mass function and the [Mg/Fe] results in a shallower [Mg/Fe]-$σ$ relation. Therefore, our results suggest that star formation in massive galaxies lasted longer than what has been reported previously, although it still occurred significantly faster than in the solar neighbourhood. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 7 pages, 5 figures, accepted for publication in A&A

Feedback mechanisms stopping the star formation in a pair of massive galaxies in the early Universe

Authors: Pablo G. Pérez-González, Francesco D`Eugenio, Bruno Rodríguez del Pino, Hannah Übler, Roberto Maiolino, Santiago Arribas, Giovanni Cresci, Isabella Lamperti, Andrew J. Bunker, Stefano Carniani, Stephane Charlot, Christopher J. Willott, Torsten Böker, Eleonora Parlanti, Jan Scholtz, Giacomo Venturi, Guillermo Barro, Luca Costantin, Ignacio Martín-Navarro, James S. Dunlop, Daniel Magee

Abstract: Feedback is the key physical mechanism regulating galaxy formation. Stars in galaxies form when baryons radiatively cool down and fall into gravitational wells. Eventually, star formation quenches as gas is depleted and/or perturbed by feedback processes, no longer being able to collapse and condense. For massive galaxies, astronomers identify feedback from accreting supermassive black holes (acti… ▽ More Feedback is the key physical mechanism regulating galaxy formation. Stars in galaxies form when baryons radiatively cool down and fall into gravitational wells. Eventually, star formation quenches as gas is depleted and/or perturbed by feedback processes, no longer being able to collapse and condense. For massive galaxies, astronomers identify feedback from accreting supermassive black holes (active galactic nuclei, AGN) as the main agent responsible for quenching. We report the first spatially resolved spectroscopic observations of a massive, completely quiescent galaxy at $z=3.714$ (Jekyll) and its neighborhood. Jekyll is part of a galaxy pair with a compact, dusty, massive star-forming companion (Hyde). We find large amounts of ionized and neutral gas in the intergalactic medium around the pair, yet Jekyll has remained quiescent for more than 500~Myr. The emitting gas is consistent with AGN photoionization, but no AGN is observed in Jekyll. We find that, in contrast to standard scenarios, AGN in satellite galaxies can be critical contributors for keeping massive galaxies quiescent in the early Universe. After the accelerated formation and quenching of the massive central galaxy, tidally induced gas stripping additionally contributes to the star-formation regulation on subsequent satellite galaxy generations. △ Less

Submitted 8 May, 2024; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted. Comments welcome

Observational constraints on the stellar recycled gas in active galactic nuclei feeding

Authors: Rogerio Riffel, Luis G. Dahmer-Hahn, Alexandre Vazdekis, Richard Davies, David Rosario, Cristina Ramos Almeida, Anelise Audibert, Ignacio Martin-Navarro, Lucimara Pires Martins, Alberto Rodriguez-Ardila, Rogemar A. Riffel, Thaisa Storchi-Bergmann, Michele Bertoldo-Coelho, Marina Trevisan, Erin Hicks, Allan Schnorr Muller, Lais Nery Marinho, Sylvain Veilleux

Abstract: 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… ▽ More 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, $<t>_L\lesssim$ 4.5 Gyr) when compared with their analogues ($<t>_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 ($ <t>_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. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS

Baryonic properties of nearby galaxies across the stellar-to-total dynamical mass relation

Authors: Laura Scholz-Diaz, Ignacio Martin-Navarro, Jesus Falcon-Barroso, Mariya Lyubenova, Glenn van de Ven

Abstract: In the standard cosmological model, the assembly of galaxies is primarily driven by the growth of their host dark matter halos. At the center of these halos, however, baryonic processes take over, leading to the plethora of observed galaxy properties. The coupling between baryonic and dark matter physics is central to our understanding of galaxies and yet, it remains a challenge for theoretical mo… ▽ More In the standard cosmological model, the assembly of galaxies is primarily driven by the growth of their host dark matter halos. At the center of these halos, however, baryonic processes take over, leading to the plethora of observed galaxy properties. The coupling between baryonic and dark matter physics is central to our understanding of galaxies and yet, it remains a challenge for theoretical models and observations. Here, we demonstrate that measured ages, metallicities, stellar angular momentum, morphology and star formation rates, correlate with both stellar and halo mass. Using dynamical modeling, we find that at fixed stellar mass, CALIFA galaxies become younger, more metal-poor and rotationally supported, have higher star formation rates and later-type morphologies as their total mass increases, with independent stellar and total masses measurements. These results indicate that the formation of galaxies and thus their baryonic properties do not vary with stellar mass alone, with halo mass also playing an important role. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 14 pages, 5 Figures. Accepted in Nature Astronomy

INSPIRE: INvestigating Stellar Population In RElics VI -- The low-mass end slope of the stellar Initial Mass Function and chemical composition

Authors: Michalina Maksymowicz-Maciata, Chiara Spiniello, Ignacio Martín-Navarro, Anna Ferré-Mateu, Davide Bevacqua, Michele Cappellari, Giuseppe D'Ago, Crescenzo Tortora, Magda Arnaboldi, Johanna Hartke, Paolo Saracco, Diana Scognamiglio

Abstract: The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.1<z<0.4 and obtained their star formation histories (SFHs). Due to their preserved very old stellar populations, relics are the perfect systems to constrain the earliest epochs of mass assembly in the Universe and the formation of massive early-type galaxies. The goal of this work is to investigate whet… ▽ More The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.1<z<0.4 and obtained their star formation histories (SFHs). Due to their preserved very old stellar populations, relics are the perfect systems to constrain the earliest epochs of mass assembly in the Universe and the formation of massive early-type galaxies. The goal of this work is to investigate whether a correlation exists between the degree of relicness (DoR), quantifying the fraction of stellar mass formed at z>2, and the other stellar population parameters.We use the Full-Index-Fitting method to fit the INSPIRE spectra to single stellar population (SSP) models. This allows us to measure, for the first time, the low-mass end slope of the IMF, as well as stellar metallicity [M/H], [Mg/Fe], [Ti/Fe] and [Na/Fe] ratios, and study correlations between them and the DoR. Similarly to normal-sized galaxies, UCMGs with larger stellar masses have overall higher metallicities. We found a correlation between the low-mass end of the IMF slope and the DoR, that, however, breaks down for systems with a more extended SFH. An even stronger dependency is found between the IMF and the fraction of mass formed at high-z. At equal velocity dispersion and metallicity, galaxies with a higher DoR have a dwarf-richer IMF than that of low-DoR counterparts. This might indicate that the cosmic epoch and formation mechanisms influence the fragmentation of the star formation cloud and hence might be the explanation for IMF variations detected in massive ETGs. △ Less

Submitted 28 January, 2024; originally announced January 2024.

Comments: 16 pages, 11 figures, submitted to MNRAS

The universal variability of the stellar initial mass function probed by the TIMER survey

Authors: Ignacio Martín-Navarro, Adriana de Lorenzo-Cáceres, Dimitri A. Gadotti, Jairo Méndez-Abreu, Jesús Falcón-Barroso, Patricia Sánchez-Blázquez, Paula Coelho, Justus Neumann, Glenn van de Ven, Isabel Pérez

Abstract: The debate about the universality of the stellar initial mass function (IMF) revolves around two competing lines of evidence. While measurements in the Milky Way, an archetypal spiral galaxy, seem to support an invariant IMF, the observed properties of massive early-type galaxies (ETGs) favor an IMF somehow sensitive to the local star formation conditions. The fundamental methodological and physic… ▽ More The debate about the universality of the stellar initial mass function (IMF) revolves around two competing lines of evidence. While measurements in the Milky Way, an archetypal spiral galaxy, seem to support an invariant IMF, the observed properties of massive early-type galaxies (ETGs) favor an IMF somehow sensitive to the local star formation conditions. The fundamental methodological and physical differences between both approaches have hampered, however, a comprehensive understanding of IMF variations. We describe here an improved modelling scheme that allows for the first time consistent IMF measurements across stellar populations with different ages and complex star formation histories. Making use of the exquisite MUSE optical data from the TIMER survey and powered by the MILES stellar population models, we show the age, metallicity, [Mg/Fe], and IMF slope maps of the inner regions of NGC 3351, a spiral galaxy with a mass similar to that of the Milky Way. The measured IMF values in NGC3351 follow the expectations from a Milky Way-like IMF, although they simultaneously show systematic and spatially coherent variations, particularly for low-mass stars. In addition, our stellar population analysis reveals the presence of metal-poor and Mg-enhanced star-forming regions that appear to be predominantly enriched by the stellar ejecta of core-collapse supernovae. Our findings showcase therefore the potential of detailed studies of young stellar populations to better understand the early stages of galaxy evolution and, in particular, the origin of the observed IMF variations beyond and within the Milky Way. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 16 pages, 11 figures. Accepted for publication in Astronomy and Astrophysics

Initial mass function variability from the integrated light of diverse stellar systems

Authors: Chloe M. Cheng, Alexa Villaume, Michael L. Balogh, Jean P. Brodie, Ignacio Martín-Navarro, Aaron J. Romanowsky, Pieter G. van Dokkum

Abstract: We present a uniform analysis of the stellar initial mass function (IMF) from integrated light spectroscopy of 15 compact stellar systems (11 globular clusters in M31 and 4 ultra compact dwarfs in the Virgo cluster, UCDs) and two brightest Coma cluster galaxies (BCGs), covering a wide range of metallicities ($-$1.7 $<$ [Fe/H] $<$ 0.01) and velocity dispersions (7.4 km~s$^{-1}$ $< σ<$ 275 km~s… ▽ More We present a uniform analysis of the stellar initial mass function (IMF) from integrated light spectroscopy of 15 compact stellar systems (11 globular clusters in M31 and 4 ultra compact dwarfs in the Virgo cluster, UCDs) and two brightest Coma cluster galaxies (BCGs), covering a wide range of metallicities ($-$1.7 $<$ [Fe/H] $<$ 0.01) and velocity dispersions (7.4 km~s$^{-1}$ $< σ<$ 275 km~s$^{-1}$). The S/N $\sim 100$ Å$^{-1}$ Keck LRIS spectra are fitted over the range $4000<λ/\mboxÅ<10,000$ with flexible, full-spectrum stellar population synthesis models. We use the models to fit simultaneously for ages, metallicities, and individual elemental abundances of the population, allowing us to decouple abundance variations from variations in IMF slope. We show that compact stellar systems do not follow the same trends with physical parameters that have been found for early-type galaxies. Most globular clusters in our sample have an IMF consistent with that of the Milky Way, over a wide range of [Fe/H] and [Mg/Fe]. There is more diversity among the UCDs, with some showing evidence for a bottom-heavy IMF, but with no clear correlation with metallicity, abundance, or velocity dispersion. The two Coma BCGs have similar velocity dispersion and metallicity, but we find the IMF of NGC~4874 is consistent with that of the Milky Way while NGC~4889 presents evidence for a significantly bottom-heavy IMF. For this sample, the IMF appears to vary between objects in a way that is not explained by a single metallicity-dependent prescription. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Report number: MN-23-2545-MJ

INSPIRE: INvestigating Stellar Population In RElics V. A catalogue of ultra-compact massive galaxies outside the local Universe and their degree of relicness

Authors: C. Spiniello, G. D'Ago, L. Coccato, J. Hartke, C. Tortora, A. Ferré-Mateu, C. Pulsoni, M. Cappellari, M. Maksymowicz-Maciata, M. Arnaboldi, D. Bevacqua, A. Gallazzi, L. K. Hunt, F. La Barbera, I. Martín-Navarro, N. R. Napolitano, M. Radovich, P. Saracco, D. Scognamiglio, M. Spavone, S. Zibetti

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

Submitted 14 November, 2023; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: accepted for publication on MNRAS, 20 pages, 15 figures, 3 tables

Disc galaxies are still settling: The discovery of the smallest nuclear discs and their young stellar bars

Authors: Camila de Sá-Freitas, Dimitri A. Gadotti, Francesca Fragkoudi, Lodovico Coccato, Paula Coelho, Adriana de Lorenzo-Cáceres, Jesús Falcón-Barroso, Tutku Kolcu, Ignacio Martín-Navarro, Jairo Mendez-Abreu, Justus Neumann, Patricia Sanchez Blazquez, Miguel Querejeta, Glenn van de Ven

Abstract: When galactic discs settle and become massive enough, they are able to form stellar bars. These non-axisymmetric structures induce shocks in the gas, causing it to flow to the centre where nuclear structures, such as nuclear discs and rings, are formed. Previous theoretical and observational studies have hinted at the co-evolution of bars and nuclear discs, suggesting that nuclear discs grow "insi… ▽ More When galactic discs settle and become massive enough, they are able to form stellar bars. These non-axisymmetric structures induce shocks in the gas, causing it to flow to the centre where nuclear structures, such as nuclear discs and rings, are formed. Previous theoretical and observational studies have hinted at the co-evolution of bars and nuclear discs, suggesting that nuclear discs grow "inside-out", thereby proposing that smaller discs live in younger bars. Nevertheless, it remains unclear how the bar and the nuclear structures form and evolve with time. The smallest nuclear discs discovered to date tend to be larger than $\sim200~\rm{pc}$, even though some theoretical studies find that when nuclear discs form they can be much smaller. Using MUSE archival data, we report for the first time two extragalactic nuclear discs with radius sizes below $100~\rm{pc}$. Additionally, our estimations reveal the youngest bars found to date. We estimate that the bars in these galaxies formed $4.50^{+1.60}_{-1.10}\rm{(sys)}^{+1.00}_{-0.75}\rm{(stat)}$ and $0.7^{+2.60}\rm{(sys)}^{+0.05}_{-0.05}\rm{(stat)}~\rm{Gyr}$ ago, for NGC\,289 and NGC\,1566, respectively. This suggests that at least some disc galaxies in the Local Universe may still be dynamically settling. By adding these results to previous findings in the literature, we retrieve a stronger correlation between nuclear disc size and bar length and we derive a tentative exponential growth scenario for nuclear discs. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: Accepted in A&A (in press), 13 pages, 5 figures, 2 tabels

Journal ref: A&A 678, A202 (2023)

Mapping the stellar population and gas excitation of MaNGA galaxies with MEGACUBES. Results for AGN versus control sample

Authors: Rogério Riffel, Nicolas D. Mallmann, Sandro B. Rembold, Gabriele S. Ilha, Rogemar A. Riffel, Thaisa Storchi-Bergmann, Daniel Ruschel-Dutra, Alexandre Vazdekis, Ignacio Martín-Navarro, Jaderson S. Schimoia, Cristina Ramos Almeida, Luiz N. da Costa, Glauber C. Vila-Verde, Lara Gatto

Abstract: 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… ▽ More 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$α$ equivalent width for proper activity classification. We present a composite ``BPT+WHAN" diagram that produces a more comprehensive mapping of the gas excitation. △ Less

Submitted 4 August, 2023; v1 submitted 21 July, 2023; originally announced July 2023.

Comments: Accepted for publication in MNRAS; small typos corrected

Effect of the initial mass function on the dynamical SMBH mass estimate in the nucleated early-type galaxy FCC 47

Authors: Sabine Thater, Mariya Lyubenova, Katja Fahrion, Ignacio Martín-Navarro, Prashin Jethwa, Dieu D. Nguyen, Glenn van de Ven

Abstract: Supermassive black holes (SMBHs) and nuclear star clusters (NSCs) co-exist in many galaxies. While the formation history of the black hole is essentially lost, NSCs preserve their evolutionary history imprinted onto their stellar populations and kinematics. Studying SMBHs and NSCs in tandem might help us to ultimately reveal the build-up of galaxy centres. In this study, we combine large-scale VLT… ▽ More Supermassive black holes (SMBHs) and nuclear star clusters (NSCs) co-exist in many galaxies. While the formation history of the black hole is essentially lost, NSCs preserve their evolutionary history imprinted onto their stellar populations and kinematics. Studying SMBHs and NSCs in tandem might help us to ultimately reveal the build-up of galaxy centres. In this study, we combine large-scale VLT/MUSE and high-resolution adaptive-optics-assisted VLT/SINFONI observations of the early-type galaxy FCC 47 with the goal being to assess the effect of a spatially (non-)variable initial mass function (IMF) on the determination of the mass of the putative SMBH in this galaxy. We achieve this by performing DYNAMITE Schwarzschild orbit-superposition modelling of the galaxy and its NSC. In order to properly take account of the stellar mass contribution to the galaxy potential, we create mass maps using a varying stellar mass-to-light ratio derived from single stellar population models with fixed and with spatially varying IMFs. Using the two mass maps, we estimate black hole masses of $(7.1^{+0.8}_{-1.1})\times 10^7\,M_{\odot}$ and $(4.4^{+1.2}_{-2.1}) \times 10^7\,M_{\odot}$ at $3σ$ signifance, respectively. Compared to models with constant stellar-mass-to-light ratio, the black hole masses decrease by 15% and 48%, respectively. Therefore, a varying IMF, both in its functional form and spatially across the galaxy, has a non-negligible effect on the SMBH mass estimate. Furthermore, we find that the SMBH in FCC 47 has probably not grown over-massive compared to its very over-massive NSC. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 23 pages 19 Figures, accepted for publication in A&A

Journal ref: A&A 675, A18 (2023)

INSPIRE: INvestigating Stellar Population In RElics IV. The Initial Mass Function slope in relics

Authors: I. Martin-Navarro, C. Spiniello, C. Tortora, L. Coccato, G. D'Ago, A. Ferre-Mateu, C. Pulsoni, J. Hartke, M. Arnaboldi, L. Hunt, N. R. Napolitano, D. Scognamiglio, M. Spavone

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

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS. 8 pages, 5 figures

INSPIRE: INvestigating Stellar Population In RElics III. Second data release (DR2): testing the systematics on the stellar velocity dispersion

Authors: G. D'Ago, C. Spiniello, L. Coccato, C. Tortora, F. La Barbera, M. Arnaboldi, D. Bevacqua, A. Ferré-Mateu, A. Gallazzi, J. Hartke, L. K. Hunt, I. Martín-Navarro, N. R. Napolitano, C. Pulsoni, M. Radovich, P. Saracco, D. Scognamiglio, S. Zibetti

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

Submitted 10 February, 2023; originally announced February 2023.

Comments: Accepted for publication, and in press on A&A

Journal ref: A&A 672, A17 (2023)

The Fornax3D project: Environmental effects on the assembly of dynamically cold disks in Fornax cluster galaxies

Authors: Y. Ding, L. Zhu, G. van de Ven, L. Coccato, E. M. Corsini, L. Costantin, K. Fahrion, J. Falcón-Barroso, D. A. Gadotti, E. Iodice, M. Lyubenova, I. Martín-Navarro, R. M. McDermid, F. Pinna, M. Sarzi

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

Submitted 15 February, 2023; v1 submitted 13 January, 2023; originally announced January 2023.

Comments: 31 pages, 32 figures, Accepted to A&A

Journal ref: A&A 672, A84 (2023)

The dark side of galaxy stellar populations II: The dependence of star formation histories on halo mass and on the scatter of the main sequence

Authors: Laura Scholz-Diaz, Ignacio Martin-Navarro, Jesus Falcon-Barroso

Abstract: Nearby galaxies are the end result of their cosmological evolution, which is predicted to be influenced by the growth of their host dark matter halos. This co-evolution potentially leaves signatures in present-day observed galaxy properties, which might be essential to further understand how the growth and properties of galaxies are connected to those of their host halos. In this work, we study th… ▽ More Nearby galaxies are the end result of their cosmological evolution, which is predicted to be influenced by the growth of their host dark matter halos. This co-evolution potentially leaves signatures in present-day observed galaxy properties, which might be essential to further understand how the growth and properties of galaxies are connected to those of their host halos. In this work, we study the evolutionary histories of nearby galaxies both in terms of their host halos and the scatter of the star-forming main sequence by investigating their time-resolved stellar populations using absorption optical spectra drawn from the Sloan Digital Sky Survey. We find that galaxy star formation histories depend on the masses of their host halos, and hence they shape the evolution of the star-forming main sequence over cosmic time. Additionally, we also find that the scatter around the z=0 star-forming main sequence is not (entirely) stochastic, as galaxies with currently different star formation rates have experienced, on average, different star formation histories. Our findings suggest that dark matter halos might play a key role in modulating the evolution of star formation in galaxies, and thus of the main sequence, and further demonstrate that galaxies at different evolutionary stages contribute to the observed scatter of this relation. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS, 13 pages, 9 Figures

A new method for age-dating the formation of bars in disc galaxies: The TIMER view on NGC1433's old bar and the inside-out growth of its nuclear disc

Authors: Camila de Sá-Freitas, Francesca Fragkoudi, Dimitri A. Gadotti, Jesús Falcón-Barroso, Adrian Bittner, Patricia Sánchez-Blázquez, Glenn van de Ven, Rebekka Bieri, Lodovico Coccato, Paula Coelho, Katja Fahrion, Geraldo Gonçalves, Taehyun Kim, Adriana de Lorenzo-Cáceres, Marie Martig, Ignacio Martín-Navarro, Jairo Mendez-Abreu, Justus Neumann, Miguel Querejeta

Abstract: The epoch in which galactic discs settle is a major benchmark to test models of galaxy formation and evolution but is as yet largely unknown. Once discs settle and become self-gravitating enough, stellar bars are able to form; therefore, determining the ages of bars can shed light on the epoch of disc settling, and on the onset of secular evolution. Nevertheless, until now, timing when the bar for… ▽ More The epoch in which galactic discs settle is a major benchmark to test models of galaxy formation and evolution but is as yet largely unknown. Once discs settle and become self-gravitating enough, stellar bars are able to form; therefore, determining the ages of bars can shed light on the epoch of disc settling, and on the onset of secular evolution. Nevertheless, until now, timing when the bar formed has proven challenging. In this work, we present a new methodology for obtaining the bar age, using the star formation history of nuclear discs. Nuclear discs are rotation-supported structures, built by gas pushed to the centre via bar-induced torques, and their formation is thus coincident with bar formation. In particular, we use integral field spectroscopic (IFS) data from the TIMER survey to disentangle the star formation history of the nuclear disc from that of the underlying main disc, which enables us to more accurately determine when the nuclear disc forms. We demonstrate the methodology on the galaxy NGC 1433 -- which we find to host an old bar that is $8.0^{+1.6}_{-1.1}\rm{(sys)}^{+0.2}_{-0.5}\rm{(stat)}$ Gyr old -- and describe a number of tests carried out on both the observational data and numerical simulations. In addition, we present evidence that the nuclear disc of NGC 1433 grows in accordance with an inside-out formation scenario. This methodology is applicable to high-resolution IFS data of barred galaxies with nuclear discs, making it ideally suited for the TIMER survey sample. In the future we will thus be able to determine the bar age for a large sample of galaxies, shedding light on the epoch of disc settling and bar formation. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted for Publication by A&A, 18 pages, 14 figures

The survival of stellar discs in Fornax-like environments, from TNG50 to real galaxies

Authors: Pablo M. Galán-de Anta, M. Sarzi, A. Pillepich, Y. Ding, L. Zhu, L. Coccato, E. M. Corsini, K. Fahrion, J. Falcón-Barroso, D. A. Gadotti, E. Iodice, M. Lyubenova, I. Martín-Navarro, R. M. McDermid, F. Pinna, G. van de Ven, P. T. de Zeeuw

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

Submitted 2 November, 2022; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: 12 pages, 7 figures. Accepted on 21 of October, published on 25 of October

Low-density star cluster formation: discovery of a young faint fuzzy on the outskirts of the low-mass spiral galaxy NGC 247

Authors: Aaron J. Romanowsky, Søren S. Larsen, Alexa Villaume, Jeffrey L. Carlin, Joachim Janz, David J. Sand, Jay Strader, Jean P. Brodie, Sukanya Chakrabarti, Chloe M. Cheng, Denija Crnojević, Duncan A. Forbes, Christopher T. Garling, Jonathan R. Hargis, Ananthan Karunakaran, Ignacio Martín-Navarro, Knut A. G. Olsen, Nicole Rider, Bitha Salimkumar, Vakini Santhanakrishnan, Kristine Spekkens, Yimeng Tang, Pieter G. van Dokkum, Beth Willman

Abstract: The classical globular clusters found in all galaxy types have half-light radii of $r_{\rm h} \sim$ 2-4 pc, which have been tied to formation in the dense cores of giant molecular clouds. Some old star clusters have larger sizes, and it is unclear if these represent a fundamentally different mode of low-density star cluster formation. We report the discovery of a rare, young "faint fuzzy" star clu… ▽ More The classical globular clusters found in all galaxy types have half-light radii of $r_{\rm h} \sim$ 2-4 pc, which have been tied to formation in the dense cores of giant molecular clouds. Some old star clusters have larger sizes, and it is unclear if these represent a fundamentally different mode of low-density star cluster formation. We report the discovery of a rare, young "faint fuzzy" star cluster, NGC 247-SC1, on the outskirts of the low-mass spiral galaxy NGC 247 in the nearby Sculptor group, and measure its radial velocity using Keck spectroscopy. We use Hubble Space Telescope imaging to measure the cluster half-light radius of $r_{\rm h} \simeq 12$ pc and a luminosity of $L_V \simeq 4\times10^5 \mathrm{L}_\odot$. We produce a colour-magnitude diagram of cluster stars and compare to theoretical isochrones, finding an age of $\simeq$ 300 Myr, a metallicity of [$Z$/H] $\sim -0.6$ and an inferred mass of $M_\star \simeq 9\times10^4 \mathrm{M}_\odot$. The narrow width of blue-loop star magnitudes implies an age spread of $\lesssim$ 50 Myr, while no old red-giant branch stars are found, so SC1 is consistent with hosting a single stellar population, modulo several unexplained bright "red straggler" stars. SC1 appears to be surrounded by tidal debris, at the end of a $\sim$ 2 kpc long stellar filament that also hosts two low-mass, low-density clusters of a similar age. We explore a link between the formation of these unusual clusters and an external perturbation of their host galaxy, illuminating a possible channel by which some clusters are born with large sizes. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 19 pages, 10 figures, Mon. Not. R. Astron. Soc., in press

The Fornax3D project: The assembly history of massive early-type galaxies in the Fornax cluster from deep imaging and integral field spectroscopy

Authors: M. Spavone, E. Iodice, G. D'Ago, G. van de Ven, L. Morelli, E. M. Corsini, M. Sarzi, L. Coccato, K. Fahrion, J. Falcón-Barroso, D. A. Gadotti, M. Lyubenova, I. Martín-Navarro, R. M. McDermid, F. Pinna, A. Pizzella, A. Poci, P. T. de Zeeuw, L. Zhu

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

Submitted 30 May, 2022; originally announced May 2022.

Comments: 23 pages, 7 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 663, A135 (2022)

The Fornax3D project: Intrinsic Correlations between Orbital Properties and the Stellar Initial Mass Function

Authors: A. Poci, R. M. McDermid, M. Lyubenova, I. Martin-Navarro, G. van de Ven, L. Coccato, E. M. Corsini, K. Fahrion, J. Falcon-Barroso, D. A. Gadotti, E. Iodice, F. Pinna, M. Sarzi, P. T. de Zeeuw, L. Zhu

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

Submitted 25 May, 2022; originally announced May 2022.

Comments: 11 pages. Accepted to MNRAS

The Fornax3D project: The environmental impact on gas metallicity gradients in Fornax cluster galaxies

Authors: M. A. Lara-Lopez, P. M. Galan-de Anta, M. Sarzi, E. Iodice, T. A. Davis, N. Zabel, E. M. Corsini, P. T. de Zeeuw, K. Fahrion, J. Falcon-Barroso, D. A. Gadotti, R. M. McDermid, F. Pinna, V. Rodriguez-Gomez, G. van de Ven, L. Zhu, L. Coccato, M. Lyubenova, I. Martin-Navarro

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

Submitted 8 February, 2022; originally announced February 2022.

Comments: Accepted by A&A

Journal ref: A&A 660, A105 (2022)

The dark side of galaxy stellar populations I: The stellar-to-halo mass relation and the velocity dispersion - halo mass relation

Authors: Laura Scholz-Diaz, Ignacio Martin-Navarro, Jesus Falcon-Barroso

Abstract: The growth and properties of galaxies are thought to be closely connected to the ones of their host dark matter halos. Despite the importance of this so-called galaxy-halo connection, the potential role of dark matter halos in regulating observed galaxy properties remains yet to be fully understood. In this work, we derive the ages, metallicites and [Mg/Fe] abundances from optical spectra from the… ▽ More The growth and properties of galaxies are thought to be closely connected to the ones of their host dark matter halos. Despite the importance of this so-called galaxy-halo connection, the potential role of dark matter halos in regulating observed galaxy properties remains yet to be fully understood. In this work, we derive the ages, metallicites and [Mg/Fe] abundances from optical spectra from the Sloan Digital Sky Survey of nearby central galaxies, and study them in terms of their host halos. We investigate how the scatter in the stellar-to-halo mass relation and the velocity dispersion - halo mass relation correlates with these stellar population parameters. In addition, we also study the differences when distinguishing between different galaxy morphologies and environments. We find that the ages and chemical enrichment of galaxies are not fully determined by their stellar masses or velocity dispersion, but also depend on the mass of the host halos. Our findings suggest that the velocity dispersion is the best proxy of the stellar population parameters with halo mass playing a secondary yet noticeable role. We interpret that the origin of the correlation between the scatter of these relations and the ages and metallicities might be related to different halo formation times. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Comments: 22 pages, 15 figures. Accepted for publication in MNRAS

The black hole population in low-mass galaxies in large-scale cosmological simulations

Authors: Houda Haidar, Melanie Habouzit, Marta Volonteri, Mar Mezcua, Jenny Greene, Nadine Neumayer, Daniel Angles-Alcazar, Ignacio Martin-Navarro, Nils Hoyer, Yohan Dubois, Romeel Dave

Abstract: Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that ar… ▽ More Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that are too massive, and that the BH occupation fraction at z=0 is not inherited from the simulation seeding modeling. The ability of BHs and their host galaxies to power an AGN depends on BH and galaxy subgrid modeling. The fraction of AGN in low-mass galaxies is not used to calibrate the simulations, and thus can be used to differentiate galaxy formation models. AGN fractions at z=0 span two orders of magnitude at fixed galaxy stellar mass in simulations, similarly to observational constraints, but uncertainties and degeneracies affect both observations and simulations. The agreement is difficult to interpret due to differences in the masses of simulated and observed BHs, BH occupation fraction affected by numerical choices, and an unknown fraction of obscured AGN. Our work advocates for more thorough comparisons with observations to improve the modeling of cosmological simulations, and our understanding of BH and galaxy physics in the low-mass regime. The mass of BHs, their ability to efficiently accrete gas, and the AGN fraction in low-mass galaxies have important implications for the build-up of the entire BH and galaxy populations with time. △ Less

Submitted 14 June, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted in MNRAS, 21 pages, 11 figures, 1 table

Rejuvenation triggers nuclear activity in nearby galaxies

Authors: Ignacio Martin-Navarro, Francesco Shankar, Mar Mezcua

Abstract: Feedback, in particular from active galactic nuclei (AGN), is believed to play a crucial role in the evolution of galaxies. In the local Universe, many galaxies with an AGN are indeed observed to reside in the so-called green valley, usually interpreted as a transition phase from a blue star-forming to a red quenched state. We use data from the Sloan Digital Sky Survey to show that such an interpr… ▽ More Feedback, in particular from active galactic nuclei (AGN), is believed to play a crucial role in the evolution of galaxies. In the local Universe, many galaxies with an AGN are indeed observed to reside in the so-called green valley, usually interpreted as a transition phase from a blue star-forming to a red quenched state. We use data from the Sloan Digital Sky Survey to show that such an interpretation requires substantial revision. Optically-selected nearby AGN galaxies follow exponentially declining star formation histories, as normal galaxies of similar stellar and dark matter halo mass, reaching in the recent past ($\sim$0.1 Gyr ago) star formation rate levels consistent with a quiescent population. However, we find that local AGN galaxies have experienced a sudden increase in their star formation rate, unfolding on timescales similar to those typical of AGN activity, suggesting that both star formation and AGN activity were triggered simultaneously. We find that this quenching followed by an enhancement in the star formation rate is common to AGN galaxies and more pronounced in early type galaxies. Our results demonstrate that local AGN galaxies are not just a simple transition type between star-forming and quiescent galaxies as previously postulated. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS Letters. 5 pages, 5 figures, and 1 table containing the sample properties

Fornax 3D project: assessing the diversity of IMF and stellar population maps within the Fornax Cluster

Authors: I. Martín-Navarro, F. Pinna, L. Coccato, J. Falcón-Barroso, G. van de Ven, M. Lyubenova, E. M. Corsini, K. Fahrion, D. A. Gadotti, E. Iodice, R. M. McDermid, A. Poci, M. Sarzi, T. W. Spriggs, S. Viaene, P. T. de Zeeuw, L. Zhu

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

Submitted 29 July, 2021; originally announced July 2021.

Comments: 14 pages and 10 figures (plus Appendix including all the stellar population maps). Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 654, A59 (2021)

The Fornax3D project: Planetary nebulae catalogue and independent distance measurements to Fornax cluster galaxies

Authors: T. W. Spriggs, M. Sarz, P. M. Galán-de Anta, R. Napiwotzki, S. Viaene, B. Nedelchev, L. Coccato, E. M. Corsini, K. Fahrion, J. Falcón-Barroso, D. A. Gadotti, E. Iodice, M. Lyubenova, I. Martín-Navarro, R. M. McDermid, L. Morelli, F. Pinna, G. van de Ven, P. T. de Zeeuw, L. Zhu

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

Submitted 22 July, 2021; v1 submitted 20 July, 2021; originally announced July 2021.

Comments: 26 pages, 16 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics Journal. Corrected Fig. 7 legend

Journal ref: A&A 653, A167 (2021)

The Fornax 3D project: PNe populations and stellar metallicity in edge-on galaxies

Authors: P. M. Galán-de Anta, M. Sarzi, T. W. Spriggs, B. Nedelchev, F. Pinna, I. Martín-Navarro, L. Coccato, E. M. Corsini, P. T. de Zeeuw, J. Falcón-Barroso, D. A. Gadotti, E. Iodice, R. J. J. Grand, K. Fahrion, M. Lyubenova, R. M. McDermid, L. Morelli, G. van de Ven, S. Viaene, L. Zhu

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

Submitted 22 June, 2021; originally announced June 2021.

Journal ref: A&A 652, A109 (2021)

Anisotropic satellite galaxy quenching modulated by supermassive black hole activity

Authors: Ignacio Martín-Navarro, Annalisa Pillepich, Dylan Nelson, Vicente Rodriguez-Gomez, Martina Donnari, Lars Hernquist, Volker Springel

Abstract: The evolution of satellite galaxies is shaped by their constant interaction with the circum galactic medium surrounding central galaxies, which in turn may be affected by gas and energy ejected from the central supermassive black hole. However, the nature of this coupling between black holes and galaxies is highly debated and observational evidence remains scarce. Here we report an analysis of arc… ▽ More The evolution of satellite galaxies is shaped by their constant interaction with the circum galactic medium surrounding central galaxies, which in turn may be affected by gas and energy ejected from the central supermassive black hole. However, the nature of this coupling between black holes and galaxies is highly debated and observational evidence remains scarce. Here we report an analysis of archival data on 124,163 satellite galaxies in the potential wells of 29,631 dark matter halos with masses between 10$^{12}$ and $10^{14}$ solar masses. We find that quiescent satellites are relatively less frequent along the minor axis of their central galaxies. This observation might appear counterintuitive as black hole activity is expected to eject mass and energy preferentially in the direction of the minor axis of the host galaxy. However, we show that the observed signal results precisely from the ejective nature of black hole feedback in massive halos, as active galactic nuclei-powered outflows clear out the circumgalactic medium, reducing the ram pressure and thus preserving star formation. This interpretation is supported by the IllustrisTNG suite of cosmological numerical simulations, where a similar modulation is observed even though the sub-grid implementation of black hole feedback is effectively isotropic. Our results provide compelling observational evidence for the role of black holes in regulating galaxy evolution over spatial scales differing by several orders of magnitude. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: To appear in the 10 June issue of Nature

X-ray bubbles in the circumgalactic medium of TNG50 Milky Way- and M31-like galaxies: signposts of supermassive black hole activity

Authors: Annalisa Pillepich, Dylan Nelson, Nhut Truong, Rainer Weinberger, Ignacio Martin-Navarro, Volker Springel, Sandy M. Faber, Lars Hernquist

Abstract: The TNG50 cosmological simulation produces X-ray emitting bubbles, shells, and cavities in the circumgalactic gas above and below the stellar disks of Milky Way- and Andromeda-like galaxies with morphological features reminiscent of the eROSITA and Fermi bubbles in the Galaxy. Two-thirds of the 198 MW/M31 analogues inspected in TNG50 at z=0 show one or more large-scale, coherent features of over-p… ▽ More The TNG50 cosmological simulation produces X-ray emitting bubbles, shells, and cavities in the circumgalactic gas above and below the stellar disks of Milky Way- and Andromeda-like galaxies with morphological features reminiscent of the eROSITA and Fermi bubbles in the Galaxy. Two-thirds of the 198 MW/M31 analogues inspected in TNG50 at z=0 show one or more large-scale, coherent features of over-pressurized gas that impinge into the gaseous halo. Some of the galaxies include a succession of bubbles or shells of increasing size, ranging from a few to many tens of kpc. These are prominent in gas pressure, X-ray emission and gas temperature, and often exhibit sharp boundaries with typical shock Mach numbers of 2-4. The gas in the bubbles outflows with maximum (95th pctl) radial velocities of 100-1500 km/s. TNG50 bubbles expand with speeds as high as 1000-2000 km/s (about 1-2 kpc/Myr), but with a great diversity and with larger bubbles expanding at slower speeds. The bubble gas is at 10^6.4-7.2 K temperatures and is enriched to metallicities of 0.5-2 solar. In TNG50, the bubbles are a manifestation of episodic, kinetic, wind-like energy injections from the supermassive black holes at the galaxy centers that accrete at low Eddington ratios. According to TNG50, X-ray, and possibly gamma-ray, bubbles similar to those observed in the Milky Way should be a frequent feature of disk-like galaxies prior to, or on the verge of, being quenched. They should be within the grasp of eROSITA in the local Universe. △ Less

Submitted 26 October, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: Accepted on MNRAS. Added details on the SMBH feedback model and a new panel in Fig. 11 with the X-ray luminosities of TNG50 SMBHs. Updated discussion and observability opportunities

Diversity of nuclear star cluster formation mechanisms revealed by their star formation histories

Authors: K. Fahrion, M. Lyubenova, G. van de Ven, M. Hilker, R. Leaman, J. Falcón-Barroso, A. Bittner, L. Coccato, E. M. Corsini, D. A. Gadotti, E. Iodice, R. M. McDermid, I. Martín-Navarro, F. Pinna, A. Poci, M. Sarzi, P. T. de Zeeuw, L. Zhu

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

Submitted 13 April, 2021; originally announced April 2021.

Comments: accepted for publication in Astronomy & Astrophysics, 18 pages, 10 figures (including appendix)

Journal ref: A&A 650, A137 (2021)

The Fornax3D project: Assembly histories of lenticular galaxies from a combined dynamical and population orbital analysis

Authors: A. Poci, R. M. McDermid, M. Lyubenova, L. Zhu, G. van de ven, E. Iodice, L. Coccato, F. Pinna, E. M. Corsini, J. Falcón-Barroso, D. A. Gadotti, R. J. J. Grand, K. Fahrion, I. Martín-Navarro, M. Sarzi, S. Viaene, P. T. de Zeeuw

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

Submitted 7 February, 2021; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 647, A145 (2021)

Galaxies within galaxies in the TIMER survey: stellar populations of inner bars are scaled replicas of main bars

Authors: Adrian Bittner, Adriana de Lorenzo-Cáceres, Dimitri A. Gadotti, Patricia Sánchez-Blázquez, Justus Neumann, Paula Coelho, Jesús Falcón-Barroso, Francesca Fragkoudi, Taehyun Kim, Ignacio Martín-Navarro, Jairo Méndez-Abreu, Isabel Pérez, Miguel Querejeta, Glenn van de Ven

Abstract: Inner bars are frequent structures in the local Universe and thought to substantially influence the nuclear regions of disc galaxies. In this study we explore the structure and dynamics of inner bars by deriving maps and radial profiles of their mean stellar population content and comparing them to previous findings in the context of main bars. To this end, we exploit observations obtained with th… ▽ More Inner bars are frequent structures in the local Universe and thought to substantially influence the nuclear regions of disc galaxies. In this study we explore the structure and dynamics of inner bars by deriving maps and radial profiles of their mean stellar population content and comparing them to previous findings in the context of main bars. To this end, we exploit observations obtained with the integral-field spectrograph MUSE of three double-barred galaxies in the TIMER sample. The results indicate that inner bars can be distinguished based on their stellar population properties alone. More precisely, inner bars show elevated metallicities and depleted [$α$/Fe] abundances. Although they exhibit slightly younger stellar ages compared to the nuclear disc, the typical age differences are small, except at their outer ends. These ends of the inner bars are clearly younger compared to their inner parts, an effect known from main bars as orbital age separation. In particular, the youngest stars (i.e. those with the lowest radial velocity dispersion) seem to occupy the most elongated orbits along the (inner) bar major axis. We speculate that these distinct ends of bars could be connected to the morphological feature of ansae. Radial profiles of metallicity and [$α$/Fe] enhancements are flat along the inner bar major axis, but show significantly steeper slopes along the minor axis. This radial mixing in the inner bar is also known from main bars and indicates that inner bars significantly affect the radial distribution of stars. In summary, based on maps and radial profiles of the mean stellar population content and in line with previous TIMER results, inner bars appear to be scaled down versions of the main bars seen in galaxies. This suggests the picture of a "galaxy within a galaxy", with inner bars in nuclear discs being dynamically equivalent to main bars in main galaxy discs. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: Accepted for publication in A&A; 11 pages, 4 figures, 2 tables

Journal ref: A&A 646, A42 (2021)

The kinematics of young and old stellar populations in nuclear rings of MUSE TIMER galaxies

Authors: D. Rosado-Belza, J. Falcón-Barroso, J. H. Knapen, A. Bittner, D. A. Gadotti, J. Neumann, A. de Lorenzo-Cáceres, J. Méndez-Abreu, M. Querejeta, I. Martín-Navarro, P. Sánchez-Blázquez, P. R. T. Coelho, M. Martig, G. van de Ven, T. Kim

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

Submitted 27 October, 2020; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 22 pages, 24 figures

Journal ref: A&A 644, A116 (2020)

The inside-out formation of nuclear discs and the absence of old central spheroids in barred galaxies of the TIMER survey

Authors: Adrian Bittner, Patricia Sánchez-Blázquez, Dimitri A. Gadotti, Justus Neumann, Francesca Fragkoudi, Paula Coelho, Adriana de Lorenzo-Cáceres, Jesús Falcón-Barroso, Taehyun Kim, Ryan Leaman, Ignacio Martín-Navarro, Jairo Méndez-Abreu, Isabel Pérez, Miguel Querejeta, Marja K. Seidel, Glenn van de Ven

Abstract: The centres of disc galaxies host a variety of structures built via both internal and external processes. In this study, we constrain the formation and evolution of these central structures, in particular nuclear rings and nuclear discs, by deriving maps of mean stellar ages, metallicities and [$α$/Fe] abundances. We use observations obtained with the MUSE integral-field spectrograph for the TIMER… ▽ More The centres of disc galaxies host a variety of structures built via both internal and external processes. In this study, we constrain the formation and evolution of these central structures, in particular nuclear rings and nuclear discs, by deriving maps of mean stellar ages, metallicities and [$α$/Fe] abundances. We use observations obtained with the MUSE integral-field spectrograph for the TIMER sample of 21 massive barred galaxies. Our results indicate that nuclear discs and nuclear rings are part of the same physical component, with nuclear rings constituting the outer edge of nuclear discs. All nuclear discs in the sample are clearly distinguished based on their stellar population properties. As expected in the picture of bar-driven secular evolution, nuclear discs are younger, more metal-rich, and show lower [$α$/Fe] enhancements, as compared to their immediate surroundings. Moreover, nuclear discs exhibit well-defined radial gradients, with ages and metallicities decreasing, and [$α$/Fe] abundances increasing with radius out to the nuclear ring. Often, these gradients show no breaks from the edge of the nuclear disc until the centre, suggesting that these structures extend to the very centres of the galaxies. We argue that continuous (stellar) nuclear discs may form from a series of bar-built (initially gas-rich) nuclear rings that grow in radius, as the bar evolves. In this picture, nuclear rings are simply the (often) star-forming outer edge of nuclear discs. Finally, by combining our results with those from a accompanying kinematic study, we do not find evidence for the presence of large, dispersion-dominated components in the centres of these galaxies. This could be a result of quiet merger histories, despite the large galaxy masses, or perhaps high angular momentum and strong feedback processes preventing the formation of these kinematically hot components. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A; 29 pages, 12 figures, 1 table, including 3 appendices; FITS files with the derived (spatially-resolved) mean stellar population properties (age, metallicity, and [$α$/Fe] abundance) for all TIMER galaxies are available at https://www.muse-timer.org/data

Journal ref: A&A 643, A65 (2020)

Kinematic signatures of nuclear discs and bar-driven secular evolution in nearby galaxies of the MUSE TIMER project

Authors: Dimitri A. Gadotti, Adrian Bittner, Jesus Falcon-Barroso, Jairo Mendez-Abreu, Taehyun Kim, Francesca Fragkoudi, Adriana de Lorenzo-Caceres, Ryan Leaman, Justus Neumann, Miguel Querejeta, Patricia Sanchez-Blazquez, Marie Martig, Ignacio Martin-Navarro, Isabel Perez, Marja K. Seidel, Glenn van de Ven

Abstract: The central regions of disc galaxies hold clues to the processes that dominate their formation and evolution. The TIMER project has obtained high signal-to-noise and spatial resolution integral-field spectroscopy data of the inner few kpc of 21 nearby massive barred galaxies, allowing studies of the stellar kinematics with unprecedented spatial resolution. We confirm theoretical predictions of the… ▽ More The central regions of disc galaxies hold clues to the processes that dominate their formation and evolution. The TIMER project has obtained high signal-to-noise and spatial resolution integral-field spectroscopy data of the inner few kpc of 21 nearby massive barred galaxies, allowing studies of the stellar kinematics with unprecedented spatial resolution. We confirm theoretical predictions of the effects of bars on stellar kinematics, and identify box/peanuts through kinematic signatures in mildly and moderately inclined galaxies, finding a lower limit to the fraction of massive barred galaxies with box/peanuts at ~62%. Further, we provide kinematic evidence of the connection between barlenses, box/peanuts and bars. We establish the presence of nuclear discs in 19 galaxies and show that their kinematics are characterised by near-circular orbits with low pressure support, and are consistent with the bar-driven secular evolution picture for their formation. In fact, we show that these nuclear discs have, in the region where they dominate, larger rotational support than the underlying main galaxy disc. We define a kinematic radius for the nuclear discs and show that it relates to bar radius, ellipticity and strength, and bar-to-total ratio. Comparing our results with photometric studies, we find that state-of-the-art galaxy image decompositions are able to discern nuclear discs from classical bulges, if the images employed have enough physical spatial resolution. In fact, we show that nuclear discs are typically identified in such image decompositions as photometric bulges with (near-)exponential profiles. However, we find that the presence of composite bulges (galaxies hosting both a classical bulge and a nuclear disc) can often be unnoticed in studies based on photometry alone, and suggest a more stringent threshold to the Sersic index to identify galaxies with pure classical bulges. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A; 24 pages, 9 figures, 3 tables, including 2 appendices. A full resolution PDF is available at http://www.sc.eso.org/~dgadotti/kinematics.pdf and FITS files with the derived (spatially-resolved) stellar kinematics (v, sigma, h3 and h4) for all TIMER galaxies are available at https://www.muse-timer.org/data

Journal ref: A&A 643, A14 (2020)

Surface Brightness Fluctuations for constraining the chemical enrichment of massive galaxies

Authors: A. Vazdekis, P. Rodríguez-Beltrán, M. Cerviño, M. Montes, I. Martín-Navarro, M. B. Beasley

Abstract: 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 bas… ▽ More 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. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: Proceedings of the IAU Symposium No. 359: Galaxy evolution and feedback across different environments. Editors: T. Storchi-Bergmann, R. Overzier, W. Forman & R. Riffel

Reconstructing the mass accretion histories of nearby Red Nuggets with their globular cluster systems

Authors: Michael A. Beasley, Ryan Leaman, Ignacio Trujillo, Mireia Montes, Alejandro Vazdekis, Núria Salvador Rusiñol, Elham Eftekhari, Anna Ferré-Mateu, Ignacio Martin-Navarro

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: Proceedings of the IAU Symposium No. 359: Galaxy evolution and feedback across different environments. Editors: T. Storchi-Bergmann, R. Overzier, W. Forman \& R. Riffel

Recovering age-metallicity distributions from integrated spectra: validation with MUSE data of a nearby nuclear star cluster

Authors: A. Boecker, M. Alfaro-Cuello, N. Neumayer, I. Martín-Navarro, R. Leaman

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

Submitted 7 May, 2020; originally announced May 2020.

Comments: 28 pages, 15 figures incl. appendix, accepted for publication in ApJ

The Fornax 3D project: Non-linear colour-metallicity relation of globular clusters

Authors: K. Fahrion, M. Lyubenova, M. Hilker, G. van de Ven, J. Falcón-Barroso, R. Leaman, I. Martín-Navarro, A. Bittner, L. Coccato, E. M. Corsini, D. A. Gadotti, E. Iodice, R. M. McDermid, F. Pinna, M. Sarzi, S. Viaene, P. T. de Zeeuw, L. Zhu

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

Submitted 13 May, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

Comments: 12 pages + appendix, updated to published version

Journal ref: A&A 637, A27 (2020)

The Fornax 3D project: Globular clusters tracing kinematics and metallicities

Authors: K. Fahrion, M. Lyubenova, M. Hilker, G. van de Ven, J. Falcón-Barroso, R. Leaman, I. Martín-Navarro, A. Bittner, L. Coccato, E. M. Corsini, D. A. Gadotti, E. Iodice, R. M. McDermid, F. Pinna, M. Sarzi, S. Viaene, P. T. de Zeeuw, L. Zhu

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

Submitted 13 May, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

Comments: 13 pages + appendix, accepted for publication in A&A, updated to published version

Journal ref: A&A 637, A26 (2020)

Fornax 3D project: automated detection of planetary nebulae in the centres of early-type galaxies and first results

Authors: T. W. Spriggs, M. Sarzi, R. Napiwotzki, P. M. Galan-de Anta, S. Viaene, B. Nedelchev, L. Coccato, E. M. Corsini, P. T. de Zeeuw, J. Falcon-Barroso, D. A. Gadotti, E. Iodice, M. Lyubenova, I. Martin-Navarro, R. M. McDermid, F. Pinna, G. van de Ven, L. Zhu

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

Submitted 18 November, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 16 pages, 17 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics Journal. Replaced with published version

Journal ref: A&A 637, A62 (2020)

Surface Brightness Fluctuation spectra to constrain stellar population properties

Authors: A. Vazdekis, M. Cerviño, M. Montes, I. Martín-Navarro, M. Beasley

Abstract: 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… ▽ More 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. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: MNRAS, accepted. 25 pages, 13 figures. Models can be downloaded from http://miles.iac.es/ (OTHER PREDICTIONS/DATA)

Black hole feedback and the evolution of massive early-type galaxies

Authors: Ignacio Martin-Navarro, Joseph N. Burchett, Mar Mezcua

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

Submitted 28 October, 2019; originally announced October 2019.

Comments: 10 pages, 5 figure, accepted for publication in MNRAS,

Quantifying the effect of black hole feedback from the central galaxy on the satellite populations of groups and clusters

Authors: Ignacio Martin-Navarro, Joseph N. Burchett, Mar Mezcua

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

Submitted 27 September, 2019; originally announced September 2019.

Comments: Accepted for publication in ApJL. 7 pages, 4 figures

IMF radial gradients in most massive early-type galaxies

Authors: F. La Barbera, A. Vazdekis, I. Ferreras, A. Pasquali, C. Allende Prieto, I. Martín-Navarro, D. S. Aguado, R. R. de Carvalho, S. Rembold, J. Falcón-Barroso, G. van de Ven

Abstract: 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<σ<350$km/s. Using… ▽ More 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<σ<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. △ Less

Submitted 6 September, 2019; v1 submitted 3 September, 2019; originally announced September 2019.

Comments: 23 pages, 13 figures, accepted for publication in MNRAS

Constraining nuclear star cluster formation using MUSE-AO observations of the early-type galaxy FCC47

Authors: Katja Fahrion, Mariya Lyubenova, Glenn van de Ven, Ryan Leaman, Michael Hilker, Ignacio Martín-Navarro, Ling Zhu, Mayte Alfaro-Cuello, Lodovico Coccato, Enrico M. Corsini, Jesús Falcón-Barroso, Enrichetta Iodice, Richard M. McDermid, Marc Sarzi, Tim de Zeeuw

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

Submitted 1 July, 2019; originally announced July 2019.

Comments: 21 pages, accepted in A&A

Journal ref: A&A 628, A92 (2019)

The Fornax3D project: Tracing the assembly history of the cluster from the kinematic and line-strength maps

Authors: E. Iodice, M. Sarzi, A. Bittner, L. Coccato, L. Costantin, E. M. Corsini, G. van de Ven, P. T. de Zeeuw, J. Falcòn-Barroso, D. A. Gadotti, M. Lyubenova, I. Martìn-Navarro, R. M. McDermid, B. Nedelchev, F. Pinna, A. Pizzella, M. Spavone, S. Viaene

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

Submitted 19 June, 2019; originally announced June 2019.

Comments: Accepted for publication in Astronomy and Astrophysics Journal. The high-resolution pdf file, including all maps in the Appendix, is available at the following link: https://www.dropbox.com/s/oukfd56dogf6v9t/F3D_sample_paper.pdf?dl=0

Journal ref: A&A 627, A136 (2019)

Star formation quenching imprinted on the internal structure of naked red nuggets

Authors: Ignacio Martín-Navarro, Glenn van de Ven, Akın Yıldırım

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

Submitted 31 May, 2019; originally announced June 2019.

Comments: Accepted for publication in MNRAS; 13 pages, 8 figures. Figure 8 summarizes the conclusions of the paper