Search

We analyze the Near Infrared ($\sim0.8-1\mu$m) rest-frame morphologies of galaxies with $\log M_*/M_\odot>9$ in the redshift range $010.5$) at $z\sim5$, and bulge-dominated galaxies also exist at these early epochs, confirming a rich and evolved morphological diversity of galaxies $\sim1$ Gyr after the Big Bang. Finally, we find that the morphology-quenching relation is already in place for massive galaxies at $z>3$, with massive quiescent galaxies ($\log M_*/M_\odot>10.5$) being predominantly bulge-dominated., Comment: Submitted to A&A, comments welcome

Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN population in the distant universe. We estimate the source properties of MIRI-selected objects by utilizing spectral energy distribution (SED) modelling, and classify them into star-forming galaxies (SF), SF-AGN mixed objects, and AGN. The source numbers of these types are 418, 111, and 31, respectively, from 4 MIRI pointings covering $\sim 9$ arcmin$^2$. The sample spans a redshift range of $\approx 0$--5. We derive the median SEDs for all three source types, respectively, and publicly release them. The median MIRI SED of AGN is similar to the typical SEDs of hot dust-obscured galaxies and Seyfert 2s, for which the mid-IR SEDs are dominantly from AGN-heated hot dust. Based on our SED-fit results, we estimate the black-hole accretion density (BHAD; i.e., total BH growth rate per comoving volume) as a function of redshift. At $z<3$, the resulting BHAD agrees with the X-ray measurements in general. At $z>3$, we identify a total of 27 AGN and SF-AGN mixed objects, leading to that our high-$z$ BHAD is substantially higher than the X-ray results ($\sim 0.5$ dex at $z \approx 3$--5). This difference indicates MIRI can identify a large population of heavily obscured AGN missed by X-ray surveys at high redshifts., Comment: ApJL in press, as part of the CEERS Focus Issue. 13 pages and 7 figures (median SEDs in Fig. 5 are available upon request)

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches., Comment: 16 pages, 7 figures, 1 table. Accepted to ApJL

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches., Comment: 16 pages, 7 figures, 1 table. Accepted to ApJL

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST., Comment: 21 Pages

We present the public release of the MULTIDARK-GALAXIES: three distinct galaxy catalogues derived from one of the Planck cosmology MULTIDARK simulations (i.e. MDPL2, with a volume of (1 h Gpc) and mass resolution of 1.5 × 109 hM) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite-mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by ongoing and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available. © 2018 The Author(s).

