Showing 1–50 of 138 results for author: Bernardi, M

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

Identification of tidal features in deep optical galaxy images with Convolutional Neural Networks

Authors: H. Domínguez Sánchez, G. Martin, I. Damjanov, F. Buitrago, M. Huertas-Company, C. Bottrell, M. Bernardi, J. H. Knapen, J. Vega-Ferrero, R. Hausen, E. Kado-Fong, D. Población-Criado, H. Souchereau, O. K. Leste, B. Robertson, B. Sahelices, K. V. Johnston

Abstract: Interactions between galaxies leave distinguishable imprints in the form of tidal features which hold important clues about their mass assembly. Unfortunately, these structures are difficult to detect because they are low surface brightness features so deep observations are needed. Upcoming surveys promise several orders of magnitude increase in depth and sky coverage, for which automated methods… ▽ More Interactions between galaxies leave distinguishable imprints in the form of tidal features which hold important clues about their mass assembly. Unfortunately, these structures are difficult to detect because they are low surface brightness features so deep observations are needed. Upcoming surveys promise several orders of magnitude increase in depth and sky coverage, for which automated methods for tidal feature detection will become mandatory. We test the ability of a convolutional neural network to reproduce human visual classifications for tidal detections. We use as training $\sim$6000 simulated images classified by professional astronomers. The mock Hyper Suprime Cam Subaru (HSC) images include variations with redshift, projection angle and surface brightness ($μ_{lim}$ =26-35 mag arcsec$^{-2}$). We obtain satisfactory results with accuracy, precision and recall values of Acc=0.84, P=0.72 and R=0.85, respectively, for the test sample. While the accuracy and precision values are roughly constant for all surface brightness, the recall (completeness) is significantly affected by image depth. The recovery rate shows strong dependence on the type of tidal features: we recover all the images showing shell features and 87% of the tidal streams; these fractions are below 75% for mergers, tidal tails and bridges. When applied to real HSC images, the performance of the model worsens significantly. We speculate that this is due to the lack of realism of the simulations and take it as a warning on applying deep learning models to different data domains without prior testing on the actual data. △ Less

Submitted 6 March, 2023; originally announced March 2023.

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

Revisiting the SFR-Mass relation at z=0 with detailed deep learning based morphologies

Authors: Helena Domínguez Sánchez, Mariangela Bernardi, Marc Huertas-Company

Abstract: Galaxy morphology is a key parameter in galaxy evolution studies. The enormous number of galaxies which current and future surveys will observe demand of automated methods for morphological classification. Supervised learning techniques have been successfully used for the morphological classification of galaxies from different datasets, including Sloan Digital Sky Survey (SDSS), Mapping Galaxies w… ▽ More Galaxy morphology is a key parameter in galaxy evolution studies. The enormous number of galaxies which current and future surveys will observe demand of automated methods for morphological classification. Supervised learning techniques have been successfully used for the morphological classification of galaxies from different datasets, including Sloan Digital Sky Survey (SDSS), Mapping Galaxies with Apache Point Observatory (MaNGA) or Dark Energy Survey (DES). With these proceedings, we release the morphological catalogue for a sample of 670,000 SDSS galaxies based on the deep learning models trained on SDSS RGB images with morphological labels from human-based classification catalogues. The released catalogue includes binary classifications (early-type versus late-type, elliptical versus lenticular, identification of edge-on and barred galaxies) plus a T-Type. The classifications also include k-fold based uncertainties. This is, as of today, the largest catalogue including a T-Type classification. As an example of the scientific potential of this classification, we show how the location of the galaxies in the star formation - stellar mass plane (SFR-M$^{*}$) depends on morphology. This is the first time the SFR-M$^{*}$ relation is combined with T-Type information for such a large sample of galaxies. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 7 pages, 2 figures, proceedings for EAS 2022 S11, to be published in Memorie della SAIt. Morphological catalogue can be downloaded from https://archive.cefca.es/ancillary_data/sdss_morphological_catalogues/sdss_morphological_catalogues.tar.gz

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

Authors: M. Bernardi, H. Dominguez Sanchez, R. K. Sheth, J. R. Brownstein, R. R. Lane

