Search

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs., Comment: 17 pages, 7 figures

We explore empirical constraints on the statistical relationship between the radial size of galaxies and the radius of their host dark matter halos from $z\sim 0.1$--3 using the GAMA and CANDELS surveys. We map dark matter halo mass to galaxy stellar mass using relationships from abundance matching, applied to the Bolshoi-Planck dissipationless N-body simulation. We define SRHR$\equiv r_e/R_h$ as the ratio of galaxy radius to halo virial radius, and SRHR$\lambda \equiv r_e/(\lambda R_h)$ as the ratio of galaxy radius to halo spin parameter times halo radius. At $z\sim 0.1$, we find an average value of SRHR $\simeq 0.018$ and SRHR$\lambda \simeq 0.5$ with very little dependence on stellar mass. SRHR and SRHR$\lambda$ have a weak dependence on cosmic time since $z\sim 3$. SRHR shows a mild decrease over cosmic time for low mass galaxies, but increases slightly or does not evolve for more massive galaxies. We find hints that at high redshift ($z\sim 2$--3), SRHR$\lambda$ is lower for more massive galaxies, while it shows no significant dependence on stellar mass at $z\lesssim 0.5$. We find that for both the GAMA and CANDELS samples, at all redshifts from $z\sim 0.1$--3, the observed conditional size distribution in stellar mass bins is remarkably similar to the conditional distribution of $\lambda R_h$. We discuss the physical interpretation and implications of these results., Comment: 27 pages, 17 figures. submitted to MNRAS. little-h related error in v1 corrected, and updated in response to referee's comments

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$., Comment: Re-submitted to MNRAS after referee revisions, the main change is that the paper has been significantly streamlined and is shorter (25 pages and 8 figures, without appendices)

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies., Comment: 16 pages, 11 figures, accepted for publication in ApJ

We study the dust content of galaxies from $z=0$ to $z=9$ in semi-analytic models of galaxy formation that include new recipes to track the production and destruction of dust. We include condensation of dust in stellar ejecta, the growth of dust in the interstellar medium (ISM), the destruction of dust by supernovae and in the hot halo, and dusty winds and inflows. The rate of dust growth in the ISM depends on the metallicity and density of molecular clouds. Our fiducial model reproduces the relation between dust mass and stellar mass from $z=0$ to $z=7$, the number density of galaxies with dust masses less than $10^{8.3}\,\rm{M}_\odot$, and the cosmic density of dust at $z=0$. The model accounts for the double power-law trend between dust-to-gas (DTG) ratio and gas-phase metallicity of local galaxies and the relation between DTG ratio and stellar mass. The dominant mode of dust formation is dust growth in the ISM, except for galaxies with $M_*<10^7\,\rm{M}_\odot$, where condensation of dust in supernova ejecta dominates. The dust-to-metal ratio of galaxies depends on the gas-phase metallicity, unlike what is typically assumed in cosmological simulations. Model variants including higher condensation efficiencies, a fixed timescale for dust growth in the ISM, or no growth at all reproduce some of the observed constraints, but fail to simultaneously reproduce the shape of dust scaling relations and the dust mass of high-redshift galaxies., Comment: Accepted for publication in MNRAS

We study the correlation of galaxy structural properties with their location relative to the SFR-M* correlation, also known as the star formation "main sequence" (SFMS), in the CANDELS and GAMA surveys and in a semi-analytic model (SAM) of galaxy formation. We first study the distribution of median Sersic index, effective radius, star formation rate (SFR) density and stellar mass density in the SFR-M* plane. We then define a redshift dependent main sequence and examine the medians of these quantities as a function of distance from this main sequence, both above (higher SFRs) and below (lower SFRs). Finally, we examine the distributions of distance from the main sequence in bins of these quantities. We find strong correlations between all of these galaxy structural properties and the distance from the SFMS, such that as we move from galaxies above the SFMS to those below it, we see a nearly monotonic trend towards higher median Sersic index, smaller radius, lower SFR density, and higher stellar density. In the semi-analytic model, bulge growth is driven by mergers and disk instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via Active Galactic Nucleus (AGN) feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and AGN feedback plays a critical role in moving galaxies off of the SFMS., Comment: 25 pages, 17 figures. Accepted for publication in MNRAS

