Showing 1–50 of 52 results for author: Rodríguez-Puebla, A

The Co-Evolution Between Galaxies and Dark Matter Halos

Authors: Aldo Rodriguez-Puebla

Abstract: The current cosmological paradigm asserts that dark matter halos provide the gravitational scaffolding for galaxy formation through a combination of hierarchical structure formation and non-linear local (g)astrophysical processes. This close relationship, known as the galaxy-halo connection, suggests that the growth and assembly of dark matter halos impact the properties of galaxies. While the ste… ▽ More The current cosmological paradigm asserts that dark matter halos provide the gravitational scaffolding for galaxy formation through a combination of hierarchical structure formation and non-linear local (g)astrophysical processes. This close relationship, known as the galaxy-halo connection, suggests that the growth and assembly of dark matter halos impact the properties of galaxies. While the stellar mass of galaxies correlates strongly with the mass of their halos, it is important to note that the galaxy-halo connection encompasses a broader distribution of galaxy and halo properties. This distribution can be constrained using data from astronomical observations and cosmological $N$-body simulations, a technique known as semi-empirical modeling. By operating at the intersection of observational data and the cosmological structure formation model, the semi-empirical modeling provides valuable insights into galaxy formation and evolution from a cosmological perspective. In this proceeding, we utilize a new sEM-emPIRical modEl, EMPIRE, to explore the star formation, SF, history of central galaxies across cosmic epochs, spanning from dwarfs to massive ellipticals. EMPIRE aims to constrain the multivariate distribution that links galaxy and halo properties. Our findings reveal distinct growth stages for progenitors of central massive galaxies. Evidence suggests that cold streams played a significant role in sustaining SF at higher $z$, while virial shock heating became more prominent at lower $z$. The maximum star formation efficiency, SFE, occurs at a factor of $\sim1.5-2$ below $M_{\rm vir \; shocks}$ for $z\lesssim1$. Furthermore, at higher redshifts, $z>1$, this peak tends towards higher masses, $M_{\rm vir}\sim 2\times 10^{12} M_{\odot}$. Notably, at redshifts higher than $z\sim2$, the peak of SFE aligns comfortably within the region characterized by cold streams. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 16 pages, 7 figures. Extended contribution to the original proceeding for the XVII Latin American Regional IAU Meeting

Non-Monotonic Relations of Galaxy Star Formation, Radius, and Structure at Fixed Stellar Mass

Authors: Jimena Stephenson, Aldo Rodriguez-Puebla, S. M. Faber, Joel R. Primack, Vladimir Avila-Reese, A. R. Calette, Carlo Cannarozzo, James Kakos, Mariana Cano-Díaz, David C. Koo, Francesco Shankar, D. F. Morell

Abstract: We investigate the relation between galaxy structure and star formation rate (SFR) in a sample of $\sim2.9\times10^{4}$ central galaxies with $z<0.0674$ and axial ratios $b/a>0.5$. The star-forming main sequence (SFMS) shows a bend around the stellar mass of $M_\ast\leq{}M_c=2\times10^{10}{}M_{\odot}$. At $M_\ast\leq{}M_c$ the SFMS follows a power-law $\text{SFR}\propto{}M_\ast^{0.85}$, while at h… ▽ More We investigate the relation between galaxy structure and star formation rate (SFR) in a sample of $\sim2.9\times10^{4}$ central galaxies with $z<0.0674$ and axial ratios $b/a>0.5$. The star-forming main sequence (SFMS) shows a bend around the stellar mass of $M_\ast\leq{}M_c=2\times10^{10}{}M_{\odot}$. At $M_\ast\leq{}M_c$ the SFMS follows a power-law $\text{SFR}\propto{}M_\ast^{0.85}$, while at higher masses it flattens. $M_c$ corresponds to a dark matter halo mass of $M_\text{vir}\sim{}10^{11.8}M_{\odot}$ where virial shocks occurs. Some galaxy structure (e.g., half-light radius, $R_e$) exhibits a non-monotonic dependence across the SFMS at a fixed $M_\ast$. We find $\text{SFR}\propto{R_e^{-0.28}}$ at fixed $M_\ast$, consistent with the global Kennicutt-Schmidt (KS) law. This finding suggests that galaxy sizes contribute to the scatter of the SFMS. However, at $M_\ast>M_c$ the relationship between SFR and $R_e$ diminishes. Low-mass galaxies above the mean of the SFMS have smaller radii, exhibit compact and centrally concentrated profiles resembling green valley (GV) and quiescent galaxies at the same mass, and have higher $M_{\text{H}_2}/M_\text{HI}$. Conversely, those below the SFMS exhibit larger radii, lower densities, have no GV or quiescent counterparts at their mass and have lower $M_{\text{H}_2}/M_\text{HI}$. The above data suggest two pathways for quenching low-mass galaxies, $M_\ast\leq{}M_c$: a fast one that changes the morphology on the SFMS and a slow one that does not. Above $M_c$, galaxies below the SFMS resemble GV and quiescent galaxies structurally, implying that they undergo a structural transformation already within the SFMS. For these massive galaxies, CG are strongly bimodal, with SFMS galaxies exhibiting negative color gradients, suggesting most star formation occurs in their outskirts, maintaining them within the SFMS. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 30 pages, 28 figures, submitted to MNRAS

A mock redshift catalogue of the dusty star-forming galaxy population with intrinsic clustering and lensing for deep millimetre surveys

Authors: Norma Araceli Nava-Moreno, Alfredo Montaña, Itziar Aretxaga, Aldo Rodríguez-Puebla, Vladimir Avila-Reese, Edgar Peralta Sánchez

Abstract: We present a new cosmologically motivated mock redshift survey of the Dusty Star-Forming Galaxy population. Our mock survey is based on the Bolshoi-Planck dark-matter halo simulation and covers an area of 5.3 sq. degree. Using a semi-empirical approach, we generate a light cone and populate the dark-matter haloes with galaxies. Infrared properties are assigned to the galaxies based on theoretical… ▽ More We present a new cosmologically motivated mock redshift survey of the Dusty Star-Forming Galaxy population. Our mock survey is based on the Bolshoi-Planck dark-matter halo simulation and covers an area of 5.3 sq. degree. Using a semi-empirical approach, we generate a light cone and populate the dark-matter haloes with galaxies. Infrared properties are assigned to the galaxies based on theoretical and empirical relations from the literature. Additionally, background galaxies are gravitationally lensed by dark-matter haloes along the line-of-sight assuming a point-mass model approximation. We characterize the mock survey by measuring the star formation rate density, integrated number counts, redshift distribution, and infrared luminosity function. When compared with single-dish and interferometric observations, the predictions from our mock survey closely follow the compiled results from the literature. We have also directed this study towards characterizing one of the extragalactic legacy surveys to be observed with the TolTEC camera at the Large Millimeter Telescope: the 0.8 sq. degree Ultra Deep Survey, with expected depths of 0.025, 0.018 and 0.012 mJy beam$^{-1}$ at 1.1, 1.4 and 2.0 mm. Exploiting the clustering information in our mock survey, we investigate its impact on the effect of flux boosting by the fainter population of dusty galaxies, finding that clustering can increase the median boosting by 0.5 per cent at 1.1 mm, 0.8 per cent at 1.4 mm and, 2.0 per cent at 2.0 mm, and with higher dispersion. △ Less

Submitted 31 May, 2024; v1 submitted 22 March, 2024; originally announced March 2024.

Comments: 17 pages, 12 figures

Star-forming and Quiescent Central Galaxies Cluster Similarly: Implications for the Galaxy-Halo Connection

Authors: James Kakos, Aldo Rodriguez-Puebla, Joel R. Primack, Sandra M. Faber, David C. Koo, Peter Behroozi, Vladimir Avila-Reese

Abstract: We measure the clustering of low-redshift SDSS galaxies as a function of stellar mass ($10.0<\log(M_*/M_\odot)<11.5$) and specific star formation rate (sSFR) and compare the results to models of the galaxy--halo connection. We find that the auto-correlation functions of central galaxies exhibit little dependence on sSFR, with the well-known stronger clustering of quiescent galaxies mainly attribut… ▽ More We measure the clustering of low-redshift SDSS galaxies as a function of stellar mass ($10.0<\log(M_*/M_\odot)<11.5$) and specific star formation rate (sSFR) and compare the results to models of the galaxy--halo connection. We find that the auto-correlation functions of central galaxies exhibit little dependence on sSFR, with the well-known stronger clustering of quiescent galaxies mainly attributable to satellites. Because halo assembly history is known to affect distinct halo clustering, this result implies that there is little net correlation between halo assembly history and central galaxy sSFR. However, cross-correlations with satellites are stronger for quiescent centrals than star-forming centrals, consistent with quiescent centrals having more satellites in their haloes at fixed $M_*$, as found in SDSS group catalogues. We model the galaxy--halo connection in an $N$-body simulation by assigning sSFRs to central galaxies in three different ways. Two of the models depend on halo assembly history (being based on halo accretion rate or concentration), while the third is independent of halo assembly history (being based on peak halo circular velocity, $V_\text{peak}$, a proxy for halo mass). All three models replicate the observed auto-correlations of central galaxies, while only the $V_\text{peak}$ model reproduces the observed cross-correlations with satellites. This further suggests that the effects of halo assembly history may not be easily seen in auto-correlations of centrals and implies that a more complete understanding of central galaxy clustering may require more than auto-correlations of centrals alone. Additionally, the good agreement with the $V_\text{peak}$ model supports the idea that quiescent galaxies reside in more massive haloes than star-forming galaxies at fixed $M_*$. △ Less

Submitted 2 March, 2024; originally announced March 2024.

Comments: MNRAS submitted, 23 pages

The evolution of radial gradients of MaNGA quiescent elliptical galaxies: inside-out quenching or outer mass growth?

Authors: V. Avila-Reese, H. Ibarra-Medel, I. Lacerna, A. Rodríguez-Puebla, J. A. Vázquez-Mata, S. F. Sánchez, H. M. Hernández-Toledo, C. Cannarozzo

Abstract: Using spatially-resolved fossil record analysis on a large sample of 'red and dead' elliptical galaxies (classical ellipticals, CLEs) from the MaNGA/SDSS-IV DR15 survey, we reconstruct the archaeological evolution of their radial gradients in mass-to-luminosity ratio ($M/L$), $g-r$ color, and specific star formation (SF) rate. We also calculate other metrics that quantify the inside-out SF quenchi… ▽ More Using spatially-resolved fossil record analysis on a large sample of 'red and dead' elliptical galaxies (classical ellipticals, CLEs) from the MaNGA/SDSS-IV DR15 survey, we reconstruct the archaeological evolution of their radial gradients in mass-to-luminosity ratio ($M/L$), $g-r$ color, and specific star formation (SF) rate. We also calculate other metrics that quantify the inside-out SF quenching and external mass growth processes. The $M/L$ gradients, $\nablaΥ_{\star}$, are approximately flat at high look-back times ($t_{\rm lb}$), but then they become negative and steeper until an epoch, when this trend reverses. These trends are shifted to later epochs the less massive the galaxies are. Color gradients follow qualitatively similar trends. We find that these trends are mainly driven by strong inside-out quenching, without significant outer growth or structural changes overall. Our results suggest a scenario where the main progenitors of local CLE galaxies evolved quasi-passively after an early dissipative phase, but underwent radial photometric changes due to the inside-out quenching that led to the systematic decrease of $\nablaΥ_{\star}$ and to an increase of the light-weighted radius. The late reversing of $\nablaΥ_{\star}$, $t_{\rm lb}\approx2-4$ Gyr, roughly coincides with the global quenching of the CLE galaxies. We have pushed archaeological inferences to the limit, but thanks to the large number of objects and an understanding of how the caveats and assumptions affect our results, we conclude that they offer an average description of evolutionary behaviors of CLE progenitors that is valid at least qualitatively. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 17 pages, 14 figures. Submitted to MNRAS, after Referee Report