Abstract: We study systematics associated with estimating simple stellar population (SSP) parameters -- age, metallicity [M/H], $α$-enhancement [$α$/Fe] and IMF shape -- and associated $M_*/L$ gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs) and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) h… ▽ More We study systematics associated with estimating simple stellar population (SSP) parameters -- age, metallicity [M/H], $α$-enhancement [$α$/Fe] and IMF shape -- and associated $M_*/L$ gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs) and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H$_β$, $\langle$Fe$\rangle$, TiO$_{\rm 2SDSS}$ and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius $R_e$: (a) ages are younger and [$α$/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about $R_e/2$; (c) this makes $M_*/L$ about $2\times$ larger in the central regions than beyond $R_e/2$. While the models of Conroy et al. (2018) return similar [M/H] and [$α$/Fe] profiles, the age and (hence) $M_*/L$ profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs and S0s: younger ages and less enhanced [$α$/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred $M_*/L$ gradients significantly. △ Less

Submitted 10 November, 2022; originally announced November 2022.

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

Lessons Learned from the Two Largest Galaxy Morphological Classification Catalogues built by Convolutional Neural Networks

Authors: Ting-Yun Cheng, H. Domínguez Sánchez, J. Vega-Ferrero, C. J. Conselice, M. Siudek, A. Aragón-Salamanca, M. Bernardi, R. Cooke, L. Ferreira, M. Huertas-Company, J. Krywult, A. Palmese, A. Pieres, A. A. Plazas Malagón, A. Carnero Rosell, D. Gruen, D. Thomas, D. Bacon, D. Brooks, D. J. James, D. L. Hollowood, D. Friedel, E. Suchyta, E. Sanchez, F. Menanteau , et al. (32 additional authors not shown)

Abstract: We compare the two largest galaxy morphology catalogues, which separate early and late type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of $\sim$21 mag. The methodologies used for the construction of the catalogues include differences such as the… ▽ More We compare the two largest galaxy morphology catalogues, which separate early and late type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of $\sim$21 mag. The methodologies used for the construction of the catalogues include differences such as the cutout sizes, the labels used for training, and the input to the CNN - monochromatic images versus $gri$-band normalized images. In addition, one catalogue is trained using bright galaxies observed with DES ($i<18$), while the other is trained with bright galaxies ($r<17.5$) and `emulated' galaxies up to $r$-band magnitude $22.5$. Despite the different approaches, the agreement between the two catalogues is excellent up to $i<19$, demonstrating that CNN predictions are reliable for samples at least one magnitude fainter than the training sample limit. It also shows that morphological classifications based on monochromatic images are comparable to those based on $gri$-band images, at least in the bright regime. At fainter magnitudes, $i>19$, the overall agreement is good ($\sim$95\%), but is mostly driven by the large spiral fraction in the two catalogues. In contrast, the agreement within the elliptical population is not as good, especially at faint magnitudes. By studying the mismatched cases we are able to identify lenticular galaxies (at least up to $i<19$), which are difficult to distinguish using standard classification approaches. The synergy of both catalogues provides an unique opportunity to select a population of unusual galaxies. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 17 pages, 14 figures (1 appendix for galaxy examples including 3 figures)

Testing the key role of the stellar mass-halo mass relation in galaxy merger rates and morphologies via DECODE, a novel Discrete statistical sEmi-empiriCal mODEl

Authors: Hao Fu, Francesco Shankar, Mohammadreza Ayromlou, Max Dickson, Ioanna Koutsouridou, Yetli Rosas-Guevara, Christopher Marsden, Kristina Brocklebank, Mariangela Bernardi, Nikolaos Shiamtanis, Joseph Williams, Lorenzo Zanisi, Viola Allevato, Lumen Boco, Silvia Bonoli, Andrea Cattaneo, Paola Dimauro, Fangzhou Jiang, Andrea Lapi, Nicola Menci, Stefani Petropoulou, Carolin Villforth

Abstract: The relative roles of mergers and star formation in regulating galaxy growth are still a matter of intense debate. We here present our DECODE, a new Discrete statistical sEmi-empiriCal mODEl specifically designed to predict rapidly and efficiently, in a full cosmological context, galaxy assembly and merger histories for any given input stellar mass-halo mass (SMHM) relation. DECODE generates objec… ▽ More The relative roles of mergers and star formation in regulating galaxy growth are still a matter of intense debate. We here present our DECODE, a new Discrete statistical sEmi-empiriCal mODEl specifically designed to predict rapidly and efficiently, in a full cosmological context, galaxy assembly and merger histories for any given input stellar mass-halo mass (SMHM) relation. DECODE generates object-by-object dark matter merger trees (hence discrete) from accurate subhalo mass and infall redshift probability functions (hence statistical) for all subhaloes, including those residing within other subhaloes, with virtually no resolution limits on mass or volume. Merger trees are then converted into galaxy assembly histories via an input, redshift dependent SMHM relation, which is highly sensitive to the significant systematics in the galaxy stellar mass function and on its evolution with cosmic time. DECODE can accurately reproduce the predicted mean galaxy merger rates and assembly histories of hydrodynamic simulations and semi-analytic models, when adopting in input their SMHM relations. In the present work we use DECODE to prove that only SMHM relations implied by stellar mass functions characterized by large abundances of massive galaxies and significant redshift evolution, at least at $M_\star \gtrsim 10^{11} \, M_\odot$, can simultaneously reproduce the local abundances of satellite galaxies, the galaxy (major merger) pairs since $z \sim 3$, and the growth of Brightest Cluster Galaxies. The same models can also reproduce the local fraction of elliptical galaxies, on the assumption that these are strictly formed by major mergers, but not the full bulge-to-disc ratio distributions, which require additional processes. △ Less

Submitted 9 August, 2022; v1 submitted 29 July, 2022; originally announced August 2022.

Comments: MNRAS, accepted, 29 pages, 25 figures

Coincidence between morphology and star-formation activity through cosmic time: the impact of the bulge growth

Authors: Paola Dimauro, Emanuele Daddi, Francesco Shankar, Andrea Cattaneo, Marc Huertas-Company, Mariangela Bernardi, Fernando Caro, Renato Dupke, Boris Häußler, Johnston Evelyn, Arianna Cortesi, Simona Mei, Reynier Peletier

Abstract: The origin of the quenching in galaxies is still highly debated. Different scenarios and processes are proposed. We use multi-band (400-1600 nm) bulge-disc decompositions of massive galaxies in the redshift range 0<z<2 to explore the distribution and the evolution of galaxies in the log SFR-logM* plane as a function of the stellar mass-weighted bulge-to-total ratio (BTM) and also for internal gala… ▽ More The origin of the quenching in galaxies is still highly debated. Different scenarios and processes are proposed. We use multi-band (400-1600 nm) bulge-disc decompositions of massive galaxies in the redshift range 0<z<2 to explore the distribution and the evolution of galaxies in the log SFR-logM* plane as a function of the stellar mass-weighted bulge-to-total ratio (BTM) and also for internal galaxy components (bulge/disc) separately. We find evidence of a clear link between the presence of a bulge and the flattening of the Main Sequence in the high-mass end. All bulgeless galaxies (BTM<0.2) lie on the main-sequence, and there is little evidence of a quenching channel without bulge growth. Galaxies with a significant bulge component (BTM>0.2) are equally distributed in number between star forming and passive regions. The vast majority of bulges in the Main Sequence galaxies are quiescent, while star formation is localized in the disc component. Our current findings underline a strong correlation between the presence of the bulge and the star formation state of the galaxy. A bulge, if present, is often quiescent, independently of the morphology or the star formation activity of the host galaxy. Additionally, if a galaxy is quiescent, with a large probability, is hosting a bulge. Conversely, if the galaxy has a disky shape is highly probable to be star forming. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 23 pages, 11 figures. Accepted for publication on MNRAS

The half mass radius of MaNGA galaxies: Effect of IMF gradients

Authors: M. Bernardi, R. K. Sheth, H. Dominguez Sanchez, B. Margalef-Bentabol, D. Bizyaev, R. R. Lane

Abstract: Gradients in the stellar populations (SP) of galaxies -- e.g., in age, metallicity, stellar Initial Mass Function (IMF) -- can result in gradients in the stellar mass to light ratio, $M_*/L$. Such gradients imply that the distribution of the stellar mass and light are different. For old SPs, e.g., in early-type galaxies at $z\sim 0$, the $M_*/L$ gradients are weak if driven by variations in age an… ▽ More Gradients in the stellar populations (SP) of galaxies -- e.g., in age, metallicity, stellar Initial Mass Function (IMF) -- can result in gradients in the stellar mass to light ratio, $M_*/L$. Such gradients imply that the distribution of the stellar mass and light are different. For old SPs, e.g., in early-type galaxies at $z\sim 0$, the $M_*/L$ gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger $M_*/L$ in the center increases the estimated total stellar mass ($M_*$) and reduces the scale which contains half this mass ($R_{e,*}$), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H$_β$, $\langle$Fe$\rangle$, [MgFe] and TiO$_{\rm 2SDSS}$ absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in $R_{e,*}$ can be significantly larger than that in $M_*$, especially when the light is more centrally concentrated. The $R_{e,*}-M_*$ correlation which results is offset by 0.3 dex to smaller sizes compared to when these gradients are ignored. Comparisons with `quiescent' galaxies at higher-$z$ must account for evolution in SP gradients (especially age and IMF) and the light profile before drawing conclusions about how $R_{e,*}$ and $M_*$ evolve. The implied merging between higher-$z$ and the present is less contrived if $R_{e,*}/R_e$ at $z\sim 0$ is closer to our IMF-driven gradient calibration than to unity. △ Less

Submitted 15 November, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 16 pages, 15 figures, accepted for publication in MNRAS

The weak dependence of velocity dispersion on disk fractions, mass-to-light ratio and redshift: Implications for galaxy and black hole evolution

Authors: Christopher Marsden, Francesco Shankar, Mariangela Bernardi, Ravi Sheth, Hao Fu, Andrea Lapi

Abstract: Velocity dispersion ($σ$) is a key driver for galaxy structure and evolution. We here present a comprehensive semi-empirical approach to compute $σ$ via detailed Jeans modelling assuming both a constant and scale-dependent mass-to-light ratio $M^*/L$. We compare with a large sample of local galaxies from MaNGA and find that both models can reproduce the Faber-Jackson (FJ) relation and the weak dep… ▽ More Velocity dispersion ($σ$) is a key driver for galaxy structure and evolution. We here present a comprehensive semi-empirical approach to compute $σ$ via detailed Jeans modelling assuming both a constant and scale-dependent mass-to-light ratio $M^*/L$. We compare with a large sample of local galaxies from MaNGA and find that both models can reproduce the Faber-Jackson (FJ) relation and the weak dependence of $σ$ on bulge-to-total ratio $B/T$ (for $B/T\gtrsim 0.25$). The dynamical-to-stellar mass ratio within $R\lesssim R_e$ can be fully accounted for by a gradient in $M^*/L$. We then build velocity dispersion evolutionary tracks $σ_{ap}[M^*,z]$ (within an aperture) along the main progenitor dark matter haloes assigning stellar masses, effective radii and Sersic indices via a variety of abundance matching and empirically motivated relations. We find: 1) clear evidence for downsizing in $σ_{ap}[M^*,z]$ along the progenitor tracks; 2) at fixed stellar mass $σ\propto(1+z)^{0.2-0.3}$ depending on the presence or not of a gradient in $M^*/L$. We extract $σ_{ap}[M^*,z]$ from the TNG50 hydrodynamic simulation and find very similar results to our models with constant $M^*/L$. The increasing dark matter fraction within $R_e$ tends to flatten the $σ_{ap}[M^*,z]$ along the progenitors at $z \gtrsim 1$ in constant $M^*/L$ models, while $σ_{ap}[M^*,z]$ have a steeper evolution in the presence of a stellar gradient. We then show that a combination of mergers and gas accretion are likely responsible for the constant or increasing $σ_{ap}[M^*,z]$ with time. Finally, our $σ_{ap}[M^*,z]$ are consistent with a nearly constant and steep $M_{bh}-σ$ relation at $z\lesssim 2$, with black hole masses derived from the $L_X-M^*$ relation. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: MNRAS accepted, 22 pages, 17 figures

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies… ▽ More This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

SDSS-IV DR17: Final Release of MaNGA PyMorph Photometric and Deep Learning Morphological Catalogs

Authors: H. Domínguez Sánchez, B. Margalef, M. Bernardi, M. Huertas-Company

Abstract: We present the MaNGA PyMorph photometric Value Added Catalogue (MPP-VAC-DR17) and the MaNGA Deep Learning Morphological VAC (MDLM-VAC-DR17) for the final data release of the MaNGA survey, which is part of the SDSS Data Release 17 (DR17). The MPP-VAC-DR17 provides photometric parameters from Sèrsic and Sèrsic+Exponential fits to the 2D surface brightness profiles of the MaNGA DR17 galaxy sample in… ▽ More We present the MaNGA PyMorph photometric Value Added Catalogue (MPP-VAC-DR17) and the MaNGA Deep Learning Morphological VAC (MDLM-VAC-DR17) for the final data release of the MaNGA survey, which is part of the SDSS Data Release 17 (DR17). The MPP-VAC-DR17 provides photometric parameters from Sèrsic and Sèrsic+Exponential fits to the 2D surface brightness profiles of the MaNGA DR17 galaxy sample in the $g$, $r$, and $i$ bands (e.g. total fluxes, half light radii, bulge-disk fractions, ellipticities, position angles, etc.). The MDLM-VAC-DR17 provides Deep Learning-based morphological classifications for the same galaxies. The MDLM-VAC-DR17 includes a number of morphological properties: e.g., a T-Type, a finer separation between elliptical and S0, as well as the identification of edge-on and barred galaxies. While the MPP-VAC-DR17 simply extends the MaNGA PyMorph photometric VAC published in the SDSS Data Release 15 (MPP-VAC-DR15) to now include galaxies which were added to make the final DR17, the MDLM-VAC-DR17 implements some changes and improvements compared to the previous release (MDLM-VAC-DR15): namely, the low-end of the T-Types are better recovered in this new version. The catalogue also includes a separation between Early- or Late-type Galaxies (ETG, LTG), which classifies the two populations in a complementary way to the T-Type, especially at the intermediate types (-1 < T-Type < 2), where the T-Type values show a large scatter. In addition, $k-$fold based uncertainties on the classifications are also provided. To ensure robustness and reliability, we have also visually inspected all the images. We describe the content of the catalogues and show some interesting ways in which they can be combined. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS. Catalogues will be released with SDSS Data Release 17