We present a sample of accreting supermassive black holes (SMBHs) in dwarf galaxies at $z<1$. We identify dwarf galaxies in the NEWFIRM Medium Band Survey with stellar masses $M_{\star}<3\times 10^{9} M_{\odot}$ that have spectroscopic redshifts from the DEEP2 survey and lie within the region covered by deep (flux limit of $\sim 5\times 10^{-17} - 6\times 10^{-16} \rm{erg \ cm}^{-2} \ \rm{s}^{-1}$) archival Chandra X-ray data. From our sample of $605$ dwarf galaxies, $10$ exhibit X-ray emission consistent with that arising from AGN activity. If black hole mass scales roughly with stellar mass, then we expect that these AGN are powered by SMBHs with masses of $\sim 10^5-10^6 \ M_{\odot}$ and typical Eddington ratios $\sim 5\%$. Furthermore, we find an AGN fraction consistent with extrapolations of other searches of $\sim 0.6-3\%$ for $10^9 \ M_{\odot} \leq M_{\star} \leq 3\times 10^{9} \ M_{\odot}$ and $0.1

Now that ALMA is reaching its full capabilities, observations of sub-mm emission line deep fields become feasible. We couple a semi-analytic model of galaxy formation with a radiative transfer code to make predictions for the luminosity function of CO J=1-0 up to CO J=6-5 and [CII] at redshifts z=0-6. We find that: 1) our model correctly reproduces the CO and [CII] emission of low- and high-redshift galaxies and reproduces the available constraints on the CO luminosity function at z<2.75; 2) we find that the CO and [CII] luminosity functions of galaxies increase from z = 6 to z = 4, remain relatively constant till z = 1 and rapidly decrease towards z = 0. The galaxies that are brightest in CO and [CII] are found at z~2; 3) the CO J=3-2 emission line is most favourable to study the CO luminosity and global H2 mass content of galaxies, because of its brightness and observability with currently available sub-mm and radio instruments; 4) the luminosity functions of high-J CO lines show stronger evolution than the luminosity functions of low-J CO lines; 5) our model barely reproduces the available constraints on the CO and [CII] luminosity function of galaxies at z>1.5 and the CO luminosity of individual galaxies at intermediate redshifts. We argue that this is driven by a lack of cold gas in galaxies at intermediate redshifts as predicted by cosmological simulations of galaxy formation., Comment: 20 pages, 12 figures, 3 tables, accepted for publication in MNRAS

To explore the connection between the global physical properties of galaxies and their far-infrared (FIR) spectral energy distributions (SEDs), we study the variation in the FIR SEDs of a set of hydrodynamically simulated galaxies that are generated by performing dust radiative transfer in post-processing. Our sample includes both isolated and merging systems at various stages of the merging process and covers infrared (IR) luminosities and dust masses that are representative of both low- and high-redshift galaxies. We study the FIR SEDs using principle component analysis (PCA) and find that 97\% of the variance in the sample can be explained by two principle components (PCs). The first PC characterizes the wavelength of the peak of the FIR SED, and the second encodes the breadth of the SED. We find that the coefficients of both PCs can be predicted well using a double power law in terms of the IR luminosity and dust mass, which suggests that these two physical properties are the primary determinants of galaxies' FIR SED shapes. Incorporating galaxy sizes does not significantly improve our ability to predict the FIR SEDs. Our results suggest that the observed redshift evolution in the effective dust temperature at fixed IR luminosity is not driven by geometry: the SEDs of $z \sim 2-3$ ultraluminous IR galaxies (ULIRGs) are cooler than those of local ULIRGs not because the high-redshift galaxies are more extended but rather because they have higher dust masses at fixed IR luminosity. Finally, based on our simulations, we introduce a two-parameter set of SED templates that depend on both IR luminosity and dust mass., Comment: Submitted to ApJ, comments welcome

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 \AA\ ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $\alpha=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $\alpha=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density., Comment: 27 pages, 17 figures, ApJ, in press