In current and future surveys, quasars play a key role. The new data will extend our knowledge of the Universe as it will be used to better constrain the cosmological model at redshift z > 1 via baryon acoustic oscillation and redshift space distortion measurements. Here, we present the first clustering study of quasars observed by the extended Baryon Oscillation Spectroscopic Survey. We measure the clustering of ~70 000 quasars located in the redshift range 0.9 < z < 2.2 that cover 1168 deg. We model the clustering and produce highfidelity quasar mock catalogues based on the BigMultiDark Planck simulation. Thus, we use a modified (sub)halo abundance matching model to account for the specificities of the halo population hosting quasars. We find that quasars are hosted by haloes with masses~1012.7M⊙ and their bias evolves from 1.54 (z = 1.06) to 3.15 (z = 1.98). Using the current extended Baryon Oscillation Spectroscopic Survey data, we cannot distinguish between models with different fractions of satellites. The high-fidelity mock light-cones, including properties of haloes hosting quasars, are made publicly available.© 2017 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass at . We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, Ultra Deep Survey, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with . Satellites for both star-forming and quiescent central galaxies ("centrals") have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at for , whereas it is only weakly significant at and . Therefore, conformity (and thus satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at . This is evidence that satellite quenching is connected to the star formation properties of the central galaxy as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star formation in the central galaxy even as the satellites are quenched.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication., Comment: Accepted for publication in ApJS. 30 pages, emulateapj format. Many figures. All data products, including images and catalogs available http://irsa.ipac.caltech.edu/data/SPITZER/SHELA

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass at . We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, Ultra Deep Survey, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with . Satellites for both star-forming and quiescent central galaxies ("centrals") have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at for , whereas it is only weakly significant at and . Therefore, conformity (and thus satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at . This is evidence that satellite quenching is connected to the star formation properties of the central galaxy as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star formation in the central galaxy even as the satellites are quenched.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication., Comment: Accepted for publication in ApJS. 30 pages, emulateapj format. Many figures. All data products, including images and catalogs available http://irsa.ipac.caltech.edu/data/SPITZER/SHELA

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used., Comment: 25 pages, 19 figures, 7 tables. Accepted for publication in MNRAS

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used., Comment: 25 pages, 19 figures, 7 tables. Accepted for publication in MNRAS

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency., Comment: Accepted for publication in ApJ; 25 pages, emulateapj format

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.54. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts., Comment: 26 pages, 23 figures. Accepted for publication on A&A

We present a semi-empirical model to infer the atomic and molecular hydrogen content of galaxies as a function of halo mass and time. Our model combines the SFR-halo mass-redshift relation (constrained by galaxy abundances) with inverted SFR-surface density relations to infer galaxy H I and H2 masses. We present gas scaling relations, gas fractions, and mass functions from z = 0 to z = 3 and the gas properties of galaxies as a function of their host halo masses. Predictions of our work include: 1) there is a ~ 0.2 dex decrease in the H I mass of galaxies as a function of their stellar mass since z = 1.5, whereas the H2 mass of galaxies decreases by > 1 dex over the same period. 2) galaxy cold gas fractions and H2 fractions decrease with increasing stellar mass and time. Galaxies with M* > 10^10 Msun are dominated by their stellar content at z < 1, whereas less-massive galaxies only reach these gas fractions at z = 0. We find the strongest evolution in relative gas content at z < 1.5. 3) the SFR to gas mass ratio decreases by an order of magnitude from z = 3 to z = 0. This is consistent with lower H2 fractions; these lower fractions in combination with smaller gas reservoirs correspond to decreased present-day galaxy SFRs. 4) an H2-based star- formation relation can simultaneously fuel the evolution of the cosmic star-formation and reproduce the observed weak evolution in the cosmic HI density. 5) galaxies residing in haloes with masses near 10^12 Msun are most efficient at obtaining large gas reservoirs and forming H2 at all redshifts. These two effects lie at the origin of the high star-formation efficiencies in haloes with the same mass., Comment: accepted for publication in MNRAS, 20 pages, 16 figures (+ 1 figure in appendix), data files are accessible through http://www.eso.org/~gpopping/Gergo_Poppings_Homepage/Data.html

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 03x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times., Comment: Published in ApJ. Structural parameter measurements publicly available at http://www.mpia-hd.mpg.de/homes/vdwel/3dhstcandels.html

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency., Comment: Accepted for publication in ApJ; 25 pages, emulateapj format

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.54. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts., Comment: 26 pages, 23 figures. Accepted for publication on A&A