SDSS-IV MaNGA: drivers of stellar metallicity in nearby galaxies

Authors: Justus Neumann, Daniel Thomas, Claudia Maraston, Daniel Goddard, Jianhui Lian, Lewis Hill, Helena Domínguez Sánchez, Mariangela Bernardi, Berta Margalef-Bentabol, Jorge K. Barrera-Ballesteros, Dmitry Bizyaev, Nicholas F. Boardman, Niv Drory, Joseé G. Fernández-Trincado, Richard Lane

Abstract: The distribution of stellar metallicities within and across galaxies is an excellent relic of the chemical evolution across cosmic time. We present a detailed analysis of spatially resolved stellar populations based on $>2.6$ million spatial bins from 7439 nearby galaxies in the SDSS-IV MaNGA survey. To account for accurate inclination corrections, we derive an equation for morphology dependent de… ▽ More The distribution of stellar metallicities within and across galaxies is an excellent relic of the chemical evolution across cosmic time. We present a detailed analysis of spatially resolved stellar populations based on $>2.6$ million spatial bins from 7439 nearby galaxies in the SDSS-IV MaNGA survey. To account for accurate inclination corrections, we derive an equation for morphology dependent determination of galaxy inclinations. Our study goes beyond the well-known global mass-metallicity relation and radial metallicity gradients by providing a statistically sound exploration of local relations between stellar metallicity $[Z/H]$, stellar surface mass density $Σ_\star$ and galactocentric distance in the global mass-morphology plane. We find a significant resolved mass density-metallicity relation $\rm rΣ_\star ZR$ for galaxies of all types and masses above $10^{9.8}\,\mathrm{M_\odot}$. Different radial distances make an important contribution to the spread of the relation. Particularly, in low and intermediate mass galaxies, we find that at fixed $Σ_\star$ metallicity increases with radius independently of morphology. For high masses, this radial dependence is only observed in high $Σ_\star$ regions of spiral galaxies. This result calls for a driver of metallicity, in addition to $Σ_\star$ that promotes chemical enrichment in the outer parts of galaxies more strongly than in the inner parts. We discuss gas accretion, outflows, recycling and radial migration as possible scenarios. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 16 pages, 11 figures. MNRAS accepted, pending a few minor revisions

The size function of massive satellites from the $R_e-R_h$ and $M_{star}-M_h$ relations: constraining the role of environment

Authors: Lorenzo Zanisi, Francesco Shankar, Mariangela Bernardi, Simona Mei, Marc Huertas-Company

Abstract: In previous work we showed that a semi-empirical model in which galaxies in host dark matter haloes are assigned stellar masses via a stellar mass-halo mass (SMHM) relation and sizes ($R_e$) via a linear and tight $R_e-R_h$ relation, can faithfully reproduce the size function of local SDSS central galaxies and the strong size evolution of massive galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$). In… ▽ More In previous work we showed that a semi-empirical model in which galaxies in host dark matter haloes are assigned stellar masses via a stellar mass-halo mass (SMHM) relation and sizes ($R_e$) via a linear and tight $R_e-R_h$ relation, can faithfully reproduce the size function of local SDSS central galaxies and the strong size evolution of massive galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$). In this third paper of the series, we focus on the population of satellite MGs. We find that without any additional calibration and irrespective of the exact SMHM relation, fraction of quenched galaxies or level of stellar stripping, the same model is able to reproduce the local size function of quiescent satellite MGs in SDSS. In addition, the same model can reproduce the puzzling weak dependence of mean size on host halo mass for both central and satellite galaxies. The model also matches the size function of starforming satellite MGs, after assuming that some of them transform into massive lenticulars in a few Gyr after infalling in the group/cluster environment. However, the vast majority of satellite lenticulars is predicted to form before infall. The $R_e-R_h$ appears to be fundamental to connect galaxies and their host haloes. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. 5 pages, 3 figures. Comments and suggestions still welcome

Identification of high-energy astrophysical point sources via hierarchical Bayesian nonparametric clustering

Authors: Andrea Sottosanti, Mauro Bernardi, Alessandra R. Brazzale, Alex Geringer-Sameth, David C. Stenning, Roberto Trotta, David A. van Dyk

Abstract: The light we receive from distant astrophysical objects carries information about their origins and the physical mechanisms that power them. The study of these signals, however, is complicated by the fact that observations are often a mixture of the light emitted by multiple localized sources situated in a spatially-varying background. A general algorithm to achieve robust and accurate source iden… ▽ More The light we receive from distant astrophysical objects carries information about their origins and the physical mechanisms that power them. The study of these signals, however, is complicated by the fact that observations are often a mixture of the light emitted by multiple localized sources situated in a spatially-varying background. A general algorithm to achieve robust and accurate source identification in this case remains an open question in astrophysics. This paper focuses on high-energy light (such as X-rays and gamma-rays), for which observatories can detect individual photons (quanta of light), measuring their incoming direction, arrival time, and energy. Our proposed Bayesian methodology uses both the spatial and energy information to identify point sources, that is, separate them from the spatially-varying background, to estimate their number, and to compute the posterior probabilities that each photon originated from each identified source. This is accomplished via a Dirichlet process mixture while the background is simultaneously reconstructed via a flexible Bayesian nonparametric model based on B-splines. Our proposed method is validated with a suite of simulation studies and illustrated with an application to a complex region of the sky observed by the \emph{Fermi} Gamma-ray Space Telescope. △ Less

Submitted 26 April, 2021; v1 submitted 23 April, 2021; originally announced April 2021.

Pushing automated morphological classifications to their limits with the Dark Energy Survey

Authors: J. Vega-Ferrero, H. Domínguez Sánchez, M. Bernardi, M. Huertas-Company, R. Morgan, B. Margalef, M. Aguena, S. Allam, J. Annis, S. Avila, D. Bacon, E. Bertin, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, A. Choi, C. Conselice, M. Costanzi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai, I. Ferrero, P. Fosalba , et al. (38 additional authors not shown)

Abstract: We present morphological classifications of $\sim$27 million galaxies from the Dark Energy Survey (DES) Data Release 1 (DR1) using a supervised deep learning algorithm. The classification scheme separates: (a) early-type galaxies (ETGs) from late-types (LTGs); and (b) face-on galaxies from edge-on. Our Convolutional Neural Networks (CNNs) are trained on a small subset of DES objects with previousl… ▽ More We present morphological classifications of $\sim$27 million galaxies from the Dark Energy Survey (DES) Data Release 1 (DR1) using a supervised deep learning algorithm. The classification scheme separates: (a) early-type galaxies (ETGs) from late-types (LTGs); and (b) face-on galaxies from edge-on. Our Convolutional Neural Networks (CNNs) are trained on a small subset of DES objects with previously known classifications. These typically have $\mathrm{m}_r \lesssim 17.7~\mathrm{mag}$; we model fainter objects to $\mathrm{m}_r < 21.5$ mag by simulating what the brighter objects with well determined classifications would look like if they were at higher redshifts. The CNNs reach 97\% accuracy to $\mathrm{m}_r<21.5$ on their training sets, suggesting that they are able to recover features more accurately than the human eye. We then used the trained CNNs to classify the vast majority of the other DES images. The final catalog comprises five independent CNN predictions for each classification scheme, helping to determine if the CNN predictions are robust or not. We obtain secure classifications for $\sim$ 87\% and 73\% of the catalog for the ETG vs. LTG and edge-on vs. face-on models, respectively. Combining the two classifications (a) and (b) helps to increase the purity of the ETG sample and to identify edge-on lenticular galaxies (as ETGs with high ellipticity). Where a comparison is possible, our classifications correlate very well with Sérsic index (\textit{n}), ellipticity ($ε$) and spectral type, even for the fainter galaxies. This is the largest multi-band catalog of automated galaxy morphologies to date. △ Less

Submitted 9 March, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: Accepted for publication in MNRAS (2021 February 22); 17 pages, 16 figures

On the Presence of a Universal Acceleration Scale in Elliptical Galaxies

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Helena Dominguez Sanchez, Ravi K. Sheth