The distribution of galaxy morphological types is a key test for models of galaxy formation and evolution, providing strong constraints on the relative contribution of different physical processes responsible for the growth of the spheroidal components. In this paper, we make use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupting galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., based on results from high-resolution hydrodynamical simulations, and we show that these new implementations reduce the efficiency of bulge formation through mergers. In addition, we compare our model results with a variety of observational measurements of the fraction of spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present uncertainties in the data represent an important limitation to our understanding of spheroid formation. Our results indicate that the main tension between theoretical models and observations does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather from the distribution of galaxies of different morphological types, as a function of their stellar mass., Comment: MNRAS in press, 11 pages, 5 figures

We study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localised to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation., Comment: MNRAS Accepted, 17 pages, 11 figures, 3 Tables

We use semi-analytic models and cosmological merger trees to provide the initial conditions for multi-merger numerical hydrodynamic simulations, and exploit these simulations to explore the effect of galaxy interaction and merging on star formation (SF). We compute numerical realisations of twelve merger trees from z=1.5 to z=0. We include the effects of the large hot gaseous halo around all galaxies, following recent obervations and predictions of galaxy formation models. We find that including the hot gaseous halo has a number of important effects. Firstly, as expected, the star formation rate on long timescales is increased due to cooling of the hot halo and refuelling of the cold gas reservoir. Secondly, we find that interactions do not always increase the SF in the long term. This is partially due to the orbiting galaxies transferring gravitational energy to the hot gaseous haloes and raising their temperature. Finally we find that the relative size of the starburst, when including the hot halo, is much smaller than previous studies showed. Our simulations also show that the order and timing of interactions are important for the evolution of a galaxy. When multiple galaxies interact at the same time, the SF enhancement is less than when galaxies interact in series. All these effects show the importance of including hot gas and cosmologically motivated merger trees in galaxy evolution models., Comment: 19 pages, 15 figures, 6 tables. Accepted for publication in MNRAS

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here., Comment: 35 pages, 20 figures, accepted for publication in MNRAS

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied., Comment: 33 pages, 20 Figures, Accepted for publication in Astrophysical Journal

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts., Comment: Accepted to ApJ. 15 pages, 5 figures, 6 tables

We implement physically motivated recipes for partitioning cold gas into different phases (atomic, molecular, and ionized) in galaxies within semi-analytic models of galaxy formation based on cosmological merger trees. We then model the conversion of molecular gas into stars using empirical recipes motivated by recent observations. We explore the impact of these new recipes on the evolution of fundamental galaxy properties such as stellar mass, star formation rate (SFR), and gas and stellar phase metallicity. We present predictions for stellar mass functions, stellar mass vs. SFR relations, and cold gas phase and stellar mass-metallicity relations for our fiducial models, from redshift $z\sim 6$ to the present day. In addition we present predictions for the global SFR, mass assembly history, and cosmic enrichment history. We find that the predicted stellar properties of galaxies (stellar mass, SFR, metallicity) are remarkably insensitive to the details of the recipes used for partitioning gas into HI and H$_2$. We see significant sensitivity to the recipes for H$_2$ formation only in very low mass halos, which host galaxies that are not detectable with current observational facilities except very nearby. The properties of low-mass galaxies are also quite insensitive to the details of the recipe used for converting H$_2$ into stars, while the formation epoch of massive galaxies does depend on this significantly. (Abridged), Comment: 35 pages, 23 figures, submitted to MNRAS

Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M_* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies., Comment: Minor updates to match published version