We present a semi-empirical model to infer the atomic and molecular hydrogen content of galaxies as a function of halo mass and time. Our model combines the SFR-halo mass-redshift relation (constrained by galaxy abundances) with inverted SFR-surface density relations to infer galaxy H I and H2 masses. We present gas scaling relations, gas fractions, and mass functions from z = 0 to z = 3 and the gas properties of galaxies as a function of their host halo masses. Predictions of our work include: 1) there is a ~ 0.2 dex decrease in the H I mass of galaxies as a function of their stellar mass since z = 1.5, whereas the H2 mass of galaxies decreases by > 1 dex over the same period. 2) galaxy cold gas fractions and H2 fractions decrease with increasing stellar mass and time. Galaxies with M* > 10^10 Msun are dominated by their stellar content at z < 1, whereas less-massive galaxies only reach these gas fractions at z = 0. We find the strongest evolution in relative gas content at z < 1.5. 3) the SFR to gas mass ratio decreases by an order of magnitude from z = 3 to z = 0. This is consistent with lower H2 fractions; these lower fractions in combination with smaller gas reservoirs correspond to decreased present-day galaxy SFRs. 4) an H2-based star- formation relation can simultaneously fuel the evolution of the cosmic star-formation and reproduce the observed weak evolution in the cosmic HI density. 5) galaxies residing in haloes with masses near 10^12 Msun are most efficient at obtaining large gas reservoirs and forming H2 at all redshifts. These two effects lie at the origin of the high star-formation efficiencies in haloes with the same mass., Comment: accepted for publication in MNRAS, 20 pages, 16 figures (+ 1 figure in appendix), data files are accessible through http://www.eso.org/~gpopping/Gergo_Poppings_Homepage/Data.html

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 03x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times., Comment: Published in ApJ. Structural parameter measurements publicly available at http://www.mpia-hd.mpg.de/homes/vdwel/3dhstcandels.html

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and study their differences in the density profile, mass fraction and two-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10 per cent level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5 per cent when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ 10 different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6D phase space analysis of particles. We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructures the nearer they are to the host haloes, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally, subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We develop a theoretical framework that combines measurements of galaxy-galaxy lensing, galaxy clustering, and the galaxy stellar mass function in a self-consistent manner. While considerable effort has been invested in exploring each of these probes individually, attempts to combine them are still in their infancy despite the potential of such combinations to elucidate the galaxy-dark matter connection, to constrain cosmological parameters, and to test the nature of gravity. In this paper, we focus on a theoretical model that describes the galaxy-dark matter connection based on standard halo occupation distribution techniques. Several key modifications enable us to extract additional parameters that determine the stellar-to-halo mass relation and to simultaneously fit data from multiple probes while allowing for independent binning schemes for each probe. In a companion paper, we demonstrate that the model presented here provides an excellent fit to galaxy-galaxy lensing, galaxy clustering, and stellar mass functions measured in the COSMOS survey from z=0.2 to z=1.0. We construct mock catalogs from numerical simulations to investigate the effects of sample variance and covariance on each of the three probes. Finally, we analyze and discuss how trends in each of the three observables impact the derived parameters of the model. In particular, we investigate the various features of the observed galaxy stellar mass function (low-mass slope, plateau, knee, and high-mass cut-off) and show how each feature is related to the underlying relationship between stellar and halo mass. We demonstrate that the observed plateau feature in the stellar mass function at Mstellar~2x10^10 Msun is due to the transition that occurs in the stellar-to-halo mass relation at Mhalo ~ 10^12 Msun from a low-mass power-law regime to a sub-exponential function at higher stellar mass., Comment: 21 pages. Accepted to ApJ

We develop a theoretical framework that combines measurements of galaxy-galaxy lensing, galaxy clustering, and the galaxy stellar mass function in a self-consistent manner. While considerable effort has been invested in exploring each of these probes individually, attempts to combine them are still in their infancy despite the potential of such combinations to elucidate the galaxy-dark matter connection, to constrain cosmological parameters, and to test the nature of gravity. In this paper, we focus on a theoretical model that describes the galaxy-dark matter connection based on standard halo occupation distribution techniques. Several key modifications enable us to extract additional parameters that determine the stellar-to-halo mass relation and to simultaneously fit data from multiple probes while allowing for independent binning schemes for each probe. In a companion paper, we demonstrate that the model presented here provides an excellent fit to galaxy-galaxy lensing, galaxy clustering, and stellar mass functions measured in the COSMOS survey from z=0.2 to z=1.0. We construct mock catalogs from numerical simulations to investigate the effects of sample variance and covariance on each of the three probes. Finally, we analyze and discuss how trends in each of the three observables impact the derived parameters of the model. In particular, we investigate the various features of the observed galaxy stellar mass function (low-mass slope, plateau, knee, and high-mass cut-off) and show how each feature is related to the underlying relationship between stellar and halo mass. We demonstrate that the observed plateau feature in the stellar mass function at Mstellar~2x10^10 Msun is due to the transition that occurs in the stellar-to-halo mass relation at Mhalo ~ 10^12 Msun from a low-mass power-law regime to a sub-exponential function at higher stellar mass., Comment: 21 pages. Accepted to ApJ