Abstract: Dark matter phenomena in rotationally supported galaxies exhibit a characteristic acceleration scale of $g_\dagger \approx 1.2\times 10^{-10}$ m s$^{-2}$. Whether this acceleration is a manifestation of a universal scale, or merely an emergent property with an intrinsic scatter, has been debated in the literature. Here we investigate whether a universal acceleration scale exists in dispersion-supp… ▽ More Dark matter phenomena in rotationally supported galaxies exhibit a characteristic acceleration scale of $g_\dagger \approx 1.2\times 10^{-10}$ m s$^{-2}$. Whether this acceleration is a manifestation of a universal scale, or merely an emergent property with an intrinsic scatter, has been debated in the literature. Here we investigate whether a universal acceleration scale exists in dispersion-supported galaxies using two uniform sets of integral field spectroscopy (IFS) data from SDSS-IV MaNGA and ATLAS$^{\rm 3D}$. We apply the spherical Jeans equation to 15 MaNGA and 4 ATLAS$^{\rm 3D}$ slow-rotator E0 (i.e., nearly spherical) galaxies. Velocity dispersion profiles for these galaxies are well determined with observational errors under control. Bayesian inference indicates that all 19 galaxies are consistent with a universal acceleration of $g_\dagger=1.5_{-0.6}^{+0.9}\times 10^{-10}$ m s$^{-2}$. Moreover, all 387 data points from the radial bins of the velocity dispersion profiles are consistent with a universal relation between the radial acceleration traced by dynamics and that predicted by the observed distribution of baryons. This universality remains if we include 12 additional non-E0 slow-rotator elliptical galaxies from ATLAS$^{\rm 3D}$. Finally, the universal acceleration from MaNGA and ATLAS$^{\rm 3D}$ is consistent with that for rotationally supported galaxies, so our results support the view that dark matter phenomenology in galaxies involves a universal acceleration scale. △ Less

Submitted 17 November, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: Published in ApJ Letters

Journal ref: 2020 November 10, 903, L31

Galaxy properties as revealed by MaNGA. III. Kinematic profiles and stellar population gradients in S0s

Authors: H. Domínguez Sánchez, M. Bernardi, F. Nikakhtar, B. Margalef-Bentabol, R. K. Sheth

