Showing papers by "Jesús Falcón-Barroso published in 2021"

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

Abstract: In order to assess the impact of the environment on the formation and evolution of galaxies, accurate assembly histories of such galaxies are needed. However, these measurements are observationally difficult owing to the diversity of formation paths that lead to the same present-day state of a galaxy. 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 ≲ 50 km 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.

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

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 relative contribution of either channel and their correlation with galaxy properties are not yet constrained observationally. Our aim was 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 gal ∼ 108 and 1010.5 M ⊙ and NSC masses between M NSC ∼ 105 and 108.5 M ⊙ . Using Multi-Unit Spectroscopic Explorer data from the Fornax 3D survey and the ESO archive, we derived star formation histories, mean ages, and metallicities of NSCs, and compared them to the host galaxies. In many low-mass galaxies, the NSCs are significantly more metal poor than their 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 (> 107 M ⊙ ) 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 gal ∼ 109 M ⊙ and NSC masses M NSC ∼ 107 M ⊙ .

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

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

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

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 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 $\sigma_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 $\lambda_z > 08$) 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

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

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 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 ([M/H] ≲ −0.1), 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.

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

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 the integral-field spectrograph MUSE of three double-barred galaxies in the TIMER sample. The results indicate that inner bars can be clearly distinguished based on their stellar population properties alone. More precisely, inner bars show significantly 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.

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

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 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 [O III] 5007 A 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, [O III] 5007 A 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 ± 0.32 Mpc (μ PNLF = 31.49 ± 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.

BAYES-LOSVD: A Bayesian framework for non-parametric extraction of the line-of-sight velocity distribution of galaxies

Abstract: We introduce BAYES-LOSVD, a novel implementation of the non-parametric extraction of line-of-sight velocity distributions (LOSVDs) in galaxies. We employed Bayesian inference to obtain robust LOSVDs and associated uncertainties. Our method relies on a principal component analysis to reduce the dimensionality on the set of templates required for the extraction and thus increase the performance of the code. In addition, we implemented several options to regularise the output solutions. Our tests, conducted on mock spectra, confirm the ability of our approach to model a wide range of LOSVD shapes, overcoming limitations of the most widely used parametric methods (e.g., Gauss-Hermite expansion). We present examples of LOSVD extractions for real galaxies with known peculiar LOSVD shapes, including NGC 4371, IC 0719, and NGC 4550, using MUSE and SAURON integral-field unit (IFU) data. Our implementation can also handle data from other popular IFU surveys (e.g., ATLAS3D , CALIFA, MaNGA, SAMI).

NGC 5746: Formation history of a massive disc-dominated galaxy

Abstract: The existence of massive galaxies lacking a classical bulge has often been proposed as a challenge to $\Lambda$CDM. However, recent simulations propose that a fraction of massive disc galaxies might have had very quiescent merger histories, and also that mergers do not necessarily build classical bulges. We test these ideas with deep MUSE observations of NGC 5746, a massive ($\sim 10^{11}$ M$_\odot$) edge-on disc galaxy with no classical bulge. We analyse its stellar kinematics and stellar populations, and infer that a massive and extended disc formed very early: 80% of the galaxy's stellar mass formed more than 10 Gyr ago. Most of the thick disc and the bar formed during that early phase. The bar drove gas towards the center and triggered the formation of the nuclear disc followed by the growth of a boxy/peanut-shaped bulge. Around $\sim$ 8 Gyr ago, a $\sim$1:10 merger happened, possibly on a low-inclination orbit. The satellite did not cause significant vertical heating, did not contribute to the growth of a classical bulge, and did not destroy the bar and the nuclear disc. It was however an important event for the galaxy: by depositing its stars throughout the whole galaxy it contributed $\sim 30$% of accreted stars to the thick disc. NGC 5746 thus did not completely escape mergers, but the only relatively recent significant merger did not damage the galaxy and did not create a classical bulge. Future observations will reveal if this is representative of the formation histories of massive disc galaxies.

Capturing the Physics of MaNGA Galaxies with Self-supervised Machine Learning

Abstract: As available data sets grow in size and complexity, advanced visualization tools enabling their exploration and analysis become more important. In modern astronomy, integral field spectroscopic galaxy surveys are a clear example of increasing dimensionality and complexity of datasets, which challenge the traditional methods used to extract the physical information they contain. We present the use of a novel self-supervised Machine Learning method to visualize the multi-dimensional information on stellar population and kinematics in the MaNGA survey in a two dimensional plane. Our framework is insensitive to non-physical properties such as the size of integral field unit (IFU) and is therefore able to order galaxies according to their resolved physical properties. Using the extracted representations, we show that galaxies naturally cluster into three well-known categories from a purely data driven perspective: rotating main-sequence disks, massive slow rotators and low-mass rotation-dominated quenched galaxies. The framework for data exploration is publicly released with this publication, ready to be used with the MaNGA or other integral field data sets.

A comparison between X-shooter spectra and PHOENIX models across the HR-diagram

Abstract: Aims. The path towards robust near-infrared extensions of stellar population models involves the confrontation between empirical and synthetic stellar spectral libraries across the wavelength ranges of photospheric emission. Indeed, the theory of stellar emission enters all population synthesis models, even when this is only implicit in the association of fundamental stellar parameters with empirical spectral library stars. With its near-ultraviolet to near-infrared coverage, the X-shooter Spectral Library (XSL) allows us to examine to what extent models succeed in reproducing stellar energy distributions (SEDs) and stellar absorption line spectra simultaneously.Methods. As a first example, this study compares the stellar spectra of XSL with those of the Gottingen Spectral Library, which are based on the PHOENIX synthesis code. The comparison was carried out both separately in the three arms of the X-shooter spectrograph known as UVB, VIS and NIR, and jointly across the whole spectrum. We did not discard the continuum in these comparisons; only reddening was allowed to modify the SEDs of the models.Results. When adopting the stellar parameters published with data release DR2 of XSL, we find that the SEDs of the models are consistent with those of the data at temperatures above 5000 K. Below 5000 K, there are significant discrepancies in the SEDs. When leaving the stellar parameters free to adjust, satisfactory representations of the SEDs are obtained down to about 4000 K. However, in particular below 5000 K and in the UVB spectral range, strong local residuals associated with intermediate resolution spectral features are then seen; the necessity of a compromise between reproducing the line spectra and reproducing the SEDs leads to dispersion between the parameters favored by various spectral ranges. We describe the main trends observed and we point out localized offsets between the parameters preferred in this global fit to the SEDs and the parameters in DR2. These depend in a complex way on the position in the Hertzsprung–Russell diagram (HRD). We estimate the effect of the offsets on bolometric corrections as a function of position in the HRD and use this for a brief discussion of their impact on the studies of stellar populations. A review of the literature shows that comparable discrepancies are mentioned in studies using other theoretical and empirical libraries.

An extension of the MILES library with derived Teff, log g, [Fe/H], and [α/Fe]

Abstract: Extragalactic astronomy and stellar astrophysics are intrinsically related. In fact, the determination of important galaxy properties such as stellar masses, star formation histories, or chemical abundances relies on the ability to model their stellar populations. One important ingredient of these models is stellar libraries. Empirical libraries must have a good coverage of Teff, [Z/H], and surface gravity, and have these parameters reliably determined. MILES is one of the most widely used empirical libraries. Here, we present an extension of this library with 205 new stars especially selected to cover important regions of the parameter space, including metal-poor stars down to [Fe/H] ∼-1.0. We describe the observations and data reductions as well as a new determination of the stellar parameters, including [α/Fe] ratio. The new MILES library contains 1070 stars with homogeneous and reliable determination of [Fe/H], Teff, log g, and [α/Fe] ratio.

Local variations of the Stellar Velocity Ellipsoid-I: the disc of galaxies in the Auriga simulations

Abstract: The connection between the Stellar Velocity Ellipsoid (SVE) and the dynamical evolution of galaxies has been a matter of debate in the last years and there is no clear consensus whether different heating agents (e.g. spiral arms, giant molecular clouds, bars and mergers) leave clear detectable signatures in the present day kinematics. Most of these results are based on a single and global SVE and have not taken into account that these agents do not necessarily equally affect all regions of the stellar disc.We study the 2D spatial distribution of the SVE across the stellar discs of Auriga galaxies, a set of high resolution magneto-hydrodynamical cosmological zoom-in simulations, to unveil the connection between local and global kinematic properties in the disc region. We find very similar, global, $\sigma_{z}/\sigma_{r}$= 0.80$\pm$ 0.08 values for galaxies of different Hubble types. This shows that the global properties of the SVE at z=0 are not a good indicator of the heating and cooling events experienced by galaxies. We also find that similar $\sigma_{z}/\sigma_{r}$radial profiles are obtained through different combinations of $\sigma_{z}$ and $\sigma_{r}$ trends: at a local level, the vertical and radial components can evolve differently, leading to similar $\sigma_{z}/\sigma_{r}$ profiles at z=0. By contrast, the 2D spatial distribution of the SVE varies a lot more from galaxy to galaxy. Present day features in the SVE spatial distribution may be associated with specific interactions such as fly-by encounters or the accretion of low mass satellites even in the cases when the global SVE is not affected. The stellar populations decomposition reveals that young stellar populations present colder and less isotropic SVEs and more complex 2D distributions than their older and hotter counterparts.

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

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 understand the later stages of stellar evolution. Previous works have indicated that a 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 very large telescope-multi unit spectroscopic explorer (VLT-MUSE) integral-field observations and dedicated PNe detection procedures, we constructed the PNe luminosity function and computed the luminosity-specific number of PNe α in both in- and off-plane regions of our edge-on systems.Results. Comparing these α 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.

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

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

Capturing the physics of MaNGA galaxies with self-supervised Machine Learning

Abstract: As available data sets grow in size and complexity, advanced visualization tools enabling their exploration and analysis become more important. In modern astronomy, integral field spectroscopic galaxy surveys are a clear example of increasing dimensionality and complexity of datasets, which challenge the traditional methods used to extract the physical information they contain. We present the use of a novel self-supervised Machine Learning method to visualize the multi-dimensional information on stellar population and kinematics in the MaNGA survey in a two dimensional plane. Our framework is insensitive to non-physical properties such as the size of integral field unit (IFU) and is therefore able to order galaxies according to their resolved physical properties. Using the extracted representations, we show that galaxies naturally cluster into three well-known categories from a purely data driven perspective: rotating main-sequence disks, massive slow rotators and low-mass rotation-dominated quenched galaxies. The framework for data exploration is publicly released with this publication, ready to be used with the MaNGA or other integral field data sets.