Understanding the mechanisms that lead dense environments to host galaxies with redder colors, more spheroidal morphologies, and lower star formation rates than field populations remains an important problem. As most candidate processes ultimately depend on host halo mass, accurate characterizations of the local environment, ideally tied to halo mass estimates and spanning a range in halo mass and redshift, are needed. In this work, we present and test a rigorous, probabilistic method for assigning galaxies to groups based on precise photometric redshifts and X-ray-selected groups drawn from the COSMOS field. The groups have masses in the range 1013 <~ M 200c/M ☉ <~ 1014 and span redshifts 0 < z < 1. We characterize our selection algorithm via tests on spectroscopic subsamples, including new data obtained at the Very Large Telescope, and by applying our method to detailed mock catalogs. We find that our group member galaxy sample has a purity of 84% and completeness of 92% within 0.5{R_{200c}}. We measure the impact of uncertainties in redshifts and group centering on the quality of the member selection with simulations based on current data as well as future imaging and spectroscopic surveys. As a first application of our new group member catalog which will be made publicly available, we show that member galaxies exhibit a higher quenched fraction compared to the field at fixed stellar mass out to z ~ 1, indicating a significant relationship between star formation and environment at group scales. We also address the suggestion that dusty star-forming galaxies in such groups may impact the high-l power spectrum of the cosmic microwave background and find that such a population cannot explain the low power seen in recent Sunyaev-Zel'dovich measurements.

We develop a theoretical framework that combines measurements of galaxy-galaxy lensing, galaxy clustering, and the galaxy stellar mass function in a self-consistent manner. While considerable effort has been invested in exploring each of these probes individually, attempts to combine them are still in their infancy. These combinations have the potential to elucidate the galaxy-dark matter connection and the galaxy formation physics responsible for it, as well as to constrain cosmological parameters and to test the nature of gravity. In this paper, we focus on a theoretical model that describes the galaxy-dark matter connection based on standard halo occupation distribution techniques. Several key modifications enable us to extract additional parameters that determine the stellar-to-halo mass relation and to simultaneously fit data from multiple probes while allowing for independent binning schemes for each probe. We construct mock catalogs from numerical simulations to investigate the effects of sample variance and covariance for each probe. Finally, we analyze how trends in each of the three observables impact the derived parameters of the model. In particular, we investigate various features of the observed galaxy stellar mass function (low-mass slope, "plateau," knee, and high-mass cutoff) and show how each feature is related to the underlying relationship between stellar and halo mass. We demonstrate that the observed "plateau" feature in the stellar mass function at M * ~ 2 × 1010 M sun is due to the transition that occurs in the stellar-to-halo mass relation at Mh ~ 1012 M sun from a low-mass power-law regime to a sub-exponential function at higher stellar mass.

We calculate the probability that a Milky-Way (MW)-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the lambda cold dark matter model are in good agreement with observations of a large sample of Sloan Digital Sky Survey (SDSS) galaxies. Under the subhalo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and halo internal velocities v max, we make detailed comparisons to similar measurements using SDSS Data Release 7 data by Liu et al. Models and observational data give very similar probabilities for having zero, one, and two MC-like satellites. In both cases, MW luminosity hosts have just a ~10% chance of hosting two satellites similar to the MCs. In addition, we present a prediction for the probability for a host galaxy to have N sats satellite galaxies as a function of the magnitudes of both the host and satellite. This probability and its scaling with host properties is significantly different from that of mass-selected objects because of scatter in the mass-luminosity relation and because of variations in the star formation efficiency with halo mass.