Abstract: This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $α$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by m… ▽ More This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $α$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (L$_r$) and central velocity dispersion ($σ_0$). We find that: 1) There is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}M_\odot$ show stronger velocity dispersion and age gradients (age and $σ_r$ decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times10^{11}M_\odot$ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed $σ_0$, more luminous S0s and E-FRs are younger, more metal rich and less $α$-enhanced. Evidently for these galaxies, the usual statement that 'massive galaxies are older' is not true if $σ_0$ is held fixed. △ Less

Submitted 15 May, 2020; originally announced May 2020.

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

The Stellar Mass Fundamental Plane: The virial relation and a very thin plane for slow-rotators

Authors: M. Bernardi, H. Domínguez Sanchez, B. Margalef-Bentabol, F. Nikakhtar, R. K. Sheth

Abstract: Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $σ_e$. Since $I_e$ and $σ_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $σ_e$ can be… ▽ More Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $σ_e$. Since $I_e$ and $σ_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $σ_e$ can be used to estimate sizes $R_e$. We show that: 1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in $R_e$ is $\sim 16\%$, of which only $\sim 9\%$ is intrinsic. Removing galaxies with $M_*<10^{11}M_\odot$ further reduces the observed scatter to $\sim 13\%$ ($\sim 4\%$ intrinsic). The observed scatter increases to the $\sim 25\%$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only $5\%$ orthogonal to the plane. 2) The structure within the FP is most easily understood as arising from the fact that $I_e$ and $σ_e$ are nearly independent, whereas the $R_e-I_e$ and $R_e-σ_e$ correlations are nearly equal and opposite. 3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be `tilted'. If we multiply $I_e$ by the global stellar mass-to-light ratio $M_*/L$ and we account for non-homology across the population by using Sérsic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for $M_*/L$ gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 13 pages, 9 figures, 3 tables, accepted for publication in MNRAS

Detecting outliers in astronomical images with deep generative networks

Authors: Berta Margalef-Bentabol, Marc Huertas-Company, Tom Charnock, Carla Margalef-Bentabol, Mariangela Bernardi, Yohan Dubois, Kate Storey-Fisher, Lorenzo Zanis

Abstract: With the advent of future big-data surveys, automated tools for unsupervised discovery are becoming ever more necessary. In this work, we explore the ability of deep generative networks for detecting outliers in astronomical imaging datasets. The main advantage of such generative models is that they are able to learn complex representations directly from the pixel space. Therefore, these methods e… ▽ More With the advent of future big-data surveys, automated tools for unsupervised discovery are becoming ever more necessary. In this work, we explore the ability of deep generative networks for detecting outliers in astronomical imaging datasets. The main advantage of such generative models is that they are able to learn complex representations directly from the pixel space. Therefore, these methods enable us to look for subtle morphological deviations which are typically missed by more traditional moment-based approaches. We use a generative model to learn a representation of expected data defined by the training set and then look for deviations from the learned representation by looking for the best reconstruction of a given object. In this first proof-of-concept work, we apply our method to two different test cases. We first show that from a set of simulated galaxies, we are able to detect $\sim90\%$ of merging galaxies if we train our network only with a sample of isolated ones. We then explore how the presented approach can be used to compare observations and hydrodynamic simulations by identifying observed galaxies not well represented in the models. △ Less

Submitted 20 March, 2020; v1 submitted 18 March, 2020; originally announced March 2020.

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

Authors: Francesco Shankar, David H. Weinberg, Christopher Marsden, Philip J. Grylls, Mariangela Bernardi, Guang Yang, Benjamin Moster, Rosamaria Carraro, David M. Alexander, Viola Allevato, Tonima T. Ananna, Angela Bongiorno, Giorgio Calderone, Francesca Civano, Emanuele Daddi, Ivan Delvecchio, Federica Duras, Fabio La Franca, Andrea Lapi, Youjun Lu, Nicola Menci, Mar Mezcua, Federica Ricci, Giulia Rodighiero, Ravi K. Sheth , et al. (3 additional authors not shown)

Abstract: The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises f… ▽ More The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 12 pages, 8 Figures. MNRAS, accepted. A discussion around the points raised by arXiv:1909.10821 is included in the Appendix

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

Authors: Romina Ahumada, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Riccardo Arcodia, Eric Armengaud, Marie Aubert, Santiago Avila, Vladimir Avila-Reese, Carles Badenes, Christophe Balland, Kat Barger, Jorge K. Barrera-Ballesteros, Sarbani Basu, Julian Bautista, Rachael L. Beaton, Timothy C. Beers, B. Izamar T. Benavides, Chad F. Bender, Mariangela Bernardi, Matthew Bershady, Florian Beutler , et al. (289 additional authors not shown)

Abstract: This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the… ▽ More This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17). △ Less

Submitted 11 May, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: DR16 release: Monday Dec 9th 2019. This is the alphabetical order SDSS-IV collaboration data release paper. 25 pages, 6 figures, accepted by ApJS on 11th May 2020. Minor changes clarify or improve text and figures relative to v1

Galaxy sizes and the galaxy-halo connection -- I: the remarkable tightness of the size distributions

Authors: Lorenzo Zanisi, Francesco Shankar, Andrea Lapi, Nicola Menci, Mariangela Bernardi, Christopher Duckworth, Marc Huertas-Company, Philip Grylls, Paolo Salucci

Abstract: The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the n… ▽ More The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 22 pages, 14 figures

Constraining black hole-galaxy scaling relations from the large-scale clustering of Active Galactic Nuclei and implied mean radiative efficiency

Authors: Francesco Shankar, Viola Allevato, Mariangela Bernardi, Christopher Marsden, Andrea Lapi, Nicola Menci, Philip J. Grylls, Mirko Krumpe, Lorenzo Zanisi, Federica Ricci, Fabio La Franca, Ranieri D. Baldi, Jorge Moreno, Ravi K. Sheth

Abstract: A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demo… ▽ More A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 31 pages, 7 Figures. Accepted

Galaxy properties as revealed by MaNGA II. Differences in stellar populations of slow and fast rotator ellipticals and dependence on environment

Authors: M. Bernardi, H. Domínguez Sánchez, J. R. Brownstein, N. Drory, R. K. Sheth

Abstract: We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $α$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $σ_0$. This explains why when one combines SRs and FRs, objects which are… ▽ More We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $α$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $σ_0$. This explains why when one combines SRs and FRs, objects which are small for their $L_r$ and $σ_0$ tend to be younger. Their SP gradients are also different. 2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. 3) SRs with $σ_0\le 250$ km s$^{-1}$ are remarkably homogeneous within $\sim R_e$: they are old, $α$-enhanced and only slightly super-solar in metallicity. These SRs show no gradients in age and $M_*/L_r$, negative gradients in metallicity, and slightly positive gradients in [$α$/Fe] (the latter are model dependent). SRs with $σ_0\ge 250$ km $s^{-1}$ are slightly younger and more metal rich, contradicting previous work suggesting that age increases with $σ_0$. They also show larger $M_*/L_r$ gradients. 4) Self-consistently accounting for $M_*/L$ gradients yields $M_{\rm dyn}\approx M_*$ because gradients reduce $M_{\rm dyn}$ by $\sim 0.2$ dex while only slightly increasing the $M_*$ inferred using a Kroupa (not Salpeter) IMF. 5) The FR population all but disappears above $M_*\ge 3\times 10^{11}M_\odot$; this is the same scale at which the size-mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale which correlates with the environment and above which mergers matter. △ Less

Submitted 27 August, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

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

Galaxy properties as revealed by MaNGA. I. Constraints on Initial Mass Function and M$_{*}$/L gradients in ellipticals

Authors: H. Domínguez Sánchez, M. Bernardi, J. R. Brownstein, N. Drory, R. K. Sheth

Abstract: We estimate ages, metallicities, $α$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distanc… ▽ More We estimate ages, metallicities, $α$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distance. This paper is the first of a series and is based on a sample of pure elliptical galaxies at z $\le$ 0.08. We show that the properties of the inner regions of Es with the largest luminosity (L$_r$) and central velocity dispersion ($σ_0$) are consistent with those associated with the commonly used Salpeter IMF, whereas a Kroupa-like IMF is a better description at $\sim$ 0.8R/Re (assuming [Ti/Fe] variations are limited). For these galaxies the stellar mass-to-light ratio decreases at most by a factor of 2 from the central regions to Re. In contrast, for lower L$_r$ and $σ_0$ galaxies, the IMF is shallower and M$_{*}$/L$_r$ in the central regions is similar to the outskirts. Although a factor of 2 is smaller than previous reports based on a handful of galaxies, it is still large enough to matter for dynamical mass estimates. Accounting self-consistently for these gradients when estimating both M$_{*}$ and M$_{dyn}$ brings the two into good agreement: gradients reduce M$_{dyn}$ by $\sim$ 0.2 dex while only slightly increasing the M$_{*}$ inferred using a Kroupa IMF. This is a different resolution of the M$_{*}$-M$_{dyn}$ discrepancy than has been followed in the recent literature where M$_{*}$ of massive galaxies is increased by adopting a Salpeter IMF while leaving Mdyn unchanged. A companion paper discusses how stellar population differences are even more pronounced if one separates slow from fast rotators. △ Less

Submitted 27 August, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 22 pages, 29 figures, accepted for publication in MNRAS

The Hubble Sequence at $z\sim0$ in the IllustrisTNG simulation with deep learning

Authors: M. Huertas-Company, V. Rodriguez-Gomez, D. Nelson, A. Pillepich, M. Bernardi, H. Domínguez-Sánchez, S. Genel, R. Pakmor, G. F. Snyder, M. Vogelsberger

Abstract: We analyze the optical morphologies of galaxies in the IllustrisTNG simulation at $z\sim0$ with a Convolutional Neural Network trained on visual morphologies in the Sloan Digital Sky Survey. We generate mock SDSS images of a mass complete sample of $\sim12,000$ galaxies in the simulation using the radiative transfer code SKIRT and include PSF and noise to match the SDSS r-band properties. The imag… ▽ More We analyze the optical morphologies of galaxies in the IllustrisTNG simulation at $z\sim0$ with a Convolutional Neural Network trained on visual morphologies in the Sloan Digital Sky Survey. We generate mock SDSS images of a mass complete sample of $\sim12,000$ galaxies in the simulation using the radiative transfer code SKIRT and include PSF and noise to match the SDSS r-band properties. The images are then processed through the exact same neural network used to estimate SDSS morphologies to classify simulated galaxies in four morphological classes (E, S0/a, Sab, Scd). The CNN model finds that $\sim95\%$ of the simulated galaxies fall in one the four main classes with high confidence. The mass-size relations of the simulated galaxies divided by morphological type also reproduce well the slope and the normalization of observed relations which confirms the realism of optical morphologies in the TNG suite. However, the Stellar Mass Functions decomposed into different morphologies still show significant discrepancies with observations both at the low and high mass end. We find that the high mass end of the SMF is dominated in TNG by massive disk galaxies while early-type galaxies dominate in the observations according to the CNN classifications. The present work highlights the importance of detailed comparisons between observations and simulations in comparable conditions. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: submitted to MNRAS, comments welcome

Modeling Nearly Spherical Pure-Bulge Galaxies with a Stellar Mass-to-Light Ratio Gradient under the $Λ$CDM and MOND Paradigms: II. The Orbital Anisotropy of Slow Rotators within the Effective Radius

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Ravi K. Sheth

Abstract: We investigate the anisotropy of the stellar velocity dispersions within the effective radius, $R_{\rm e}$, in 24 ATLAS$^{\rm 3D}$ pure-bulge galaxies, 16 of which are kinematic slow rotators (SRs). We allow the spherical anisotropy parameter $β$ to be radially varying and allow a radial gradient in the stellar mass-to-light ratio ($M_\star/L$) through the parameter $K$ introduced earlier. The med… ▽ More We investigate the anisotropy of the stellar velocity dispersions within the effective radius, $R_{\rm e}$, in 24 ATLAS$^{\rm 3D}$ pure-bulge galaxies, 16 of which are kinematic slow rotators (SRs). We allow the spherical anisotropy parameter $β$ to be radially varying and allow a radial gradient in the stellar mass-to-light ratio ($M_\star/L$) through the parameter $K$ introduced earlier. The median anisotropy for SRs depends on $K$ as follows: $\langleβ_{\rm m}\rangle = a + b K$ with $a=0.19\pm 0.05$, $b=-0.13\pm 0.07$ ($Λ$CDM) or $a=0.21\pm 0.05$, $b=-0.26\pm 0.08$ (MOND), where $β_{\rm m}$ refers to the radially averaged quantity. Under the $Λ$CDM paradigm this scaling is tied to a scaling of $\langle f_{\rm DM}\rangle = (0.16\pm 0.03) +(0.31\pm 0.06) K$, where $f_{\rm DM}$ refers to the DM fraction within a sphere of $r=R_{\rm e}$. For $K=0$ (constant $M_\star/L$), we obtain radially biased results with $\langleβ_{\rm m}\rangle \approx 0.2$ consistent with previous results. However, marginalizing over $0 < K < 1.5$ yields $\langleβ_{\rm m}\rangle = 0.06 ^{+0.11}_{-0.14}$ with $\langle f_{\rm DM}\rangle = 0.35 \pm 0.08$: isotropy is preferred. This isotropy hides the fact that $β_{\rm m}$ is correlated with kinematic features such as counter rotating cores (CRCs), kinematically distinct cores (KDCs), and low-level velocities (LVs): SRs with LVs are likely to be radially biased while SRs with CRCs are likely to be tangentially biased, and SRs with KDCs are intermediate. Existing cosmological simulations allow us to understand these results qualitatively in terms of their dynamical structures and formation histories although there exist quantitative tensions. More realistic cosmological simulations, particularly allowing for $M_\star/L$ gradients, may be required to better understand SRs. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 11 figures, 3 tables, ApJ, accepted for publication

The structural properties of classical bulges and discs from z~2

Authors: Paola Dimauro, Marc Huertas-Company, Emanuele Daddi, Pablo G. Pérez-González, Mariangela Bernardi, Fernando Caro, Andrea Cattaneo, Boris Häußler, Ulrike Kuchner, Francesco Shankar, Guillermo Barro, Fernando Buitrago, Sandra M. Faber, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Simona Mei, Reynier Peletier, Joel Primack, Aldo Rodriguez-Puebla, Mara Salvato, Diego Tuccillo

Abstract: We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component Sérsic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxi… ▽ More We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component Sérsic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers. △ Less

Submitted 13 February, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 20 pages, 15 figures, Submitted to MNRAS

Black hole scaling relations of active and quiescent galaxies: Addressing selection effects and constraining virial factors

Authors: Francesco Shankar, Mariangela Bernardi, Kayleigh Richardson, Christopher Marsden, Ravi K. Sheth, Viola Allevato, Luca Graziani, Mar Mezcua, Federica Ricci, Samantha J. Penny, Fabio La Franca, Fabio Pacucci

Abstract: Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-… ▽ More Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 15 pages, 5 Figures. MNRAS, accepted

The Fifteenth Data Release of the Sloan Digital Sky Surveys: First Release of MaNGA Derived Quantities, Data Visualization Tools and Stellar Library

Authors: D. S. Aguado, Romina Ahumada, Andres Almeida, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Sandro Barboza Rembold, Kat Barger, Jorge Barrera-Ballesteros, Dominic Bates, Julian Bautista, Rachael L. Beaton, Timothy C. Beers, Francesco Belfiore, Mariangela Bernardi, Matthew Bershady, Florian Beutler, Jonathan Bird, Dmitry Bizyaev , et al. (209 additional authors not shown)

Abstract: Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar… ▽ More Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data. △ Less

Submitted 10 December, 2018; v1 submitted 6 December, 2018; originally announced December 2018.

Comments: Paper to accompany DR15. 25 pages, 6 figures. Accepted for publication in ApJSS. The two papers on the MaNGA Data Analysis Pipeline (DAP, Westfall et al. and Belfiore et al., see Section 4.1.2), and the paper on Marvin (Cherinka et al., see Section 4.2) have been submitted for collaboration review and will be posted to arXiv in due course. v2 fixes some broken URLs in the PDF

A Statistical Semi-Empirical Model: Satellite galaxies in Groups and Clusters

Authors: P. J. Grylls, F. Shankar, L. Zanisi, M. Bernardi

Abstract: We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass func… ▽ More We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass functions the satellite mass function is progressively built in time via abundance matching techniques constrained by number densities of centrals in the local Universe. By enforcing dynamical merging timescales as predicted by high-resolution N-body simulations, we obtain satellite distributions as a function of stellar mass and halo mass consistent with current data. We show that stellar stripping, star formation, and quenching play all a secondary role in setting the number densities of massive satellites above $M_*\gtrsim 3\times 10^{10}\, M_{\odot}$. We further show that observed star formation rates used in our empirical model over predict low-mass satellites below $M_*\lesssim 3\times 10^{10}\, M_{\odot}$, whereas, star formation rates derived from a continuity equation approach yield the correct abundances similar to previous results for centrals. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: 21 pages, 17 Figures. MNRAS, in press

SDSS-IV MaNGA PyMorph Photometric and Deep Learning Morphological Catalogs and implications for bulge properties and stellar angular momentum

Authors: J. -L. Fischer, H. Domínguez Sánchez, M. Bernardi

Abstract: We describe the SDSS-IV MaNGA PyMorph Photometric (MPP-VAC) and MaNGA Deep Learning Morphology (MDLM-VAC) Value Added Catalogs. The MPP-VAC provides photometric parameters from Sérsic and Sérsic+Exponential fits to the 2D surface brightness profiles of the MaNGA DR15 galaxy sample. Compared to previous PyMorph analyses of SDSS imaging, our analysis of the MaNGA DR15 incorporates three improvements… ▽ More We describe the SDSS-IV MaNGA PyMorph Photometric (MPP-VAC) and MaNGA Deep Learning Morphology (MDLM-VAC) Value Added Catalogs. The MPP-VAC provides photometric parameters from Sérsic and Sérsic+Exponential fits to the 2D surface brightness profiles of the MaNGA DR15 galaxy sample. Compared to previous PyMorph analyses of SDSS imaging, our analysis of the MaNGA DR15 incorporates three improvements: the most recent SDSS images; modified criteria for determining bulge-to-disk decompositions; and the fits in MPP-VAC have been eye-balled, and re-fit if necessary, for additional reliability. A companion catalog, the MDLM-VAC, provides Deep Learning-based morphological classifications for the same galaxies. The MDLM-VAC includes a number of morphological properties (e.g., a TType, and a finer separation between elliptical and S0 galaxies). Combining the MPP- and MDLM-VACs allows to show that the MDLM morphological classifications are more reliable than previous work. It also shows that single-Sérsic fits to late- and early-type galaxies are likely to return Sérsic indices of $n \le 2$ and $\ge 4$, respectively, and this correlation between $n$ and morphology extends to the bulge component as well. While the former is well-known, the latter contradicts some recent work suggesting little correlation between $n$-bulge and morphology. Combining both VACs with MaNGA's spatially resolved spectroscopy allows us to study how the stellar angular momentum depends on morphological type. We find correlations between stellar kinematics, photometric properties, and morphological type even though the spectroscopic data played no role in the construction of the MPP- and MDLM-VACs. △ Less

Submitted 11 December, 2018; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: 22 pages, 30 figures and 4 tables. Accepted for publication in MNRAS. The MPP-VAC and MDLM-VAC are available online at https://www.sdss.org/dr15/data_access/value-added-catalogs/manga-pymorph-dr15-photometric-catalog and https://www.sdss.org/dr15/data_access/value-added-catalogs/manga-morphology-deep-learning-dr15-catalog

SDSS-IV MaNGA: The Formation Sequence of S0 Galaxies

Authors: Amelia Fraser-McKelvie, Alfonso Aragón-Salamanca, Michael Merrifield, Martha Tabor, Mariangela Bernardi, Niv Drory, Taniya Parikh, Maria Argudo-Fernández

Abstract: Gas stripping of spiral galaxies or mergers are thought to be the formation mechanisms of lenticular galaxies. In order to determine the conditions in which each scenario dominates, we derive stellar populations of both the bulge and disk regions of 279 lenticular galaxies in the MaNGA survey. We find a clear bimodality in stellar age and metallicity within the population of S0s and this is strong… ▽ More Gas stripping of spiral galaxies or mergers are thought to be the formation mechanisms of lenticular galaxies. In order to determine the conditions in which each scenario dominates, we derive stellar populations of both the bulge and disk regions of 279 lenticular galaxies in the MaNGA survey. We find a clear bimodality in stellar age and metallicity within the population of S0s and this is strongly correlated with stellar mass. Old and metal-rich bulges and disks belong to massive galaxies, and young and metal-poor bulges and disks are hosted by low-mass galaxies. From this we conclude that the bulges and disks are co-evolving. When the bulge and disk stellar ages are compared, we find that the bulge is almost always older than the disk for massive galaxies ($\textrm{M}_{\star} > 10^{10}~\textrm{M}_{\odot}$). The opposite is true for lower mass galaxies. We conclude that we see two separate populations of lenticular galaxies. The old, massive, and metal-rich population possess bulges that are predominantly older than their disks, which we speculate may have been caused by morphological or inside-out quenching. In contrast, the less massive and more metal-poor population have bulges with more recent star formation than their disks. We postulate they may be undergoing bulge rejuvenation (or disk fading), or compaction. Environment doesn't play a distinct role in the properties of either population. Our findings give weight to the notion that while the faded spiral scenario likely formed low-mass S0s, other processes, such as mergers, may be responsible for high-mass S0s. △ Less

Submitted 12 September, 2018; originally announced September 2018.

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

Transfer learning for galaxy morphology from one survey to another

Authors: H. Domínguez Sánchez, M. Huertas-Company, M. Bernardi, S. Kaviraj, J. L. Fischer, T. M. C. Abbott, F. B. Abdalla, J. Annis, S. Avila, D. Brooks, E. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, C. E. Cunha, C. B. D'Andrea, L. N. da Costa, C. Davis, J. De Vicente, P. Doel, A. E. Evrard, P. Fosalba, J. Frieman, J. García-Bellido, E. Gaztanaga , et al. (34 additional authors not shown)

Abstract: Deep Learning (DL) algorithms for morphological classification of galaxies have proven very successful, mimicking (or even improving) visual classifications. However, these algorithms rely on large training samples of labelled galaxies (typically thousands of them). A key question for using DL classifications in future Big Data surveys is how much of the knowledge acquired from an existing survey… ▽ More Deep Learning (DL) algorithms for morphological classification of galaxies have proven very successful, mimicking (or even improving) visual classifications. However, these algorithms rely on large training samples of labelled galaxies (typically thousands of them). A key question for using DL classifications in future Big Data surveys is how much of the knowledge acquired from an existing survey can be exported to a new dataset, i.e. if the features learned by the machines are meaningful for different data. We test the performance of DL models, trained with Sloan Digital Sky Survey (SDSS) data, on Dark Energy survey (DES) using images for a sample of $\sim$5000 galaxies with a similar redshift distribution to SDSS. Applying the models directly to DES data provides a reasonable global accuracy ($\sim$ 90%), but small completeness and purity values. A fast domain adaptation step, consisting in a further training with a small DES sample of galaxies ($\sim$500-300), is enough for obtaining an accuracy > 95% and a significant improvement in the completeness and purity values. This demonstrates that, once trained with a particular dataset, machines can quickly adapt to new instrument characteristics (e.g., PSF, seeing, depth), reducing by almost one order of magnitude the necessary training sample for morphological classification. Redshift evolution effects or significant depth differences are not taken into account in this study. △ Less

Submitted 16 January, 2019; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS

Deep Learning Identifies High-z Galaxies in a Central Blue Nugget Phase in a Characteristic Mass Range

Authors: M. Huertas-Company, J. R. Primack, A. Dekel, D. C. Koo, S. Lapiner, D. Ceverino, R. C. Simons, G. F. Snyder, M. Bernardi, Z. Chen, H. Domínguez-Sánchez, Z. Chen, C. T. Lee, B. Margalef-Bentabol, D. Tuccillo

Abstract: We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Ne… ▽ More We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: Accepted for publication in ApJ

Modeling Nearly Spherical Pure-bulge Galaxies with a Stellar Mass-to-Light Ratio Gradient under the $Λ$CDM and MOND Paradigms: I. Methodology, Dynamical Stellar Mass, and Fundamental Mass Plane

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Ravi K. Sheth

Abstract: We carry out spherical Jeans modeling of nearly round pure-bulge galaxies selected from the ATLAS$^{\rm 3D}$ sample. Our modeling allows for gradients in the stellar mass-to-light ratio ($M_\star/L$) through analytic prescriptions parameterized with a `gradient strength' $K$ introduced to accommodate any viable gradient. We use a generalized Osipkov-Merritt model for the velocity dispersion (VD) a… ▽ More We carry out spherical Jeans modeling of nearly round pure-bulge galaxies selected from the ATLAS$^{\rm 3D}$ sample. Our modeling allows for gradients in the stellar mass-to-light ratio ($M_\star/L$) through analytic prescriptions parameterized with a `gradient strength' $K$ introduced to accommodate any viable gradient. We use a generalized Osipkov-Merritt model for the velocity dispersion (VD) anisotropy. We produce Monte Carlo sets of models based on the stellar VD profiles under both the $Λ$CDM and MOND paradigms. Here, we describe the galaxy data, the empirical inputs, and the modeling procedures of obtaining the Monte Carlo sets. We then present the projected dynamical stellar mass, $M_{\rm \star e}$, within the effective radius $R_{\rm e}$, and the fundamental mass plane (FMP) as a function of $K$. We find the scaling of the $K$-dependent mass with respect to the ATLAS$^{\rm 3D}$ reported mass as: $\log_{10} \left[M_{\star{\rm e}}(K)/M_{\star{\rm e}}^{\rm A3D} \right] = a' + b' K$ with $a'=-0.019\pm 0.012$ and $b'=-0.18\pm 0.02$ ($Λ$CDM), or $a'=-0.023\pm 0.014$ and $b'=-0.23\pm 0.03$ (MOND), for $0\le K < 1.5$. The FMP has coefficients consistent with the virial expectation and only the zero-point scales with $K$. The median value of $K$ for the ATLAS$^{\rm 3D}$ galaxies is $\langle K\rangle =0.53^{+0.05}_{-0.04}$. We perform a similar analysis of the much larger SDSS DR7 spectroscopic sample. In this case, only the VD within a single aperture is available, so we impose the additional requirement that the VD slope be similar to that in the ATLAS$^{\rm 3D}$ galaxies. Our analysis of the SDSS galaxies suggests a positive correlation of $K$ with stellar mass. △ Less