We examine the spheroid growth and star formation quenching experienced by galaxies from z~3 to the present by studying the evolution with redshift of the quiescent and spheroid-dominated fractions of galaxies from the CANDELS and GAMA surveys. We compare the observed fractions with predictions from a semi-analytic model which includes prescriptions for bulge growth and AGN feedback due to mergers and disk instabilities. We facilitate direct morphological comparison by converting our model bulge-to-total stellar mass ratios to Sersic indices. We then subdivide our population into the four quadrants of the sSFR-Sersic index plane and study the buildup of each of these subpopulations. We find that the fraction of star forming disks declines steadily, while the fraction of quiescent spheroids builds up over cosmic time. The fractions of star forming spheroids and quiescent disks are both non-negligible, and stay nearly constant over the period we have studied, at about 10% and 15-20% respectively. Our model is qualitatively successful at reproducing the evolution of the two "main" populations (star forming disk-dominated galaxies and quiescent spheroid-dominated galaxies), and approximately reproduces the relative fractions of all four types, but predicts a stronger decline in star forming spheroids, and increase in quiescent disks, than seen in the observations. A model with an additional channel for bulge growth via disk instabilities agrees better overall with the observations than a model in which bulges may grow only through mergers. We study evolutionary tracks of some individual galaxies as they experience morphological transformation and quenching, and examine the importance of different physical drivers of this transformation (major and minor mergers and disk instabilities). We find that complex histories with multiple transformative events are the norm., Comment: 32 pages, 21 figures; Tables fixed, figure labels changed, references added. Accepted for publication in MNRAS

Modeling galaxy formation in a cosmological context presents one of the greatest challenges in astrophysics today, due to the vast range of scales and numerous physical processes involved. Here we review the current status of models that employ two leading techniques to simulate the physics of galaxy formation: semi-analytic models and numerical hydrodynamic simulations. We focus on a set of observational targets that describe the evolution of the global and structural properties of galaxies from roughly Cosmic High Noon ($z\sim 2-3$) to the present. Although minor discrepancies remain, overall, models show remarkable convergence between different methods and make predictions that are in qualitative agreement with observations. Modelers seem to have converged on a core set of physical processes that are critical for shaping galaxy properties. This core set includes cosmological accretion, strong stellar-driven winds that are more efficient at low masses, black hole feedback that preferentially suppresses star formation at high masses, and structural and morphological evolution through merging and environmental processes. However, all cosmological models currently adopt phenomenological implementations of many of these core processes, which must be tuned to observations. Many details of how these diverse processes interact within a hierarchical structure formation setting remain poorly understood. Emerging multi-scale simulations are helping to bridge the gap between stellar and cosmological scales, placing models on a firmer, more physically grounded footing. Concurrently, upcoming telescope facilities will provide new challenges and constraints for models, particularly by directly constraining inflows and outflows through observations of gas in and around galaxies., Comment: Submitted version of review article to appear in Vol. 53 of Annual Reviews of Astronomy and Astrophysics. 67 pages, 8 figures

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift., Comment: 28 pages, 23 figures, 2 appendices. Accepted for publication in MNRAS, August 7 2014