We introduce a probabilistic approach to the problem of counting dwarf satellites around host galaxies in databases with limited redshift information. This technique is used to investigate the occurrence of satellites with luminosities similar to the Magellanic Clouds around hosts with properties similar to the Milky Way (MW) in the object catalog of the Sloan Digital Sky Survey (SDSS). Our analysis uses data from SDSS Data Release 7, selecting candidate MW-like hosts from the spectroscopic catalog and candidate analogs of the Magellanic Clouds from the photometric catalog. Our principal result is the probability for an MW-like galaxy to host N sat close satellites with luminosities similar to the Magellanic Clouds. We find that 81% of galaxies like the MW have no such satellites within a radius of 150 kpc, 11% have one, and only 3.5% of hosts have two. The probabilities are robust to changes in host and satellite selection criteria, background-estimation technique, and survey depth. These results demonstrate that the MW has significantly more satellites than a typical galaxy of its luminosity; this fact is useful for understanding the larger cosmological context of our home galaxy.

We use high-resolution Hydro+N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ∼106−7M⊙, total mass 1010M⊙) in Λ-dominated cold dark matter (CDM) and 2keV warm dark matter (WDM) cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the dark matter (DM) model, but proportionally to the SF efficiency, gas outflows lower the central mass density through ‘dynamical heating', such that all realizations have circular velocities <20 km s−1 at 500pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial colour-magnitude diagrams and star formation histories (SFHs) in order to directly compare to available observations. The simulated galaxies formed most of their stars in many ∼10Myr long bursts. The CDM galaxy has a global SFH, H i abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including ‘baryon physics' in simulations when (1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and (2) differentiating between CDM and non-standard models with different DM or power spectra

Spectroscopic+photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 0 < z < 3. Separating early- and late-type galaxies on the basis of star-formation activity, we confirm that early-type galaxies are on average smaller than late-type galaxies at all redshifts, and we find a significantly different rate of average size evolution at fixed galaxy mass, with fast evolution for the early-type population, $R_{\rm {eff}}\propto (1 + z)^{-1.48}$, and moderate evolution for the late-type population, $R_{\rm {eff}}\propto (1 + z)^{-0.75}$. The large sample size and dynamic range in both galaxy mass and redshift, in combination with the high fidelity of our measurements due to the extensive use of spectroscopic data, not only fortify previous results but also enable us to probe beyond simple average galaxy size measurements. At all redshifts the slope of the size-mass relation is shallow, $R_{\rm {eff}}\propto M_*^{0.22}$, for late-type galaxies with stellar mass >3 × 10$^{9} M_\odot$, and steep, $R_{\rm {eff}}\propto M_*^{0.75}$, for early-type galaxies with stellar mass >2 × 10$^{10} M_\odot$. The intrinsic scatter is ≲ 0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric but is skewed toward small sizes: at all redshifts and masses, a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10$^{11} M_\odot$), compact ($R_{\rm {eff}}$ < 2 kpc) early-type galaxies increases from z = 3 to z = 1.5-2 and then strongly decreases at later cosmic times.

Spectroscopic+photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 0 < z < 3. Separating early- and late-type galaxies on the basis of star-formation activity, we confirm that early-type galaxies are on average smaller than late-type galaxies at all redshifts, and we find a significantly different rate of average size evolution at fixed galaxy mass, with fast evolution for the early-type population, $R_{\rm {eff}}\propto (1 + z)^{-1.48}$, and moderate evolution for the late-type population, $R_{\rm {eff}}\propto (1 + z)^{-0.75}$. The large sample size and dynamic range in both galaxy mass and redshift, in combination with the high fidelity of our measurements due to the extensive use of spectroscopic data, not only fortify previous results but also enable us to probe beyond simple average galaxy size measurements. At all redshifts the slope of the size-mass relation is shallow, $R_{\rm {eff}}\propto M_*^{0.22}$, for late-type galaxies with stellar mass >3 × 10$^{9} M_\odot$, and steep, $R_{\rm {eff}}\propto M_*^{0.75}$, for early-type galaxies with stellar mass >2 × 10$^{10} M_\odot$. The intrinsic scatter is ≲ 0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric but is skewed toward small sizes: at all redshifts and masses, a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10$^{11} M_\odot$), compact ($R_{\rm {eff}}$ < 2 kpc) early-type galaxies increases from z = 3 to z = 1.5-2 and then strongly decreases at later cosmic times.