Submitted 31 May, 2018; v1 submitted 30 March, 2018; originally announced April 2018.

Comments: ApJ, in press (32 pages in preprint style including 22 figures; Fig 21 corrected)

A catalog of polychromatic bulge-disk decompositions of ~ 17.600 galaxies in CANDELS

Authors: Paola Dimauro, Marc Huertas-Company, Emanuele Daddi, Pablo G. Pérez-González, Mariangela Bernardi, Guillermo Barro, Fernando Buitrago, Fernando Caro, Andrea Cattaneo, Helena Dominguez-Sánchez, Sandra M. Faber, Boris Häußler, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Christoph T. Lee, Simona Mei, Berta Margalef-Bentabol, Joel Primack, Aldo Rodriguez-Puebla, Mara Salvato, Francesco Shankar, Diego Tuccillo

Abstract: Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge… ▽ More Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works. △ Less

Submitted 30 March, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: 18 pages, 18 figures, 3 table, submitted to MNRAS after the first referee's report

$M_*/L$ gradients driven by IMF variation: Large impact on dynamical stellar mass estimates

Authors: M. Bernardi, R. K. Sheth, H. Dominguez-Sanchez, J. -L. Fischer, K. -H. Chae, M. Huertas-Company, F. Shankar

Abstract: Within a galaxy the stellar mass-to-light ratio $Υ_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $Υ_*$ gradients than if the IMF is held fixed. If $Υ_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as… ▽ More Within a galaxy the stellar mass-to-light ratio $Υ_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $Υ_*$ gradients than if the IMF is held fixed. If $Υ_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $Υ_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $Υ_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $Υ_*$ gradients changes the high-mass slope of the stellar mass function $φ(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $Υ_*$ gradients in larger samples. △ Less

Submitted 21 March, 2018; v1 submitted 14 December, 2017; originally announced December 2017.

Comments: 13 pages, 7 figures, MNRAS in press

Revisiting the bulge-halo conspiracy II: Towards explaining its puzzling dependence on redshift

Authors: Francesco Shankar, Alessandro Sonnenfeld, Philip Grylls, Lorenzo Zanisi, Carlo Nipoti, Kyu-Hyun Chae, Mariangela Bernardi, Carlo Enrico Petrillo, Marc Huertas-Company, Gary A. Mamon, Stewart Buchan

Abstract: We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and Sérsic. We assign dark matter haloes to galaxies via a… ▽ More We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and Sérsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be Sérsic and the input Sérsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing Sérsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution. △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: 14 pages, 8 figures. MNRAS, accepted. Main result of the paper in Figure 2

Improving galaxy morphologies for SDSS with Deep Learning

Authors: H. Domínguez Sánchez, M. Huertas-Company, M. Bernardi, D. Tuccillo, J. L. Fischer

Abstract: We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are ob… ▽ More We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are obtained with Deep Learning algorithms using Convolutional Neural Networks (CNNs). We use two visual classification catalogues, GZ2 and Nair & Abraham (2010), for training CNNs with colour images in order to obtain T-Types and a series of GZ2 type questions (disk/features, edge-on galaxies, bar signature, bulge prominence, roundness and mergers). We also provide an additional probability enabling a separation between pure elliptical (E) from S0, where the T-Type model is not so efficient. For the T-Type, our results show smaller offset and scatter than previous models trained with support vector machines. For the GZ2 type questions, our models have large accuracy (> 97\%), precision and recall values (> 90\%) when applied to a test sample with the same characteristics as the one used for training. The catalogue is publicly released with the paper. △ Less

Submitted 7 February, 2018; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 18 pages, 21 figures; Accepted for publication in MNRAS

Stellar mass functions and implications for a variable IMF

Authors: M. Bernardi, R. K. Sheth, J. -L. Fischer, A. Meert, K. -H. Chae, H. Dominguez-Sanchez, M. Huertas-Company, F. Shankar, V. Vikram

Abstract: Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $σ_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with… ▽ More Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $σ_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $σ$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm α_{JAM}}$, from $M_{*}$ and $σ_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,σ_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(σ_R/σ_e) = (R/R_e)^{-γ(n)}$. The resulting stellar mass functions, based on $M_*^{\rm α_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm α_{JAM}}/M_* - σ_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $φ(M_*^{\rm α_{JAM}})$ is an estimate of the stellar mass function in which $σ_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $σ_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form. △ Less