We present a new technique for overcoming confusion noise in deep far-infrared \Herschel space telescope images making use of prior information from shorter $\lambda<2$\micron wavelengths. For the deepest images obtained by \Herschels, the flux limit due to source confusion is about a factor of three brighter than the flux limit due to instrumental noise and (smooth) sky background. We have investigated the possibility of de-confusing simulated \Herschel PACS-160\micron images by using strong Bayesian priors on the positions and weak priors on the flux of sources. We find the blended sources and group them together and simultaneously fit their fluxes. We derive the posterior probability distribution function of fluxes subject to these priors through Monte Carlo Markov Chain (MCMC) sampling by fitting the image. Assuming we can predict FIR flux of sources based on ultraviolet-optical part of their SEDs to within an order of magnitude, the simulations show that we can obtain reliable fluxes and uncertainties at least a factor of three fainter than the confusion noise limit of $3\sigma_{c} $=2.7 mJy in our simulated PACS-160 image. This technique could in principle be used to mitigate the effects of source confusion in any situation where one has prior information of positions and plausible fluxes of blended sources. For \Herschel, application of this technique will improve our ability to constrain the dust content in normal galaxies at high redshift.

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $\sigma(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation., Comment: 31 pages, 24 figures, accepted for publication in ApJ

Both numerical hydrodynamic and semi-analytic cosmological models of galaxy formation struggle to match observed star formation histories of galaxies in low mass halos (M$_{\rm{H}} \lesssim 10^{11} M_\odot$), predicting more star formation at high redshift and less star formation at low redshift than observed. The fundamental problem is that galaxies' gas accretion and star formation rates are too closely coupled in the models: the accretion rate largely drives the star formation rate. Observations point to gas accretion rates that outpace star formation at high redshift, resulting in a buildup of gas and a delay in star formation until lower redshifts. We present three empirical adjustments of standard recipes in a semi-analytic model motivated by three physical scenarios that could cause this decoupling: 1) the mass-loading factors of outflows driven by stellar feedback may have a steeper dependence on halo mass at earlier times, 2) the efficiency of star formation may be lower in low mass halos at high redshift, and 3) gas may not be able to accrete efficiently onto the disk in low mass halos at high redshift. These new recipes, once tuned, better reproduce the evolution of $f_\star \equiv M_\star/M_{\rm{H}}$ as a function of halo mass as derived from abundance matching over redshifts $z=0$ to $3$, though they have different effects on cold gas fractions, star formation rates, and metallicities. Changes to gas accretion and stellar-driven winds are promising, while direct modification of the star formation timescale requires drastic measures that are not physically well-motivated., Comment: 27 pages, 29 figures

We use a large suite of hydrodynamical simulations of binary galaxy mergers to construct and calibrate a physical prescription for computing the effective radii and velocity dispersions of spheroids. We implement this prescription within a semi-analytic model embedded in merger trees extracted from the Bolshoi Lambda-CDM N-body simulation, accounting for spheroid growth via major and minor mergers as well as disk instabilities. We find that without disk instabilities, our model does not predict sufficient numbers of intermediate mass early-type galaxies in the local universe. Spheroids also form earlier in models with spheroid growth via disk instabilities. Our model correctly predicts the normalization, slope, and scatter of the low-redshift size-mass and Fundamental Plane relations for early type galaxies. It predicts a degree of curvature in the Faber-Jackson relation that is not seen in local observations, but this could be alleviated if higher mass spheroids have more bottom-heavy initial mass functions. The model also correctly predicts the observed strong evolution of the size-mass relation for spheroids out to higher redshifts, as well as the slower evolution in the normalization of the Faber-Jackson relation. We emphasize that these are genuine predictions of the model since it was tuned to match hydrodynamical simulations and not these observations., Comment: Submitted to MNRAS

We investigate the star formation histories (SFHs) of high redshift (3 <~ z <~ 5) star-forming galaxies selected based on their rest-frame ultraviolet (UV) colors in the CANDELS/GOODS-S field. By comparing the results from the spectral-energy-distribution-fitting analysis with two different assumptions about the SFHs --- i.e., exponentially declining SFHs as well as increasing ones, we conclude that the SFHs of high-redshift star-forming galaxies increase with time rather than exponentially decline. We also examine the correlations between the star formation rates (SFRs) and the stellar masses. When the galaxies are fit with rising SFRs, we find that the trend seen in the data qualitatively matches the expectations from a semi-analytic model of galaxy formation. The mean specific SFR is shown to increase with redshift, also in agreement with the theoretical prediction. From the derived tight correlation between stellar masses and SFRs, we derive the mean SFH of star-forming galaxies in the redshift range of 3 <~ z <~ 5, which shows a steep power-law (with power alpha = 5.85) increase with time. We also investigate the formation timescales and the mean stellar population ages of these star-forming galaxies. Our analysis reveals that UV-selected star-forming galaxies have a broad range of the formation redshift. The derived stellar masses and the stellar population ages show positive correlation in a sense that more massive galaxies are on average older, but with significant scatter. This large scatter implies that the galaxies' mass is not the only factor which affects the growth or star formation of high-redshift galaxies., Comment: 31 pages, 8 figures, 2 tables

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth., Comment: Accepted by The Astrophysical Journal, 17 pages, 8 figures

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^111 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts., Comment: 23 pages, 16 figures, and 1 table. Accepted for publication in MNRAS

We study the evolution of atomic and molecular gas in galaxies in semi-analytic models of galaxy formation that include new modeling of the partitioning of cold gas in galactic discs into atomic, molecular, and ionised phases. We adopt two scenarios for the formation of molecules: one pressure-based and one metallicity-based. We find that both recipes successfully reproduce the gas fractions and gas-to-stellar mass ratios of HI and H2 in local galaxies,as well as the HI and H2 disc sizes up to z<2. We reach good agreement with the locally observed HI and H2 mass function, although both recipes slightly overpredict the low-mass end of the HI mass function. Both of our models predict that the high-mass end of the HI mass function remains nearly constant at redshifts z < 2.0. The metallicity-based recipe yields a higher cosmic density of cold gas and much lower cosmic H2 fraction over the entire redshift range probed than the pressure based recipe. These strong differences in HI mass function and cosmic density between the two recipes are driven by low mass galaxies (log (M*/Msun) < 7) residing in low mass halos (log (Mvir/Msun) < 10). Both recipes predict that galaxy gas fractions remain high from z ~ 6-3 and drop rapidly at lower redshift. The galaxy H2 fractions show a similar trend, but drop even more rapidly. We provide predictions for the CO J = 1-0 luminosity of galaxies, which will be directly comparable with observations with sub-mm and radio instruments., Comment: 24 pages, 17 figures, accepted for publication in MNRAS. For a high resolution version see http://www.astro.rug.nl/~gpopp/popping.samgas.pdf

We investigate the properties of damped Ly{\alpha} absorption systems (DLAs) in semi-analytic models of galaxy formation, including partitioning of cold gas in galactic discs into atomic, molecular, and ionized phases with a molecular gas-based star formation recipe. We investigate two approaches for partitioning gas into these constituents: a pressure-based and a metallicity-based recipe. We identify DLAs by passing lines of sight through our simulations to compute HI column densities. We find that models with "standard" gas radial profiles - where the average specific angular momentum of the gas disc is equal to that of the host dark matter halo - fail to reproduce the observed column density distribution of DLAs. These models also fail to reproduce the distribution of velocity widths {\Delta}v, overproducing low {\Delta}v relative to high {\Delta}v systems. Models with "extended" radial gas profiles - corresponding to gas discs with higher specific angular momentum - are able to reproduce quite well the column density distribution of absorbers over the column density range 19 < log NHI < 22.5 in the redshift range 2 < z < 3.5. The model with pressure-based gas partitioning also reproduces the observed line density of DLAs, HI gas density, and {\Delta}v distribution at z < 3 remarkably well. However all of the models investigated here underproduce DLAs and the HI gas density at z > 3. If this is the case, the flatness in the number of DLAs and HI gas density over the redshift interval 0 < z < 5 may be due to a cosmic coincidence where the majority of DLAs at z > 3 arise from intergalactic gas in filaments while those at z < 3 arise predominantly in galactic discs. We further investigate the dependence of DLA metallicity on redshift and {\Delta}v, and find reasonably good agreement with the observations, particularly when including the effects of metallicity gradients (abbrv.)., Comment: 27 pages, 15 figures, submitted to MNRAS

There is anecdotal evidence that spatially and physically unassociated galaxies blended into a single submillimetre (submm) source contribute to the submm galaxy (SMG) population. This work is the first to theoretically predict the number counts of such sources. We generate mock SMG catalogues using lightcones derived from the Bolshoi cosmological simulation; to assign submm flux densities to the mock galaxies, we use a fitting function previously derived from the results of dust radiative transfer performed on hydrodynamical simulations of isolated disc and merging galaxies. We then calculate submm number counts for different beam sizes and without blending. We predict that > ~50 per cent of blended SMGs have at least one spatially unassociated component with S_850 > 1 mJy. For a 15-arcsec beam, blends of >2 galaxies in which at least one component is spatially unassociated dominate the blended sources with total S_850 > ~3 mJy. The distribution of the redshift separations amongst the components is strongly bimodal. The typical redshift separation of spatially unassociated blended sources is ~1. Our predictions for the contributions of spatially unassociated components and the distribution of redshift separations are not testable with currently available data, but they will be easily tested once sufficiently accurate redshifts for the individual subcomponents (resolved by, e.g., ALMA) of a sufficient number of single-dish-detected blended SMGs are available., Comment: Accepted for publication in MNRAS. 11 pages, 4 figures, 2 tables

We present a new approach to study galaxy evolution in a cosmological context. We combine cosmological merger trees and semi-analytic models of galaxy formation to provide the initial conditions for multi-merger hydrodynamic simulations. In this way we exploit the advantages of merger simulations (high resolution and inclusion of the gas physics) and semi-analytic models (cosmological background and low computational cost), and integrate them to create a novel tool. This approach allows us to study the evolution of various galaxy properties, including the treatment of the hot gaseous halo from which gas cools and accretes onto the central disc, which has been neglected in many previous studies. This method shows several advantages over other methods. As only the particles in the regions of interest are included, the run time is much shorter than in traditional cosmological simulations, leading to greater computational efficiency. Using cosmological simulations, we show that multiple mergers are expected to be more common than sequences of isolated mergers, and therefore studies of galaxy mergers should take this into account. In this pilot study, we present our method and illustrate the results of simulating ten Milky Way-like galaxies since z=1. We find good agreement with observations for the total stellar masses, star formation rates, cold gas fractions and disc scale length parameters. We expect that this novel numerical approach will be very useful for pursuing a number of questions pertaining to the transformation of galaxy internal structure through cosmic time., Comment: 20 pages, 11 figures, 2 tables, submitted to MNRAS

We combine star-formation histories derived from observations of high redshift galaxies with measurements of the z~0 relation between gas-phase metallicity, stellar mass, and star formation rate to make an explicit and completely empirical connection between near-field and distant galaxy observations. Our approach relies on two basic assumptions: 1) galaxies' average paths through time in stellar mass vs. star formation rate space are represented by a family of smooth functions that are determined by the galaxies' final stellar mass, and 2) galaxies grow and become enriched with heavy elements such that they always evolve along the mass--metallicity--star formation rate relation. By integrating over these paths, we can track the chemical evolution of stars in galaxies in a model independent way, without the need for explicit assumptions about gas inflow, outflow, or star formation efficiency. Using this approach, we present predictions of stellar metallicity (i.e., O/H) distribution functions for present day star-forming galaxies of different stellar masses and the evolution of the alpha-element stellar metallicity-mass relation since z~1. The metallicity distribution functions are fairly well described as Gaussians, truncated at high metallicity, with power-law tails to low metallicity. We find that the stellar metallicity distribution for Milky Way mass galaxies is in reasonable agreement with observations for our Galaxy, and that the predicted stellar mass vs. mean stellar metallicity relation at z=0 agrees quite well with results derived from galaxy surveys. This validates the assumptions that are implicit in our simple approach. Upcoming observations will further test these assumptions and their range of validity, by measuring the mean stellar mass-metallicity relation up to z~1, and by measuring the stellar metallicity distributions over a range of galaxy masses., Comment: 8 pages, 2 figures, 2 tables. MNRAS, in press