The growth of brightest cluster galaxies in the TNG300 simulation:dissecting the contributions from mergers and in situ star formation

Authors: Daniel Montenegro-Taborda, Vicente Rodriguez-Gomez, Annalisa Pillepich, Vladimir Avila-Reese, Laura V. Sales, Aldo Rodríguez-Puebla, Lars Hernquist

Abstract: We investigate the formation of brightest cluster galaxies (BCGs) in the TNG300 cosmological simulation of the IllustrisTNG project. Our cluster sample consists of 700 haloes with $M_{200} \geq 5 \times 10^{13} \, \mathrm{M}_{\odot}$ at $z=0$, along with their progenitors at earlier epochs. This includes 280 systems with $M_{200} \geq 10^{14} \, \mathrm{M}_{\odot}$ at $z=0$, as well as three haloe… ▽ More We investigate the formation of brightest cluster galaxies (BCGs) in the TNG300 cosmological simulation of the IllustrisTNG project. Our cluster sample consists of 700 haloes with $M_{200} \geq 5 \times 10^{13} \, \mathrm{M}_{\odot}$ at $z=0$, along with their progenitors at earlier epochs. This includes 280 systems with $M_{200} \geq 10^{14} \, \mathrm{M}_{\odot}$ at $z=0$, as well as three haloes with $M_{200} \geq 10^{15} \, \mathrm{M}_{\odot}$. We find that the stellar masses and star formation rates of our simulated BCGs are in good agreement with observations at $z \lesssim 0.4$, and that they have experienced, on average, $\sim$2 ($\sim$3) major mergers since $z=1$ ($z=2$). Separating the BCG from the intracluster light (ICL) by means of a fixed 30 kpc aperture, we find that the fraction of stellar mass contributed by ex situ (i.e. accreted) stars at $z=0$ is approximately 70, 80, and 90 per cent for the BCG, BCG+ICL, and ICL, respectively. Tracking our simulated BCGs back in time using the merger trees, we find that they became dominated by ex situ stars at $z \sim $1-2, and that half of the stars that are part of the BCG at $z=0$ formed early ($z \sim 3$) in other galaxies, but `assembled' onto the BCG until later times ($z \approx 0.8$ for the whole sample, $z \approx 0.5$ for BCGs in $M_{200} \geq 5 \times 10^{14} \, \mathrm{M}_{\odot}$ haloes). Finally, we show that the stellar mass profiles of BCGs are often dominated by ex situ stars at all radii, with stars from major mergers being found closer to the centre, while stars that were tidally stripped from other galaxies dominate the outer regions. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 19 pages, 11 figures. Accepted for publication in MNRAS

SDSS-IV MaNGA: The MaNGA Dwarf Galaxy Sample

Authors: M. Cano-Díaz, H. M. Hernández-Toledo, A. Rodríguez-Puebla, H. J. Ibarra-Medel, V. Ávila-Reese, O. Valenzuela, A. E. Medellin-Hurtado, J. A. Vázquez-Mata, A. Weijmans, J. J. González, E. Aquino-Ortiz, L. A. Martínez-Vázquez, Richard R. Lane

Abstract: We present the MaNGA Dwarf galaxy, MaNDala, Value-Added-Catalog, VAC, from the final release of the Sloan Digital Sky Survey-IV program. MaNDala consists of 136 randomly selected bright dwarf galaxies with $M_{*}<10^{9.1}\odot$ and $M_{g}>-18.5$ making it the largest Integral Field Spectroscopy homogeneous sample of dwarf galaxies. We release a photometric analysis of the $g,r$ and $z$ broadband i… ▽ More We present the MaNGA Dwarf galaxy, MaNDala, Value-Added-Catalog, VAC, from the final release of the Sloan Digital Sky Survey-IV program. MaNDala consists of 136 randomly selected bright dwarf galaxies with $M_{*}<10^{9.1}\odot$ and $M_{g}>-18.5$ making it the largest Integral Field Spectroscopy homogeneous sample of dwarf galaxies. We release a photometric analysis of the $g,r$ and $z$ broadband imaging based on the DESI Legacy Imaging Surveys as well as the spectroscopic analysis based on the Pipe3D SDSS-IV VAC. Our release includes the surface brightness (SB), geometric parameters and color profiles, Sérsic fits as well as stellar population properties (such as, stellar ages, metallicities, star formation histories), and emission lines fluxes within the FOV and the effective radii of the galaxies. We find that the majority of the MaNDala galaxies are star forming late-type galaxies with $\langle{}n_{\text{Sersic,r}}\rangle\sim1.6$ that are centrals (central/satellite dichotomy). MaNDala covers a large range of SB values (we find 11 candidates of ultra diffuse galaxies and 3 compact ones), filling the gap between classical dwarfs and low-mass galaxies in the Kormendy Diagram and in the size-mass/luminosity relation, whichseems to flatten at $10^8<M_{*}/\odot<10^{9}$ with $\langle{}R_{e,r}\rangle\sim2.7$ kpc. A large fraction of MaNDala galaxies formed from an early low-metallicity burst of SF but also of late SF events from more metal-enriched gas: half of the MaNDala galaxies assembled $50\%$ of their mass at $\langle{}z\rangle>2$, while the last $20\%$ was at $\langle{}z\rangle<0.3$. Finally, a bending of the sSFR-$M_{*}$ relation at $M_{*}\sim10^{9}\odot$ for the main sequence galaxies seems to be supported by MaNDala. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 29 pages, 17 figures. Accepted to appear in The Astronomical Journal

SDSS-IV MaNGA: pyPipe3D analysis release for 10,000 galaxies

Authors: S. F. Sánchez, J. K. Barrera-Ballesteros, E. Lacerda, A. Mejıa-Narvaez, A. Camps-Fariña, G. Bruzual, C. Espinosa-Ponce, A. Rodrıguez-Puebla, A. R. Calette, H. Ibarra-Medel, V. Avila-Reese, H. Hernandez-Toledo, M. A. Bershady, M. Cano-Diaz, A. M. Munguia-Cordova

Abstract: We present here the analysis performed using the pyPipe3D pipeline for the final MaNGA dataset included in the SDSS seventeenth data-release. This dataset comprises more than 10,000 individual datacubes, being the integral field spectroscopy galaxy survey with the largest number of galaxies. pyPipe3D processes the IFS datacubes to extract spatially-resolved spectroscopic properties of both the ste… ▽ More We present here the analysis performed using the pyPipe3D pipeline for the final MaNGA dataset included in the SDSS seventeenth data-release. This dataset comprises more than 10,000 individual datacubes, being the integral field spectroscopy galaxy survey with the largest number of galaxies. pyPipe3D processes the IFS datacubes to extract spatially-resolved spectroscopic properties of both the stellar population and the ionized-gas emission lines. A brief summary of the properties of the sample and the characteristics of the analyzed data are included. The article provides details on (i) the performed analysis, (ii) a description of the pipeline, (iii) the adopted stellar population library, (iv) the morphological and photometric analysis, (v) the adopted datamodel for the derived spatially resolved properties and (vi) the individual integrated and characteristic galaxy properties included in a final catalog. Comparisons with results from a previous version of the pipeline for earlier data releases and from other tools using this dataset are included. A practical example on how to use of the full dataset, and final catalog illustrates how to handle the delivered product. Our full analysis can be accessed and downloaded from the webpage http://ifs.astroscu.unam.mx/MaNGA/Pipe3D_v3_1_1/. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 66 pages, 29 figures, 18 Tables, submitted to ApJS

SDSS IV MaNGA: Visual Morphological and Statistical Characterization of the DR15 sample

Authors: J. A. Vázquez-Mata, H. M. Hernández-Toledo, V. Avila-Reese, M. Herrera-Endoqui, A. Rodríguez-Puebla, M. Cano-Díaz, I. Lacerna, L. A. Martínez-Vázquez, R. Lane

Abstract: We present a detailed visual morphological classification for the 4614 MaNGA galaxies in SDSS Data Release 15, using image mosaics generated from a combination of r-band (SDSS and deeper DESI Legacy Surveys) images and their digital post-processing. We distinguish 13 Hubble types and identify the presence of bars and bright tidal debris. After correcting the MaNGA sample for volume completeness, w… ▽ More We present a detailed visual morphological classification for the 4614 MaNGA galaxies in SDSS Data Release 15, using image mosaics generated from a combination of r-band (SDSS and deeper DESI Legacy Surveys) images and their digital post-processing. We distinguish 13 Hubble types and identify the presence of bars and bright tidal debris. After correcting the MaNGA sample for volume completeness, we calculate the morphological fractions, the bi-variate distribution of type and stellar mass M*-where we recognise a morphological transition "valley" around S0a-Sa types- and the variations of the g-i colour and luminosity-weighted age over this distribution. We identified bars in 46.8% of galaxies, present in all Hubble types later than S0. This fraction amounts to a factor ~2 larger when compared with other works for samples in common. We detected 14% of galaxies with tidal features, with the fraction changing with M* and morphology. For 355 galaxies, the classification was uncertain; they are visually faint, mostly of low/intermediate masses, low concentrations, and disky in nature. Our morphological classification agrees well with other works for samples in common, though some particular differences emerge, showing that our image procedures allow us to identify a wealth of added value information as compared to SDSS-based previous estimates. Based on our classification, we also propose an alternative criteria for the E-S0 separation, in the structural semi-major to semi-minor axis versus bulge to total light ratio (b/a-B/T) and concentration versus semi-major to semi-minor axis (C-b/a) space. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 23 pages, 20 figures, 4 tables, accepted for publication in MNRAS

Galaxy Correlation Function and Local Density from Photometric Redshifts Using the Stochastic Order Redshift Technique (SORT)

Authors: James Kakos, Joel R. Primack, Aldo Rodriguez-Puebla, Nicolas Tejos, L. Y. Aaron Yung, Rachel S. Somerville

Abstract: The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small ($\sim$10 per cent) reference sample of high-quality redshifts. Within pencil-beam-like sub-volumes surrounding each galaxy, we use the precise dN/d$z$ distribution of the reference sample to recover new redshifts and assign t… ▽ More The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small ($\sim$10 per cent) reference sample of high-quality redshifts. Within pencil-beam-like sub-volumes surrounding each galaxy, we use the precise dN/d$z$ distribution of the reference sample to recover new redshifts and assign them one-to-one to galaxies such that the original rank order of redshifts is preserved. Preserving the rank order is motivated by the fact that random variables drawn from Gaussian probability density functions with different means but equal standard deviations satisfy stochastic ordering. The process is repeated for sub-volumes surrounding each galaxy in the survey. This results in every galaxy with an uncertain redshift being assigned multiple "recovered" redshifts from which a new redshift estimate can be determined. An earlier paper applied SORT to a mock Sloan Digital Sky Survey at $z \lesssim$ 0.2 and accurately recovered the two-point correlation function on scales $\gtrsim$4 $h^{-1}$Mpc. In this paper, we test the performance of SORT in surveys spanning the redshift range 0.75$<z<$2.25. We used two mock surveys extracted from the Small MultiDark-Planck and Bolshoi-Planck N-body simulations with dark matter haloes that were populated by the Santa Cruz semi-analytic model. We find that SORT is able to improve redshift estimates and recover distinctive large-scale features of the cosmic web. Further, it provides unbiased estimates of the redshift-space two-point correlation function $ξ(s)$ on scales $\gtrsim$2.5 $h^{-1}$Mpc, as well as local densities in regions of average or higher density. This may allow improved understanding of how galaxy properties relate to their local environments. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: MNRAS submitted, 24 pages

The differences between mass- and light-derived structural parameters over time for MaNGA Elliptical galaxies

Authors: H. Ibarra-Medel, V. Avila-Reese, I. Lacerna, A. Rodríguez-Puebla, J. A. Vázquez-Mata, H. M. Hernández-Toledo, S. F. Sánchez

Abstract: We apply stellar population synthesis analysis to obtain spatially-resolved archaeological inferences for a large sample of "red and dead" Elliptical galaxies (Classical Ellipticals; CLEs) from the MaNGA/SDSS-IV DR15 survey. From their 2D stellar light and mass maps, we explore the differences between the radial mass and light distributions in the rest-frame bands $g,$ $r,$ and $i$ as functions of… ▽ More We apply stellar population synthesis analysis to obtain spatially-resolved archaeological inferences for a large sample of "red and dead" Elliptical galaxies (Classical Ellipticals; CLEs) from the MaNGA/SDSS-IV DR15 survey. From their 2D stellar light and mass maps, we explore the differences between the radial mass and light distributions in the rest-frame bands $g,$ $r,$ and $i$ as functions of look-back time, $t_{\rm lb}$, or redshift, $z$. We characterize these differences through the ratios between the following mass- and light-derived global properties: sizes, concentrations, and effective surface densities. We find that the mass-to-light ratios of these properties change with $t_{\rm lb}$, more the more massive the galaxies are. The CLE galaxy archaeological progenitors are, on average, less compact, concentrated, and dense in light than in mass as $z$ decreases. However, at later times, when also the evolution of the progenitors becomes passive at all radii, there is an upturn in these trends and the differences between mass and light in compactness/concentration decrease towards $z\sim 0$. The trends in the ratios of mass to light sizes agree qualitatively with results from direct observations in galaxy surveys at different redshifts. We discuss the caveats and interpretations of our results, and speculate that the strong structural evolution found in some previous studies for early-type galaxies could be explained partially by photometric changes rather than by intrinsic structural changes. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS. 19 pages, 11 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)

The galaxy HI-(sub)halo connection and the HI spatial clustering of local galaxies

Authors: A. R. Calette, Aldo Rodríguez-Puebla, Vladimir Avila-Reese, Claudia del P. Lagos

Abstract: We extend the local stellar galaxy-(sub)halo connection to the atomic hydrogen (HI) component by seeding semi-empirically galaxies into a large N-body dark matter (DM) simulation. The main input to construct the mock galaxy catalogue are: our constrained stellar mass-to-(sub)halo circular velocity ($M_{\ast}$-$V_{\rm DM}$) relation, assuming a scatter independent of any galaxy property, and the em… ▽ More We extend the local stellar galaxy-(sub)halo connection to the atomic hydrogen (HI) component by seeding semi-empirically galaxies into a large N-body dark matter (DM) simulation. The main input to construct the mock galaxy catalogue are: our constrained stellar mass-to-(sub)halo circular velocity ($M_{\ast}$-$V_{\rm DM}$) relation, assuming a scatter independent of any galaxy property, and the empirical $M_{\rm HI}$ conditional probability distributions given $M_{\ast}$ for central and satellite galaxies. We find that the $\langle\log M_{\rm HI}\rangle-\log M_{\rm DM}$ relation is not a monotonic increasing function. It increases with mass up to $M_{\rm DM}\sim 10^{12}$ $M_{\odot}$, attaining a maximum of $\langle\log(M_{\rm HI}/M_{\odot})\rangle \sim 9.2$, and at higher (sub)halo masses, $\langle\log(M_{\rm HI})\rangle$ decreases slightly with $M_{\rm DM}$. The scatter around it is also large and mass dependent. The bivariate $M_{\rm HI}$ and $M_{\rm DM}$ distribution is broad and bimodal, specially at $M_{\rm DM}\gtrsim 10^{12}$ $M_\odot$, which is inherited from the input $M_{\rm HI}$ conditional distributions. We also report the total (central+satellites) HI gas mass within halos, $\langle M^{\rm tot}_{\rm HI}(M_{\rm DM})\rangle$, as a function of $M_{\rm DM}$. The mean $M^{\rm tot}_{\rm HI}-M_{\rm DM}$ relation is an increasing monotonic function. The galaxy spatial clustering increases weakly as the $M_{\rm HI}$ threshold increases. Our HI mock galaxies cluster more in comparison to the blind HI ALFALFA (Arecibo Fast Legacy ALFA) survey but we show that it is mainly due to the selection effects. We discuss the implications of our results in the light of predictions from semi-analytical models and hydrodynamics simulations of galaxy evolution. △ Less

Submitted 7 July, 2021; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: In press, 20 pages, 19 figures. Some corrections in text, references added. Results unchanged

The evolution of compact massive quiescent and starforming galaxies derived from the $R_e-R_h$ and $M_{\rm star}-M_h$ relations

Authors: L. Zanisi, F. Shankar, H. Fu, A. Rodriguez-Puebla, V. Avila-Reese, A. Faisst, E. Daddi, L. Boco, A. Lapi, M. Giavalisco, P. Saracco, F. Buitrago, M. Huertas-Company, A. Puglisi, A. Dekel

Abstract: The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a… ▽ More The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a data-driven, flexible theoretical framework to probe the structural evolution of MGs. We assign galaxies to dark matter haloes via stellar mass-halo mass (SMHM) relations with varying high-mass slopes and scatters $σ_{\rm SMHM}$ in stellar mass at fixed halo mass, and assign sizes to galaxies using an empirically-motivated, constant and linear relationship between $R_e$ and the host dark matter halo radius $R_h$. We find that: 1) the fast mean size growth of MGs is well reproduced independently of the shape of the input SMHM relation; 2) the numbers of compact MGs grow steadily until $z\gtrsim2$ and fall off at lower redshifts, suggesting a lesser role of progenitor bias at later epochs; 3) a time-independent scatter $σ_{\rm SMHM}$ is consistent with a scenario in which compact starforming MGs transition into quiescent MGs in a few $10^8$yr with a negligible structural evolution during the compact phase, while a scatter increasing at high redshift implies significant size growth during the starforming phase. A robust measurement of the size function of MGs at high redshift can set strong constraints on the scatter of the SMHM relation and, by extension, on models of galaxy evolution. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. Comments still welcome

The HI and stellar mass bivariate distribution of centrals and satellites for all, late- and early-type local galaxies

Authors: A. R. Calette, Vladimir Avila-Reese, Aldo Rodríguez-Puebla, Claudia del P. Lagos, Barbara Catinella

Abstract: We characterize the conditional distributions of the HI gas-to-stellar mass ratio, $R_{HI}\equiv M_{HI}/M_{\ast}$, given the stellar mass, $M_{\ast}$, of local galaxies from $M_{\ast}\sim 10^7$ to $10^{12}$ $M_{\odot}$ separated into centrals and satellites as well as into late- and early-type galaxies (LTGs and ETGs, respectively). To do so, we use (1) the homogeneous "eXtended GALEX Arecibo SDSS… ▽ More We characterize the conditional distributions of the HI gas-to-stellar mass ratio, $R_{HI}\equiv M_{HI}/M_{\ast}$, given the stellar mass, $M_{\ast}$, of local galaxies from $M_{\ast}\sim 10^7$ to $10^{12}$ $M_{\odot}$ separated into centrals and satellites as well as into late- and early-type galaxies (LTGs and ETGs, respectively). To do so, we use (1) the homogeneous "eXtended GALEX Arecibo SDSS Survey", xGASS (Catinella et al. 2018), by re-estimating their upper limits and taking into account them in our statistical analysis; and (2) the results from a large compilation of HI data reported in Calette et al. (2018). We use the $R_{HI}$ conditional distributions combined with the Galaxy Stellar Mass Function to infer the bivariate $M_{HI}$ and $M_{\ast}$ distribution of all galaxies as well of the late/early-type and central/satellite subsamples and their combinations. Satellites are on average less HI gas-rich than centrals at low and intermediate masses, with differences being larger for ETGs than LTGs; at $M_{\ast}>3-5\times 10^{10}$ $M_{\odot}$ the differences are negligible. The differences in the HI gas content are much larger between LTGs and ETGs than between centrals and satellites. Our empirical HI Mass Function is strongly dominated by central galaxies at all masses. The empirically constrained bivariate $M_{HI}$ and $M_{\ast}$ distributions presented here can be used to compare and constrain theoretical predictions as well as to generate galaxy mock catalogues. △ Less

Submitted 26 May, 2021; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: In press, 21 pages, 15 figures. Minor corrections in text, results unchanged. Supplementary material available at ancillary files

The star-forming main sequence and the contribution of dust-obscured star formation since $z\sim4$ from the FUV+IR luminosity functions

Authors: Aldo Rodriguez-Puebla, Vladimir Avila-Reese, Mariana Cano-Diaz, S. M. Faber, Joel R. Primack, Jose Franco, I. Aretxaga, Eder Santiago-Mayoral

Abstract: An analytical approach is proposed to study the evolution of the star-forming galaxy (SFG) main sequence (MS) and the fraction of dust-obscured SF up to $z\sim4$. Far-ultraviolet (FUV) and infrared (IR) star formation rates, SFRs, are described as conditional probability functions of $M_{\ast}$. We convolve them with the galaxy stellar mass function (GSMF) of SFGs to derive the FUV and IR LFs. The… ▽ More An analytical approach is proposed to study the evolution of the star-forming galaxy (SFG) main sequence (MS) and the fraction of dust-obscured SF up to $z\sim4$. Far-ultraviolet (FUV) and infrared (IR) star formation rates, SFRs, are described as conditional probability functions of $M_{\ast}$. We convolve them with the galaxy stellar mass function (GSMF) of SFGs to derive the FUV and IR LFs. The 2 SF modes formalism is used to describe starburst galaxies. By fitting observed FUV and IR LFs, the parametrization of SFR$_{\rm FUV}-M_{\ast}$ and SFR$_{\rm IR}-M_{\ast}$ are constrained. Our derived SFR$_{\rm FUV+IR}-M_{\ast}$ reproduces the evolution of the MS as compared to other observational inferences. At any redshift, we find that the sSFR$_{\rm FUV+IR}-M_{\ast}$ relation for MS SFGs approaches to a power law at the high-mass end. At lower masses, it bends and eventually the slope sign changes from negative to positive at very low masses. At $z\sim0$, this change of sign is at $M_{\ast}\sim5\times10^{8}{\rm M}_{\odot}$ close to dust-obscured SF regime, $M_{\ast}\sim6\times10^{8}{\rm M}_{\odot}$. The slope sign change is related to the knee of the FUV LF. Our derived dust-obscured fractions agree with previous determinations at $0\leq z\leq2.5$. Dust-obscured fractions depend strongly on mass with almost no dependence with redshift at $z\gtrsim1.2$. At $z\lesssim0.75$ high-mass galaxies become more "transparent" compared to their high redshift counterparts. On the opposite, low- and intermediate-mass galaxies have become more obscured by dust. The joint evolution of the GSMF and the FUV and IR LFs is a promising approach to study mass growth and dust formation/destruction mechanisms. △ Less

Submitted 21 October, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Submitted to ApJ, 20 pages, 12 figures. Comments are welcome. This is an improved version after the first Referee report. Results remained unchanged

Clustering and Halo Abundances in Early Dark Energy Cosmological Models

Authors: Anatoly Klypin, Vivian Poulin, Francisco Prada, Joel Primack, Marc Kamionkowski, Vladimir Avila-Reese, Aldo Rodriguez-Puebla, Peter Behroozi, Doug Hellinger, Tristan L. Smith

Abstract: LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% o… ▽ More LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: 13 pages, 12 figures, submitted to MNRAS

A Universal fundamental plane and the $M_{dyn}-M_{\star}$ relation for galaxies with CALIFA and MaNGA

Authors: E. Aquino-Ortíz, S. F. Sánchez, O. Valenzuela, H. Hernández-Toledo, Yunpeng Jin, Ling Zhu, Glenn van de Ven, J. K. Barrera-Ballesteros, V. Avila-Reese, A. Rodríguez-Puebla, Patricia B. Tissera

Abstract: We use the stellar kinematics for $2458$ galaxies from the MaNGA survey to explore dynamical scaling relations between the stellar mass $M_{\star}$ and the total velocity parameter at the effective radius, $R_e$, defined as $S_{K}^{2}=KV_{R_e}^{2}+σ_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $σ_{\star_e}$. We confirm that spheroidal and spiral galaxies foll… ▽ More We use the stellar kinematics for $2458$ galaxies from the MaNGA survey to explore dynamical scaling relations between the stellar mass $M_{\star}$ and the total velocity parameter at the effective radius, $R_e$, defined as $S_{K}^{2}=KV_{R_e}^{2}+σ_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $σ_{\star_e}$. We confirm that spheroidal and spiral galaxies follow the same $M_{\star}-S_{0.5}$ scaling relation with lower scatter than the $M_{\star}-V_{R_e}$ and $M_{\star}-σ_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(Υ_{e}) = log (S_{0.5}^{2}) - log (I_{e}) - log (R_{e}) + C$, which in addition to kinematics, $S_{0.5}$, and effective radius, $R_e$, it includes surface brightness, $I_e$, and dynamical mass-to-light ratio, $Υ_e$. We use sophisticated Schwarzschild dynamical models for a sub-sample of 300 galaxies from the CALIFA survey to calibrate the so called Universal Fundamental Plane. That calibration allows us to propose both: (i) a parametrization to estimate the difficult-to-measure dynamical mass-to-light ratio at the effective radius; and (ii) a new dynamical mass proxy consistent with dynamical models within $0.09\ dex$. We reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. We use the estimated dynamical mass-to-light ratio from our analysis, $Υ_{e}^{fit}$, to explore the Universal Fundamental Plane with the MaNGA data set. We find that all classes of galaxies, from spheroids to disks, follow this Universal Fundamental Plane with a scatter significantly smaller $(0.05\ dex)$ than the one reported for the $M_{\star}-S_{0.5}$ relation $(0.1\ dex)$, the Fundamental Plane $(\sim 0.09\ dex)$ and comparable with Tully-Fisher studies $(\sim 0.05\ dex)$, but for a wider range of galaxy types. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Comments: Submitted to ApJ. 21 pages, 12 figures

The bivariate gas-stellar mass distributions and the mass functions of early- and late-type galaxies at $z\sim0$

Authors: Aldo Rodriguez-Puebla, A. R. Calette, Vladimir Avila-Reese, Vicente Rodriguez-Gomez, Marc Huertas-Company

Abstract: We report the bivariate HI- and H$_2$-stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for HI, H$_2$, cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the HI and H$_2$ conditional distributions and the galaxy stellar mass function, GSMF. For the conditional distributions we use the compil… ▽ More We report the bivariate HI- and H$_2$-stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for HI, H$_2$, cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the HI and H$_2$ conditional distributions and the galaxy stellar mass function, GSMF. For the conditional distributions we use the compilation presented in Calette et al. 2018. For determining the GSMF from $M_{\ast}\sim3\times10^{7}$ to $3\times10^{12}$ $M_{\odot}$, we combine two spectroscopic samples from the SDSS at the redshift range $0.0033<z<0.2$. We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is $α\approx-1.4$, consistent with previous determinations. The obtained HI MFs agree with radio blind surveys. Similarly, the H$_2$ MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while $\sim96$% of the HI and H$_2$ mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of H$_2$ and total neutral H in late-type galaxies of $\sim 1.3$ and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting H$_2$ into stars and in transforming HI gas into H$_2$. Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas--stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models. △ Less

Submitted 20 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: 37 pages, 17 figures. Accepted for publication in PASA. A code that displays tables and figures with all the relevant statistical distributions and correlations discussed in this paper is available here https://github.com/arcalette/Python-code-to-generate-Rodriguez-Puebla-2020-results

The SFR-radius connection: data and implications for wind strength and halo concentration

Authors: Lin Lin, S. M. Faber, David C. Koo, Samir Salim, Aaron A. Dutton, Jerome J. Fang, Fangzhou Jiang, Cristoph T. Lee, Aldo Rodríguez-Puebla, A. van der Wel, Yicheng Guo, Guillermo Barro, Joel R. Primack, Avishai Dekel, Zhu Chen, Yifei Luo, Viraj Pandya, Rachel S. Somerville, Henry C. Ferguson, Susan Kassin, Anton M. Koekemoer, Norman A. Grogin, Audrey Galametz, P. Santini, Hooshang Nayyeri , et al. (4 additional authors not shown)

Abstract: This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend… ▽ More This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR. △ Less

Submitted 20 July, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 21 pages, 9 figures, ApJ accepted

Quenching as a Contest between Galaxy Halos and their Central Black Holes

Authors: Zhu Chen, S. M. Faber, David. C. Koo, Rachel S. Somerville, Joel R. Primack, Avishai Dekel, Aldo Rodríguez-Puebla, Yicheng Guo, Guillermo Barro, Dale D. Kocevski, A. van der Wel, Joanna Woo, Eric F. Bell, Jerome J. Fang, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company, Fangzhou Jiang, Susan Kassin, Lin Lin, Fengshan Liu, Yifei Luo, Zhijian Luo, Camilla Pacifici, Viraj Pandya , et al. (5 additional authors not shown)

Abstract: Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-formi… ▽ More Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws. △ Less

Submitted 26 May, 2020; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

Structural and Stellar Population Properties vs. Bulge Types in Sloan Digital Sky Survey Central Galaxies

Authors: Yifei Luo, S. M. Faber, Aldo Rodriguez-Puebla, Joanna Woo, Yicheng Guo, David C. Koo, Joel R. Primack, Zhu Chen, Hassen M. Yesuf, Lin Lin, Guillermo Barro, Jerome J. Fang, Viraj Pandya, M. Huertas-Company, Shude Mao

Abstract: This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with… ▽ More This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with $Δ\langleμ_e\rangle$. $ΔΣ_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $ΔΣ_1$ and $Δ\langleμ_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench. △ Less

Submitted 9 March, 2020; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: 23 pages, 17 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 2, 2020, p.1686-1707

SDSS-IV MaNGA: Effects of Morphology in the global and local Star Formation Main Sequences

Authors: M. Cano-Díaz, V. Ávila-Reese, S. F. Sánchez, H. M. Hernández-Toledo, A. Rodríguez-Puebla, M. Boquien, H. Ibarra-Medel

Abstract: We study the global star-formation rate (SFR) vs. stellar mass (M$_*$) correlation, and the spatially-resolved SFR surface density ($Σ_{SFR}$) vs. stellar mass surface density (\Sm) correlation, in a sample of $\sim2,000$ galaxies from the MaNGA MPL-5 survey. We classify galaxies and spatially-resolved areas into star-forming and retired according to their ionization processes. We confirm the exis… ▽ More We study the global star-formation rate (SFR) vs. stellar mass (M$_*$) correlation, and the spatially-resolved SFR surface density ($Σ_{SFR}$) vs. stellar mass surface density (\Sm) correlation, in a sample of $\sim2,000$ galaxies from the MaNGA MPL-5 survey. We classify galaxies and spatially-resolved areas into star-forming and retired according to their ionization processes. We confirm the existence of a Star-Forming Main Sequence (SFMS) for galaxies and spatially-resolved areas, and show that they have the same nature, with the global as a consequence of the local one. The latter presents a bend below a limit \Sm value, $\approx 3\times 10^7$ M$_\odot$kpc$^{-2}$, which is not physical. Using only star-forming areas (SFAs) above this limit, a slope and a scatter of $\approx1$ and $\approx0.27$ dex are determined. The retired galaxies/areas strongly segregate from their respective SFMS's, by $\sim -1.5$ dex on average. We explore how the global/local SFMS's depend on galaxy morphology, finding that for star-forming galaxies and SFAs, there is a trend to lower values of star-formation activity with earlier morphological types, which is more pronounced for the local SFMS. The morphology not only affects the global SFR due to the diminish of SFAs with earlier types, but also affects the local SF process. Our results suggest that the local SF at all radii is established by some universal mechanism partially modulated by morphology. Morphology seems to be connected to the slow aging and sharp decline of the SF process, and on its own it may depend on other properties as the environment. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 24 pages, 16 figues. Accepted to appear in MNRAS

Astro2020: Empirically Constraining Galaxy Evolution

Authors: Peter Behroozi, Matthew R. Becker, Frank C. van den Bosch, Jarle Brinchmann, Charlie Conroy, Mark Dickinson, Christopher M. Hirata, Andrew Hearin, Alexie Leauthaud, Chun Ly, Yao-Yuan Mao, Benjamin P. Moster, Christine O'Donnell, Casey Papovich, Aldo Rodríguez-Puebla, Rachel Somerville, Erik Tollerud, Jeremy Tinker, Yun Wang, Risa H. Wechsler, Charity Woodrum, Ann Zabludoff, Dennis Zaritsky, Andrew R. Zentner, Huanian Zhang

Abstract: Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive bl… ▽ More Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

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

The Galaxy Halo Connection in Modified Gravity Cosmologies: Environment Dependence of Galaxy Luminosity function

Authors: N. Chandrachani Devi, Aldo Rodríguez-Puebla, O. Valenzuela, Vladimir Avila-Reese, César Hernández-Aguayo, Baojiu Li

Abstract: We investigate the dependence of the galaxy-halo connection and galaxy density field in modified gravity models using the $N-$body simulations for $f(R)$ and nDGP models at $z=0$. Because of the screening mechanisms employed by these models, chameleon and Vainshtein, halos are clustered differently in the non-linear regime of structure formation. We quantify their deviations in the galaxy density… ▽ More We investigate the dependence of the galaxy-halo connection and galaxy density field in modified gravity models using the $N-$body simulations for $f(R)$ and nDGP models at $z=0$. Because of the screening mechanisms employed by these models, chameleon and Vainshtein, halos are clustered differently in the non-linear regime of structure formation. We quantify their deviations in the galaxy density field from the standard $Λ$CDM model under different environments. We populate galaxies in halos via the (Sub)Halo Abundance Matching. Our main results are: 1) The galaxy-halo connection {\it strongly} depends on the gravity model; a maximum variation of $\sim40\%$ is observed between Halo Occupational Distribution (HOD) parameters; 2) $f(R)$ gravity models predict an excess of galaxies in low density environments of $\sim10\%$ but predict a deficit of $\sim10\%$ at high density environments for $|f_{R0}| = 10^{-4}$ and $10^{-6}$ while $|f_{R0}| = 10^{-5}$ predicts more high density structures; nDGP models are consistent with $Λ$CDM; 3) Different gravity models predict different dependences of the galaxy luminosity function (GLF) with the environment, especially in void-like regions we find differences around $\sim10\%$ for the $f(R)$ models while nDPG models remain closer to $Λ$CDM for low-luminosity galaxies but there is a deficit of $\sim11\%$ for high-luminosity galaxies in all environments. We conclude that the dependence of the GLF with environment might provide a test to distinguish between gravity models and their screening mechanisms from the $Λ$CDM. We provide HOD parameters for the gravity models analyzed in this paper. △ Less

Submitted 25 June, 2019; v1 submitted 7 January, 2019; originally announced January 2019.

Comments: 22 pages, 14 figures, 4 tables

Optical Integral Field Spectroscopy observations applied to simulated galaxies: Testing the fossil record method

Authors: Héctor J. Ibarra-Medel, Vladimir Avila-Reese, Sebastián F. Sánchez, Alejandro González-Samaniego, Aldo Rodríguez-Puebla

Abstract: By means of post-processed cosmological Hydrodynamics simulations we explore the ability of the fossil record method to recover the stellar mass, age gradients, and global/radial star formation and mass growth histories of galaxies observed with an optical Integral Field Unity. We use two simulations of representative Milky Way-sized galaxies, one disk dominated and another bulge-dominated. We gen… ▽ More By means of post-processed cosmological Hydrodynamics simulations we explore the ability of the fossil record method to recover the stellar mass, age gradients, and global/radial star formation and mass growth histories of galaxies observed with an optical Integral Field Unity. We use two simulations of representative Milky Way-sized galaxies, one disk dominated and another bulge-dominated. We generate data cubes emulating MaNGA/SDSS-IV instrumental and observational conditions, and analyze them using the Pipe3D pipeline. For optimal MaNGA instrumental/observational setups, the global masses, ages, and dust extinctions of both galaxies are well recovered within the uncertainties, while for the radial distributions, the trend is to get slightly flatter age gradients and weaker inside-out growth modes than the true ones. The main bias is in recovering younger stellar populations in the inner, older regions, and slightly older stellar populations in the outer, younger regions, as well as lower global stellar masses. These trends are enhanced as lower is the spatial sampling/resolution (less number of fibers), more inclined and dust-attenuated is the observed galaxy, lower is the signal-to-noise ratio, and worse is the seeing. The recovered properties and histories are more affected for the disk-dominated galaxy with a prominent inside-out growth mode than for the bulge-dominated galaxy assembled earlier and more coherently. We make publicly available the code and the mock IFU datacubes used in this paper. △ Less

Submitted 27 November, 2018; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: Accepted for publication to MNRAS

SDSS IV MaNGA - An archaeological view of the Cosmic Star Formation History

Authors: S. F. Sánchez, V. Avila-Reese, A. Rodríguez-Puebla, H. Ibarra-Medel, R. Calette, M. Bershady, H. Hernández-Toledo, K. Pan, D. Bizyaev

Abstract: We present the results of the archaeological analysis of the stellar populations of a sample of ~4,000 galaxies observed by the SDSS-IV-MaNGA survey using Pipe3D. Based on this analysis we extract a sample of ~150,000 SFRs and stellar masses that mimic a single cosmological survey covering the redshift range between z~0 to z~7. We confirm that the Star-Forming Main Sequence holds as a tight relati… ▽ More We present the results of the archaeological analysis of the stellar populations of a sample of ~4,000 galaxies observed by the SDSS-IV-MaNGA survey using Pipe3D. Based on this analysis we extract a sample of ~150,000 SFRs and stellar masses that mimic a single cosmological survey covering the redshift range between z~0 to z~7. We confirm that the Star-Forming Main Sequence holds as a tight relation in this range of redshifts, and evolves strongly in both the zero-point and slope. This evolution is different for the population of local star-forming (SFGs) and retired (RGs) galaxies, with the latter presenting a stronger evolution in the zero-point and a weaker evolution in the slope. The fraction of RGs decreases rapidly with z, particularly for those classified as RGs at z~0. Contrary to previous studies we detect RGs well above z>1, although not all of them are progenitors of local RGs. Finally, adopting the required corrections to make the survey complete in mass in a limited volume, we recover the cosmic star-formation rate (SFR), stellar mass density, and average specific SFR histories of the Universe in this wide range of look-back times, with a remarkable agreement with the values reported by various cosmological surveys. We demonstrate that the progenitors of local RGs were more actively forming stars in the past, contributing to most of the cosmic SFR density at z>0.5, and to most of the cosmic stellar mass density at any redshift. They suffer a general quenching in the SFR at z~0.35. Below this redshift the progenitors of local SFGs dominate the SFR density of the Universe. △ Less

Submitted 2 November, 2018; v1 submitted 30 July, 2018; originally announced July 2018.

Comments: 33 pages, 23 figures, accepted to publish to MNRAS

Kinematic scaling relations of CALIFA galaxies: A dynamical mass proxy for galaxies across the Hubble sequence

Authors: E. Aquino-Ortíz, O. Valenzuela, S. F. Sánchez, H. Hernández-Toledo, V. Ávila-Reese, G. van de Ven, A. Rodríguez-Puebla, L. Zhu, B. Mancillas, M. Cano-Díaz, R. García-Benito

Abstract: We used ionized gas and stellar kinematics for 667 spatially resolved galaxies publicly available from the Calar Alto Legacy Integral Field Area survey (CALIFA) 3rd Data Release with the aim of studying kinematic scaling relations as the Tully $\&$ Fisher (TF) relation using rotation velocity, $V_{rot}$, the Faber $\&$ Jackson (FJ) relation using velocity dispersion, $σ$, and also a combination of… ▽ More We used ionized gas and stellar kinematics for 667 spatially resolved galaxies publicly available from the Calar Alto Legacy Integral Field Area survey (CALIFA) 3rd Data Release with the aim of studying kinematic scaling relations as the Tully $\&$ Fisher (TF) relation using rotation velocity, $V_{rot}$, the Faber $\&$ Jackson (FJ) relation using velocity dispersion, $σ$, and also a combination of $V_{rot}$ and $σ$ through the $S_{K}$ parameter defined as $S_{K}^2 = KV_{rot}^2 + σ^2$ with constant $K$. Late-type and early-type galaxies reproduce the TF and FJ relations. Some early-type galaxies also follow the TF relation and some late-type galaxies the FJ relation, but always with larger scatter. On the contrary, when we use the $S_{K}$ parameter, all galaxies, regardless of the morphological type, lie on the same scaling relation, showing a tight correlation with the total stellar mass, $M_\star$. Indeed, we find that the scatter in this relation is smaller or equal to that of the TF and FJ relations. We explore different values of the $K$ parameter without significant differences (slope and scatter) in our final results with respect the case $K=0.5$ besides than a small change in the zero point. We calibrate the kinematic $S_{K}^2$ dynamical mass proxy in order to make it consistent with sophisticated published dynamical models within $0.15\ dex$. We show that the $S_{K}$ proxy is able to reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. Our result may be useful in order to produce fast estimations of the central dynamical mass in galaxies and to study correlations in large galaxy surveys. △ Less

Submitted 23 June, 2018; originally announced June 2018.

Comments: accepted by MNRAS

Dark Matter Halo Properties vs. Local Density and Cosmic Web Location

Authors: Tze Goh, Joel Primack, Christoph T. Lee, Miguel Aragon-Calvo, Doug Hellinger, Peter Behroozi, Aldo Rodriguez-Puebla, Elliot Eckholm, Kathryn Johnston

Abstract: We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo prope… ▽ More We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo properties that we study are the specific mass accretion rate, spin parameter, concentration, prolateness, scale factor of the last major merger, and scale factor when the halo had half of its z = 0 mass. For all these properties, we find that there is no discernible difference between the halo properties in filaments, walls, or voids when compared at the same environmental density. As a result, we conclude that environmental density is the core attribute that affects these properties. This conclusion is in line with recent findings that properties of galaxies in redshift surveys are independent of their cosmic web environment at the same environmental density at z ~ 0. We also find that the local web environment of the Milky Way and the Andromeda galaxies near the centre of a cosmic wall does not appear to have any effect on the properties of these galaxies' dark matter halos except for their orientation, although we find that it is rather rare to have such massive halos near the centre of a relatively small cosmic wall. △ Less

Submitted 22 October, 2019; v1 submitted 13 May, 2018; originally announced May 2018.

Comments: 23 pages

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

The $HI$- and $H_{2}$-to-stellar mass correlations of late- and early-type galaxies and their consistency with the observational mass functions

Authors: A. R. Calette, Vladimir Avila-Reese, Aldo Rodríguez-Puebla, Héctor Hernández-Toledo, Emmanouil Papastergis

Abstract: We compile and carrefully homogenize local galaxy samples with available information on stellar, $\rm HI$ and/or $\rm H_{2}$ masses, and morphology. After processing the information on upper limits in the case of non gas detections, we determine the $\rm HI$- and $\rm H_{2}$-to-stellar mass relations and their $1σ$ scatter for both late- and early-type galaxies. The obtained relations are fitted t… ▽ More We compile and carrefully homogenize local galaxy samples with available information on stellar, $\rm HI$ and/or $\rm H_{2}$ masses, and morphology. After processing the information on upper limits in the case of non gas detections, we determine the $\rm HI$- and $\rm H_{2}$-to-stellar mass relations and their $1σ$ scatter for both late- and early-type galaxies. The obtained relations are fitted to single or double power laws. Late-type galaxies are significantly gas richer than early-type ones, specially at high masses. The respective $\rm H_{2}$-to-$\rm HI$ mass ratios as a function of $M_{\ast}$ are discussed. Further, we constrain the full mass-dependent distribution functions of the $\rm HI$- and $\rm H_{2}$-to-stellar mass ratios. We find that they can be described by a Schechter function for late types and a (broken) Schechter + uniform function for early types. By using the observed galaxy stellar mass function and the volume-complete late-to-early-type galaxy ratio as a function of $M_{\ast}$, these empirical distribution functions are mapped into $\rm HI$ and $\rm H_{2}$ mass functions. The obtained mass functions are consistent with those inferred from large surveys. The empirical gas-to-stellar mass relations and their distributions for local late- and early-type galaxies presented here can be used to constrain models and simulations of galaxy evolution. △ Less

Submitted 29 June, 2018; v1 submitted 20 March, 2018; originally announced March 2018.

Comments: 43 pages, 18 figures, to appear in RMxAA. Minor corrections introduced. The presented results are optimal for comparisons with theoretical predictions. Py-code to generate the HI- and H2-to-stellar mass relations and their 1sigma scatter, as well as the full mass-dependent distribution functions of the MHI/Ms and MH2/Ms ratios are available upon request to A.R. Calette

The global and radial stellar mass assembly of Milky Way-sized galaxies

Authors: V. Avila-Reese, A. González-Samaniego, P. Colín, H. Ibarra-Medel, A. Rodríguez-Puebla

Abstract: We study the global and radial stellar mass assembly of eight zoomed-in MW-sized galaxies produced in Hydrodynamics cosmological simulations. The disk-dominated galaxies (4) show a fast initial stellar mass growth in the innermost parts, driven mostly by in-situ SF, but since $z\sim2-1$ the SF enters in a long-term quenching phase. The outer regions follow this trend but more gentle as more extern… ▽ More We study the global and radial stellar mass assembly of eight zoomed-in MW-sized galaxies produced in Hydrodynamics cosmological simulations. The disk-dominated galaxies (4) show a fast initial stellar mass growth in the innermost parts, driven mostly by in-situ SF, but since $z\sim2-1$ the SF enters in a long-term quenching phase. The outer regions follow this trend but more gentle as more external they are. As the result, the radial stellar mass growth is highly inside-out due to both the inside-out structural growth and inside-out SF quenching. The half-mass radius evolves fast; for instance, $R_{0.5}$($z=1$)$<0.5$$R_{0.5}$($z=0$). Two other runs resemble lenticular galaxies. One shows also a pronounced inside-out growth and the other one presents a nearly uniform radial mass assembly. The other two galaxies suffered late major mergers. Their normalized radial mass growth histories (MGHs) are nearly close among them but with periods of outside-in assembly during or after the mergers. For all the simulations, the archaeological radial MGHs calculated from the $z = 0$ stellar-particles age distribution are similar to the current MGHs, which evidences that the mass assembly by ex-situ stars and the radial mass transport do not change significantly their radial mass distributions. Our results agree qualitatively with observational inferences from the fossil record method applied to a survey of local galaxies and from look-back observations of progenitors of MW-sized galaxies. However, the inside-out growth mode is more pronounced and the $R_{0.5}$ growth is faster in simulations than in observational inferences. △ Less

Submitted 31 January, 2018; v1 submitted 5 January, 2018; originally announced January 2018.

Comments: 23 pages, 15 figures. Accepted for publication in ApJ. Minor corrections and more references added. The discussion was extended in section 6.1

Tidal Stripping and Post-Merger Relaxation of Dark Matter Halos: Causes and Consequences of Mass Loss

Authors: Christoph T. Lee, Joel R. Primack, Peter Behroozi, Aldo Rodríguez-Puebla, Doug Hellinger, Avishai Dekel

Abstract: We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initial… ▽ More We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initially boost $M_{\mathrm{vir}}$ and typically cause the final halo to become more prolate and less relaxed and to have higher spin and lower NFW concentration. As the halo relaxes, high energy material from the recent merger gradually escapes beyond the virial radius, temporarily resulting in a net negative accretion rate that reduces the halo mass by $5-15\%$ on average. Halos that experience a major merger around $z = 0.4$ typically reach a minimum mass near $z = 0$. Tidal stripping mainly occurs in dense regions, and it causes halos to become less prolate and have lower spins and higher NFW concentrations. Tidally stripped halos often lose a large fraction of their peak mass ($> 20\%$) and most never recover (or even reattain a positive accretion rate). Low mass halos can be strongly affected by both post-merger mass loss and tidal stripping, while high mass halos are predominantly influenced by post-merger mass loss and show few signs of significant tidal stripping. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Does the Galaxy-Halo Connection Vary with Environment?

Authors: Radu Dragomir, Aldo Rodriguez-Puebla, Joel R. Primack, Christoph T. Lee

Abstract: SubHalo Abundance Matching (SHAM) assumes that one (sub)halo property, such as mass Mvir or peak circular velocity Vpeak, determines properties of the galaxy hosted in each (sub)halo such as its luminosity or stellar mass. This assumption implies that the dependence of Galaxy Luminosity Functions (GLFs) and the Galaxy Stellar Mass Function (GSMF) on environmental density is determined by the corre… ▽ More SubHalo Abundance Matching (SHAM) assumes that one (sub)halo property, such as mass Mvir or peak circular velocity Vpeak, determines properties of the galaxy hosted in each (sub)halo such as its luminosity or stellar mass. This assumption implies that the dependence of Galaxy Luminosity Functions (GLFs) and the Galaxy Stellar Mass Function (GSMF) on environmental density is determined by the corresponding halo density dependence. In this paper, we test this by determining from an SDSS sample the observed dependence with environmental density of the ugriz GLFs and GSMF for all galaxies, and for central and satellite galaxies separately. We then show that the SHAM predictions are in remarkable agreement with these observations, even when the galaxy population is divided between central and satellite galaxies. However, we show that SHAM fails to reproduce the correct dependence between environmental density and g-r color for all galaxies and central galaxies, although it better reproduces the color dependence on environmental density of satellite galaxies. △ Less

Submitted 8 February, 2018; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: 21 pages, 11 figures. Accepted for publication in MNRAS

Demographics of Star-forming Galaxies since $z\sim2.5$. I. The $UVJ$ Diagram in CANDELS

Authors: Jerome J. Fang, S. M. Faber, David C. Koo, Aldo Rodriguez-Puebla, Yicheng Guo, Guillermo Barro, Peter Behroozi, Gabriel Brammer, Zhu Chen, Avishai Dekel, Henry C. Ferguson, Eric Gawiser, Mauro Giavalisco, Jeyhan Kartaltepe, Dale D. Kocevski, Anton M. Koekemoer, Elizabeth J. McGrath, Daniel McIntosh, Jeffrey A. Newman, Camilla Pacifici, Viraj Pandya, Pablo G. Perez-Gonzalez, Joel R. Primack, Brett Salmon, Jonathan R. Trump , et al. (19 additional authors not shown)

Abstract: This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correla… ▽ More This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $μ\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift. △ Less

Submitted 6 April, 2018; v1 submitted 15 October, 2017; originally announced October 2017.

Comments: 36 pages, 26 figures, ApJ accepted

SSDSS IV MaNGA - Properties of AGN host galaxies

Authors: S. F. Sanchez, V. Avila-Reese, H. Hernandez-Toledo, E. Cortes-Suarez, A. Rodriguez-Puebla, H. Ibarra-Medel, M. Cano-Diaz, J. K. Barrera-Ballesteros, C. A. Negrete, A. R. Calette, A. de Lorenzo-Caceres, R. A. Ortega-Minakata, E. Aquino, O. Valenzuela, J. C. Clemente, T. Storchi-Bergmann, R. Riffel, J. Schimoia, R. A. Riffel, S. B. Rembold, J. R. Brownstein, K. Pan, R. Yates, N. Mallmann, T. Bitsakis

Abstract: We present here the characterization of the main properties of a sample of 98 AGN host galaxies, both type-II and type-I, in comparison with those of about 2700 non-active galaxies observed by the MaNGA survey. We found that AGN hosts are morphologically early-type or early-spirals. For a given morphology AGN hosts are, in average, more massive, more compact, more central peaked and rather pressur… ▽ More We present here the characterization of the main properties of a sample of 98 AGN host galaxies, both type-II and type-I, in comparison with those of about 2700 non-active galaxies observed by the MaNGA survey. We found that AGN hosts are morphologically early-type or early-spirals. For a given morphology AGN hosts are, in average, more massive, more compact, more central peaked and rather pressurethan rotational-supported systems. We confirm previous results indicating that AGN hosts are located in the intermediate/transition region between star-forming and non-star-forming galaxies (i.e., the so-called green valley), both in the ColorMagnitude and the star formation main sequence diagrams. Taking into account their relative distribution in terms of the stellar metallicity and oxygen gas abundance and a rough estimation of their molecular gas content, we consider that these galaxies are in the process of halting/quenching the star formation, in an actual transition between both groups. The analysis of the radial distributions of the starformation rate, specific star-formation rate, and molecular gas density shows that the quenching happens from inside-out involving both a decrease of the efficiency of the star formation and a deficit of molecular gas. All the intermediate data-products used to derive the results of our analysis are distributed in a database including the spatial distribution and average properties of the stellar populations and ionized gas, published as a Sloan Digital Sky Survey Value Added Catalog being part of the 14th Data Release: http://www.sdss.org/dr14/manga/manga-data/manga-pipe3d-value-added-catalog/ △ Less

Submitted 19 February, 2018; v1 submitted 15 September, 2017; originally announced September 2017.

Comments: 48 pages, 14 figures, in press in RMxAA

Stochastic Order Redshift Technique (SORT): a simple, efficient and robust method to improve cosmological redshift measurements

Authors: Nicolas Tejos, Aldo Rodriguez-Puebla, Joel R. Primack

Abstract: We present a simple, efficient and robust approach to improve cosmological redshift measurements. The method is based on the presence of a reference sample for which a precise redshift number distribution (dN/dz) can be obtained for different pencil-beam-like sub-volumes within the original survey. For each sub-volume we then impose: (i) that the redshift number distribution of the uncertain redsh… ▽ More We present a simple, efficient and robust approach to improve cosmological redshift measurements. The method is based on the presence of a reference sample for which a precise redshift number distribution (dN/dz) can be obtained for different pencil-beam-like sub-volumes within the original survey. For each sub-volume we then impose: (i) that the redshift number distribution of the uncertain redshift measurements matches the reference dN/dz corrected by their selection functions; and (ii) the rank order in redshift of the original ensemble of uncertain measurements is preserved. The latter step is motivated by the fact that random variables drawn from Gaussian probability density functions (PDFs) of different means and arbitrarily large standard deviations satisfy stochastic ordering. We then repeat this simple algorithm for multiple arbitrary pencil-beam-like overlapping sub-volumes; in this manner, each uncertain measurement has multiple (non-independent) "recovered" redshifts which can be used to estimate a new redshift PDF. We refer to this method as the Stochastic Order Redshift Technique (SORT). We have used a state-of-the-art N-body simulation to test the performance of SORT under simple assumptions and found that it can improve the quality of cosmological redshifts in an efficient and robust manner. Particularly, SORT redshifts are able to recover the distinctive features of the 'cosmic web' and can provide unbiased measurement of the two-point correlation function on scales > 4 Mpc/h. Given its simplicity, we envision that a method like SORT can be incorporated into more sophisticated algorithms aimed to exploit the full potential of large extragalactic photometric surveys. △ Less

Submitted 11 September, 2017; originally announced September 2017.

Comments: MNRAS accepted; 15 pages, 9 figures

CANDELS Sheds Light on the Environmental Quenching of Low-mass Galaxies

Authors: Yicheng Guo, Eric F. Bell, Yu Lu, David C. Koo, S. M. Faber, Anton M. Koekemoer, Peter Kurczynski, Seong-Kook Lee, Casey Papovich, Zhu Chen, Avishai Dekel, Henry C. Ferguson, Adriano Fontana, Mauro Giavalisco, Dale D. Kocevski, Hooshang Nayyeri, Pablo G. Pérez-González, Janine Pforr, Aldo Rodríguez-Puebla, Paola Santini

Abstract: We investigate the environmental quenching of galaxies, especially those with stellar masses (M*)$<10^{9.5} M_\odot$, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test {\it whether or not} such a dwarf QG--massive central galaxy co… ▽ More We investigate the environmental quenching of galaxies, especially those with stellar masses (M*)$<10^{9.5} M_\odot$, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test {\it whether or not} such a dwarf QG--massive central galaxy connection exists beyond the local universe. To this purpose, we only need a statistically representative, rather than a complete, sample of low-mass galaxies, which enables our study to $z\gtrsim1.5$. For each low-mass galaxy, we measure the projected distance ($d_{proj}$) to its nearest massive neighbor (M*$>10^{10.5} M_\odot$) within a redshift range. At a given redshift and M*, the environmental quenching effect is considered to be observed if the $d_{proj}$ distribution of QGs ($d_{proj}^Q$) is significantly skewed toward lower values than that of star-forming galaxies ($d_{proj}^{SF}$). For galaxies with $10^{8} M_\odot < M* < 10^{10} M_\odot$, such a difference between $d_{proj}^Q$ and $d_{proj}^{SF}$ is detected up to $z\sim1$. Also, about 10\% of the quenched galaxies in our sample are located between two and four virial radii ($R_{Vir}$) of the massive halos. The median projected distance from low-mass QGs to their massive neighbors, $d_{proj}^Q / R_{Vir}$, decreases with satellite M* at $M* \lesssim 10^{9.5} M_\odot$, but increases with satellite M* at $M* \gtrsim 10^{9.5} M_\odot$. This trend suggests a smooth, if any, transition of the quenching timescale around $M* \sim 10^{9.5} M_\odot$ at $0.5<z<1.0$. △ Less

Submitted 16 May, 2017; v1 submitted 4 May, 2017; originally announced May 2017.

Comments: 8 pages, 5 figures. ApJL accepted. Typos corrected

The Galaxy-Halo Connection Over The Last 13.3 Gyrs

Authors: Aldo Rodriguez-Puebla, Joel R. Primack, Vladimir Avila-Reese, S. M. Faber

Abstract: We present new determinations of the stellar-to-halo mass relation (SHMR) at $z=0-10$ that match the evolution of the galaxy stellar mass function, the SFR$-M_*$ relation,and the cosmic star formation rate. We utilize a compilation of 40 observational studies from the literature and correct them for potential biases. Using our robust determinations of halo mass assembly and the SHMR, we infer star… ▽ More We present new determinations of the stellar-to-halo mass relation (SHMR) at $z=0-10$ that match the evolution of the galaxy stellar mass function, the SFR$-M_*$ relation,and the cosmic star formation rate. We utilize a compilation of 40 observational studies from the literature and correct them for potential biases. Using our robust determinations of halo mass assembly and the SHMR, we infer star formation histories, merger rates, and structural properties for average galaxies, combining star-forming and quenched galaxies. Our main findings: (1) The halo mass $M_{50}$ above which 50\% of galaxies are quenched coincides with sSFR/sMAR$\sim1$, where sMAR is the specific halo mass accretion rate. (2) $M_{50}$ increases with redshift, presumably due to cold streams being more efficient at high redshift while virial shocks and AGN feedback become more relevant at lower redshifts. (3) The ratio sSFR/sMAR has a peak value, which occurs around $M_{\rm vir}\sim2\times10^{11}M_{\odot}$. (4) The stellar mass density within 1 kpc, $Σ_1$, is a good indicator of the galactic global sSFR. (5) Galaxies are statistically quenched after they reach a maximum in $Σ_1$, consistent with theoretical expectations of the gas compaction model; this maximum depends on redshift. (6) In-situ star formation is responsible for most galactic stellar mass growth, especially for lower-mass galaxies. (7) Galaxies grow inside out. The marked change in the slope of the size--mass relation when galaxies became quenched, from $d\log R_{\rm eff}/d\log M_*\sim0.35$ to $\sim2.5$, could be the result of dry minor mergers. △ Less

Submitted 15 May, 2017; v1 submitted 13 March, 2017; originally announced March 2017.

Comments: 42 pages, 31 Figures. Accepted for publication in MNRAS. Typos corrected. The results presented in this paper are available via an interactive website: https://132.248.1.39/galaxy/galaxy_halo.html

Journal ref: 2017, MNRAS, 470, 651

Properties of Dark Matter Halos as a Function of Local Environment Density

Authors: Christoph T Lee, Joel R. Primack, Peter Behroozi, Aldo Rodriguez-Puebla, Doug Hellinger, Avishai Dekel

Abstract: We study how properties of discrete dark matter halos depend on halo environment, characterized by the mass density around the halos on scales from 0.5 to 16 $h^{-1}{\rm Mpc}$. We find that low mass halos (those less massive than the characteristic mass $M_{\rm C}$ of halos collapsing at a given epoch) in high-density environments have lower accretion rates, lower spins, higher concentrations, and… ▽ More We study how properties of discrete dark matter halos depend on halo environment, characterized by the mass density around the halos on scales from 0.5 to 16 $h^{-1}{\rm Mpc}$. We find that low mass halos (those less massive than the characteristic mass $M_{\rm C}$ of halos collapsing at a given epoch) in high-density environments have lower accretion rates, lower spins, higher concentrations, and rounder shapes than halos in median density environments. Halos in median and low-density environments have similar accretion rates and concentrations, but halos in low density environments have lower spins and are more elongated. Halos of a given mass in high-density regions accrete material earlier than halos of the same mass in lower-density regions. All but the most massive halos in high-density regions are losing mass (i.e., being stripped) at low redshifts, which causes artificially lowered NFW scale radii and increased concentrations. Tidal effects are also responsible for the decreasing spins of low mass halos in high density regions at low redshifts $z < 1$, by preferentially removing higher angular momentum material from halos. Halos in low-density regions have lower than average spins because they lack nearby halos whose tidal fields can spin them up. We also show that the simulation density distribution is well fit by an Extreme Value Distribution, and that the density distribution becomes broader with cosmic time. △ Less

Submitted 6 January, 2017; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: 27 pages, 20 figures

Halo and Subhalo Demographics with Planck Cosmological Parameters: Bolshoi-Planck and MultiDark-Planck Simulations

Authors: Aldo Rodriguez-Puebla, Peter Behroozi, Joel Primack, Anatoly Klypin, Christoph Lee, Doug Hellinger

Abstract: We report and provide fitting functions for the abundance of dark matter halos and subhalos as a function of mass, circular velocity, and redshift from the new Bolshoi-Planck and MultiDark-Planck $Λ$CDM cosmological simulations, based on the Planck cosmological parameters. We also report the halo mass accretion rates, which may be connected with galaxy star formation rates. We show that the higher… ▽ More We report and provide fitting functions for the abundance of dark matter halos and subhalos as a function of mass, circular velocity, and redshift from the new Bolshoi-Planck and MultiDark-Planck $Λ$CDM cosmological simulations, based on the Planck cosmological parameters. We also report the halo mass accretion rates, which may be connected with galaxy star formation rates. We show that the higher cosmological matter density of the Planck parameters compared with the WMAP parameters leads to higher abundance of massive halos at high redshifts. We find that the median halo spin parameter $λ_{\rm B} = J(2M_{\rm vir}R_{\rm vir}V_{\rm vir})^{-1}$ is nearly independent of redshift, leading to predicted evolution of galaxy sizes that is consistent with observations, while the significant decrease with redshift in median $λ_{\rm P} = J|E|^{-1/2}G^{-1}M^{-5/2}$ predicts more decrease in galaxy sizes than is observed. Using the Tully-Fisher and Faber-Jackson relations between galaxy velocity and mass, we show that a simple model of how galaxy velocity is related to halo maximum circular velocity leads to increasing overprediction of cosmic stellar mass density as redshift increases beyond redshifts $z\sim1$, implying that such velocity-mass relations must change at redshifts $z>1$. By making a realistic model of how observed galaxy velocities are related to halo circular velocity, we show that recent optical and radio observations of the abundance of galaxies are in good agreement with our $Λ$CDM simulations. Our halo demographics are based on updated versions of the \rockstar\ and \ctrees\ codes, and this paper includes appendices explaining all of their outputs. This paper is an introduction to a series of related papers presenting other analyses of the Bolshoi-Planck and MultiDark-Planck simulations. △ Less

Submitted 14 March, 2017; v1 submitted 15 February, 2016; originally announced February 2016.

Comments: 27 pages, 26 figures. To match the published version at MNRAS. Rockstar and Consistent Trees outputs are available at http://hipacc.ucsc.edu/Bolshoi/MergerTrees.html

Journal ref: 2016MNRAS, 462, 893

Is Main Sequence Galaxy Star Formation Controlled by Halo Mass Accretion?

Authors: Aldo Rodriguez-Puebla, Joel R. Primack, Peter Behroozi, S. M. Faber

Abstract: The galaxy stellar-to-halo mass relation (SHMR) is nearly time-independent for $z < 4$. We therefore construct a time-independent SHMR model for central galaxies, wherein the in-situ star formation rate (SFR) is determined by the halo mass accretion rate (MAR), which we call Stellar-Halo Accretion Rate Coevolution (SHARC). We show that the $\sim0.3$ dex dispersion of the halo MAR matches the obser… ▽ More The galaxy stellar-to-halo mass relation (SHMR) is nearly time-independent for $z < 4$. We therefore construct a time-independent SHMR model for central galaxies, wherein the in-situ star formation rate (SFR) is determined by the halo mass accretion rate (MAR), which we call Stellar-Halo Accretion Rate Coevolution (SHARC). We show that the $\sim0.3$ dex dispersion of the halo MAR matches the observed dispersion of the SFR on the star-formation main sequence (MS). In the context of "bathtub"-type models of galaxy formation, SHARC leads to mass-dependent constraints on the relation between SFR and MAR. Despite its simplicity and the simplified treatment of mass growth from mergers, the SHARC model is likely to be a good approximation for central galaxies with $M_*=10^9- 10^{10.5}M_\odot$ that are on the MS, representing most of the star formation in the Universe. SHARC predictions agree with observed SFRs for galaxies on the MS at low redshifts, agree fairly well at $z\sim4$, but exceed observations at $z>4$. Assuming that the interstellar gas mass is constant for each galaxy (the "equilibrium condition" in bathtub models), the SHARC model allows calculation of net mass loading factors for inflowing and outflowing gas. With assumptions about preventive feedback based on simulations, SHARC allows calculation of galaxy metallicity evolution. If galaxy SFRs indeed track halo MARs, especially at low redshifts, that may help explain the success of models linking galaxy properties to halos (including age-matching) and the similarities between two-halo galaxy conformity and halo mass accretion conformity. △ Less

Submitted 27 October, 2015; v1 submitted 19 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS. No substantial changes. Abstract slightly modified in order to match the MNRAS Abstract length limit of 250 words. 18 pages, 12 figures

The stellar-to-halo mass relations of local galaxies segregates by color

Authors: A. Rodriguez-Puebla, V. Avila-Reese, X. Yang, S. Foucaud, N. Drory, Y. P. Jing

Abstract: We derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the… ▽ More We derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the value of red centrals. 2) The intrinsic scatters of the SHMRs of red and blue centrals are ~0.14 and ~0.11dex, respectively. The intrinsic scatter of the average SHMR of all central galaxies changes from ~0.20dex to ~0.14dex in the 11.3<logMh<15 range. 3) The fraction of halos hosting blue centrals at Mh=1E11Msun is 87%, but at 2x1E12Msun decays to ~20%, approaching to a few per cents at higher masses. The characteristic mass at which this fraction is the same for blue and red galaxies is Mh~7x1E11Msun. Our results suggest that the SHMR of central galaxies at large masses is shaped by halo mass quenching (likely through shock virial heating and AGN feedback), but group richness also plays an important role: central galaxies living in less dense environments quenched their star formation later or did not quench it yet. At low masses, processes that delay star formation without invoking too strong supernova-driven outflows could explain the high Ms/Mh ratios of blue centrals as compared to those of the scarce red centrals. △ Less

Submitted 6 January, 2015; v1 submitted 22 August, 2014; originally announced August 2014.

Comments: Accepted for publication in ApJ (20/Nov/2014). 23 pages, 16 figures. This version has substantial changes after Referee's report. Parts of Section 5 from our previous version will be publish in a follow-up paper

Central galaxies in different environments: Do they have similar properties?

Authors: Ivan Lacerna, Aldo Rodriguez-Puebla, Vladimir Avila-Reese, Hector M. Hernandez-Toledo

Abstract: We perform an exhaustive comparison among central galaxies from SDSS catalogs in different local environments at 0.01<=z<=0.08. The central galaxies are separated into two categories: group centrals (host halos containing satellites) and field centrals (host halos without satellites). From the latter, we select other two subsamples: isolated centrals and bright field centrals, both with the same m… ▽ More We perform an exhaustive comparison among central galaxies from SDSS catalogs in different local environments at 0.01<=z<=0.08. The central galaxies are separated into two categories: group centrals (host halos containing satellites) and field centrals (host halos without satellites). From the latter, we select other two subsamples: isolated centrals and bright field centrals, both with the same magnitude limit. The stellar mass (Ms) distributions of the field and group central galaxies are different, which explains why in general the field central galaxies are mainly located in the blue cloud/star forming regions, whereas the group central galaxies are strongly biased to the red sequence/passive regions. The isolated centrals occupy the same regions as the bright field centrals since both populations have similar Ms distributions. At parity of Ms, the color and specific star formation rate (sSFR) distributions of the samples are similar, specially between field and group centrals. Furthermore, we find that the stellar-to-halo mass (Ms-Mh) relation of isolated galaxies does not depend on the color, sSFR and morphological type. For systems without satellites, the Ms-Mh relation steepens at high halo masses compared to group centrals, which is a consequence of assuming a one-to-one relation between group total stellar mass and halo mass. Under the same assumption, the scatter around the Ms-Mh relation of centrals with satellites increases with halo mass. Our results suggest that the mass growth of central galaxies is mostly driven by the halo mass, with environment and mergers playing a secondary role. △ Less

Submitted 13 April, 2014; v1 submitted 25 November, 2013; originally announced November 2013.

Comments: 17 pages, 11 figures after last Referee's report. Accepted for publication in ApJ

Journal ref: 2014, ApJ, 788, 29

Simulations of isolated dwarf galaxies formed in dark matter halos with different mass assembly histories

Authors: A. González-Samaniego, P. Colín, V. Avila-Reese, A. Rodríguez-Puebla, O. Valenzuela

Abstract: We present high-resolution N-body/hydrodynamics simulations of dwarf galaxies formed in isolated CDM halos with the same virial mass, Mv~2.5x10^10 Msun at z=0, in order to (1) study the mass assembly histories (MAHs) of the halo, stars, and gas components, and (2) explore the effects of the halo MAHs on the stellar/baryonic assembly of the simulated dwarfs and on their z~0 properties. Overall, the… ▽ More We present high-resolution N-body/hydrodynamics simulations of dwarf galaxies formed in isolated CDM halos with the same virial mass, Mv~2.5x10^10 Msun at z=0, in order to (1) study the mass assembly histories (MAHs) of the halo, stars, and gas components, and (2) explore the effects of the halo MAHs on the stellar/baryonic assembly of the simulated dwarfs and on their z~0 properties. Overall, the simulated dwarfs are roughly consistent with observations. Our main results are: a) The stellar-to-halo mass ratio is ~0.01 and remains roughly constant since z~1 (the stellar MAHs follow closely the halo MAHs), with a smaller value at higher z's for those halos that assemble their mass later. b) The evolution of the galaxy gas fraction, fg, is episodic and higher, most of the time, than the stellar fraction. When fg decreases (increases), the gas fraction in the halo typically increases (decreases), showing that the SN driven outflows play an important role in regulating the gas fractions -and hence the SFR- of the dwarfs. However, in most cases, an important fraction of the gas escapes the virial radius, Rv; at z=4 the total baryon fraction inside Rv is 1.5-2 times smaller than the universal one, while at z=0 is 2-6 times smaller, with the earlier assembled halos ejecting more gas. c) The SF histories are episodic with changes in the SFRs of factors 2-10 on average. d) Although the dwarfs formed in late assembled halos show more extended SF histories, their z~0 SFRs are still below the ones measured for local isolated dwarfs. e) The effects of baryons on Mv are such that at almost any time Mv is 10-20% smaller than the corresponding Mv obtained in pure N-body simulations. Our results suggest that rather than increasing the strength of the SN-driven outflows, processes that reduce the SF efficiency even more will help to solve the potential issues faced by the CDM-based simulations of dwarfs. △ Less

Submitted 27 March, 2014; v1 submitted 21 August, 2013; originally announced August 2013.

Comments: 14 pages, 12 figures. ApJ, published. Minor changes after final Referee's report

Journal ref: ApJ, 785, 58 (2014)

The clustering of ALFALFA galaxies: dependence on HI mass, relationship to optical samples & clues on host halo properties

Authors: Emmanouil Papastergis, Riccardo Giovanelli, Martha P. Haynes, Aldo Rodríguez-Puebla, Michael G. Jones

Abstract: We use a sample of ~6000 galaxies detected by the Arecibo Legacy Fast ALFA (ALFALFA) 21cm survey, to measure the clustering properties of HI-selected galaxies. We find no convincing evidence for a dependence of clustering on the galactic atomic hydrogen (HI) mass, over the range M_HI ~ 10^{8.5} - 10^{10.5} M_sun. We show that previously reported results of weaker clustering for low-HI mass galaxie… ▽ More We use a sample of ~6000 galaxies detected by the Arecibo Legacy Fast ALFA (ALFALFA) 21cm survey, to measure the clustering properties of HI-selected galaxies. We find no convincing evidence for a dependence of clustering on the galactic atomic hydrogen (HI) mass, over the range M_HI ~ 10^{8.5} - 10^{10.5} M_sun. We show that previously reported results of weaker clustering for low-HI mass galaxies are probably due to finite-volume effects. In addition, we compare the clustering of ALFALFA galaxies with optically selected samples drawn from the Sloan Digital Sky Survey (SDSS). We find that HI-selected galaxies cluster more weakly than even relatively optically faint galaxies, when no color selection is applied. Conversely, when SDSS galaxies are split based on their color, we find that the correlation function of blue optical galaxies is practically indistinguishable from that of HI-selected galaxies. At the same time, SDSS galaxies with red colors are found to cluster significantly more than HI-selected galaxies, a fact that is evident in both the projected as well as the full two-dimensional correlation function. A cross-correlation analysis further reveals that gas-rich galaxies "avoid" being located within ~3 Mpc of optical galaxies with red colors. Next, we consider the clustering properties of halo samples selected from the Bolshoi LambdaCDM simulation. A comparison with the clustering of ALFALFA galaxies suggests that galactic HI mass is not tightly related to host halo mass, and that a sizable fraction of subhalos do not host HI galaxies. Lastly, we find that we can recover fairly well the correlation function of HI galaxies by just excluding halos with low spin parameter. This finding lends support to the hypothesis that halo spin plays a key role in determining the gas content of galaxies. △ Less

Submitted 12 August, 2013; originally announced August 2013.

Comments: Accepted to ApJ. 27 figures, 1 table; 50 pages in aastex (equivalent to 19 pages of emulateapj)

The Massive Satellite Population of Milky-Way Sized Galaxies

Authors: A. Rodriguez-Puebla, V. Avila-Reese, N. Drory

Abstract: Several occupational distributions for satellite galaxies more massive than ms~4E7 Msun around MW-sized hosts are presented and used to predict their internal dynamics. For this, a galaxy group mock catalog is constructed on the basis of (sub)halo-to-stellar mass relations constrained with the observed Galaxy Stellar Mass Function decomposed into centrals and satellites, and the 2-point correlatio… ▽ More Several occupational distributions for satellite galaxies more massive than ms~4E7 Msun around MW-sized hosts are presented and used to predict their internal dynamics. For this, a galaxy group mock catalog is constructed on the basis of (sub)halo-to-stellar mass relations constrained with the observed Galaxy Stellar Mass Function decomposed into centrals and satellites, and the 2-point correlation functions. For the MW-sized hosts we find that: 6.6% have 2 satellites in the SMC-LMC mass range; their probabilities to have 1, 2 or 3 satellites equal or larger than the LMC, the SMC or Sgr, respectively, are 0. 26, 0. 14, and 0. 14; those hosts with 3 satellites >= Sgr (as the MW) are among the most common cases, though the most and 2nd most massive satellites in these systems are SMALLER than the LMC and SMC by 0.7 and 0.8 dex, respectively. The cumulative satellite mass function is quite broad, the particular configuration of the MW being of low frequency but NOT an outlier (it is within the upper 1sig yet). For the MW satellite configuration, we constrain the halo mass to be Mh>1.4E12 Msun (1sig level). We find that the abundance of massive subhalos should agree with the abundance of massive satellites in all MW-sized hosts (there is not a massive satellite missing problem for LCDM). However, the vmax of the subhalos of satellites with ms<1E8 Msun is systematically larger than the vmax inferred from current observational studies of the MW bright dwarfs; at difference of previous works, this conclusion is based on an analysis of the overall population of MW-sized galaxies. Some pieces of evidence suggest that the issue could refer only to satellite dwarfs but not to central dwarfs; then, environmental processes associated to dwarfs inside host halos combined with SN-driven core expansion could be at the basis of the lowering of vmax.[Abridged] △ Less

Submitted 25 June, 2013; v1 submitted 18 June, 2013; originally announced June 2013.

Comments: 10 pages, 6 figures. Accepted for publication in ApJ. arXiv abstract improved

Journal ref: 2013ApJ,773,172

The Galaxy-Halo/Subhalo connection: Mass relations and implications for some satellite occupational distributions

Authors: A. Rodriguez-Puebla, V. Avila-Reese, N. Drory

Abstract: We infer the local stellar-to-halo/subhalo mass relations (MRs) for central and satellite galaxies separately. We constraint this relations by using several combinations of observational data, consisting of the total galaxy stellar mass function, its decomposition into centrals and satellites, and the projected two-point correlation functions measured in different stellar mass bins. The difference… ▽ More We infer the local stellar-to-halo/subhalo mass relations (MRs) for central and satellite galaxies separately. We constraint this relations by using several combinations of observational data, consisting of the total galaxy stellar mass function, its decomposition into centrals and satellites, and the projected two-point correlation functions measured in different stellar mass bins. The differences among the resulting MRs are within the model-fit uncertainties (which are very small, smaller than the intrinsic scatter between galaxy and halo mass), no matter what combination of data are used. We also constrain the values of the intrinsic scatter around the central-halo (CH) and satellite-subhalo (SS) MRs assuming them to be constant: sigma_c= 0.168+-0.051 dex and sigma_s = 0.172+-0.057 dex, respectively. The CH and SS MRs are actually different, in particular when we take the subhalo mass at the present-day epoch instead of at their accretion time. We calculate several population statistics. We find that the central galaxy is not on average within the mass distribution of the most-massive satellite, even for cluster-sized halos, i.e., centrals are not a mere realization of the high-end of the satellite mass function; however for > 3x10^13 msun halos, ~15% of centrals could be. We also find that the probabilities of MW-sized halos of having N Magellanic-Clouds (MCs)-sized satellites agree well with observational measures; for a halo mass of 2x10^12 msun, the probability to have 2 MCs is 5.4%, but if we exclude those systems with satellites larger than the MCs, then the probability decreases to <2.2%. △ Less

Submitted 26 March, 2013; v1 submitted 31 January, 2013; originally announced February 2013.

Comments: 18 pages, 7 figures. Accepted for publication in ApJ. Some references were added

Journal ref: 2013ApJ,767,92

The stellar-subhalo mass relation of satellite galaxies

Authors: A. Rodriguez-Puebla, N. Drory, V. Avila-Reese

Abstract: We extend the abundance matching technique (AMT) to infer the satellite-subhalo and central-halo mass relations (MRs) of galaxies, as well as the corresponding satellite conditional mass functions (CMFs). We use the observed galaxy stellar mass function (GSMF) decomposed into centrals and satellites and the LCDM halo/subhalo mass functions as inputs. We explore the effects of defining the subhalo… ▽ More We extend the abundance matching technique (AMT) to infer the satellite-subhalo and central-halo mass relations (MRs) of galaxies, as well as the corresponding satellite conditional mass functions (CMFs). We use the observed galaxy stellar mass function (GSMF) decomposed into centrals and satellites and the LCDM halo/subhalo mass functions as inputs. We explore the effects of defining the subhalo mass at the time of accretion (m_acc) vs. at the time of observation (m_obs). We test the standard assumption that centrals and satellites follow the same MRs, showing that this assumption leads to predictions in disagreement with observations, specially for m_obs. Instead, when the satellite-subhalo MRs are constrained following our AMT, they are always different from the central-halo MR: the smaller the stellar mass (Ms), the less massive is the subhalo of satellites as compared to the halo of centrals of the same Ms. On average, for Ms<2x10^11Msol, the dark mass of satellites decreased by 60-65% with respect to their masses at accretion time. The resulting MRs for both definitions of subhalo mass yield satellite CMFs in agreement with observations. Also, when these MRs are used in a HOD model, the predicted correlation functions agree with observations. We show that the use of m_obs leads to less uncertain MRs than m_acc, and discuss implications of the obtained satellite-subhalo MR. For example, we show that the tension between abundance and dynamics of MW satellites in LCDM gives if the slope of the GSMF faint-end slope upturns to -1.6. △ Less

Submitted 3 July, 2012; v1 submitted 3 April, 2012; originally announced April 2012.

Comments: 13, pages, 4 figures. Accepted for publication in ApJ. Minor changes to previous version

Journal ref: 2012ApJ,756,2