Submitted 15 December, 2017; v1 submitted 19 October, 2017; originally announced October 2017.

Comments: 17 pages, 19 figures, 4 tables. Accepted for publication by MNRAS. Tables 1, C1 and C2 are provided as ancillary files

The Fourteenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the extended Baryon Oscillation Spectroscopic Survey and from the second phase of the Apache Point Observatory Galactic Evolution Experiment

Authors: Bela Abolfathi, D. S. Aguado, Gabriela Aguilar, Carlos Allende Prieto, Andres Almeida, Tonima Tasnim Ananna, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Eric Armengaud, Metin Ata, Eric Aubourg, Vladimir Avila-Reese, Carles Badenes, Stephen Bailey, Christophe Balland, Kathleen A. Barger, Jorge Barrera-Ballesteros, Curtis Bartosz, Fabienne Bastien, Dominic Bates, Falk Baumgarten , et al. (323 additional authors not shown)

Abstract: The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulativ… ▽ More The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V. △ Less

Submitted 6 May, 2018; v1 submitted 28 July, 2017; originally announced July 2017.

Comments: SDSS-IV collaboration alphabetical author data release paper. DR14 happened on 31st July 2017. 19 pages, 5 figures. Accepted by ApJS on 28th Nov 2017 (this is the "post-print" and "post-proofs" version; minor corrections only from v1, and most of errors found in proofs corrected)

Radial Acceleration Relation between Baryons and Dark or Phantom Matter in the Super-critical Acceleration Regime of Nearly Spherical Galaxies

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Ravi K. Sheth, In-Taek Gong

