Search

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., Comment: 16 pages, 7 figures. Extended contribution to the original proceeding for the XVII Latin American Regional IAU Meeting

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., Comment: 30 pages, 28 figures, submitted to MNRAS

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., Comment: 17 pages, 12 figures

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_*$., Comment: MNRAS accepted, 26 pages

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\Upsilon_{\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\Upsilon_{\star}$ and to an increase of the light-weighted radius. The late reversing of $\nabla\Upsilon_{\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., Comment: 17 pages, 14 figures. Submitted to MNRAS, after Referee Report

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

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\'ersic 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^82$, 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., Comment: 29 pages, 17 figures. Accepted to appear in The Astronomical Journal

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/., Comment: 66 pages, 29 figures, 18 Tables, submitted to ApJS

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., Comment: 23 pages, 20 figures, 4 tables, accepted for publication in MNRAS

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

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., Comment: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

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., Comment: In press, 20 pages, 19 figures. Some corrections in text, references added. Results unchanged

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 $\sigma_{\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 $\sigma_{\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., Comment: Accepted for publication in MNRAS. Comments still welcome

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., Comment: In press, 21 pages, 15 figures. Minor corrections in text, results unchanged. Supplementary material available at ancillary files

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., Comment: Submitted to ApJ, 20 pages, 12 figures. Comments are welcome. This is an improved version after the first Referee report. Results remained unchanged

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., Comment: 13 pages, 12 figures, submitted to MNRAS

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}+\sigma_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $\sigma_{\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}-\sigma_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(\Upsilon_{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, $\Upsilon_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, $\Upsilon_{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., Comment: Submitted to ApJ. 21 pages, 12 figures

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

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., Comment: 21 pages, 9 figures, ApJ accepted

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., Comment: 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_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. $\Delta\Sigma_1$ and $\Delta\langle\mu_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., Comment: 23 pages, 17 figures, published in MNRAS

We study the global star-formation rate (SFR) vs. stellar mass (M$_*$) correlation, and the spatially-resolved SFR surface density ($\Sigma_{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., Comment: 24 pages, 16 figues. Accepted to appear in MNRAS

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., Comment: Science white paper submitted to the Astro2020 Decadal Survey

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\'ersic 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., Comment: 20 pages, 15 figures, Submitted to MNRAS

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

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., Comment: 33 pages, 23 figures, accepted to publish to MNRAS

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, $\sigma$, and also a combination of $V_{rot}$ and $\sigma$ through the $S_{K}$ parameter defined as $S_{K}^2 = KV_{rot}^2 + \sigma^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., Comment: accepted by MNRAS

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., Comment: 23 pages

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

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\sigma$ 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., Comment: 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

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., Comment: 23 pages, 15 figures. Accepted for publication in ApJ. Minor corrections and more references added. The discussion was extended in section 6.1

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.

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

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2

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/, Comment: 48 pages, 14 figures, in press in RMxAA

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

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, $\Sigma_1$, is a good indicator of the galactic global sSFR. (5) Galaxies are statistically quenched after they reach a maximum in $\Sigma_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., Comment: 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

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., Comment: 27 pages, 20 figures

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 $\Lambda$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 $\lambda_{\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 $\lambda_{\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 $\Lambda$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., Comment: 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

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

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