Abstract: The central regions of nearby elliptical galaxies are dominated by baryons (stars) and provide interesting laboratories for studying the radial acceleration relation (RAR). We carry out exploratory analyses and discuss the possibility of constraining the RAR in the super-critical acceleration range $(10^{-9.5},\hspace{1ex}10^{-8})$~${\rm m}~{\rm s}^{-2}$ using a sample of nearly round pure-bulge (… ▽ More The central regions of nearby elliptical galaxies are dominated by baryons (stars) and provide interesting laboratories for studying the radial acceleration relation (RAR). We carry out exploratory analyses and discuss the possibility of constraining the RAR in the super-critical acceleration range $(10^{-9.5},\hspace{1ex}10^{-8})$~${\rm m}~{\rm s}^{-2}$ using a sample of nearly round pure-bulge (spheroidal, dispersion-dominated) galaxies including 24 ATLAS$^{\rm 3D}$ galaxies and 4201 SDSS galaxies covering a wide range of masses, sizes and luminosity density profiles. We consider a range of current possibilities for the stellar mass-to-light ratio ($M_\star/L$), its gradient and dark or phantom matter (DM/PM) halo profiles. We obtain the probability density functions (PDFs) of the parameters of the considered models via Bayesian inference based on spherical Jeans Monte Carlo modeling of the observed velocity dispersions. We then constrain the DM/PM-to-baryon acceleration ratio $a_{\rm X}/a_{\rm B}$ from the PDFs. Unless we ignore observed radial gradients in $M_\star/L$, or assume unreasonably strong gradients, marginalization over nuisance factors suggests $a_{\rm X}/a_{\rm B} = 10^{p} (a_{\rm B}/a_{+1})^q$ with $p = -1.00 \pm 0.03$ (stat) $^{+0.11}_{-0.06}$ (sys) and $q=-1.02 \pm 0.09$ (stat) $^{+0.16}_{-0.00}$ (sys) around a super-critical acceleration $a_{+1}\equiv 1.2\times 10^{-9}~{\rm m}~{\rm s}^{-2}$. In the context of the $Λ$CDM paradigm, this RAR suggests that the NFW DM halo profile is a reasonable description of galactic halos even after the processes of galaxy formation and evolution. In the context of the MOND paradigm, this RAR favors the Simple interpolating function but is inconsistent with the vast majority of other theoretical proposals and fitting functions motivated mainly from sub-critical acceleration data. △ Less

Submitted 12 April, 2019; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: ApJ, accepted (22 pages including 19 figures)

Revisiting the bulge-halo conspiracy I: Dependence on galaxy properties and halo mass

Authors: Francesco Shankar, Alessandro Sonnenfeld, Gary A. Mamon, Kyu-Hyun Chae, Raphael Gavazzi, Tommaso Treu, Benedikt Diemer, Carlo Nipoti, Stewart Buchan, Mariangela Bernardi, Ravi K. Sheth, Marc Huertas-Company

Abstract: We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with m… ▽ More We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and Sérsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a Sérsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations. △ Less

Submitted 17 March, 2017; originally announced March 2017.

Comments: 23 pages, 10 figures, 3 Appendices (with an extra 7 figures). ApJ, accepted. Main results in Figures 3, 5, 6, 8

Sloan Digital Sky Survey IV: Mapping the Milky Way, Nearby Galaxies, and the Distant Universe

Authors: Michael R. Blanton, Matthew A. Bershady, Bela Abolfathi, Franco D. Albareti, Carlos Allende Prieto, Andres Almeida, Javier Alonso-García, Friedrich Anders, Scott F. Anderson, Brett Andrews, Erik Aquino-Ortíz, Alfonso Aragón-Salamanca, Maria Argudo-Fernández, Eric Armengaud, Eric Aubourg, Vladimir Avila-Reese, Carles Badenes, Stephen Bailey, Kathleen A. Barger, Jorge Barrera-Ballesteros, Curtis Bartosz, Dominic Bates, Falk Baumgarten, Julian Bautista, Rachael Beaton , et al. (328 additional authors not shown)

Abstract: We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spat… ▽ More We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016. △ Less

Submitted 29 June, 2017; v1 submitted 28 February, 2017; originally announced March 2017.

Comments: Published in Astronomical Journal

Journal ref: The Astronomical Journal, Volume 154, Number 1, pp. 28-62 (2017)

Comparing PyMorph and SDSS photometry. II. The differences are more than semantics and are not dominated by intracluster light

Authors: M. Bernardi, J. -L. Fischer, R. K. Sheth, A. Meert, M. Huertas-Company, F. Shankar, V. Vikram

Abstract: The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMor… ▽ More The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects. △ Less

Submitted 20 March, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: 15 pages, 17 figures, accepted for publication in MNRAS. The PyMorph luminosities and stellar masses are available at https://www.physics.upenn.edu/~ameert/SDSS_PhotDec/

Comparing PyMorph and SDSS photometry. I. Background sky and model fitting effects

Authors: J. -L. Fischer, M. Bernardi, A. Meert

Abstract: A number of recent estimates of the total luminosities of galaxies in the SDSS are significantly larger than those reported by the SDSS pipeline. This is because of a combination of three effects: one is simply a matter of defining the scale out to which one integrates the fit when defining the total luminosity, and amounts on average to < 0.1 mags even for the most luminous galaxies. The other tw… ▽ More A number of recent estimates of the total luminosities of galaxies in the SDSS are significantly larger than those reported by the SDSS pipeline. This is because of a combination of three effects: one is simply a matter of defining the scale out to which one integrates the fit when defining the total luminosity, and amounts on average to < 0.1 mags even for the most luminous galaxies. The other two are less trivial and tend to be larger; they are due to differences in how the background sky is estimated and what model is fit to the surface brightness profile. We show that PyMorph sky estimates are fainter than those of the SDSS DR7 or DR9 pipelines, but are in excellent agreement with the estimates of Blanton et al. (2011). Using the SDSS sky biases luminosities by more than a few tenths of a magnitude for objects with half-light radii > 7 arcseconds. In the SDSS main galaxy sample these are typically luminous galaxies, so they are not necessarily nearby. This bias becomes worse when allowing the model more freedom to fit the surface brightness profile. When PyMorph sky values are used, then two component Sersic-Exponential fits to E+S0s return more light than single component deVaucouleurs fits (up to ~0.2 mag), but less light than single Sersic fits (0.1 mag). Finally, we show that PyMorph fits of Meert et al. (2015) to DR7 data remain valid for DR9 images. Our findings show that, especially at large luminosities, these PyMorph estimates should be preferred to the SDSS pipeline values. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: 12 pages, 17 figures, accepted for publication in MNRAS

Selection bias in dynamically-measured supermassive black hole samples: Scaling relations and correlations between residuals in semi-analytic galaxy formation models

Authors: Enrico Barausse, Francesco Shankar, Mariangela Bernardi, Yohan Dubois, Ravi K. Sheth

Abstract: Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by usi… ▽ More Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by using a comprehensive galaxy evolution semi-analytic model, which self-consistently evolves supermassive black holes from high-redshift seeds via gas accretion and mergers, and also includes AGN feedback. Once tuned to reproduce the (mean) correlation of black-hole mass with velocity dispersion, the model is unable to also account for the correlation with stellar mass. This behaviour is independent of the model's parameters, thus suggesting an internal inconsistency in the data. The predicted distributions, especially at the low-mass end, are also much broader than observed. However, if selection effects are included, the model's predictions tend to align with the observations. We also demonstrate that the correlations between the residuals of the local scaling relations are more effective than the scaling relations themselves at constraining AGN feedback models. In fact, we find that our semi-analytic model, while in apparent broad agreement with the scaling relations when accounting for selection biases, yields very weak correlations between their residuals at fixed stellar mass, in stark contrast with observations. This problem persists when changing the AGN feedback strength, and is also present in the $z\sim 0$ outputs of the hydrodynamic cosmological simulation Horizon-AGN, which includes state-of-the-art treatments of AGN feedback. This suggests that current AGN feedback models may be too weak or are simply not capturing the effect of the black hole on the stellar velocity dispersion. △ Less

Submitted 18 April, 2017; v1 submitted 6 February, 2017; originally announced February 2017.

Comments: 11 pages, 7 figures. Minor changes to match version to be published in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 468 (4): 4782-4791

Selection bias in dynamically-measured super-massive black hole samples: dynamical masses and dependence on Sérsic index

Authors: Francesco Shankar, Mariangela Bernardi, Ravi K. Sheth

Abstract: We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and Sérsic indices of the bulges of their host galaxies, have dynamical mass and Sérsic index distributions… ▽ More We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and Sérsic indices of the bulges of their host galaxies, have dynamical mass and Sérsic index distributions consistent with those of SDSS early-type galaxies of similar bulge stellar mass. The host galaxies of local black hole samples thus do not appear structurally different from SDSS galaxies, sharing similar dynamical masses, light profiles and light distributions. Analysis of the residuals reveals that velocity dispersion is more fundamental than Sérsic index n in the scaling relations between black holes and galaxies. Indeed, residuals with Sérsic index could be ascribed to the (weak) correlation with bulge mass or even velocity dispersion. Finally, targetted Monte Carlo simulations that include the effects of the sphere of influence of the black hole, and tuned to reproduce the observed residuals and scaling relations in terms of velocity dispersion and stellar mass, show that, at least for galaxies with Mbulge > 1e10 Msun and n>5, the observed mean black hole mass at fixed Sérsic index is biased significantly higher than the intrinsic value. △ Less

Submitted 6 January, 2017; originally announced January 2017.

Comments: 12 pages, 8 Figures. MNRAS, in press

Mass assembly and morphological transformations since $z\sim3$ from CANDELS

Authors: M. Huertas-Company, M. Bernardi, P. G. Pérez-González, M. L. N. Ashby, G. Barro, C. Conselice, E. Daddi, A. Dekel, P. Dimauro, S. M. Faber, N. A. Grogin, J. S. Kartaltepe, D. D. Kocevski, A. M. Koekemoer, D. C. Koo, S. Mei, F. Shankar

Abstract: [abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+dis… ▽ More [abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: resubmitted to MNRAS after addressing minor comments from the referee