The James Webb Space Telescope (JWST) is expected to observe galaxies at $z>10$ that are presently inaccessible. Here, we use a self-consistent empirical model, the UniverseMachine, to generate mock galaxy catalogues and lightcones over the redshift range $z=0-15$. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy-halo relationships, and galaxy-galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at $z<10$. We predict that Cycle 1 JWST surveys will very likely detect galaxies with $M_*>10^7 M_\odot$ and/or $M_{1500}<-17$ out to at least $z\sim 13.5$. Number density uncertainties at $z>12$ expand dramatically, so efforts to detect $z>12$ galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits ($M_\ast < 10^7M_\odot$ or $M_\mathrm{1500}>-17$). For reionization models, extrapolating luminosity functions with a constant faint-end slope from $M_{1500}=-17$ down to $M_{1500}=-12$ gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to $z\sim 12$. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and lightcones for common fields are available online.

Physical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation of the IllustrisTNG project to explore the physical parameter space comprised of 8 resolved ISM properties in star-forming regions across a wide range of host galaxy mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$). We deconstruct our simulated galaxies into regions of 1 kpc in size, and measure the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region, and study them in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate weighted distributions of these properties, we show that stars primarily form in two distinct regimes, which are brought about by an underlying bi-component radial star formation rate profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time, and compare our findings with observations from the MaNGA IFU survey. Further, using principal component analysis, we characterise the aforementioned parameter space to reveal a high-degree of multicollinearity in relationships between ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we find that a reduced three dimensional representation underpinned in essence by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

In the coming decade, a new generation of telescopes, including JWST and WFIRST, will probe the period of the formation of first galaxies and quasars, and open up the last frontier for structure formation. Recent simulations as well as observations have suggested that these galaxies are strongly clustered (with large scale bias $\gtrsim6$), and therefore have significant cosmic variance. In this work, we use \texttt{BlueTides}, the largest volume cosmological simulation of galaxy formation, to directly estimate the cosmic variance for current and upcoming surveys. Given its resolution and volume, \texttt{BlueTides} can probe the bias and cosmic variance of $z>7$ galaxies between magnitude $M_{UV}\sim-16$ to $M_{UV}\sim-22$ over survey areas $\sim0.1\ \mathrm{arcmin}^2$ to $\sim 10~\mathrm{deg}^2$. Within this regime, the cosmic variance decreases with survey area/ volume as a power law with exponents between $\sim-0.25$ to $\sim-0.45$. For the planned $10~\mathrm{deg}^2$ field of WFIRST, the cosmic variance is between $3\%$ to $10\%$. Upcoming JWST medium/ deep surveys with areas up to $A\sim100\ \mathrm{arcmin}^2$ will have cosmic variance ranging from $\sim 20-50\%$. Lensed surveys have the highest cosmic variance $\gtrsim 40\%$; the cosmic variance of $M_{UV}\lesssim-16$ galaxies is $\lesssim100\%$ up to $z\sim11$. At higher redshifts such as $z\sim12~(14)$, effective volumes of $\gtrsim(8~\mathrm{Mpc}/h)^3$ ($\gtrsim(12\ \mathrm{Mpc}/h)^3$) are required to limit the cosmic variance to within $100\%$. Finally, we find that cosmic variance is larger than Poisson variance and forms the dominant component of the overall uncertainty in all current and upcoming surveys. We present our calculations in the form of simple fitting functions and an online cosmic variance calculator (CV_AT_COSMIC_DAWN) which we publicly release.

Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level' galaxy formation physics to `top-level' cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.

Using the CANDELS photometric catalogs for the HST/ACS and WFC3, we identified massive evolved galaxies at $3 < z < 4.5$, employing three different selection methods. We find the comoving number density of these objects to be $\sim 2 \times 10^{-5}$ and $8 \times 10^{-6}Mpc^{-3}$ after correction for completeness for two redshift bins centered at $z=3.4, 4.7$. We quantify a measure of how much confidence we should have for each candidate galaxy from different selections and what are the conservative error estimates propagated into our selection. Then we compare the evolution of the corresponding number densities and their stellar mass density with numerical simulations, semi-analytical models, and previous observational estimates, which shows slight tension at higher redshifts as the models tend to underestimate the number and mass densities. By estimating the average halo masses of the candidates ($M_h \approx 4.2, 1.9, 1.3 \times 10^{12} M_\odot$ for redshift bins centered at $z=3.4, 4.1, 4.7$), we find them to be consistent with halos that were efficient in turning baryons to stars and were relatively immune to the feedback effects and on the verge of transition into hot-mode accretion. This can suggest the relative cosmological starvation of the cold gas followed by an overconsumption phase in which the galaxy consumes the available cold gas rapidly as one of the possible drivers for the quenching of the massive evolved population at high redshift.

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

We present the first predictions for the $L_{\rm [CII]}$ - SFR relation and [CII] luminosity function (LF) in the Epoch of Reionization (EoR) based on cosmological hydrodynamics simulations using the SIMBA suite plus radiative transfer calculations via S\'IGAME. The sample consists of 11,137 galaxies covering halo mass $\log M_{\rm halo}\in$[9, 12.4] $M_\odot$, star formation rate SFR$\in$[0.01, 330] $M_\odot$ yr$^{-1}$, and metallicity $<Z_{\rm gas}>_{\rm SFR}\in$[0.1, 0.9] $Z_\odot$. The simulated $L_{\rm [CII]}$-SFR relation is consistent with the range observed, but with a spread of $\simeq$0.3 dex at the high end of SFR ($>$100 $M_\odot$ yr$^{-1}$) and $\simeq$0.6 dex at the lower end, and there is tension between our predictions and the values of $L_{\rm [CII]}$ above 10$^{8.5}$ $L_\odot$ observed in some galaxies reported in the literature. The scatter in the $L_{\rm [CII]}$-SFR relation is mostly driven by galaxy properties, such that at a given SFR, galaxies with higher molecular gas mass and metallicity have higher $L_{\rm [CII]}$. The [CII] LF predicted by SIMBA is consistent with the upper limits placed by the only existing untargeted flux-limited [CII] survey at the EoR (ASPECS) and those predicted by semi-analytic models. We compare our results with existing models and discuss differences responsible for the discrepant slopes in the $L_{\rm [CII]}$-SFR relatiion.

The James Webb Space Telescope (JWST) is expected to enable transformational progress in studying galaxy populations in the very early Universe, during the Epoch of Reionization (EoR). A critical parameter for understanding the sources that reionized the Universe is the Lyman-continuum production efficiency, $\xi_{\rm ion}$, defined as the rate of production of ionizing photons divided by the intrinsic UV luminosity. In this work, we combine self-consistent star formation and chemical enrichment histories predicted by semi-analytic models of galaxy formation with stellar population synthesis (SPS) models to predict the expected dependence of $\xi_{\rm ion}$ on galaxy properties and cosmic epoch from $z=4-10$. We then explore the sensitivity of the production rate of ionizing photons, $\dot{N}_{\rm ion}$, to the choice of SPS model and the treatment of stellar feedback in our galaxy formation model. We compare our results to those of other simulations, constraints from empirical models, and observations. We find that adopting SPS models that include binary stars predict about a factor of two more ionizing radiation than models that only assume single stellar populations. We find that UV-faint, low-mass galaxies have values of $\xi_{\rm ion}$ about 0.25 dex higher than those of more massive galaxies, but find weak evolution with cosmic time, about 0.2 dex from $z\sim 12-4$ at fixed rest-UV luminosity. We provide predictions of $\dot{N}_{\rm ion}$ as a function of $M_{\rm h}$ and a number of other galaxy properties. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $\mu$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $\mu$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $\mu$m number counts is almost exhausted.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

Cosmic voids, the under-dense regions of the cosmic web, are widely used to constrain cosmology. Voids contain few, isolated galaxies, presumably expected to be less evolved and preserving memory of the pristine Universe. We use the cosmological hydrodynamical simulation Horizon-AGN coupled to the void finder {\sc \texttt{VIDE}} to investigate properties of galaxies in voids at z=0. We find that, closer to void centers, low-mass galaxies are more common than their massive counterparts. At fixed dark matter halo mass, they have smaller stellar masses than in denser regions. The star formation rate of void galaxies diminishes when approaching void centers, but their sSFR slightly increases, suggesting that void galaxies form stars more efficiently with respect to their stellar mass. We find that this can not only be attributed to the prevalence of low-mass galaxies. The inner region of voids also predominantly host low-mass BHs. However, the BH mass to galaxy mass ratios resemble those of the whole simulation at z=0. Our results suggest that even if the growth channels in cosmic voids are different than in denser environments, voids grow their galaxies and BHs in a similar way. While a large fraction of the BHs have low Eddington ratios, we find that 20\% could be observed as AGN with log10 L=41.5-42.5 erg/s in hard X-ray (2-10 keV). These results pave the way to future work with larger next-generation hydro simulations, aiming to confirm our findings and prepare the application on data from upcoming large surveys such as PFS, Euclid and WFIRST.

We explore the galaxy-galaxy merger rate with the empirical model for galaxy formation, EMERGE. On average, we find that between $2$ per cent and $6$ per cent of massive galaxies ($\log_{10}(m_{*}/M_{\odot}) \geq 10.3$) will experience a major merger per Gyr. Our model predicts galaxy merger rates that do not scale as a power-law with redshift, and exhibit a clear stellar mass and mass-ratio dependence, unlike the underlying halo-halo merger rate. Specifically, major mergers are more frequent at high masses and at low redshift. We show mergers are significant for the stellar mass growth of galaxies $\log_{10}(m_{*}/M_{\odot}) \gtrsim 11.0$. For the most massive galaxies major mergers dominate the accreted mass fraction, contributing as much as $90$ per cent of the total accreted stellar mass. We argue that these phenomena are a direct result of the stellar-to-halo mass relation, which results in massive galaxies having a higher likelihood of experiencing major mergers than low mass galaxies. Furthermore, we find no statistical evidence for mergers as a driving mechanism for quenching, with more than $58$ per cent of the most massive galaxies experiencing no major merger within $2$ Gyr of quenching. Our model produces a galaxy pair fraction consistent with resent observations, exhibiting a form best described by a power-law exponential function. Translating these pair fractions into merger rates results in an over prediction compared to the model intrinsic values. We find the pair fraction can be best mapped to the intrinsic merger rate by adopting a constant observation timescale $T_{obs}=2.56$ Gyr, assuming all observed pairs merge by $z=0$.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole (BH) growth is sustained by an enhanced number of mergers and by being at the intersection of filaments bringing fresh gas. This assumption should be supported by an enhancement of the number counts of galaxies in the field of view of quasars. While the current observations of quasar environments do not lead to a consensus on a possible excess of galaxies, the future missions JWST, WFIRST, and Euclid will provide new insights on quasar environments, and will substantially increase the number of study-cases. We are in a crucial period, where we need to both understand the current observations and predict how upcoming missions will improve our understanding of BH environments. Using the large-scale simulation Horizon-AGN, we find that statistically the most massive BHs reside in environments with the largest number counts of galaxies. However, we find a large variance in galaxy number counts, and some massive BHs do not show enhanced counts in their neighborhood. Interestingly, some massive BHs have a very close galaxy companion but no further enhancement of the galaxy number counts at larger scales, in agreement with recent observations. We find that AGN feedback in the surrounding galaxies is able to decrease their luminosity and stellar mass, and therefore to make them un-observable when using restrictive galaxy selection criteria. Radiation from the quasars can spread over large distances, which could affect the formation history of surrounding galaxies, but a careful analysis of these processes requires radiative transfer simulations.

We investigate the effects of massive black hole growth on the structural evolution of dwarf galaxies within the Romulus25 cosmological hydrodynamical simulation. We study a sample of 228 central, isolated dwarf galaxies with stellar masses $M_{star} < 10^{10} M_\odot$ and a central BH. We find that the local $M_{BH} - M_{star}$ relation exhibits a high degree of scatter below $M_{star} < 10^{10} M_\odot$, which we use to classify BHs as overmassive or undermassive relative to their host $M_{star}$. Overmassive BHs grow through a mixture of BH mergers and relatively high average accretion rates, while undermassive BHs grow slowly through accretion. We find that isolated dwarf galaxies that host overmassive BHs also follow different evolutionary tracks relative to their undermassive BH counterparts, building up their stars and dark matter earlier and experiencing star formation suppression starting around $z=2$. By $z=0.05$, overmassive BH hosts above $M_{star} > 10^{9} M_\odot$ are more likely to exhibit lower central stellar mass density, lower HI gas content, and lower star formation rates than their undermassive BH counterparts. Our results suggest that overmassive BHs in isolated galaxies above $M_{star} > 10^{9} M_\odot$ are capable of driving feedback, in many cases suppressing and even quenching star formation by late times.

(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall.

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

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-\sigma$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-\sigma$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-\sigma$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $\sigma>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

The long-anticipated James Webb Space Telescope (JWST) will be able to directly detect large samples of galaxies at very high redshift. Using the well-established, computationally efficient Santa Cruz semi-analytic model, with recently implemented multiphase gas partitioning and H2-based star formation recipes, we make predictions for a wide variety of galaxy properties for galaxy populations at $z$ = 4-10. In this work, we provide forecasts for the physical properties of high-redshift galaxies and links to their photometric properties. With physical parameters calibrated only to $z\sim0$ observations, our model predictions are in good agreement with current observational constraints on stellar mass and star formation rate distribution functions up to $z \sim 8$. We also provide predictions representing wide, deep, and lensed JWST survey configurations. We study the redshift evolution of key galaxy properties and the scaling relations among them. Taking advantage of our models' high computational efficiency, we study the impact of systematically varying the model parameters. All distribution functions and scaling relations presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Str\"omgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas.

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.

We study the destruction of interstellar dust via sputtering in supernova (SN) shocks using three-dimensional hydrodynamical simulations. With a novel numerical framework, we follow both sputtering and dust dynamics governed by direct collisions, plasma drag and betatron acceleration. Grain-grain collisions are not included and the grain-size distribution is assumed to be fixed. The amount of dust destroyed per SN is quantified for a broad range of ambient densities and fitting formulae are provided. Integrated over the grain-size distribution, nonthermal (inertial) sputtering dominates over thermal sputtering for typical ambient densities. We present the first simulations that explicitly follow dust sputtering within a turbulent multiphase interstellar medium. We find that the dust destruction timescales $\tau$ are 0.35 Gyr for silicate dust and 0.44 Gyr for carbon dust in solar neighborhood conditions. The SN environment has an important impact on $\tau$. SNe that occur in preexisting bubbles destroy less dust as the destruction is limited by the amount of dust in the shocked gas. This makes $\tau$ about 2.5 times longer than the estimate based on results from a single SN explosion. We investigate the evolution of the dust-to-gas mass ratio (DGR), and find that a spatial inhomogeneity of $\sim$ 14\% develops for scales below 10 pc. It locally correlates positively with gas density but negatively with gas temperature even in the exterior of the bubbles due to incomplete gas mixing. This leads to a $\sim$ 30\% lower DGR in the volume filling warm gas compared to that in the dense clouds.

One of the key questions in understanding the formation and evolution of galaxies is how starbursts affect the assembly of stellar populations in galaxies over time. We define a burst indicator ($\eta$), which compares a galaxy's star formation rates on short ($\sim10$ Myr) and long ($\sim100$ Myr) timescales. To estimate $\eta$, we apply the detailed time-luminosity relationship for H$\alpha$ and near-ultraviolet emission to simulated star formation histories (SFHs) from semi-analytic models and the Mufasa hydrodynamical cosmological simulations. The average of $\eta$ is not a good indicator of star formation stochasticity (burstiness); indeed, we show that this average should be close to zero unless the population has an average SFH which is rising or falling rapidly. Instead, the width of the $\eta$ distribution characterizes the burstiness of a galaxy population's recent star formation. We find this width to be robust to variations in stellar initial mass function and metallicity. We apply realistic noise and selection effects to the models to generate mock HST and JWST galaxy catalogs and compare these catalogs with 3D-HST observations of 956 galaxies at $0.65<z<1.5$ detected in H$\alpha$. Measurements of $\eta$ are unaffected by dust measurement errors under the assumption that $E(B-V)_\mathrm{stars}=0.44\,E(B-V)_\mathrm{gas}$ (i.e., $Q_\mathrm{sg}=0.44$). However, setting $Q_\mathrm{sg}=0.8^{+0.1}_{-0.2}$ removes an unexpected dependence of the average value of $\eta$ upon dust attenuation and stellar mass in the 3D-HST sample while also resolving disagreements in the distribution of star formation rates. However, even varying the dust law cannot resolve all discrepancies between the simulated and the observed galaxies.

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $\beta=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) provides new constraints for galaxy formation models on the molecular gas properties of galaxies. We compare results from ASPECS to predictions from two cosmological galaxy formation models: the IllustrisTNG hydrodynamical simulations and the Santa Cruz semi-analytic model (SC SAM). We explore several recipes to model the H$_2$ content of galaxies, finding them to be consistent with one another, and take into account the sensitivity limits and survey area of ASPECS. For a canonical CO-to-H$_2$ conversion factor of $\alpha_{\rm CO} = 3.6\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$ the results of our work include: (1) the H$_2$ mass of $z>1$ galaxies predicted by the models as a function of their stellar mass is a factor of 2-3 lower than observed; (2) the models do not reproduce the number of H$_2$-rich ($M_{\rm H2} > 3\times 10^{10}\,\rm{M}_\odot$) galaxies observed by ASPECS; (3) the H$_2$ cosmic density evolution predicted by IllustrisTNG (the SC SAM) is in tension (only just agrees) with the observed cosmic density, even after accounting for the ASPECS selection function and field-to-field variance effects. The tension between models and observations at $z>1$ can be alleviated by adopting a CO-to-H$_2$ conversion factor in the range $\alpha_{\rm CO} = 2.0 - 0.8\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$. Additional work on constraining the CO-to-H$_2$ conversion factor and CO excitation conditions of galaxies through observations and theory will be necessary to more robustly test the success of galaxy formation models.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

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

Achieving a precise understanding of galaxy formation in a cosmological context is one of the great challenges in theoretical astrophysics, due to the vast range of spatial scales involved in the relevant physical processes. Observations in the millimeter bands, particularly those using the cosmic microwave background (CMB) radiation as a "backlight", provide a unique probe of the thermodynamics of these processes, with the capability to directly measure the density, pressure, and temperature of ionized gas. Moreover, these observations have uniquely high sensitivity into the outskirts of the halos of galaxies and clusters, including systems at high redshift. In the next decade, the combination of large spectroscopic and photometric optical galaxy surveys and wide-field, low-noise CMB surveys will transform our understanding of galaxy formation via these probes.

The star formation histories (SFHs) of galaxies contain imprints of the physical processes responsible for regulating star formation during galaxy growth and quenching. We improve the Dense Basis SFH reconstruction method of Iyer & Gawiser (2017), introducing a nonparametric description of the SFH based on the lookback times at which a galaxy assembles certain quantiles of its stellar mass. The method uses Gaussian Processes to create smooth SFHs that are independent of any functional form, with a flexible number of parameters that is adjusted to extract the maximum possible amount of SFH information from the SEDs being fit. We apply the method to reconstruct the SFHs of 48,791 galaxies with $H<25$ at $0.5 < z < 3.0$ across the five CANDELS fields. Using these SFHs, we study the evolution of galaxies as they grow more massive over cosmic time. We quantify the fraction of galaxies that show multiple major episodes of star formation, finding that the median time between two peaks of star formation is $\sim 0.42_{-0.10}^{+0.15}t_{univ}$ Gyr, where $t_{univ}$ is the age of the universe at a given redshift and remains roughly constant with stellar mass. Correlating SFHs with morphology, we find that studying the median SFHs of galaxies at $0.5<z<1.0$ at the same mass ($10^{10}< M_* < 10^{10.5}M_\odot$) allows us to compare the timescales on which the SFHs decline for different morphological classifications, ranging from $0.60^{-0.54}_{+1.54}$ Gyr for galaxies with spiral arms to $2.50^{-1.50}_{+2.25}$ Gyr for spheroids. The Gaussian Process-based SFH description provides a general approach to reconstruct smooth, nonparametric SFH posteriors for galaxies with a flexible number of parameters that can be incorporated into Bayesian SED fitting codes to minimize the bias in estimating physical parameters due to SFH parametrization.

We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good agreement with observations, but we note an excess of faint AGN in hard X-ray (L_x ~ 10^{43-44} erg/s), and a lower number of bright AGN (L_x>10^{44} erg/s), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the 10^{9} Msun BHs compared to observations suggest that AGN feedback may be too efficient in this regime. We study galaxy star formation activity and structural properties, and design sample-dependent criteria to identify different galaxy types (star-forming/quiescent, extended/compact) that we apply both to the simulations and observations from the candels fields. We analyze how the simulated and observed galaxies populate the specific star formation rate - stellar mass surface density diagram. A large fraction of the z=0 M_{star}>10^{11} Msun quiescent galaxies first experienced a compaction phase (i.e. reduction of galaxy size) while still forming stars, and then a quenching event. We measure the dependence of AGN fraction on galaxies' locations in this diagram. After correcting the simulations with a redshift and AGN luminosity-dependent model for AGN obscuration, we find good qualitative and quantitative agreement with observations. The AGN fraction is the highest among compact star-forming galaxies (16-20% at z~1.5-2), and the lowest among compact quiescent galaxies (6-10% at z~1.5-2).

We present the $ugriz$-band Dark Energy Camera (DECam) plus 3.6 and 4.5 $\mu$m IRAC catalogs for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim24$ deg$^{2}$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 region, with seven bandpasses spanning a wavelength range of 0.35 to 4.5 $\mu$m. SHELA falls within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which will provide spectroscopic redshifts for $\sim200{,}000$ Ly$\alpha$ emitters at $1.9<z<3.5$ and also for $\sim200{,}000$ [OII] emitters at $z<0.5$. SHELA's deep, wide-area multiwavelength images combined with HETDEX's spectroscopic information, will facilitate many extragalactic studies, including measuring the evolution of galaxy stellar mass, halo mass, and environment from $1.5<z<3.5$. Here we present $riz$-band selected $ugriz$-band DECam catalogs that reach a $5\sigma$ depth of $\sim24.5$ AB mag (for point sources with an aperture that encloses $70\%$ of the total flux) and cover $17.5$ deg$^{2}$ of the overall SHELA field. We validate our DECam catalog by comparison to the DECam Legacy Survey (DECaLS) DR5 and the Dark Energy Survey (DES) DR1. We perform IRAC forced photometry with The Tractor image modeling code to measure 3.6 and 4.5 $\mu$m fluxes for all objects within our DECam catalog. We demonstrate the utility of our catalog by computing galaxy number counts and estimating photometric redshifts. Our photometric redshifts recover the available $\left\langle z \right\rangle = 0.33 $ SDSS spectroscopic redshifts with a $1\sigma$ scatter in $\Delta z/(1 +z)$ of 0.04.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

In anticipation of the upcoming deployment of the James Webb Space Telescope (JWST), we present high-redshift predictions by the well-established Santa Cruz semi-analytic model. We update the models by re-calibrating them after adopting cosmological parameters consistent with recent constraints from Planck. We provide predictions for rest-frame UV luminosity functions for galaxy populations over a wide range of $M_\text{UV}$ from $\sim-6$ to $\sim-24$ between $z = 4 - 10$. In addition, we present the corresponding predictions for observed-frame galaxy number counts in different redshift bins in the full set of NIRCam filters. We provide predictions of the quantitative effect on these observables of varying the physical recipes implemented in the models, such as the molecular gas depletion time (star formation efficiency) scalings or the scalings of outflow rates driven by stars and supernovae with galaxy circular velocity. Based on these results, we discuss what may be learned about the physical processes that shape galaxy formation from JWST observations of galaxy number densities at different intrinsic luminosities. All data tables for the results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $\xi_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/{\sigma}$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $\xi_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content.

We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which optical- and UV-line diagnostic diagrams can best help to discern between the main ionizing sources, as traced by the ratio of BH accretion to star formation rates in model galaxies, over a wide range of redshifts. At low redshift, simulated star-forming galaxies, galaxies dominated by active galactic nuclei and composite galaxies are appropriately differentiated by standard selection criteria in the classical [OIII]$\lambda$5007/H$\beta$ versus [NII]$\lambda$6584/H$\alpha$ diagram. At redshifts $z \gt 1$, however, this optical diagram fails to discriminate between active and inactive galaxies at metallicities below $0.5\ Z_\odot$. To robustly classify the ionizing radiation of such metal-poor galaxies, which dominate in the early Universe, we confirm 3 previous, and propose 11 novel diagnostic diagrams based on equivalent widths and luminosity ratios of UV emission lines, such as EW(OIII]$\lambda$1663) versus OIII]$\lambda$1663/HeII$\lambda$1640, CIII]$\lambda$1908/HeII$\lambda$1640 versus OIII]$\lambda$1663/HeII$\lambda$1640, and CIV$\lambda$1550/CIII]$\lambda$1908 versus CIII]$\lambda$1908/CII$\lambda$2326. We formulate associated UV selection criteria and discuss some caveats of our results (e.g., uncertainties in the modelling of the HeII$\lambda$1640 line). These UV diagnostic diagrams are potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.

The advent of new sub-millimeter observational facilities has stimulated the desire to model the sub-mm line emission of galaxies within cosmological galaxy formation models. This is typically done in post-processing by applying sub-resolution recipes to describe the properties of the unresolved interstellar medium. At the same time, while there is freedom in how one implements these sorts of recipes, the impact of various choices has yet to be systematically explored in simulations. In this paper, we do just that. We combine a semi-analytic model of galaxy formation with chemical equilibrium networks and numerical radiative transfer models and explore how different choices for the sub-grid modeling affect the predicted CO, [CI], and [CII] emission of galaxies. We find that a key component for a successful model includes a molecular cloud mass-size relation and scaling for the ultraviolet and cosmic ray radiation field that varies based on local ISM properties. Our most successful model adopts a Plummer profile for the radial density distribution of gas within molecular clouds. On the other hand, different assumptions for the clumping of gas within molecular clouds and changes in the slope of the molecular cloud mass distribution function hardly affect the CO, [CI], and [CII] luminosities of galaxies. At fixed star-formation rate the [CII]-SFR ratio of galaxies scales inversely with the pressure acting on molecular clouds, increasing the molecular clouds density and hence decreasing the importance of [CII] line cooling. Overall we find that it is essential that a wide range of sub-mm emission lines arising in vastly different phases of the ISM are used as model constraints in order to limit the freedom in sub-grid choices and we present the details of the most successful model variant.

The next generation Very Large Array (ngVLA) will revolutionize our understanding of the distant Universe via the detection of cold molecular gas in the first galaxies. Its impact on studies of galaxy characterization via detailed gas dynamics will provide crucial insight on dominant physical drivers for star-formation in high redshift galaxies, including the exchange of gas from scales of the circumgalactic medium down to resolved clouds on mass scales of $\sim10^{5}\,M_\odot$. In this study, we employ a series of high-resolution, cosmological, hydrodynamic zoom simulations from the MUFASA simulation suite and a CASA simulator to generate mock ngVLA observations of a $z\sim4.5$ gas rich star-forming galaxy. Using the DESPOTIC radiative transfer code that encompasses simultaneous thermal, chemical, and statistical equilibrium in calculating the molecular and atomic level transitions of CO from ALMA for comparison. We find that observations of CO(1-0) are especially important for tracing the systemic redshift of the galaxy and the total mass of the well-shielded molecular gas reservoir, while even CO(2-1) can predominantly trace denser gas regions distinct from CO(1-0). The factor of 100 times improvement in mapping speed for the ngVLA beyond the Jansky VLA and the proposed ALMA Band 1 will make these detailed, high-resolution imaging and kinematic studies of CO(1-0) routine at $z\sim2-5$.

For the first time in human history, technologies have matured sufficiently to enable a mission capable of discovering and characterizing habitable planets like Earth orbiting sunlike stars other than the Sun. At the same time, such a platform would enable unique science not possible from ground-based facilities. This science is broad and exciting, ranging from new investigations of our own solar system to a full range of astrophysics disciplines. The Habitable Exoplanet Observatory, or HabEx, is one of four studies currently being undertaken by NASA in preparation for the 2020 Astrophysics Decadal Survey. HabEx has been designed to be the Great Observatory of the 2030s, with community involvement through a competed and funded Guest Observer (GO) program. This interim report describes the HabEx baseline concept, which is a space-based 4-meter diameter telescope mission concept with ultraviolet (UV), optical, and near-infrared (near-IR) imaging and spectroscopy capabilities. More information on HabEx can be found at https://www.jpl.nasa.gov/habex

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ correlation. The reconstructed SFHs represent trajectories in SFR-M$_*$ space, enabling us to study galaxies at epochs earlier than observed by propagating them backwards in time along these trajectories. We study the SFR-M$_*$ correlation at $z=1,2,3,4,5,6$ using both direct fits to galaxies observed at those epochs and SFR-M$_*$ trajectories of galaxies observed at lower redshifts. The SFR-M$_*$ correlations obtained using the two approaches are found to be consistent with each other through a KS test. Validation tests using SFHs from semi-analytic models and cosmological hydrodynamical simulations confirm the sensitivity of the method to changes in the slope, normalization and shape of the SFR-M$_*$ correlation. This technique allows us to further probe the low-mass regime of the correlation at high-z by $\sim 1$ dex and over an effective volume of $\sim 10\times$ larger than possible with just direct fits. We find that the SFR-M$_*$ correlation is consistent with being linear down to M$_*\sim 10^7 M_\odot$ at $z>4$. The evolution of the correlation is well described by $\log SFR= (0.80\pm 0.029 - 0.017\pm 0.010\times t_{univ})\log M_*$ $- (6.487\pm 0.282-0.039\pm 0.008\times t_{univ})$, where $t_{univ}$ is the age of the universe in Gyr.

We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nucleus (AGN) feedback and (2) with AGN feedback arising from winds and X-ray radiation. We find that AGN feedback can alter the stellar density distribution, reduce the core density within the central 1 kpc by 0.3 dex from z=1, and enhance the size growth of massive galaxies. We also find that galaxies simulated with AGN feedback evolve along similar tracks to those characterized by observations in specific star formation versus compactness. We confirm that AGN feedback plays an important role in transforming galaxies from blue compact galaxies into red extended galaxies in two ways: (1) it effectively quenches the star formation, transforming blue compact galaxies into compact quiescent galaxies and (2) it also removes and prevents new accretion of cold gas, shutting down in-situ star formation and causing subsequent mergers to be gas-poor or mixed. Gas poor minor mergers then build up an extended stellar envelope. AGN feedback also puffs up the central region through the fast AGN driven winds as well as the slow expulsion of gas while the black hole is quiescent. Without AGN feedback, large amounts of gas accumulate in the central region, triggering star formation and leading to overly massive blue galaxies with dense stellar cores.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

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<z<0.6$. Our AGN fraction is in good agreement with a semi-analytic model, suggesting that as we search larger volumes we may use comparisons between observed AGN fractions and models to understand seeding mechanisms in the early universe.

We present a joint analysis of the rest-frame ultraviolet (UV) luminosity functions of continuum-selected star-forming galaxies and galaxies dominated by active galactic nuclei (AGNs) at $z \sim$ 4. These 3,740 $z \sim$ 4 galaxies are selected from broad-band imaging in nine photometric bands over 18 deg$^2$ in the \textit{Spitzer}/HETDEX Exploratory Large Area Survey (SHELA) field. The large area and moderate depth of our survey provide a unique view of the intersection between the bright end of the galaxy UV luminosity function (M$_{AB}<-$22) and the faint end of the AGN UV luminosity function. We do not separate AGN-dominated galaxies from star-formation-dominated galaxies, but rather fit both luminosity functions simultaneously. These functions are best fit with a double power-law (DPL) for both the galaxy and AGN components, where the galaxy bright-end slope has a power-law index of $-3.80\pm0.10$, and the corresponding AGN faint-end slope is $\alpha_{AGN} = -1.49^{+0.30}_{-0.21}$. We cannot rule out a Schechter-like exponential decline for the galaxy UV luminosity function, and in this scenario the AGN luminosity function has a steeper faint-end slope of $-2.08^{+0.18}_{-0.11}$. Comparison of our galaxy luminosity function results with a representative cosmological model of galaxy formation suggests that the molecular gas depletion time must be shorter, implying that star formation is more efficient in bright galaxies at $z=4$ than at the present day. If the galaxy luminosity function does indeed have a power-law shape at the bright end, the implied ionizing emissivity from AGNs is not inconsistent with previous observations. However, if the underlying galaxy distribution is Schechter, it implies a significantly higher ionizing emissivity from AGNs at this epoch.

We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN' model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective' feedback - the outflows themselves which drive gas out of galaxies - and `preventative' feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies.

The similarity between the distributions of spins for galaxies ($\lambda_{\rm g}$) and for dark-matter haloes ($\lambda_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeq\lambda_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass radius of the galaxy and $R_{\rm v}$ the halo radius. Utilizing two different suites of zoom-in cosmological simulations, we find that $\lambda_{\rm g}$ and $\lambda_{\rm h}$ are in fact only barely correlated, especially at $z\geq 1$. A general smearing of this correlation is expected based on the different spin histories, where the more recently accreted baryons through streams gain and then lose significant angular momentum compared to the gradually accumulated dark matter. Expecting the spins of baryons and dark matter to be correlated at accretion into $R_{\rm v}$, the null correlation at the end reflects an anti-correlation between $\lambda_{\rm g}/\lambda_{\rm h}$ and $\lambda_{\rm h}$, which can partly arise from mergers and a compact star-forming phase that many galaxies undergo. On the other hand, the halo and galaxy spin vectors tend to be aligned, with a median $\cos\theta=0.6$-0.7 between galaxy and halo, consistent with instreaming within a preferred plane. The galaxy spin is better correlated with the spin of the inner halo, but this largely reflects the effect of the baryons on the halo. Following the null spin correlation, $\lambda_{\rm h}$ is not a useful proxy for $R_{\rm e}$. While our simulations reproduce a general relation of the sort $R_{\rm e}=AR_{\rm vir}$, in agreement with observational estimates, the relation becomes tighter with $A=0.02(c/10)^{-0.7}$, where $c$ is the halo concentration, which in turn introduces a dependence on mass and redshift.

We estimate the cosmic number density of the recently identified class of HI-bearing ultra-diffuse sources (HUDs) based on the completeness limits of the ALFALFA survey. These objects fall in the range $8.5 < \log M_{\rm{HI}}/\rm{M_{\odot}} < 9.5$, have average $r$-band surface brightnesses fainter than 24 mag arsec$^{-2}$, half-light radii greater than 1.5 kpc, and are separated from neighbours by at least 350 kpc. We find that HUDs contribute at most 6% of the population of HI-bearing dwarfs, have a total cosmic number density of $(1.5 \pm 0.6) \times 10^{-3}$ $\rm{Mpc^{-3}}$, and an HI mass density of $(6.0 \pm 0.8) \times 10^{5}$ $\rm{M_{\odot}\,Mpc^{-3}}$. We estimate that this is similar to the total cosmic number density of UDGs in groups and clusters, and conclude that the relation between the number of UDGs hosted in a halo and the halo mass, must have a break below $M_{200} \sim 10^{12}$ $\rm{M_{\odot}}$ in order to account for the abundance of HUDs. The distribution of the velocity widths of HUDs rises steeply towards low values, indicating a preference for slow rotation rates. These objects have been absent from measurements of the galaxy stellar mass function owing to their low surface brightness. However, we estimate that due to their low number density, their inclusion would constitute a correction of less than 1%. Comparison with the Santa Cruz SAM shows that it produces HI-rich central UDGs that have similar colours to HUDs, but these are currently produced in much great a number. While previous results from this sample have favoured formation scenarios where HUDs form in high spin parameter halos, comparisons with the results of Rong et al. 2017, which invokes that formation mechanism, reveal that this model produces an order of magnitude more field UDGs than we observe in the HUD population.(Abridged)

The \emph{Spitzer} Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) provides unparalleled depth at $3.6$ and $4.5$~$\mu$m over $\sim0.66$~deg$^2$ of the COSMOS field, allowing precise photometric determinations of redshift and stellar mass. From this unique dataset we can connect galaxy samples, selected by stellar mass, to their host dark matter halos for $1.5<z<5.0$, filling in a large hitherto unexplored region of the parameter space. To interpret the observed galaxy clustering we utilize a phenomenological halo model, combined with a novel method to account for uncertainties arising from the use of photometric redshifts. We find that the satellite fraction decreases with increasing redshift and that the clustering amplitude (e.g., comoving correlation length / large-scale bias) displays monotonic trends with redshift and stellar mass. Applying $\Lambda$CDM halo mass accretion histories and cumulative abundance arguments for the evolution of stellar mass content we propose pathways for the coevolution of dark matter and stellar mass assembly. Additionally, we are able to estimate that the halo mass at which the ratio of stellar to halo mass is maximized is $10^{12.5_{-0.08}^{+0.10}}$~M$_{\odot}$ at $z\sim2.5$. This peak halo mass is here inferred for the first time from stellar mass-selected clustering measurements at $z\gtrsim2$, and implies mild evolution of this quantity for $z\lesssim3$, consistent with constraints from abundance-matching techniques.

We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various star formation history for each galaxy. Like others, we find that the slope of the Main Sequence (MS) of star formation in the $(M_*;SFR)$ plane bends at high mass. We observe clear morphological differences among galaxies across the MS, which also correlate with stellar age. At all redshifts, galaxies that are quenching or quenched, and thus old, have high $\Sigma_1$ (the projected density within the central 1 kpc), while younger, star-forming galaxies span a much broader range of $\Sigma_1$, which includes the high values observed for quenched galaxies, but also extends to much lower values. As galaxies age and quench, the stellar age and the dispersion of $\Sigma_1$ for fixed values of $M_{*}$ shows two different regimes, one, at the low--mass end, where quenching might be driven by causes external to the galaxies; the other, at the high--mass end, where quenching is driven by internal causes, very likely the mass given the low scatter of $\Sigma_1$ (mass quenching). We suggest that the monotonic increase of central density as galaxies grow is one manifestation of a more general phenomenon of structural transformation that galaxies undergo as they evolve.

The rate of major galaxy-galaxy merging is theoretically predicted to steadily increase with redshift during the peak epoch of massive galaxy development ($1{\leq}z{\leq}3$). We use close-pair statistics to objectively study the incidence of massive galaxies (stellar $M_{1}{\geq}2{\times}10^{10}M_{\odot}$) hosting major companions ($1{\leq}M_{1}/M_{2}{\leq}4$; i.e., $<$4:1) at six epochs spanning $0{<}z{<}3$. We select companions from a nearly complete, mass-limited ($\geq5{\times}10^{9}M_{\odot}$) sample of 23,696 galaxies in the five CANDELS fields and the SDSS. Using $5-50$ kpc projected separation and close redshift proximity criteria, we find that the major companion fraction $f_{\mathrm{mc}}(z)$ based on stellar mass-ratio (MR) selection increases from 6% ($z{\sim}0$) to 16% ($z{\sim}0.8$), then turns over at $z{\sim}1$ and decreases to 7% ($z{\sim}3$). Instead, if we use a major F160W flux ratio (FR) selection, we find that $f_{\mathrm{mc}}(z)$ increases steadily until $z=3$ owing to increasing contamination from minor (MR$>$4:1) companions at $z>1$. We show that these evolutionary trends are statistically robust to changes in companion proximity. We find disagreements between published results are resolved when selection criteria are closely matched. If we compute merger rates using constant fraction-to-rate conversion factors ($C_{\mathrm{merg,pair}}{=}0.6$ and $T_{\mathrm{obs,pair}}{=}0.65\mathrm{Gyr}$), we find that MR rates disagree with theoretical predictions at $z{>}1.5$. Instead, if we use an evolving $T_{\mathrm{obs,pair}}(z){\propto}(1+z)^{-2}$ from Snyder et al., our MR-based rates agree with theory at $0{<}z{<}3$. Our analysis underscores the need for detailed calibration of $C_{\mathrm{merg,pair}}$ and $T_{\mathrm{obs,pair}}$ as a function of redshift, mass and companion selection criteria to better constrain the empirical major merger history.

The initial mass function (IMF) is a key ingredient in many studies of galaxy formation and evolution. Although the IMF is often assumed to be universal, there is continuing evidence that it is not universal. Spectroscopic studies that derive the IMF of the unresolved stellar populations of a galaxy often assume that this spectrum can be described by a single stellar population (SSP). To alleviate these limitations, in this paper we have developed a unique hierarchical Bayesian framework for modelling composite stellar populations (CSPs). Within this framework we use a parameterized IMF prior to regulate a direct inference of the IMF. We use this new framework to determine the number of SSPs that is required to fit a set of realistic CSP mock spectra. The CSP mock spectra that we use are based on semi-analytic models and have an IMF that varies as a function of stellar velocity dispersion of the galaxy. Our results suggest that using a single SSP biases the determination of the IMF slope to a higher value than the true slope, although the trend with stellar velocity dispersion is overall recovered. If we include more SSPs in the fit, the Bayesian evidence increases significantly and the inferred IMF slopes of our mock spectra converge, within the errors, to their true values. Most of the bias is already removed by using two SSPs instead of one. We show that we can reconstruct the variable IMF of our mock spectra for signal-to-noise ratios exceeding $\sim$75.

We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

Galaxies occupy different regions of the [OIII]$\lambda5007$/H$\beta$-versus-[NII]$\lambda6584$/H$\alpha$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$\beta$, [NII]$\lambda6584$/H$\alpha$, [SII]$\lambda\lambda6717,6731$/H$\alpha$ and [OI]$\lambda6300$/H$\alpha$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$\beta$ increases and [NII]/H$\alpha$, [SII]/H$\alpha$ and [OI]/H$\alpha$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$\beta$ and [NII]/H$\alpha$. For [SII]/H$\alpha$ and [OI]/H$\alpha$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.

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$.

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.

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.

At $z=1-3$, the formation of new stars is dominated by dusty galaxies whose far-IR emission indicates they contain colder dust than local galaxies of a similar luminosity. We explore the reasons for the evolving IR emission of similar galaxies over cosmic time using: 1) Local galaxies from GOALS $(L_{\rm IR}=10^{11}-10^{12}\,L_\odot)$; 2) Galaxies at $z\sim0.1-0.5$ from the 5MUSES ($L_{\rm IR}=10^{10}-10^{12}\,L_\odot$); 3) IR luminous galaxies spanning $z=0.5-3$ from GOODS and Spitzer xFLS ($L_{\rm IR}>10^{11}\,L_\odot$). All samples have Spitzer mid-IR spectra, and Herschel and ground-based submillimeter imaging covering the full IR spectral energy distribution, allowing us to robustly measure $L_{\rm IR}^{\rm\scriptscriptstyle SF}$, $T_{\rm dust}$, and $M_{\rm dust}$ for every galaxy. Despite similar infrared luminosities, $z>0.5$ dusty star forming galaxies have a factor of 5 higher dust masses and 5K colder temperatures. The increase in dust mass is linked with an increase in the gas fractions with redshift, and we do not observe a similar increase in stellar mass or star formation efficiency. $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$, a proxy for $T_{\rm dust}$, is strongly correlated with $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ independently of redshift. We measure merger classification and galaxy size for a subsample, and there is no obvious correlation between these parameters and $L_{\rm IR}^{\rm \scriptscriptstyle SF}/M_{\rm dust}$ or $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$. In dusty star forming galaxies, the change in $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ can fully account for the observed colder dust temperatures, suggesting that any change in the spatial extent of the interstellar medium is a second order effect.

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.

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.

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.

Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the CANDELS dataset. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected because the Universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$ Gyr). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms that operate on different timescales.

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.

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.

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.

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.

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.

We study the evolution of the core (r<1 kpc) and effective (r<r_e) stellar-mass surface densities, in star-forming and quiescent galaxies. Since z=3, both populations occupy distinct, linear relations in log(Sigma_e) and log(Sigma_1) vs. log(M). These structural relations exhibit slopes and scatter that remain almost constant with time while their normalizations decline. For SFGs, the normalization declines by less than a factor of 2 from z=3, in both Sigma_e and Sigma_1. Such mild declines suggest that SFGs build dense cores by growing along these relations. We define this evolution as the structural main sequence (Sigma-MS). Quiescent galaxies follow different relations (Sigma^Q_e, Sigma^Q_1) off the Sigma-MS by having higher densities than SFGs of the same mass and redshift. The normalization of Sigma^Q_e declines by a factor of 10 since z=3, but only a factor of 2 in Sigma^Q_1. Thus, the common denominator for quiescent galaxies at all redshifts is the presence of a dense stellar core, and the formation of such cores in SFGs is the main requirement for quenching. Expressed in 2D as deviations off the SFR-MS and off Sigma^Q_1 at each redshift, the distribution of massive galaxies forms a universal, L-shaped sequence that relates two fundamental physical processes: compaction and quenching. Compaction is a process of substantial core-growth in SFGs relative to that in the Sigma-MS. This process increases the core-to-total mass and Sersic index, thereby, making compact SFGs. Quenching occurs once compact SFGs reach a maximum central density above Sigma^Q_1 > 9.5 M_sun/kpc^2. This threshold provides the most effective selection criterion to identify the star-forming progenitors of quiescent galaxies at all redshifts.

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.

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.

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.

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.

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

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.

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.

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.

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.

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)

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.

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.

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.

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

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.

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.

We combine a semi-analytic model of galaxy formation, tracking atomic and molecular phases of cold gas, with a three-dimensional radiative-transfer and line tracing code to study the sub-mm emission from atomic and molecular species (CO, HCN, [CI], [CII], [OI]) in galaxies. We compare the physics that drives the formation of stars at the epoch of peak star formation (SF) in the Universe (z = 2.0) with that in local galaxies. We find that normal star-forming galaxies at high redshift have much higher CO-excitation peaks than their local counterparts and that CO cooling takes place at higher excitation levels. CO line ratios increase with redshift as a function of galaxy star-formation rate, but are well correlated with H2 surface density independent of redshift. We find an increase in the [OI]/[CII] line ratio in typical star-forming galaxies at z = 1.2 and z = 2.0 with respect to counterparts at z = 0. Our model results suggest that typical star-forming galaxies at high redshift consist of much denser and warmer star-forming clouds than their local counterparts. Galaxies belonging to the tail of the SF activity peak at z = 1.2 are already less dense and cooler than counterparts during the actual peak of SF activity (z = 2.0). We use our results to discuss how future ALMA surveys can best confront our predictions and constrain models of galaxy formation.

Recent observations have probed the formation histories of nearby elliptical galaxies by tracking correlations between the stellar population parameters, age and metallicity, and the structural parameters that enter the Fundamental Plane, radius and velocity dispersion sigma. These studies have found intriguing correlations between these four parameters. In this work, we make use of a semi-analytic model, based on halo merger trees extracted from the Bolshoi cosmological simulation, that predicts the structural properties of spheroid-dominated galaxies based on an analytic model that has been tested and calibrated against an extensive suite of hydrodynamic+N-body binary merger simulations. We predict the radius, sigma, luminosity, age, and metallicity of spheroid-dominated galaxies, enabling us to compare directly to observations. Our model predicts a strong correlation between age and sigma for early-type galaxies, and no significant correlation between age and radius, in agreement with observations. In addition we predict a strong correlation between metallicity and sigma, and a weak correlation between metallicity and radius, in qualitative agreement with observations. We find that the correlations with sigma arise as a result of the strong link between sigma and the galaxy assembly time. Minor mergers produce a large change in radius while leaving sigma nearly the same, which explains the weaker trends with radius.

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.

We present Keck-I MOSFIRE near-infrared spectroscopy for a sample of 13 compact star-forming galaxies (SFGs) at redshift $2\leq z \leq2.5$ with star formation rates of SFR$\sim$100M$_{\odot}$ y$^{-1}$ and masses of log(M/M$_{\odot}$)$\sim10.8$. Their high integrated gas velocity dispersions of $\sigma_{\rm{int}}$=230$^{+40}_{-30}$ km s$^{-1}$, as measured from emission lines of H$_{\alpha}$ and [OIII], and the resultant M$_{\star}-\sigma_{\rm{int}}$ relation and M$_{\star}$$-$M$_{\rm{dyn}}$ all match well to those of compact quiescent galaxies at $z\sim2$, as measured from stellar absorption lines. Since log(M$_{\star}$/M$_{\rm{dyn}}$)$=-0.06\pm0.2$ dex, these compact SFGs appear to be dynamically relaxed and more evolved, i.e., more depleted in gas and dark matter ($<$13$^{+17}_{-13}$\%) than their non-compact SFG counterparts at the same epoch. Without infusion of external gas, depletion timescales are short, less than $\sim$300 Myr. This discovery adds another link to our new dynamical chain of evidence that compact SFGs at $z\gtrsim2$ are already losing gas to become the immediate progenitors of compact quiescent galaxies by $z\sim2$.

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.

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.).

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.

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

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.

The processes that trigger Active Galactic Nuclei (AGN) remain poorly understood. While lower luminosity AGN may be triggered by minor disturbances to the host galaxy, stronger disturbances are likely required to trigger luminous AGN. Major wet mergers of galaxies are ideal environments for AGN triggering since they provide large gas supplies and galaxy scale torques. There is however little observational evidence for a strong connection between AGN and major mergers. We analyse the morphological properties of AGN host galaxies as a function of AGN and host galaxy luminosity and compare them to a carefully matched sample of control galaxies. AGN are X-ray selected in the redshift range 0.5 < z < 0.8 and have luminosities 41 < log(L_X [erg/s]) < 44.5. 'Fake AGN' are simulated in the control galaxies by adding point sources with the magnitude of the matched AGN. We find that AGN host and control galaxies have comparable assymetries, Sersic indices and ellipticities at restframe ~950nm. AGN host galaxies show neither higher average asymmetries nor higher fractions of very disturbed objects. There is no increase in the prevalence of merger signatures with AGN luminosity. At 95% confidence we find that major mergers are responsible for <6% of all AGN in our sample as well as <40% of the highest luminosity AGN log(L_X [erg/s]) ~ 43.5). Major mergers therefore either play only a very minor role in the triggering of AGN in the luminosity range studied or time delays are too long for merger features to remain visible.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

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<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.

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^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 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.

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

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.

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

We use multi-band spatially resolved photometry from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) in the 4 MSec Chandra Deep Field-South (CDF-S) to explore the nuclear and extended colors, color gradients and stellar populations of X-ray selected AGN host galaxies out to z=3. Based on a study of their central light, we develop X-ray based criteria to exclude objects with strong AGN contamination. We use stellar masses from the FIREWORKS database to understand and account for stellar mass selection effects, and carefully study, for the first time, the resolved host galaxy properties of AGNs at z~2 in their rest-frame optical light without substantial nuclear contamination. AGN hosts span a sizable range of stellar masses, colors and color gradients at these redshifts. Their colors, color gradients and stellar population properties are very similar to inactive galaxies of the same stellar mass. At z~1, we find a slightly narrower range in host colors compared to inactive galaxies, as well as hints of more recent star-formation. These differences are weaker or non-existent among AGN hosts at z~2. We discuss the importance of AGN driven feedback in the quenching of galaxies at z>~1 and speculate on possible evolution in the relationship between black hole accretion and the host galaxy towards high redshifts.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

We discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

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.

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.

Observational studies have revealed a "downsizing" trend in black hole (BH) growth: the number densities of luminous AGN peak at higher redshifts than those of faint AGN. This would seem to imply that massive black holes formed before low mass black holes, in apparent contradiction to hierarchical clustering scenarios. We investigate whether this observed "downsizing" in BH growth is reproduced in a semi-analytic model for the formation and evolution of galaxies and black holes, set within the hierarchical paradigm for structure formation (Somerville et al. 2008; S08). In this model, black holes evolve from light seeds (\sim100M\odot) and their growth is merger-driven. The original S08 model (baseline model) reproduces the number density of AGN at intermediate redshifts and luminosities, but underproduces luminous AGN at very high redshift (z > 3) and overproduces them at low redshift (z < 1). In addition, the baseline model underproduces low-luminosity AGN at low redshift (z < 1). To solve these problems we consider several modifications to the physical processes in the model: (1) a 'heavy' black hole seeding scenario (2) a sub-Eddington accretion rate ceiling that depends on the cold gas fraction, and (3) an additional black hole accretion mode due to disk instabilities. With these three modifications, the models can explain the observed downsizing, successfully reproduce the bolometric AGN luminosity function and simultaneously reproduce galaxy and black hole properties in the local Universe. We also perform a comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for torus-level obscuration, and reach similar conclusions. Our best-fit model suggests a scenario in which disk instabilities are the main driver for moderately luminous Seyfert galaxies at low redshift, while major mergers are the main trigger for luminous AGN.

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

We explore the effects of mergers on the evolution of massive early-type galaxies by modeling the evolution of their stellar populations in a hierarchical context. We investigate how a realistic red sequence population set up by z~1 evolves under different assumptions for the merger and star formation histories, comparing changes in color, luminosity and mass. The purely passive fading of existing red sequence galaxies, with no further mergers or star formation, results in dramatic changes at the bright end of the luminosity function and color-magnitude relation. Without mergers there is too much evolution in luminosity at a fixed space density compared to observations. The change in color and magnitude at a fixed mass resemble that of a passively evolving population that formed relatively recently, at z~2. Mergers amongst the red sequence population ("dry mergers") occurring after z=1 build up mass, counteracting the fading of the existing stellar populations to give smaller changes in both color and luminosity for massive galaxies. By allowing some galaxies to migrate from the blue cloud onto the red sequence after z=1 through gas-rich mergers, younger stellar populations are added to the red sequence. This manifestation of the progenitor bias increases the scatter in age and results in even smaller changes in color and luminosity between z=1 and z=0 at a fixed mass. The resultant evolution appears much slower, resembling the passive evolution of a population that formed at high redshift (z~3-5) and is in closer agreement with observations. Measurements of the luminosity and color evolution alone are not sufficient to distinguish between the purely passive evolution of an old population and cosmologically motivated hierarchical growth, although these scenarios have very different implications for the mass growth of early-type galaxies over the last half of cosmic history.

Attenuation of high-energy gamma rays by pair-production with UV, optical and IR extragalactic background light (EBL) photons provides a link between the history of galaxy formation and high-energy astrophysics. We present results from our latest semi-analytic models (SAMs), which employ the main ingredients thought to be important to galaxy formation and evolution, as well as an improved model for reprocessing of starlight by dust to mid- and far-IR wavelengths. These SAMs are based upon a Lambda-CDM hierarchical structural formation scenario, and are successful in reproducing a large variety of observational constraints such as number counts, luminosity and mass functions, and color bimodality. Our fiducial model is based upon a WMAP5 cosmology, and treats dust emission using empirical templates. This model predicts a background flux considerably lower than optical and near-IR measurements that rely on subtraction of zodiacal and galactic foregrounds, and near the lower bounds set by number counts of resolvable sources at a large number of wavelengths. We also show the results of varying cosmological parameters and dust attenuation model used in our SAM. For each EBL prediction, we show how the optical depth due to electron-positron pair-production is affected by redshift and gamma-ray energy, and the effect of gamma-ray absorption on the spectra of a variety of extragalactic sources. We conclude with a discussion of the implications of our work, comparisons to other models and key measurements of the EBL and a discussion of how the burgeoning science of gamma-ray astronomy will continue to help constrain cosmology. The low EBL flux predicted by our fiducial model suggests an optimistic future for further studies of distant gamma-ray sources.

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

We combine a semi-analytic model of galaxy formation with simple analytic recipes describing the absorption and re-emission of starlight by dust in the interstellar medium of galaxies. We use the resulting models to predict galaxy counts and luminosity functions from the far-ultraviolet to the sub-mm, from redshift five to the present, and compare with an extensive compilation of observations. We find that in order to reproduce the rest-UV and optical luminosity functions at high redshift, we must assume an evolving normalization in the dust-to-metal ratio, implying that galaxies of a given bolometric luminosity (or metal column density) must be less extinguished than their local counterparts. In our best-fit model, we find remarkably good agreement with observations from rest-frame 1500 Angstroms to 250 microns. At longer wavelengths, most dramatically in the sub-mm, our models underpredict the number of bright galaxies by a large factor. The models reproduce the observed total IR luminosity function fairly well. We show the results of varying several ingredients of the models, including various aspects of the dust attenuation recipe, the dust emission templates, and the cosmology. We use our models to predict the integrated Extragalactic Background Light (EBL), and compare with an observationally-motivated EBL model and with other available observational constraints.

Recent low-redshift observations have attempted to determine the star formation histories of elliptical galaxies by tracking correlations between the stellar population parameters (age and metallicity) and the structural parameters that enter the fundamental plane (size and velocity dispersion). These studies have found that velocity dispersion, rather than effective radius or dynamical mass, is the main predictor of a galaxy's stellar age and metallicity. In this work, we apply an analytic model that predicts the structural properties of remnants formed in major mergers to progenitor disk galaxies with properties taken from two different semi-analytic models. We predict the effective radius, velocity dispersion, luminosity, age, and metallicity of the merger remnants, enabling us to compare directly to observations of early-type galaxies. While we find a tight correlation between age and velocity dispersion, we find a stronger dependence of age and metallicity on effective radius than observations report. The correlations arise as a result of the dependence of gas fraction, age, and metallicity on the stellar mass in the progenitor disk galaxies. These dependences induce a rotation in the radius-velocity plane between the correlations with effective radius and circular velocity in the disk galaxy progenitors, and the correlations with effective radius and velocity dispersion in the elliptical galaxy remnants. The differences between our results and those from observations suggest that major mergers alone cannot produce the observed lack of correlation between effective radius and stellar population parameters. Simulations have suggested that subsequent minor mergers introduce scatter in the effective radius while leaving the velocity dispersion essentially unchanged. Incorporating such minor mergers into the model may, then, bring the simulations into closer agreement with observations.

One key piece of information missing from high redshift galaxy surveys is the galaxies' cold gas contents. We present a new method to indirectly determine cold gas surface densities and integrated gas masses from galaxy star formation rates and to separate the atomic and molecular gas components. Our predicted molecular and total gas surface densities and integrated masses are in very good agreement with direct measurements quoted in the literature for low and high-z galaxies. We apply this method to predict the gas content for a sample of $\sim 57000$ galaxies in the COSMOS field at $0.5 \leq z \leq 2.0$, selected to have $I_{AB} < 24$ mag. This approach allows us to investigate in detail the redshift evolution of galaxy cold and molecular gas content versus stellar mass and to provide fitting formulae for galaxy gas fractions. We find a clear trend between galaxy gas fraction, molecular gas fraction and stellar mass with redshift, suggesting that massive galaxies consume and/or expel their gas at higher redshift than less massive objects and have lower fractions of their gas in molecular form. The characteristic stellar mass separating gas- from stellar-dominated galaxies decreases with time. This indicates that massive galaxies reach a gas-poor state earlier than less massive objects. These trends can be considered to be another manifestation of downsizing in star formation activity.

[Abridged] We make use of a semi-analytic cosmological model that includes simple prescriptions for dust attenuation and emission to make predictions for the observable and physical properties of galaxies that may be detected by the recently launched Herschel Space Observatory in deep fields such as GOODS-Herschel. We compare our predictions for differential galaxy number counts in the PACS (100 & 160) and SPIRE (250, 350, and 500 micron) bands with available observations. We find very good agreement with the counts in the PACS bands, for the overall counts and for galaxies binned by redshift at z< 2. At z > 2 our model underpredicts the number of bright galaxies by a factor of ten. The agreement is much worse for all three SPIRE bands, and becomes progressively worse with increasing wavelength. We discuss a number of possible reasons for these discrepancies, and hypothesize that the effect of blending on the observational flux estimates is likely to be the dominant issue. We note that the PACS number counts are relatively robust to changes in the dust emission templates, while the predicted SPIRE number counts are more template dependent. We present quantitative predictions for the relationship between the observed PACS 160 and SPIRE 250 micron fluxes and physical quantities such as halo mass, stellar mass, cold gas mass, star formation rate, and total infrared (IR) luminosity, at different redshifts. We also present quantitative predictions for the correlation between PACS 160 micron flux and the probability that a galaxy has experienced a recent major or minor merger. Although our models predict a strong correlation between these quantities, such that more IR-luminous galaxies are more likely to be merger-driven, we find that more than half of all high redshift IR-luminous galaxies detected by Herschel are able to attain their high star formation rates without enhancement by a merger.

We used Early Release Science (ERS) observations taken with the Wide Field Camera 3 (WFC3) in the GOODS-S field to study the galaxy stellar mass function (GSMF) at 0.6<=z<4.5. Deep WFC3 near-IR data (for Y as faint as 27.3, J and H as faint as 27.4 AB mag at 5 sigma), as well as deep Ks (as faint as 25.5 at 5 sigma) Hawk-I band data, provide an exquisite data set with which determine in an unprecedented way the low-mass end of the GSMF, allowing an accurate probe of masses as low as M~7.6 10^9 Msun at z~3. Although the area used is relatively small (~33 arcmin^2), we found generally good agreement with previous studies on the entire mass range. Our results show that the slope of the faint-end increases with redshift, from alpha=-1.44+/-0.03 at z~0.8 to alpha=-1.86+/-0.16 at z~3, although indications exist that it does not steepen further between z~3 and z~4. This result is insensitive to any uncertainty in the M* parameter. The steepness of the GSMF faint-end solves the well-known disagreement between the stellar mass density (SMD) and the integrated star formation history at z>2. However, we confirm the that there appears to be an excess of integrated star formation with respect to the SMD at z<2, by a factor of ~2-3. Our comparison of the observations with theoretical predictions shows that the models forecast a greater abundance of low mass galaxies, at least up to z~3, as well as a dearth of massive galaxies at z~4 with respect to the data, and that the predicted SMD is generally overestimated at z<~2.

[abridged] We study the evolution of the Star Formation Rate Function (SFRF) of massive galaxies over the 0.4<z<1.8 redshift range and its implications for our understanding of the physical processes responsible for galaxy evolution. We use multiwavelength observations included in the GOODS-MUSIC catalogue, which provides a suitable coverage of the spectral region from 0.3 to 24 micron and either spectroscopic or photometric redshifts for each object. Individual SFRs have been obtained by combining UV and 24 micron observations, when the latter were available. For all other sources an "SED fitting" SFR estimate has been considered. We then define a stellar mass limited sample, complete in the Mstar>1.e10 Msun range and determine the SFRF using the 1/Vmax algorithm. We define simulated galaxy catalogues based on three different semi-analytical models of galaxy formation and evolution. We show that the theoretical SFRFs are well described by a double power law functional form and its redshift evolution is approximated with high accuracy by a pure evolution of the typical SFR. We find good agreement between model predictions and the high-SFR end of the SFRF, when the observational errors on the SFR are taken into account. However, the observational SFRF is characterised by a double peaked structure, which is absent in its theoretical counterparts. At z>1.0 the observed SFRF shows a relevant density evolution, which is not reproduced by SAMs, due to the well known overprediction of intermediate mass galaxies at z~2. The agreement at the low-SFR end is poor: all models overpredict the space density of SFR~1 Msun/yr and no model reproduces the double peaked shape of the observational SFRF. If confirmed by deeper IR observations, this discrepancy will provide a key constraint on theoretical modelling of star formation and stellar feedback.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present a detailed comparison between numerical cosmological hydrodynamic zoom simulations and semi-analytic models (SAMs) run within merger trees extracted from the simulations. The high-resolution simulations represent 48 individual halos with virial masses in the range 2.4*10^11M_sun < M_Halo < 3.3*10^13M_sun. They include radiative H & He cooling, photo-ionization, star formation and thermal SN feedback. We compare with different SAM versions including only this complement of physical processes, and also ones including supernova driven winds, metal cooling, and feedback from AGN. Our analysis is focused on the cosmic evolution of the baryon content in galaxies and its division into various components (stars, cold gas, and hot gas). Both the SAMs and simulations are compared with observational relations between halo mass and stellar mass, and between stellar mass and star formation rate, at low and high redshift. The simulations turn out to have much higher star formation efficiencies (by about a factor of ten) than the SAMs. Therefore the cold gas is consumed much more rapidly in the simulations and stars form much earlier. Also, simulations show a transition between stellar mass growth that is dominated by in situ formation of stars to growth that is predominantly through accretion of stars formed in external galaxies. In SAMs, stellar growth is always dominated by in situ star formation. In addition, SAMs overpredict the overall gas accretion rates relative to the simulations, and overestimate the fraction of "hot" relative to "cold" accretion. We discuss the reasons for these discrepancies, and identify several physical processes that are missing in our SAM. We also highlight physical processes that are neglected in the simulations studied here, but which appear to be crucial in order to understand the properties of real galaxies.

Using HST/WFC3 imaging taken as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), we examine the role that major galaxy mergers play in triggering active galactic nuclei (AGN) activity at z~2. Our sample consists of 72 moderate-luminosity (Lx ~ 1E42-1E44 erg/s) AGN at 1.5<z<2.5 that are selected using the 4 Msec Chandra observations in the Chandra Deep Field South, the deepest X-ray observations to date. Employing visual classifications, we have analyzed the rest-frame optical morphologies of the AGN host galaxies and compared them to a mass-matched control sample of 216 non-active galaxies at the same redshift. We find that most of the AGN reside in disk galaxies (51.4%), while a smaller percentage are found in spheroids (27.8%). Roughly 16.7% of the AGN hosts have highly disturbed morphologies and appear to be involved in a major merger or interaction, while most of the hosts (55.6%) appear relatively relaxed and undisturbed. These fractions are statistically consistent with the fraction of control galaxies that show similar morphological disturbances. These results suggest that the hosts of moderate-luminosity AGN are no more likely to be involved in an ongoing merger or interaction relative to non-active galaxies of similar mass at z~2. The high disk fraction observed among the AGN hosts also appears to be at odds with predictions that merger-driven accretion should be the dominant AGN fueling mode at z~2, even at moderate X-ray luminosities. Although we cannot rule out that minor mergers are responsible for triggering these systems, the presence of a large population of relatively undisturbed disk-like hosts suggests that secular processes play a greater role in fueling AGN activity at z~2 than previously thought.

Calculating the galaxy merger rate requires both a census of galaxies identified as merger candidates, and a cosmologically-averaged `observability' timescale T_obs(z) for identifying galaxy mergers. While many have counted galaxy mergers using a variety of techniques, T_obs(z) for these techniques have been poorly constrained. We address this problem by calibrating three merger rate estimators with a suite of hydrodynamic merger simulations and three galaxy formation models. We estimate T_obs(z) for (1) close galaxy pairs with a range of projected separations, (2) the morphology indicator G-M20, and (3) the morphology indicator asymmetry A. Then we apply these timescales to the observed merger fractions at z < 1.5 from the recent literature. When our physically-motivated timescales are adopted, the observed galaxy merger rates become largely consistent. The remaining differences between the galaxy merger rates are explained by the differences in the range of mass-ratio measured by different techniques and differing parent galaxy selection. The major merger rate per unit co-moving volume for samples selected with constant number density evolves much more strongly with redshift (~ (1+z)^(+3.0 \pm 1.1)) than samples selected with constant stellar mass or passively evolving luminosity (~ (1+z)^(+0.1 \pm 0.4)). We calculate the minor merger rate (1:4 < M_{sat}/M_{primary} <~ 1:10) by subtracting the major merger rate from close pairs from the `total' merger rate determined by G-M20. The implied minor merger rate is ~3 times the major merger rate at z ~ 0.7, and shows little evolution with redshift.

We employ hydrodynamical simulations to study the effects of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. For the first time we present a suite of simulations that includes a diffuse, rotating, cooling, hot gaseous halo, as predicted by cosmological hydrodynamical simulations as well as models of galaxy formation. We study the effect of this new gaseous component on the vertical structure of a Milky Way-like stellar disc during 1:10 and 1:5 mergers. For 1:10 mergers we find no increased final thin disc scale height compared to the isolated simulation, leading to the conclusion that thin discs can be present even after a 1:10 merger if a reasonable amount of hot gas is present. The reason for this is the accretion of new cold gas, leading to the formation of a massive new thin stellar disc that dominates the surface brightness profile. In a previous study, in which we included only cold gas in the disk, we showed that the presence of cold gas decreased the thickening by a minor merger relative to the no-gas case. Here, we show that the evolution of the scale height in the presences of a cooling hot halo is dominated by the formation of the new stellar disc. In this scenario, the thick disc is the old stellar disc that has been thickened in a minor merger at z>1, while the thin disc is the new stellar disc that reforms after this merger. In addition, we study the evolution of the scale height during a 1:5 merger and find that a thin disc can be present even after this merger, provided enough hot gas is available. The final scale height in our simulations depends on the mass of the hot gaseous halo, the efficiency of the winds and the merger mass ratio. We find post-merger values in the range 0.5<z0<1.0 kpc in good agreement with observational constraints by local galaxies.

Data from (non-) attenuation of gamma rays from active galactic nuclei (AGN) and gamma ray bursts (GRBs) give upper limits on the extragalactic background light (EBL) from the UV to the mid-IR that are only a little above the lower limits from observed galaxies. These upper limits now rule out some EBL models and purported observations, with improved data likely to provide even stronger constraints. We present EBL calculations both based on multiwavelength observations of thousands of galaxies and also based on semi-analytic models, and show that they are consistent with these lower limits from observed galaxies and with the gamma-ray upper limit constraints. Such comparisons "close the loop" on cosmological galaxy formation models, since they account for all the light, including that from galaxies too faint to see. We compare our results with those of other recent works, and discuss the implications of these new EBL calculations for gamma ray attenuation. Catching a few GRBs with groundbased atmospheric Cherenkov Telescope (ACT) arrays or water Cherenkov detectors could provide important new constraints on the high-redshift star formation history of the universe.

There are strong correlations between the three structural properties of elliptical galaxies -- stellar mass, velocity dispersion and size -- in the form of a tight "fundamental plane" and a "scaling relation" between each pair. Major mergers of disk galaxies are assumed to be a mechanism for producing ellipticals, but semi-analytic galaxy formation models (SAM) have encountered apparent difficulties in reproducing the observed slope and scatter of the size-mass relation. We study the scaling relations of merger remnants using progenitor properties from two SAMs. We apply a simple merger model that includes gas dissipation and star formation based on theoretical considerations and simulations. Combining the SAMs and the merger model allows calculation of the structural properties of the remnants of major mergers that enter the population of elliptical galaxies at a given redshift. Without tuning the merger model parameters for each SAM, the results roughly match the slope and scatter in the observed scaling relations and their evolution in the redshift range $z=0-3$. Within this model, the observed scaling relations, including the tilt of the fundamental plane relative to the virial plane, result primarily from the decrease of gas fraction with increasing progenitor mass. The scatter in the size-mass relation of the remnants is reduced from that of the progenitors because of a correlation between progenitor size and gas fraction at a given mass.

Interstellar dust grains are responsible for modifying the spectral energy distribution (SED) of galaxies, both absorbing starlight at UV and optical wavelengths and converting this energy into thermal emission in the infrared. The detailed description of these phenomena is of fundamental importance in order to compare the predictions of theoretical models of galaxy formation and evolution with the most recent observations in the infrared region. In this paper we compare the results of GRASIL, a code explicitly solving for the equation of radiative transfer in a dusty medium, with the predictions of a variety of IR template libraries. We employ star formation history samples extracted from the semi-analytical galaxy formation model MORGANA to create libraries of synthetic SEDs from the near- to the far-infrared. We consider model predictions at different redshift ranges to explore any possible influence in the shape and normalization of the SEDs due to the expected evolution of the galaxy properties. We compute the total absorbed starlight predicted by GRASIL at optical wavelengths to statistically compare the synthetic SEDs with the selected IR templates. We show that synthetic SEDs at a given total infrared luminosity are predicted to be systematically different at different redshift and for different properties of the underlying model galaxy. However, we determine spectral regions where the agreement between the results of radiative transfer and IR templates is good in a statistical sense (i.e. in terms of the luminosity functions). Moreover, we highlight some potentially relevant discrepancies between the different approaches, both in the region dominated by PAH emission and at sub-mm wavelengths. These results determine potentially critical issues in the infrared luminosity functions as predicted by semi-analytical models coupled with different IR flux estimators.

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy. However, the overall spectrum of the EBL between 0.1 and 1000 microns has never been determined directly neither from observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncertainties in our modelling based directly on the data are quantified, and their consequences for attenuation of very-high-energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level roughly matching galaxy count data. Independent efforts from IR and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

AGN feedback is believed to play an important role in shaping a variety of observed galaxy properties, as well as the evolution of their stellar masses and star formation rates. In particular, in the current theoretical paradigm of galaxy formation, AGN feedback is believed to play a crucial role in regulating the levels of activity in galaxies, in relatively massive halos at low redshift. Only in recent years, however, has detailed statistical information on the dependence of galaxy activity on stellar mass, parent halo mass and hierarchy has become available. In this paper, we compare the fractions of galaxies belonging to different activity classes (star-forming, AGN and radio active) with predictions from four different and independently developed semi-analytical models. We adopt empirical relations to convert physical properties into observables (H_alpha emission lines, OIII line strength and radio power). We demonstrate that all models used in this study reproduce the overall distributions of galaxies belonging to different activity classes as a function of stellar mass and halo mass: star forming galaxies and the strongest radio sources are preferentially associated with low-mass and high-mass galaxies/halos respectively. However, model predictions differ from observational measurements in a number of ways. All models used in our study predict that almost every >1.e12 Msun dark matter halo and/or >1.e11 Msun galaxy should host a bright radio source, while only a small fraction of galaxies belong to this class in the data. In addition, radio brightness is expected to depend strongly on the mass of the parent halo mass in the models, while strong and weak radio galaxies are found in similar environments in data. Our results highlight that the distribution of AGN as a function of stellar mass provides one of the most promising discriminants between different gas accretion schemes.

We explore methods to improve the estimates of star formation rates and mean stellar population ages from broadband photometry of high redshift star-forming galaxies. We use synthetic spectral templates with a variety of simple parametric star formation histories to fit broadband spectral energy distributions. These parametric models are used to infer ages, star formation rates and stellar masses for a mock data set drawn from a hierarchical semi-analytic model of galaxy evolution. Traditional parametric models generally assume an exponentially declining rate of star-formation after an initial instantaneous rise. Our results show that star formation histories with a much more gradual rise in the star formation rate are likely to be better templates, and are likely to give better overall estimates of the age distribution and star formation rate distribution of Lyman break galaxies. For B- and V-dropouts, we find the best simple parametric model to be one where the star formation rate increases linearly with time. The exponentially-declining model overpredicts the age by 100 % and 120 % for B- and V-dropouts, on average, while for a linearly-increasing model, the age is overpredicted by 9 % and 16 %, respectively. Similarly, the exponential model underpredicts star-formation rates by 56 % and 60 %, while the linearly-increasing model underpredicts by 15 % 22 %, respectively. For U-dropouts, the models where the star-formation rate has a peak (near z ~ 3) provide the best match for age -- overprediction is reduced from 110 % to 26 % -- and star-formation rate -- underprediction is reduced from 58 % to 22 %. We classify different types of star-formation histories in the semi-analytic models and show how the biases behave for the different classes. We also provide two-band calibration formulae for stellar mass and star formation rate estimations.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy, but the overall spectrum of the EBL between 0.1-1000 microns has never been determined directly from galaxy spectral energy distribution (SED) observations over a wide redshift range. The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. This is achieved from the observed evolution of the rest-frame K-band galaxy luminosity function up to redshift 4 (Cirasuolo et al. 2010), combined with a determination of galaxy SED-type fractions. These are based on fitting SWIRE templates to a multiwavelength sample of about 6000 galaxies in the redshift range from 0.2 to 1 from the All-wavelength Extended Groth Strip International Survey (AEGIS). The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and AGN galaxies in that redshift range are estimated, and two alternative extrapolations of SED-types to higher redshifts are considered. This allows calculation of the evolution of the luminosity densities from the UV to the IR, the evolving star formation rate density of the universe, the evolving contribution to the bolometric EBL from the different galaxy populations including AGN galaxies and the buildup of the EBL. Our EBL calculations are compared with those from a semi-analytic model, from another observationally-based model and observational data. The EBL uncertainties in our modeling based directly on the data are quantified, and their consequences for attenuation of very high energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL is well constrained from the UV to the mid-IR, but independent efforts from infrared and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

We explore the rate and impact of galaxy mergers on the massive galaxy population using the amplitude of the two-point correlation function on small scales for M > 5e10 M_sun galaxies from the COSMOS and COMBO-17 surveys. Using a pair fraction derived from the Sloan Digital Sky Survey as a low-redshift benchmark, the large survey area at intermediate redshifts allows us to determine the evolution of the close pair fraction with unprecedented accuracy for a mass-selected sample: we find that the fraction of galaxies more massive than 5e10M_sun in pairs separated by less than 30 kpc in 3D space evolves as F(z) = (0.0130+/-0.0019)x(1+z)^1.21+/-0.25 between z = 0 and z = 1.2. Assuming a merger time scale of 0.5 Gyrs, the inferred merger rate is such that galaxies with mass in excess of 1e11 M_sun have undergone, on average, 0.5 (0.7) mergers involving progenitor galaxies both more massive than 5e10 M_sun since z = 0.6 (1.2). We also study the number density evolution of massive red sequence galaxies using published luminosity functions and constraints on the M/L evolution from the fundamental plane. Moreover, we demonstrate that the measured merger rate of massive galaxies is sufficient to explain this observed number density evolution in massive red sequence galaxies since z = 1.

We use the cosmological semi-analytic model (SAM) for galaxy formation presented in Paper I to study the metallicities and abundance ratios of the intracluster medium (ICM) within the hierarchical structure formation paradigm. By requiring a slightly flat IMF (x=1.15) and a two-population delay-time-distribution (DTD) for SN Ia explosions we found previously that this model is able to reproduce the abundance ratios and supernova rates of early-type galaxies in the local Universe. Predictions for elemental abundances in the ICM pose a further test of the model. We find that with the fiducial model from Paper I the overall metal content of the ICM is too low, although the abundance ratios are in good agreement with the data. However, we find that allowing a fraction of the metal-enriched material ejected by stars to be deposited directly into the hot ICM, instead of being deposited into the cold ISM, appears to be a plausible and physically-motivated solution.

We present results on the evolution of the intrinsic scatter of black hole masses considering different implementations of a model in which black holes only grow via mergers. We demonstrate how merger driven growth affects the correlations between black hole mass and host bulge mass. The simple case of an initially log-normal distributed scatter in black hole and bulge masses combined with random merging within the galaxy population results in a decreasing scatter with merging generation/number as predicted by the Central-limit theorem. In general we find that the decrease in scatter {\sigma} is well approximated by {\sigma}merg(m) = {\sigma}ini \times (m + 1)^(-a/2) with a = 0.42 for a range of mean number of mergers m < 50. For a large mean number of mergers (m > 100) we find a convergence to a = 0.61. This is valid for a wide range of different initial distributions, refill-scenarios or merger mass-ratios. Growth scenarios based on halo merger trees of a (100 Mpc)^3 dark matter LambdaCDM-simulation show a similar behaviour with a scatter decrease of a = 0.30 with typical number of mergers m < 50 consistent with random merging (best matching model: a = 0.34). Assuming a present day scatter of 0.3 dex in black hole mass and a mean number of mergers not exceeding m = 50 our results imply a scatter of 0.6 dex at z = 3 and thus a possible scenario in which overmassive (and undermassive) black holes at high redshift are a consequence of a larger intrinsic scatter in black hole mass. A simple toy model connecting the growth of black holes to the growth of LambdaCDM dark matter halos via mergers, neglecting any contribution from accretion, yields a consistent M\cdot -MBulge relation at z = 0 - if we assume the correct initial relation.

Deep pencil beam surveys (<1 deg^2) are of fundamental importance for studying the high-redshift universe. However, inferences about galaxy population properties are in practice limited by 'cosmic variance'. This is the uncertainty in observational estimates of the number density of galaxies arising from the underlying large-scale density fluctuations. This source of uncertainty can be significant, especially for surveys which cover only small areas and for massive high-redshift galaxies. Cosmic variance for a given galaxy population can be determined using predictions from cold dark matter theory and the galaxy bias. In this paper we provide tools for experiment design and interpretation. For a given survey geometry we present the cosmic variance of dark matter as a function of mean redshift z and redshift bin size Dz. Using a halo occupation model to predict galaxy clustering, we derive the galaxy bias as a function of mean redshift for galaxy samples of a given stellar mass range. In the linear regime, the cosmic variance of these galaxy samples is the product of the galaxy bias and the dark matter cosmic variance. We present a simple recipe using a fitting function to compute cosmic variance as a function of the angular dimensions of the field, z, Dz and stellar mass m*. We also provide tabulated values and a software tool. We find that for GOODS at z=2 and with Dz=0.5 the relative cosmic variance of galaxies with m*>10^11 Msun is ~38%, while it is ~27% for GEMS and ~12% for COSMOS. For galaxies of m*~10^10 Msun the relative cosmic variance is ~19% for GOODS, ~13% for GEMS and ~6% for COSMOS. This implies that cosmic variance is a significant source of uncertainty at z=2 for small fields and massive galaxies, while for larger fields and intermediate mass galaxies cosmic variance is less serious.

We study the effect of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. We produce a suite of minor merger simulations for Milky Way-like galaxies. This suite consists of collisionless simulations as well as hydrodynamical runs including a gaseous component in the galactic disc. We find that in dissipationless simulations minor mergers cause the scale height of the disc to increase by up to a factor of ~2. When the presence of gas in the disc is taken into account this thickening is reduced by 25% (50%) for an initial disc gas fraction of 20% (40%), leading to a final scale height z0 between 0.6 and 0.7 kpc, regardless of the initial scale height. We argue that the presence of gas reduces disc heating via two mechanisms: absorption of kinetic impact energy by the gas and/or formation of a new thin stellar disc that can cause heated stars to recontract towards the disc plane. We show that in our simulations most of the gas is consumed during the merger and thus the regrowth of a new thin disc has a negligible impact on the z0 of the post merger galaxy. Final disc scale heights found in our simulations are in good agreement with studies of the vertical structure of spiral galaxies where the majority of the systems are found to have scale heights of 0.4 kpc < z0 < 0.8 kpc. We also found no tension between recent measurements of the scale height of the Milky Way thin disc and results coming from our hydrodynamical simulations. We conclude that the existence of a thin disc in the Milky Way and in external galaxies is not in obvious conflict with the predictions of the CDM model.

We use a statistical approach to determine the relationship between the stellar masses of galaxies and the masses of the dark matter halos in which they reside. We obtain a parameterized stellar-to-halo mass (SHM) relation by populating halos and subhalos in an N-body simulation with galaxies and requiring that the observed stellar mass function be reproduced. We find good agreement with constraints from galaxy-galaxy lensing and predictions of semi-analytic models. Using this mapping, and the positions of the halos and subhalos obtained from the simulation, we find that our model predictions for the galaxy two-point correlation function (CF) as a function of stellar mass are in excellent agreement with the observed clustering properties in the SDSS at z=0. We show that the clustering data do not provide additional strong constraints on the SHM function and conclude that our model can therefore predict clustering as a function of stellar mass. We compute the conditional mass function, which yields the average number of galaxies with stellar masses in the range [m, m+dm] that reside in a halo of mass M. We study the redshift dependence of the SHM relation and show that, for low mass halos, the SHM ratio is lower at higher redshift. The derived SHM relation is used to predict the stellar mass dependent galaxy CF and bias at high redshift. Our model predicts that not only are massive galaxies more biased than low mass ones at all redshifts, but the bias increases more rapidly with increasing redshift for massive galaxies than for low mass ones. We present convenient fitting functions for the SHM relation as a function of redshift, the conditional mass function, and the bias as a function of stellar mass and redshift.

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

The star formation rate (SFR) and black hole accretion rate (BHAR) functions are measured to be proportional to each other at z < ~3. This close correspondence between SF and BHA would naturally yield a BH mass-galaxy mass correlation, whereas a BH mass-bulge mass correlation is observed. To explore this apparent contradiction we study the SF in spheroid-dominated galaxies between z=1 and the present day. We use 903 galaxies from the COMBO-17 survey with M* >2x10^10M_sun, ultraviolet and infrared-derived SFRs from Spitzer and GALEX, and morphologies from GEMS HST/ACS imaging. Using stacking techniques, we find that <25% of all SF occurs in spheroid-dominated galaxies (Sersic index n>2.5), while the BHAR that we would expect if the global scalings held is three times higher. This rules out the simplest picture of co-evolution, in which SF and BHA trace each other at all times. These results could be explained if SF and BHA occur in the same events, but offset in time, for example at different stages of a merger event. However, one would then expect to see the corresponding star formation activity in early-stage mergers, in conflict with observations. We conclude that the major episodes of SF and BHA occur in different events, with the bulk of SF happening in isolated disks and most BHA occurring in major mergers. The apparent global co-evolution results from the regulation of the BH growth by the potential well of the galactic spheroid, which includes a major contribution from disrupted disk stars.

We present a study of galaxy mergers and the influence of environment in the Abell 901/902 supercluster at z~0.165. We use HST ACS F606W data from the STAGES survey, COMBO-17, Spitzer 24um, and XMM-Newton X-ray data. Our analysis utilizes both a visual classification system, and quantitative CAS parameters to identify systems which show evidence of a recent or ongoing merger of mass ratio >1/10. Our results are: (1) After visual classification and minimizing the contamination from false projection pairs, we find that the merger fraction f_merge is 0.023+/-0.007. The estimated fractions of likely major mergers, likely minor mergers, and ambiguous cases are 0.01+/-0.004, 0.006+/-0.003, and 0.007+/-0.003, respectively. (2) The mergers lie outside the cluster core of radius R < 0.25 Mpc: the lack of mergers in the core is likely due to the large galaxy velocity dispersion in the core. Mergers populate the region (0.25 Mpc < R <= 2 Mpc) between the core and outskirt. In this region, the estimated frequency of mergers is similar to those seen at typical group overdensities. This suggests ongoing growth of the clusters via accretion of group and field galaxies. (3) We compare our observed merger fraction with those reported in other clusters and groups out to z~0.4. Existing data points on the merger fraction for L<= L* galaxies in clusters allow for a range of evolutionary scenarios. (4) The fraction of mergers, which lie on the blue cloud is 80%+/-18% versus 34%+/-7% for non-interacting galaxies, implying that interacting galaxies are preferentially blue. (5) The average SFR, based on UV or UV+IR data, is enhanced by a factor of ~1.5 to 2 in mergers compared to non-interacting galaxies. However, mergers in the clusters contribute only a small fraction (between 10% and 15%) of the total SFR density.(Abridged)

We study the metallicities and abundance ratios of early-type galaxies in cosmological semi-analytic models (SAMs) within the hierarchical galaxy formation paradigm. To achieve this we implemented a detailed galactic chemical evolution (GCE) model and can now predict abundances of individual elements for the galaxies in the semi-analytic simulations. This is the first time a SAM with feedback from Active Galactic Nuclei (AGN) has included a chemical evolution prescription that relaxes the instantaneous recycling approximation. We find that the new models are able to reproduce the observed mass-metallicity (M*-[Z/H]) relation and, for the first time in a SAM, we reproduce the observed positive slope of the mass-abundance ratio (M*-[$alpha$/Fe]) relation. Our results indicate that in order to simultaneously match these observations of early-type galaxies, the use of both a very mildly top-heavy IMF (i.e., with a slope of x=1.15 as opposed to a standard x=1.3), and a lower fraction of binaries that explode as Type Ia supernovae appears to be required. We also examine the rate of supernova explosions in the simulated galaxies. In early-type (non-star forming) galaxies, our predictions are also consistent with the observed SNe rates. However, in star-forming galaxies, a higher fraction of SN Ia binaries than in our preferred model is required to match the data. If, however, we deviate from the classical model and introduce a population of SNe Ia with very short delay times, our models simultaneously produce a good match to the observed metallicities, abundance ratios and SN rates.

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in halos according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T~0.5 at <10^10 M_sun), with almost no bulgeless systems), even if they have avoided major mergers. Including the different physical behavior of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving B/T~0.1 at <10^10 M_sun, with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behavior of gas have difficulty reproducing the strong observed morphology-mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.

I investigate whether useful constraints on the evolution of the relationship between galaxy mass and black hole (BH) mass can be obtained from recent measurements of galaxy stellar mass functions and QSO bolometric luminosity functions at high redshift. I assume a simple power-law relationship between galaxy mass and BH mass, as implied by BH mass measurements at low redshift, and consider only evolution in the zero-point of the relation. I argue that one can obtain a lower limit on the zero-point evolution by assuming that every galaxy hosts a BH, shining at its Eddington rate. One can obtain an upper limit by requiring that the number of massive BH at high redshift does not exceed that observed locally. I find that, under these assumptions, and neglecting scatter in the BH-galaxy mass relation, BH must have been a factor of about 2 times larger at z=1 and 5 to 6 times more massive relative to their host galaxies at z=2. However, accounting for intrinsic scatter in the BH-galaxy mass relationship considerably relaxes these constraints. With a logarithmic scatter of 0.3 to 0.5 dex in black hole mass at fixed galaxy mass, similar to estimates of the intrinsic scatter in the observed relation today, there are enough massive BH to produce the observed population of luminous QSOs at z=2 even in the absence of any zero-point evolution. Adopting more realistic estimates for the fraction of galaxies that host active BH and the Eddington ratios of the associated quasars, I find that the zero-point of the BH-galaxy mass relation at z=2 cannot be much more than a factor of two times larger than the present-day value, as the number of luminous quasars predicted would exceed the observed population.

We present new calculations of the evolving UV background out to the epoch of cosmological reionization and make predictions for the amount of GeV gamma-ray attenuation by electron-positron pair production. Our results are based on recent semi-analytic models of galaxy formation, which provide predictions of the dust-extinguished UV radiation field due to starlight, and empirical estimates of the contribution due to quasars. We account for the reprocessing of ionizing photons by the intergalactic medium. We test whether our models can reproduce estimates of the ionizing background at high redshift from flux decrement analysis and proximity effect measurements from quasar spectra, and identify a range of models that can satisfy these constraints. Pair-production against soft diffuse photons leads to a spectral cutoff feature for gamma rays observed between 10 and 100 GeV. This cutoff varies with redshift and the assumed star formation and quasar evolution models. We find only negligible amounts of absorption for gamma rays observed below 10 GeV for any emission redshift. With observations of high-redshift sources in sufficient numbers by the Fermi Gamma-ray Space Telescope and new ground-based instruments it should be possible to constrain the extragalactic background light in the UV and optical portion of the spectrum.

Both observations and simulations show that major tidal interactions or mergers between gas-rich galaxies can lead to intense bursts of starformation. Yet, the average enhancement in star formation rate (SFR) in major mergers and the contribution of such events to the cosmic SFR are not well estimated. Here we use photometric redshifts, stellar masses and UV SFRs from COMBO-17, 24 micron SFRs from Spitzer and morphologies from two deep HST cosmological survey fields (ECDFS/GEMS and A901/STAGES) to study the enhancement in SFR as a function of projected galaxy separation. We apply two-point projected correlation function techniques, which we augment with morphologically-selected very close pairs (separation <2 arcsec) and merger remnants from the HST imaging. Our analysis confirms that the most intensely star-forming systems are indeed interacting or merging. Yet, for massive (M* > 10^10 Msun) star-forming galaxies at 0.4<z<0.8, we find that the SFRs of galaxies undergoing a major interaction (mass ratios <1:4 and separations < 40 kpc) are only 1.80 +/- 0.30 times higher than the SFRs of non-interacting galaxies when averaged over all interactions and all stages of the interaction, in good agreement with other observational works. We demonstrate that these results imply that <10% of star formation at 0.4 < z < 0.8 is triggered directly by major mergers and interactions; these events are not important factors in the build-up of stellar mass since z=1.

We investigate the effect of dry merging on the color-magnitude relation (CMR) of galaxies and find that the amount of merging predicted by a hierarchical model results in a red sequence that compares well with the observed low-redshift relation. A sample of 29,000 early-type galaxies selected from the Sloan Digital Sky Survey Data Release 6 shows that the bright end of the CMR has a shallower slope and smaller scatter than the faint end. This magnitude dependence is predicted by a simple toy model in which gas-rich mergers move galaxies onto a "creation red sequence" (CRS) by quenching their star formation, and subsequent mergers between red, gas-poor galaxies (so-called "dry" mergers) move galaxies along the relation. We use galaxy merger trees from a semianalytic model of galaxy formation to test the amplitude of this effect and find a change in slope at the bright end that brackets the observations, using gas fraction thresholds of 10-30% to separate wet and dry mergers. A more realistic model that includes scatter in the CRS shows that dry merging decreases the scatter at the bright end. Contrary to previous claims, the small scatter in the observed CMR thus cannot be used to constrain the amount of dry merging.

[abridged] It has been widely claimed that several lines of observational evidence point towards a "downsizing" (DS) of the process of galaxy formation over cosmic time. This behavior is sometimes termed "anti-hierarchical", and contrasted with the "bottom-up" assembly of the dark matter structures in Cold Dark Matter models. In this paper we address three different kinds of observational evidence that have been described as DS: the stellar mass assembly, star formation rate and the ages of the stellar populations in local galaxies. We compare a broad compilation of available data-sets with the predictions of three different semi-analytic models of galaxy formation within the Lambda-CDM framework. In the data, we see only weak evidence at best of DS in stellar mass and in star formation rate. We find that, when observational errors on stellar mass and SFR are taken into account, the models acceptably reproduce the evolution of massive galaxies, over the entire redshift range that we consider. However, lower mass galaxies are formed too early in the models and are too passive at late times. Thus, the models do not correctly reproduce the DS trend in stellar mass or the archaeological DS, while they qualitatively reproduce the mass-dependent evolution of the SFR. We demonstrate that these discrepancies are not solely due to a poor treatment of satellite galaxies but are mainly connected to the excessively efficient formation of central galaxies in high-redshift haloes with circular velocities ~100-200 km/s. [abridged]

We present a simple formalism to interpret two galaxy statistics, the UV luminosity function and two-point correlation functions for star-forming galaxies at z~4, 5, 6 in the context of LCDM cosmology. Both statistics are the result of how star formation takes place in DM halos, and thus are used to constrain how UV light depends on halo properties such as mass. The two measures were taken from the GOODS data, thus ideal for joint analysis. The two physical quantities we explore are the SF duty cycle, and the range of L_UV that a halo of mass M can have (mean and variance). The former addresses the typical duration of SF activity in halos while the latter addresses the averaged SF history and regularity of gas inflow into these systems. We explore various physical models consistent with data, and find the following: 1) the typical duration of SF observed in the data is <0.4 Gyr (1 sig), 2) the inferred scaling law between L_UV and halo mass M from the observed slope of the LFs is roughly linear at all redshifts, and 3) L_UV for a fixed halo mass decreases with time, implying that the SF efficiency (after dust extinction) is higher at earlier times. We explore several physical scenarios relating star formation to halo mass, but find that these scenarios are indistinguishable due to the limited range of halo mass probed by our data. In order to discriminate between different scenarios, we discuss constraining the bright-faint galaxy cross-correlation functions and luminosity-dependence of galaxy bias. (Abridged)

We perform a comprehensive estimate of the frequency of galaxy mergers and their impact on star formation over z~0.24--0.80 (lookback time T_b~3--7 Gyr) using 3698 (M*>=1e9 Msun) galaxies with GEMS HST, COMBO-17, and Spitzer data. Our results are: (1) Among 790 high mass (M*>=2.5e10 Msun) galaxies, the visually-based merger fraction over z~0.24--0.80, ranges from 9%+-5% to 8%+-2%. Lower limits on the major and minor merger fractions over this interval range from 1.1% to 3.5%, and 3.6% to 7.5%, respectively. This is the first approximate empirical estimate of the frequency of minor mergers at z<1. For a visibility timescale of ~0.5 Gyr, it follows that over T_b~3--7 Gyr, ~68% of high mass systems have undergone a merger of mass ratio >1/10, with ~16%, 45%, and 7% of these corresponding respectively to major, minor, and ambiguous `major or minor' mergers. The mean merger rate is a few x 1e-4 Gyr-1 Mpc-3. (2) We compare the empirical merger fraction and rate for high mass galaxies to a suite of Lambda CDM-based models: halo occupation distribution models, semi-analytic models, and hydrodynamic SPH simulations. We find qualitative agreement between observations and models such that the (major+minor) merger fraction or rate from different models bracket the observations, and show a factor of five dispersion. Near-future improvements can now start to rule out certain merger scenarios. (3) Among ~3698 M*>=1e9 Msun galaxies, we find that the mean SFR of visibly merging systems is only modestly enhanced compared to non-interacting galaxies over z~0.24--0.80. Visibly merging systems only account for less than 30% of the cosmic SFR density over T_b~3--7 Gyr. This suggests that the behavior of the cosmic SFR density over the last 7 Gyr is predominantly shaped by non-interacting galaxies.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

We investigate the effect of gas-poor (so-called "dry") mergers on the color-magnitude relation (CMR) of early-type galaxies through a simple toy model and compare with low-z observations from the Sloan Digital Sky Survey (SDSS). The observed red sequence shows a tilt towards bluer colors and a decrease in scatter at the bright end. These characteristics are predicted by our model, based on merger trees from a semi-analytic model of galaxy formation. We assume galaxies move onto a "creation red sequence" when they undergo major gas-rich mergers. Subsequent dry mergers move galaxies along the relation by increasing their mass, but also make them slightly bluer. This occurs because bright galaxies are most likely to merge with one of the more numerous fainter and consequently bluer galaxies that lie further down the relation. Bright galaxies undergo a higher fraction of dry mergers than faint galaxies, which causes a change in the slope of the CMR. A more realistic model that includes scatter in the initial relation shows that dry merging causes a tightening of the CMR towards the bright end. The small scatter in the observed CMR thus cannot be used to argue against significant mass growth from dry merging.

Stellar population parameters derived from spectral line-strengths provide a powerful probe of galaxy properties and formation histories. We implement the machinery for extracting single-stellar-population-equivalent stellar population parameters from synthetic spectra generated by a hierarchical galaxy formation model. We find that the SSP-equivalent age is related to the light-weighted age in a complicated fashion that reflects the influence of recently-formed stars and is poorly correlated with the mass-weighted age. The tendency for SSP-equivalent ages to be biased young means that archaeological downsizing overstates the mass-weighted downsizing in age with mass. We find that the SSP-equivalent metallicity closely tracks the mass- and light-weighted metallicities, so that observed mass--metallicity relations for old galaxies closely reflect the underlying trends. We construct mock catalogues of early-type galaxies in a Coma cluster-sized halo and compare them directly to observations of early-type galaxies in the Coma cluster. The similarity of the SSP-equivalent ages in the observational samples and the mock catalogues gives us confidence that the star-formation quenching implemented in the hierarchical galaxy formation model produces roughly the correct amount of recent star formation. The SSP-equivalent metallicities are however too low and have the wrong slope as a function of velocity dispersion, and the SSP-equivalent ages of the model galaxies may have an incorrect slope as a function of velocity dispersion. (Abridged)

We investigate the correlation of star formation quenching with internal galaxy properties and large scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based Group Catalog of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer (GALEX) are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of "central" and "satellite" galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with AGN feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This {\em satellite overquenching problem} is caused by the too-rapid stripping of the satellites' hot gas halos, which leads to rapid strangulation of star formation.

We investigate the biases and uncertainties in estimates of physical parameters of high-redshift Lyman break galaxies (LBGs), such as stellar mass, mean stellar population age, and star formation rate (SFR), obtained from broad-band photometry. By combining LCDM hierarchical structure formation theory, semi-analytic treatments of baryonic physics, and stellar population synthesis models, we construct model galaxy catalogs from which we select LBGs at redshifts z ~ 3.4, 4.0, and 5.0. The broad-band spectral energy distributions (SEDs) of these model LBGs are then analysed by fitting galaxy template SEDs derived from stellar population synthesis models with smoothly declining SFRs. We compare the statistical properties of LBGs' physical parameters -- such as stellar mass, SFR, and stellar population age -- as derived from the best-fit galaxy templates with the intrinsic values from the semi-analytic model. We find some trends in these distributions: first, when the redshift is known, SED-fitting methods reproduce the input distributions of LBGs' stellar masses relatively well, with a minor tendency to underestimate the masses overall, but with substantial scatter. Second, there are large systematic biases in the distributions of best-fit SFRs and mean ages, in the sense that single-component SED-fitting methods underestimate SFRs and overestimate ages. We attribute these trends to the different star formation histories predicted by the semi-analytic models and assumed in the galaxy templates used in SED-fitting procedure, and to the fact that light from the current generation of star-formation can hide older generations of stars. These biases, which arise from the SED-fitting procedure, can significantly affect inferences about galaxy evolution from broadband photometry.

Attenuation of high--energy gamma rays by pair--production with UV, optical and IR background photons provides a link between the history of galaxy formation and high--energy astrophysics. We present results from our latest semi-analytic models (SAMs), based upon a $\Lambda$CDM hierarchical structural formation scenario and employing all ingredients thought to be important to galaxy formation and evolution, as well as reprocessing of starlight by dust to mid- and far-IR wavelengths. Our models also use results from recent hydrodynamic galaxy merger simulations. These latest SAMs are successful in reproducing a large variety of observational constraints such as number counts, luminosity and mass functions, and color bimodality. We have created 2 models that bracket the likely ranges of galaxy emissivities, and for each of these we show how the optical depth from pair--production is affected by redshift and gamma-ray energy. We conclude with a discussion of the implications of our work, and how the burgeoning science of gamma-ray astronomy will continue to help constrain cosmology.

We present two models for the cosmological UV background light, and calculate the opacity of GeV gamma--rays out to redshift 9. The contributors to the background include 2 possible quasar emissivities, and output from star--forming galaxies as determined by recent a semi--analytic model (SAM) of structure formation. The SAM used in this work is based upon a hierarchical build-up of structure in a $\Lambda$CDM universe and is highly successful in reproducing a variety of observational parameters. Above 1 Rydberg energy, ionizing radiation is subject to reprocessing by the IGM, which we treat using our radiative transfer code, CUBA. The two models for quasar emissivity differing above z = 2.3 are chosen to match the ionization rates observed using flux decrement analysis and the higher values of the line-of-sight proximity effect. We also investigate the possibility of a flat star formation rate density at z $>5$. We conclude that observations of gamma--rays from 10 to 100 GeV by Fermi (GLAST) and the next generation of ground based experiments should confirm a strongly evolving opacity from $1<$ z $<4$. Observation of attenuation in the spectra of gamma--ray bursts at higher redshift could constrain emission of UV radiation at these early times, either from a flat or increasing star-formation density or an unobserved population of sources.

The treatment of dust attenuation is crucial in order to compare the predictions of galaxy formation models with multiwavelength observations. Most past studies have either used simple analytic prescriptions or else full radiative transfer (RT) calculations. Here, we couple star formation histories and morphologies predicted by the semi-analytic galaxy formation model MORGANA with RT calculations from the spectrophotometric and dust code GRASIL to create a library of galaxy SEDs from the UV/optical through the far Infrared, and compare the predictions of the RT calculations with analytic prescriptions. We consider a low and high redshift sample, as well as an additional library constructed with empirical, non-cosmological star formation histories and simple (pure bulge or disc) morphologies. Based on these libraries, we derive fitting formulae for the effective dust optical depth as a function of galaxy physical properties such as metallicity, gas mass, and radius. We show that such fitting formulae can predict the V-band optical depth with a scatter smaller than 0.4 dex for both the low and high redshift samples, but that there is a large galaxy-to-galaxy scatter in the shapes of attenuation curves, probably due to geometrical variations, which our simple recipe does not capture well. However, our new recipe provides a better approximation to the GRASIL results at optical wavelength than standard analytic prescriptions from the literature, particularly at high redshift.

We present a new semi-analytic model that self-consistently traces the growth of supermassive black holes (BH) and their host galaxies within the context of the LCDM cosmological framework. In our model, the energy emitted by accreting black holes regulates the growth of the black holes themselves, drives galactic scale winds that can remove cold gas from galaxies, and produces powerful jets that heat the hot gas atmospheres surrounding groups and clusters. We present a comprehensive comparison of our model predictions with observational measurements of key physical properties of low-redshift galaxies, such as cold gas fractions, stellar metallicities and ages, and specific star formation rates. We find that our new models successfully reproduce the exponential cutoff in the stellar mass function and the stellar and cold gas mass densities at z~0, and predict that star formation should be largely, but not entirely, quenched in massive galaxies at the present day. We also find that our model of self-regulated BH growth naturally reproduces the observed relation between BH mass and bulge mass. We explore the global formation history of galaxies in our models, presenting predictions for the cosmic histories of star formation, stellar mass assembly, cold gas, and metals. We find that models assuming the "concordance" LCDM cosmology overproduce star formation and stellar mass at high redshift (z>2). A model with less small-scale power predicts less star formation at high redshift, and excellent agreement with the observed stellar mass assembly history, but may have difficulty accounting for the cold gas in quasar absorption systems at high redshift (z~3-4).

The observational link between Supermassive Black Holes (SMBH) and galaxies at low redshift seems to be very tight, and statistically the global evolution of star formation activity and BH accretion activity also seem to trace each other closely. However, pinning down the co-evolution of galaxies and BH on an object-by-object basis remains elusive. I present results from new models for the joint evolution of galaxies, SMBH, and AGN, which may be able to help resolve some of the observational puzzles. A unique aspect of these models is our treatment of self-regulated BH growth based on hydrodynamic simulations of galaxy-galaxy mergers. Although these models do quite well at reproducing the observed evolution of galaxies, they do not reproduce the observed history of BH accretion, predicting too much early accretion and not enough at late times. I suggest two possible resolutions to this problem.

The QSO luminosity function at z>5 provides strong constraints on models of joint evolution of QSO and their hosts. However, these observations are challenging because the low space densities of these objects necessitate surveying of large areas, in order to obtain statistically meaningful samples, while at the same time cosmological redshifting and dimming means that rather deep Near Infrared (NIR) imaging must be carried out. Several upcoming and proposed facilities with wide-field NIR imaging capabilities will open up this new region of parameter space. In this paper we present predictions for the expected number counts of z>5 QSOs, based on simple empirical and semi-empirical models of QSO evolution, as a function of redshift, depth and surveyed area. We compute the evolution of observed-frame QSO magnitudes and colors in a representative photometric system covering the wavelength range 550nm<lambda<1800nm, and combine this information with different estimates for the evolution of the QSO luminosity function. We conclude that planned ground-based surveys such as Pan-STARRS and VISTA should be able to detect a large number of luminous QSOs up to z<7.5, but that space-based missions such as EUCLID (formerly SPACE/DUNE) or SNAP are probably required in order to obtain substantial samples at higher redshift. We also use our models to predict the expected number counts for future X-ray space missions (such as XEUS and Constellation-X), and show that because of their small field-of-view, these telescopes are unlikely to discover significant numbers of AGN at very high redshift. However, X-ray follow-up of objects detected at longer wavelength will be an important means of confirming their identity as AGN and constraining obscuration.

We use a SDSS galaxy group catalogue to study the dependence of galaxy activity on stellar mass, halo mass, and group hierarchy (centrals vs. satellites). We split our galaxy sample in star-forming galaxies, galaxies with optical AGN activity and radio sources. We find a smooth transition in halo mass as the activity of central galaxies changes from star formation to optical AGN activity to radio emission. Star-forming centrals preferentially reside in haloes with M<10^{12} Msun, central galaxies with optical-AGN activity typically inhabit haloes with M \sim 10^{13} Msun, and centrals emitting in the radio mainly reside in haloes more massive than 10^{14} Msun. Although this seems to suggest that the environment (halo mass) determines the type of activity of its central galaxy, we find a similar trend with stellar mass: central star formers typically have stellar masses below 10^{10} Msun, while optical-AGN hosts and central radio sources have characteristic stellar masses of 10^{10.8} Msun and 10^{11.6} Msun, respectively. Since more massive haloes typically host more massive centrals, it is unclear whether the activity of a central galaxy is causally connected to its stellar mass or to its halo mass. In general, satellite galaxies have their activity suppressed wrt central galaxies of the same stellar mass. At fixed stellar mass, we find that the activity of satellite galaxies depends only weakly on halo mass. In fact, for satellite galaxies the dependence of galaxy activity on halo mass is more than four times weaker than the dependence on stellar mass. As we discuss, all these results are consistent with a picture in which low mass haloes accrete cold gas, while massive haloes have coronae of hot gas that promote radio activity of their central galaxies. [Abridged]

We present a study of galaxies in the STAGES survey of the Abell 901/902 supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps. We characterize galaxies with strong, externally-triggered morphological distortions and normal, relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We report here our preliminary results.

We explore the history and impact of galaxy mergers and interactions over z~0.24 to 0.80, based on HST ACS, Combo-17, and Spitzer 24 mu data of ~4500 galaxies in the GEMS survey. Using visual and quantitative parameters,we identify galaxies with strong distortions indicative of recent strong interactions and mergers versus normal galaxies (E/S0, Sa, Sb-Sc, Sd/Irr). Our results are: (1) The observed fraction F of strongly disturbed systems among high mass (M>=2.5E10 Msun) galaxies is ~9% to 12% in every Gyr bin over z~0.24 to 0.80. The corresponding merger rate is a few times 10^-4 galaxies Gyr-1 Mpc-3. The fraction F shows fair agreement with the merger fraction of mass ratio >=1:10 predicted by several LCDM-based simulations. (2) For M>=1E9 Msun systems, the average star formation rate (SFR) of strongly disturbed systems is only modestly enhanced with respect to normal galaxies, in agreement with recent simulations. In fact, over z~0.24 to 0.80, strongly disturbed systems only account for a small fraction (<~30%) of the total SFR density. This suggests that the behaviour of the cosmic SFR density over the last seven billion years is predominantly shaped by normal galaxies.

We explore galaxy assembly over the last seven billion years by characterizing "normal" galaxies along the Hubble sequence, against strongly disturbed merging/interacting galaxies with the widely used CAS system of concentration (C), asymmetry (A), and 'clumpiness' (S) parameters, as well as visual classification. We analyze Hubble Space Telescope (HST) ACS images of ~4000 intermediate and high mass (> 10^9 solar masses) galaxies from the GEMS survey, one of the largest HST surveys conducted to date in two filters. We explore the effectiveness of the CAS criteria [A>S and A>~0.35] in separating normal and strongly disturbed galaxies at different redshifts, and quantify the recovery and contamination rate. We also compare the average star formation rate and the cosmic star formation rate density as a function of redshift between normal and interacting systems identified by CAS.

We present a study of galaxies in the Abell 901/902 Supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps, as part of the STAGES survey. We characterize galaxies with strong externally-triggered morphological distortions and normal relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We revisit the morphology density relation, which postulates a higher fraction of early type galaxies in dense environments, by considering separately galaxies with a low bulge-to-disk (B/D) ratio and a low gas content as these two parameters may not be correlated in clusters. We report here on our preliminary analysis.

Using hydrodynamic simulations of disc-galaxy major mergers, we investigate the star formation history and remnant properties when various parametrizations of a simple stellar feedback model are implemented. The simulations include radiative cooling, a density-dependent star formation recipe and a model for feedback from massive stars. The feedback model stores supernova feedback energy within individual gas particles and dissipates this energy on a time-scale specified by two free parameters; tau_fb, which sets the dissipative time-scale, and n, which sets the effective equation of state in star-forming regions. Using a self-consistent disc galaxy, modelled after a local Sbc spiral, in both isolated and major-merger simulations, we investigate parametrizations of the feedback model that are selected with respect to the quiescent disc stability. These models produce a range of star formation histories that are consistent with the star formation relation found by Kennicutt. All major mergers produce a population of new stars that is highly centrally concentrated, demonstrating a distinct break in the r1/4 surface density profile, consistent with previous findings. The half-mass radius and one-dimensional velocity dispersion are affected by the feedback model used. Finally, we compare our results to those of previous simulations of star formation in disc-galaxy major mergers, addressing the effects of star formation normalization, the version of smoothed particle hydrodynamics (SPH) employed and assumptions about the interstellar medium.

We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger--driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disk galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the local Universe. We find that the merger--driven star formation is a strong function of merger mass ratio, with very little, if any, induced star formation for large mass ratio mergers. We define a burst efficiency that is useful to characterize the merger--driven star formation and test that it is insensitive to uncertainties in the feedback parameterization. In accord with previous work we find that the burst efficiency depends on the structure of the primary galaxy. In particular, the presence of a massive stellar bulge stabilizes the disk and suppresses merger--driven star formation for large mass ratio mergers. Direct, co--planar merging orbits produce the largest tidal disturbance and yield that most intense burst of star formation. Contrary to naive expectations, a more compact distribution of gas or an increased gas fraction both decrease the burst efficiency. Owing to the efficient feedback model and the newer version of SPH employed here, the burst efficiencies of the mergers presented here are smaller than in previous studies.

Surveys of distant galaxies with the Hubble Space Telescope and from the ground have shown that there is only mild evolution in the relationship between radial size and stellar mass for galactic disks from z~1 to the present day. Using a sample of nearby disk-dominated galaxies from the Sloan Digital Sky Survey (SDSS), and high redshift data from the GEMS (Galaxy Evolution from Morphology and SEDs) survey, we investigate whether this result is consistent with theoretical expectations within the hierarchical paradigm of structure formation. The relationship between virial radius and mass for dark matter halos in the LCDM model evolves by about a factor of two over this interval. However, N-body simulations have shown that halos of a given mass have less centrally concentrated mass profiles at high redshift. When we compute the expected disk size-stellar mass distribution, accounting for this evolution in the internal structure of dark matter halos and the adiabatic contraction of the dark matter by the self-gravity of the collapsing baryons, we find that the predicted evolution in the mean size at fixed stellar mass since z~1 is about 15-20 percent, in good agreement with the observational constraints from GEMS. At redshift z~2, the model predicts that disks at fixed stellar mass were on average only 60% as large as they are today. Similarly, we predict that the rotation velocity at a given stellar mass (essentially the zero-point of the Tully-Fisher relation) is only about 10 percent larger at z~1 (20 percent at z~2) than at the present day.

We present the clustering properties of 151 Lyman alpha emitting galaxies at z ~ 4.5 selected from the Large Area Lyman Alpha (LALA) survey. Our catalog covers an area of 36' x 36' observed with five narrowband filters. We assume that the angular correlation function w(theta) is well represented by a power law A_w = Theta^(-beta) with slope beta = 0.8, and we find A_w = 6.73 +/- 1.80. We then calculate the correlation length r_0 of the real-space two-point correlation function xi(r) = (r/r_0)^(-1.8) from A_w through the Limber transformation, assuming a flat, Lambda-dominated universe. Neglecting contamination, we find r_0 = 3.20 +/- 0.42 Mpc/h. Taking into account a possible 28% contamination by randomly distributed sources, we find r_0 = 4.61 +/- 0.6 Mpc/h. We compare these results with the expectations for the clustering of dark matter halos at this redshift in a Cold Dark Matter model, and find that the measured clustering strength can be reproduced if these objects reside in halos with a minimum mass of 1-2 times 10^11 Solar masses/h. Our estimated correlation length implies a bias of b ~ 3.7, similar to that of Lyman-break galaxies (LBG) at z ~ 3.8-4.9. However, Lyman alpha emitters are a factor of ~ 2-16 rarer than LBGs with a similar bias value and implied host halo mass. Therefore, one plausible scenario seems to be that Lyman alpha emitters occupy host halos of roughly the same mass as LBGs, but shine with a relatively low duty cycle of 6-50%.

In the context of measuring structure and morphology of intermediate redshift galaxies with recent HST/ACS surveys, we tune, test, and compare two widely used fitting codes (GALFIT and GIM2D) for fitting single-component Sersic models to the light profiles of both simulated and real galaxy data. We find that fitting accuracy depends sensitively on galaxy profile shape. Exponential disks are well fit with Sersic models and have small measurement errors, whereas fits to de Vaucouleurs profiles show larger uncertainties owing to the large amount of light at large radii. We find that both codes provide reliable fits and little systematic error, when the effective surface brightness is above that of the sky. Moreover, both codes return errors that significantly underestimate the true fitting uncertainties, which are best estimated with simulations. We find that GIM2D suffers significant systematic errors for spheroids with close companions owing to the difficulty of effectively masking out neighboring galaxy light; there appears to be no work around to this important systematic in GIM2D's current implementation. While this crowding error affects only a small fraction of galaxies in GEMS, it must be accounted for in the analysis of deeper cosmological images or of more crowded fields with GIM2D. In contrast, GALFIT results are robust to the presence of neighbors because it can simultaneously fit the profiles of multiple companions thereby deblending their effect on the fit to the galaxy of interest. We find GALFIT's robustness to nearby companions and factor of >~20 faster runtime speed are important advantages over GIM2D for analyzing large HST/ACS datasets. Finally we include our final catalog of fit results for all 41,495 objects detected in GEMS.

We combine Spitzer 24micron observations with data from the COMBO-17 survey for ~15,000 0.2<z<1 galaxies to determine how the average star formation rates (SFR) have evolved for galaxy sub-populations of different stellar masses. In the determination of <SFR> we consider both the ultraviolet (UV) and the infrared (IR) luminosities, and account for the contributions of galaxies that are individually undetected at 24micron through image stacking. For all redshifts we find that higher-mass galaxies have substantially lower specific SFR, <SFR>/<M*>, than lower-mass ones. However, we find the striking result that the rate of decline in cosmic SFR with redshift is nearly the same for massive and low-mass galaxies, i.e. NOT a strong function of stellar mass. This analysis confirms one version of what has been referred to as `downsizing', namely that the epoch of major mass build-up in massive galaxies is substantially earlier than the epoch of mass build-up in low-mass galaxies. Yet it shows that star formation activity is NOT becoming increasingly limited to low-mass galaxies towards the present epoch. We argue that this suggests that heating by AGN-powered radio jets is not the dominant mechanism responsible for the decline in cosmic SFR since z~1, which is borne out by comparison with semi-analytic models that include this effect.

We analyze star formation (SF) as a function of stellar mass (M*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS). For 2905 field galaxies, complete to 10^10(10^10.8) Msun at z<0.7(1), with Keck spectroscopic redshifts out to z=1.1, we compile SF rates (SFR) from emission lines, GALEX, and Spitzer MIPS 24 micron photometry, optical-NIR M* measurements, and HST morphologies. Galaxies with reliable signs of SF form a distinct "main sequence (MS)", with a limited range of SFR at a given M* and z (1 sigma < +-0.3 dex), and log(SFR) approximately proportional to log(M*). The range of log(SFR) remains constant to z>1, while the MS as a whole moves to higher SFR as z increases. The range of SFR along the MS constrains the amplitude of episodic variations of SF, and the effect of mergers on SFR. Typical galaxies spend ~67(95)% of their lifetime since z=1 within a factor of <~ 2(4) of their average SFR at a given M* and z. The dominant mode of the evolution of SF since z~1 is apparently a gradual decline of the average SFR in most individual galaxies, not a decreasing frequency of starburst episodes, or a decreasing factor by which SFR are enhanced in starbursts. LIRGs at z~1 seem to mostly reflect the high SFR typical for massive galaxies at that epoch. The smooth MS may reflect that the same set of few physical processes governs star formation prior to additional quenching processes. A gradual process like gas exhaustion may play a dominant role.

We analyze star formation (SF) as a function of stellar mass (M_*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS), for star-forming field galaxies with M_* >~ 10^10 M_sun out to z=1.1. The data indicate that the high specific SF rates (SFR) of many less massive galaxies do not represent late, irregular or recurrent, starbursts in evolved galaxies. They rather seem to reflect the onset (initial burst) of the dominant SF episode of galaxies, after which SF gradually declines on Gyr timescales to z=0 and forms the bulk of a galaxy's M_*. With decreasing mass, this onset of major SF shifts to decreasing z for an increasing fraction of galaxies (staged galaxy formation). This process may be an important component of the ``downsizing'' phenomenon. We find that the predominantly gradual decline of SFR (Noeske et al. 2007, this volume) can be reproduced by exponential SF histories (tau models), if less massive galaxies have systematically longer e-folding times tau, and a later onset of SF (z_f). Our model can provide a first parametrization of SFR as a function of M_* and z, and quantify mass-dependences of tau and z_f, from direct observations of M_* and SFR up to z>1. The observed evolution of SF in galaxies can plausibly reflect the dominance of gradual gas exhaustion. The data are also consistent with the history of cosmological accretion onto Dark Matter halos.

Using a model for self-regulated growth of black holes (BHs) in mergers involving gas-rich galaxies, we study the relationship between quasars and the population of merging galaxies and predict the merger-induced star formation rate density of the Universe. Mergers drive nuclear gas inflows, fueling starbursts and 'buried quasars' until accretion feedback expels the gas, rendering a briefly visible optical quasar. Star formation is shut down and accretion declines, leaving a passively evolving remnant with properties typical of red, elliptical galaxies. Based on evolution of these events in our simulations, we demonstrate that the observed statistics of merger rates, luminosity functions (LFs) and mass functions, SFR distributions, specific SFRs, quasar and quasar host galaxy LFs, and elliptical/red galaxy LFs are self-consistent and follow from one another as predicted by the merger hypothesis. We use our simulations to de-convolve both quasar and merging galaxy LFs to determine the birthrate of black holes of a given final mass and merger rates as a function of stellar mass. We use this to predict the merging galaxy LF in several observed wavebands, color-magnitude relations, mass functions, absolute and specific SFR distributions and SFR density, and quasar host galaxy LFs, as a function of redshift from z=0-6. We invert this and predict e.g. quasar LFs from observed merger LFs or SFR distributions. Our results agree well with observations, but idealized models of quasar lightcurves are ruled out by comparison of merger and quasar observations at >99.9% confidence. Using only observations of quasars, we estimate the contribution of mergers to the SFR density of the Universe even to high redshifts z~4.

We calculate the projected two point correlation function for samples of luminous and massive galaxies in the COMBO-17 photometric redshift survey, focusing particularly on the amplitude of the correlation function at small projected radii and exploring the constraints such measurements can place on the galaxy merger rate. For nearly volume-limited samples with 0.4<z<0.8, we find that 4+/-1% of luminous M_B<-20 galaxies are in close physical pairs (with real space separation of <30 proper kpc). The corresponding fraction for massive galaxies with M_*>2.5e10 M_sun is 5+/-1%. Incorporating close pair fractions from the literature, the 2dFGRS and the SDSS, we find a fairly rapid evolution of the merger fraction of massive galaxies between z=0.8 and the present day. Assuming that the major merger timescale is of order the dynamical timescale for close massive galaxy pairs, we tentatively infer that ~50% (70%) of all galaxies with present-day masses M_*>5e10 M_sun (remnants of mergers between galaxies with M_*>2.5e10 M_sun) have undergone a major merger since z=0.8(1): major mergers between massive galaxies are a significant driver of galaxy evolution over the last eight billion years.

This paper presents the Hubble Ultra Deep Field (HUDF), a one million second exposure of an 11 square minute-of-arc region in the southern sky with the Hubble Space Telescope. The exposure time was divided among four filters, F435W (B435), F606W (V606), F775W (i775), and F850LP (z850), to give approximately uniform limiting magnitudes mAB~29 for point sources. The image contains at least 10,000 objects presented here as a catalog. Few if any galaxies at redshifts greater than ~4 resemble present day spiral or elliptical galaxies. Using the Lyman break dropout method, we find 504 B-dropouts, 204 V-dropouts, and 54 i-dropouts. Using these samples that are at different redshifts but derived from the same data, we find no evidence for a change in the characteristic luminosity of galaxies but some evidence for a decrease in their number densities between redshifts of 4 and 7. The ultraviolet luminosity density of these samples is dominated by galaxies fainter than the characteristic luminosity, and the HUDF reveals considerably more luminosity than shallower surveys. The apparent ultraviolet luminosity density of galaxies appears to decrease from redshifts of a few to redshifts greater than 6. The highest redshift samples show that star formation was already vigorous at the earliest epochs that galaxies have been observed, less than one billion years after the Big Bang.

We present deep Hubble Space Telescope Advanced Camera for Surveys / High Resolution Channel U, narrow-V, and I images of nine "ultracompact" blue dwarf galaxies (UCBDs) selected from the Sloan Digital Sky Survey. We define UCBDs as local (z < 0.01) star-forming galaxies having angular diameters < 6" and physical sizes < 1 kpc. They are also among the most metal-poor galaxies known, and are found to reside within voids. Both the HST images and the objects' optical spectra reveal that they are composites of young (~1 Myr) populations that dominate their light, and older (~10 Gyr) populations that dominate their stellar masses, which we estimate to be ~10^7 - 10^8 Msol. An intermediate-age population is also indicated in most cases. The objects are not as dynamically disturbed as the prototype UCBD, POX 186, but the structure of several of them suggests that their current starbursts have been triggered by the collisions/mergers of smaller clumps of stars. In one case, HS 0822+3542, the ACS/HRC images resolve the object into two small (~100 pc) components which appear to have recently collided, supporting this interpretation. In six of the objects much of their star formation is concentrated in Young Massive Star clusters. The evidence that the galaxies consist mainly of ~10 Gyr old stars establishes that they are not protogalaxies; their low metallicities are more likely to be the result of the escape of supernova ejecta, as opposed to youth. These results are consistent with recent galaxy formation simulations which predict that cosmic re-ionization at z ~ 6 significantly limited the subsequent star formation of dwarf galaxies in voids due to the photo-evaporation of baryons from their cold dark matter halos (Abridged).

We present constraints on the evolution of the virial to stellar mass ratio of galaxies with high stellar masses in the redshift range 0.2<z<0.8, by comparing weak lensing measurements of virial mass M_vir to estimates of stellar mass M_star from COMBO-17. For a complete sample of galaxies with log(M_star/ M_\odot) > 10.5, where the majority show an early-type morphology, we find that the virial mass to stellar mass ratio is given by M_vir/M_star = 53^{+13}_{-16}. Assuming a baryon fraction from the concordance cosmology, this corresponds to a stellar fraction of baryons in massive galaxies of Omega_b^*/\Omega_b = 0.10 +/- 0.03. Analysing the galaxy sample in different redshift slices, we find little or no evolution in the virial to stellar mass ratio, and place an upper limit of ~2.5 on the growth of massive galaxies through the conversion of gas into stars from z=0.8 to the present day.

We report on the angular correlation function of Lyman-break galaxies (LBGs) at z~4 and 5 from deep samples obtained from the Great Observatories Deep Origins Survey (GOODS). Similar to LBGs at z~3, the shape of w(theta) of the GOODS LBGs is well approximated by a power-law with slope beta~0.6 at angular separation theta > 10 arcsec. The clustering strength of z~4, 5 LBGs also depends on the rest-frame UV luminosity, with brighter galaxies more strongly clustered than fainter ones, implying a general correlation between halos' mass and LBGs' star-formation rate. At smaller separations, w(theta) of deep samples significantly exceeds the extrapolation of the large-scale power-law fit, implying enhanced spatial clustering at scales r < 1 Mpc. We also find that bright LBGs statistically have more faint companions on scales theta < 20 arcsec than fainter ones, showing that the enhanced small-scale clustering is very likely due to sub-structure, namely the fact that massive halos can host multiple galaxies. A simple model for the halo occupation distribution and the CDM halo mass function reproduce well the observed w(theta). The scaling relationship of the clustering strength with volume density and with redshift is quantitatively consistent with that of CDM halos. A comparison of the clustering strength of three samples of equal luminosity limit at z ~ 3, 4 and 5 shows that the LBGs at z~5 are hosted in halos about one order of magnitude less massive than those in the lower redshift bins, suggesting that star-formation was more efficient at higher-redshift.

We have used the 28'x 28' HST image mosaic from the GEMS (Galaxy Evolution from Morphology and SEDs) survey in conjunction with the COMBO-17 photometric redshift survey to constrain the incidence of major mergers between spheroid-dominated galaxies with little cold gas (dry mergers) since z = 0.7. A set of N-body merger simulations was used to explore the morphological signatures of such interactions: they are recognizable either as < 5kpc separation close pairs or because of broad, low surface brightness tidal features and asymmetries. Data with the depth and resolution of GEMS are sensitive to dry mergers between galaxies with M_V < -20.5 for z < 0.7; dry mergers at higher redshifts are not easily recovered in single-orbit HST imaging. Six dry mergers (12 galaxies) with luminosity ratios between 1:1 and 4:1 were found from a sample of 379 red early-type galaxies with M_V < -20.5 and 0.1 < z < 0.7. The simulations suggest that the morphological signatures of dry merging are visible for ~250Myr and we use this timescale to convert the observed merger incidence into a rate. On this basis we find that present day spheroidal galaxies with M_V < -20.5 on average have undergone between 0.5 and 2 major dry mergers since z ~ 0.7. We have compared this result with the predictions of a Cold Dark Matter based semi-analytic galaxy formation model. The model reproduces the observed declining major merger fraction of bright galaxies and the space density of luminous early-type galaxies reasonably well. The predicted dry merger fraction is consistent with our observational result. Hence, hierarchical models predict and observations now show that major dry mergers are an important driver of the evolution of massive early-type galaxies in recent epochs.

We describe the data reduction and object cataloging for the GEMS survey, a large-area (800 arcmin^2) two-band (F606W and F850LP) imaging survey with the Advanced Camera for Surveys on HST, centered on the Chandra Deep Field South.

We present deep U, narrow-V, and I-band images of the ultracompact blue dwarf galaxy HS 0822+3542, obtained with the Advanced Camera for Surveys / High Resolution Channel of the Hubble Space Telescope. This object is extremely metal-poor (12 + log(O/H) = 7.45) and resides in a nearby void. The images resolve it into two physically separate components that were previously described as star clusters in a single galaxy. The primary component is only \~100 pc in maximum extent, and consists of starburst region surrounded by a ring-like structure of relatively redder stars. The secondary component is ~50 pc in size and lies at a projected distance of ~80 pc away from the primary, and is also actively star-forming. We estimate masses ~10^7 M(sol) and ~10^6 M(sol) for the two components, based on their luminosities, with an associated dynamical timescale for the system of a few Myr. This timescale and the structure of the components suggests that a collision between them triggered their starbursts. The spectral energy distributions of both components can be fitted by the combination of recent (few Myr old) starburst and an evolved (several Gyr old) underlying stellar population, similar to larger blue compact dwarf galaxies. This indicates that despite its metal deficiency the object is not forming its first generation of stars. However, the small sizes and masses of the two components suggests that HS 0822+3542 represents a dwarf galaxy in the process of assembling from clumps of stars intermediate in size between globular clusters and objects previously classified as galaxies. Its relatively high ratio of neutral gas mass to stellar mass (~1) and high specific star formation rate, log(SFR/M(sol) = -9.2, suggests that it is still converting much of its gas to stars.

We present our cosmic shear analysis of GEMS, one of the largest wide-field surveys ever undertaken by the Hubble Space Telescope. Imaged with the Advanced Camera for Surveys (ACS), GEMS spans 795 square arcmin in the Chandra Deep Field South. We detect weak lensing by large-scale structure in high resolution F606W GEMS data from ~60 resolved galaxies per square arcminute. We measure the two-point shear correlation function, the top-hat shear variance and the shear power spectrum, performing an E/B mode decomposition for each statistic. We show that we are not limited by systematic errors and use our results to place joint constraints on the matter density parameter Omega_m and the amplitude of the matter power spectrum sigma_8. We find sigma_8(Omega_m/0.3)^{0.65}=0.68 +/- 0.13 where the 1sigma error includes both our uncertainty on the median redshift of the survey and sampling variance. Removing image and point spread function (PSF) distortions are crucial to all weak lensing analyses. We therefore include a thorough discussion on the degree of ACS PSF distortion and anisotropy which we characterise directly from GEMS data. Consecutively imaged over 20 days, GEMS data also allows us to investigate PSF instability over time. We find that, even in the relatively short GEMS observing period, the ACS PSF ellipticity varies at the level of a few percent which we account for with a semi-time dependent PSF model. Our correction for the temporal and spatial variability of the PSF is shown to be successful through a series of diagnostic tests.

Using Hubble Space Telescope and ground-based U through K- band photometry from the Great Observatories Origins Deep Survey (GOODS), we measure the evolution of the luminosity function and luminosity density in the rest-frame optical (UBR) to z ~ 2, bridging the poorly explored ``redshift desert'' between z~1 and z~2. We also use deep near-infrared observations to measure the evolution in the rest-frame J-band to z~1. Compared to local measurements from the SDSS, we find a brightening of the characteristic magnitude, (M*), by ~2.1, \~0.8 and ~0.7 mag between z=0.1 and z=1.9, in U, B, and R bands, respectively. The evolution of M* in the J-band is in the opposite sense, showing a dimming between redshifts z=0.4 and z=0.9. This is consistent with a scenario in which the mean star formation rate in galaxies was higher in the past, while the mean stellar mass was lower, in qualitative agreement with hierarchical galaxy formation models. We find that the shape of the luminosity function is strongly dependent on spectral type and that there is strong evolution with redshift in the relative contribution from the different spectral types to the luminosity density. We find good agreement in the luminosity function derived from an R-selected and a K-selected sample at z~1, suggesting that optically selected surveys of similar depth (R < 24) are not missing a significant fraction of objects at this redshift relative to a near-infrared-selected sample. We compare the rest-frame B-band luminosity functions from z~0--2 with the predictions of a semi-analytic hierarchical model of galaxy formation, and find qualitatively good agreement. In particular, the model predicts at least as many optically luminous galaxies at z~1--2 as are implied by our observations.

We combine high-resolution images from GEMS with redshifts and spectral energy distributions from COMBO-17 to explore the morphological types of galaxies that dominate the z~0.7 UV luminosity density. We analysed rest-frame 280 nm and V-band luminosities of 1483 galaxies with 0.65<z< 0.75, combining these with visual morphological classifications from F850LP images (approximately rest-frame V-band) taken with HST/ACS on the Extended Chandra Deep Field South. We derive UV luminosity functions and j_280 luminosity densities for spheroid-dominated galaxies, spiral galaxies, Magellanic irregulars, and clearly-interacting galaxies with morphologies suggestive of ongoing major mergers. We check the reliability of GEMS morphologies against the deeper GOODS images and quantify an incompleteness of the GEMS merger identification at the faint end. We derive the fractions of the global UV luminosity density j_280 originating from the galaxy types, and find that spiral galaxies and Magellanic irregulars dominate with about 40% each. Interacting and merging galaxies account for roughly 20% of j_280, while the contribution of early types is negligible. These results imply that the strong decline in the UV luminosity density of the Universe observed from z~1 until today is dominated by the decreasing UV luminosities of normal spiral galaxies, accompanied by the migration of UV-luminous star formation in irregular galaxies to systems of progressively lower mass and luminosity. (abridged)

A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak luminosities ~90% of the bolometric luminosity is absorbed by dust. In general, the detailed appearance of the merging event depends on the stage of the merger and the geometry of the encounter. The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy formation.

We combine HST imaging from the GEMS survey with photometric redshifts from COMBO-17 to explore the evolution of disk-dominated galaxies since z<1.1. The sample is comprised of all GEMS galaxies with Sersic indices n<2.5, derived from fits to the galaxy images. We account fully for selection effects through careful analysis of image simulations; we are limited by the depth of the redshift and HST data to the study of galaxies with absolute magnitudes M(V)<-20, or equivalently stellar masses log(M)>10. We find strong evolution in the magnitude-size scaling relation for galaxies with M(V)<-20, corresponding to a brightening of 1 mag per sqarcsec in rest-frame V-band by z=1. Yet, disks at a given absolute magnitude are bluer and have lower stellar mass-to-light ratios at z=1 than at the present day. As a result, our findings indicate weak or no evolution in the relation between stellar mass and effective disk size for galaxies with log(M)>10 over the same time interval. This is strongly inconsistent with the most naive theoretical expectation, in which disk size scales in proportion to the halo virial radius, which would predict that disks are a factor of two denser at fixed mass at z=1. The lack of evolution in the stellar mass-size relation is consistent with an ``inside-out'' growth of galaxy disks on average (galaxies increasing in size as they grow more massive), although we cannot rule out more complex evolutionary scenarios.

We discuss how measurements of the absorption of gamma-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. Semi-analytic models (SAMs) of galaxy formation, based on the flat LCDM hierarchical structure formation scenario, are used to make predictions of the EBL from 0.1 to 1000 microns. SAMs incorporate simplified physical treatments of the key processes of galaxy formation -- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production. We will summarize SAM successes and failures in accounting for observations at low and high redshift. New ground- and space-based gamma ray telescopes will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the stellar initial mass function, the history of star formation, and the reprocessing of light by dust. On a separate topic concerning gamma ray cosmology, we discuss a new theoretical insight into the distribution of dark matter at the center of the Milky Way, and its implications concerning the high energy gamma rays observed from the Galactic center.

In the past year, the Multimission Archive at the Space Telescope Science Institute (MAST) has taken major steps in making MAST's holdings available using VO-defined protocols and standards, and in implementing VO-based tools. For example, MAST has implemented the Simple Cone Search protocol, and all MAST mission searches may be returned in the VOTable format, allowing other archives to use MAST data for their VO applications. We have made many of our popular High Level Science Products available through Simple Image Access Protocol (SIAP), and are implementing the VO Simple Spectral Access Protocol (SSAP). The cross correlation of VizieR catalogs with MAST missions is now possible, and illustrates the integration of VO services into MAST. The user can easily display the results from searches within MAST using the plotting tool VOPlot. MAST also participates in the NVO registry service. Thus, the user can harvest MAST holdings simultaneously with data from many other surveys and missions through the VO DataScope Data Inventory Service.

We combine HST/ACS imaging from the GEMS survey with redshifts and rest-frame quantities from COMBO-17 to study the evolution of morphologically early-type galaxies with red colors since z=1. We use a new large sample of 728 galaxies with centrally-concentrated radial profiles (Sersic n>2.5) and rest-frame U-V colors on the red sequence. By appropriate comparison with the local relations from SDSS, we find that the luminosity-size (L-R) and stellar mass-size (M-R) relations evolve in a manner that is consistent with the passive aging of ancient stars. By itself, this result is consistent with a completely passive evolution of the red early-type galaxy population. If instead, as demonstrated by a number of recent surveys, the early-type galaxy population builds up in mass by a factor of 2 since z=1, our results imply that new additions to the early-type galaxy population follow similar L-R and M-R correlations, compared to the older subset of early-type galaxies. Adding early-type galaxies to the red sequence through disk fading appears to be consistent with the data. Through comparison with models, the role of dissipationless merging is limited to <1 major merger on average since z=1 for the most massive galaxies. Predictions from models of gas-rich mergers are not yet mature enough to allow a detailed comparison to our observations. We find tentative evidence that the amount of luminosity evolution depends on galaxy stellar mass, such that the least massive galaxies show stronger luminosity evolution compared to more massive early types. This could reflect a different origin of low-mass early-type galaxies and/or younger stellar populations; the present data is insufficient to discriminate between these possibilities. (abridged)

One third of present-day spirals host optically visible strong bars that drive their dynamical evolution. However, the fundamental question of how bars evolve over cosmological times has yet to be addressed, and even the frequency of bars at intermediate redshifts remains controversial. We investigate the frequency of bars out to z~1.0 drawing on a sample of 1590 galaxies from the GEMS survey, which provides morphologies from HST ACS two-color images, and highly accurate redshifts from the COMBO-17 survey. We identify spiral galaxies using the Sersic index, concentration parameter, and rest-frame color. We characterize bars and disks by fitting ellipses to F606W and F850LP images, taking advantage of the two bands to minimize bandpass shifting. We exclude highly inclined (i>60 deg) galaxies to ensure reliable morphological classifications, and apply completeness cuts of M_v <= -19.3 and -20.6. More than 40% of the bars that we detect have semi major axes a<0.5" and would be easily missed in earlier surveys without the small PSF of ACS. The bars that we can reliably detect are fairly strong (with ellipticities e>=0.4) and have a in the range ~1.2-13 kpc. We find that the optical fraction of such strong bars remains at ~(30% +- 6%) from the present-day out to look-back times of 2-6 Gyr (z~0.2-0.7) and 6-8 Gyr (z~0.7-1.0); it certainly shows no sign of a drastic decline at z>0.7. Our findings of a large and similar bar fraction at these three epochs favor scenarios in which cold gravitationally unstable disks are already in place by z~1, and where on average bars have a long lifetime (well above 2 Gyr). The distributions of structural bar properties in the two slices are, however, not statistically identical and therefore allow for the possibility that the bar strengths and sizes may evolve over time.

Bars drive the dynamical evolution of disk galaxies by redistributing mass and angular momentum, and they are ubiquitous in present-day spirals. Early studies of the Hubble Deep Field reported a dramatic decline in the rest-frame optical bar fraction f_opt to below 5% at redshifts z>0.7, implying that disks at these epochs are fundamentally different from present-day spirals. The GEMS bar project, based on ~8300 galaxies with HST-based morphologies and accurate redshifts over the range 0.2-1.1, aims at constraining the evolution and impact of bars over the last 8 Gyr. We present early results indicating that f_opt remains nearly constant at ~30% over the range z=0.2-1.1,corresponding to lookback times of ~2.5-8 Gyr. The bars detected at z>0.6 are primarily strong with ellipticities of 0.4-0.8. Remarkably, the bar fraction and range of bar sizes observed at z>0.6 appear to be comparable to the values measured in the local Universe for bars of corresponding strengths. Implications for bar evolution models are discussed.

We present the results from a study of the host galaxies of 15 optically selected AGNs with 0.5<z<1.1 from GEMS. GEMS is a Hubble Space Telescope imaging survey of a ~28"x28" contiguous field centered on the Chandra Deep Field South in the F606W and F850LP filter bands. It incorporates the SEDs and redshifts of ~10000 objects, obtained by the COMBO-17 project. We have detected the host galaxies of all 15 AGNs in the F850LP-band (and 13/15 in the F606W-band), recovering their fluxes, morphologies and structural parameters. We find that 80% of the host galaxies have early-type (bulge-dominated) morphologies, while the rest have structures characteristic of late-type (disk-dominated) galaxies. We find that 25% of the early types, and 30% of the late types, exhibit disturbances consistent with galaxy interactions. The hosts show a wide range of colors, from those of red-sequence galaxies to blue colors consistent with ongoing star formation. Roughly 70% of the morphologically early-type hosts have rest-frame blue colors, a much larger fraction than those typical of non-active morphologically early-type galaxies in this redshift and luminosity range. Yet, we find that the early-type hosts are structurally similar to red-sequence ellipticals, inasmuch as they follow an absolute magnitude versus half-light size correlation that are consistent with the mean relation for early-type galaxies at similar redshifts.

We have performed HST imaging of a sample of 23 high-redshift (1.8<z<2.75) AGN, drawn from the COMBO-17 survey. The sample contains moderately luminous quasars (M_B ~ -23). The data are part of the GEMS imaging survey that provides high resolution optical images obtained with ACS in two bands (F606W and F850LP), sampling the rest-frame UV flux of the targets. To deblend the AGN images into nuclear and resolved (host galaxy) components we use a PSF subtraction technique that is strictly conservative with respect to the flux of the host galaxy. We resolve the host galaxies in both filter bands in 9 of the 23 AGN, whereas the remaining 14 objects are considered non-detections, with upper limits of less than 5% of the nuclear flux. However, when we coadd the unresolved AGN images into a single high signal-to-noise composite image we find again an unambiguously resolved host galaxy. The recovered host galaxies have apparent magnitudes of 23.0<F606W<26.0 and 22.5<F850LP<24.5 with rest-frame UV colours in the range -0.2<(F606W-F850LP)_obs<2.3. The rest-frame absolute magnitudes at 200 nm are -20.0<M_200nm<-22.2. The photometric properties of the composite host are consistent with the individual resolved host galaxies. We find that the UV colors of all host galaxies are substantially bluer than expected from an old population of stars with formation redshift z=5, independent of the assumed metallicities. These UV colours and luminosities range up to the values found for Lyman-break galaxies (LBGs) at z=3. Our results agree with the recent discovery of enhanced blue stellar light in AGN hosts at lower redshifts. [abridged]

We have mapped the AGN luminosity function and its evolution between z=1 and z=5 down to apparent magnitudes of $R<24$. Within the GEMS project we have analysed HST-ACS images of many AGN in the Extended Chandra Deep Field South, enabling us to assess the evolution of AGN host galaxy properties with cosmic time.

We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. The non-adiabatic heating, i.e. shocks, generated in the final merger forces gas out of the central region of the merger remnant and into the dark-matter halo. We demonstrate that the amount of heating depends on the size of the progenitor disk galaxy as well as the initial orbit the galaxies are placed on. Based upon these dependencies, we motivate a possible recipe for including this effect in semi-analytic models of galaxy formation.

Pan-chromatic galaxy surveys are providing tightening constraints on the global mass assembly history, and high-resolution imaging of large fields is telling us when and where stars were formed. How well are state-of-the-art hierarchical galaxy formation models currently doing at reproducing these observations? I present results here that suggest that hierarchical models are doing quite well at reproducing the global star formation and stellar mass assembly history obtained from galaxies selected in the optical and Near IR. However, the same models fail to reproduce two very important populations at high redshift: quiescent red spheroids and vigorously star-forming, dust-enshrouded starbursts. This mismatch carries important lessons about how star formation is triggered and regulated in early galaxies, and may force us to consider new ideas about the formation of massive spheroids.

We examine the evolution of the sizes and number densities of disk galaxies using the high resolution images obtained by the Great Observatories Origins Deep Survey (GOODS) with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). The multiwavelength (B,V,i,z) images allow us to classify galaxies based on their rest-frame B-band morphologies out to redshift z = 1.25. In order to minimize the effect of selection biases, we confine our analysis to galaxies which occupy the region of magnitude-size plane where the survey is almost 90% complete at all redshifts. The observed size distribution is consistent with a log-normal distribution as seen for the disk galaxies in the local Universe and does not show any significant evolution over the redshift range 0.25 < z < 1.25. We find that the number densities of disk galaxies remains fairly constant over this redshift range, although a modest evolution by a factor of four may be possible within the 2-sigma uncertainties.

We investigate the constraints on models of supermassive black hole (SMBH) and quasar formation obtainable from two recent observational developments: the discovery of luminous quasars at z~6, and estimates of the local mass density of SMBHs. If ~90 per cent of this mass was accreted at redshifts z<3, as suggested by the observed quasar luminosity functions, these joint constraints pose a challenge for models, which must account for the observed luminous quasar population at z~6 within a very limited `mass budget'. We investigate a class of models based within the hierarchical structure formation scenario, in which major mergers lead to black hole formation and fuelling, and the resulting quasars shine at their Eddington-limited rate until their fuel is exhausted. We show that the simplest such model, in which a constant fraction of the gas within the halo is accreted in each major merger, cannot satisfy both constraints simultaneously. When this model is normalized to reproduce the number density of luminous quasars at z~6, the mass budget is grossly exceeded due to an overabundance of lower mass SMBHs. We explore a range of modifications to the simple model designed to overcome this problem. We show that both constraints can be satisfied if the gas accretion fraction scales as a function of the halo virial velocity. Similar scalings have been proposed in order to reproduce the local M-sigma relation. Successful models can also be constructed by restricting the formation of seed black holes to redshifts above z~11.5 or to haloes above a velocity threshold ~55 km/s, or assuming that only a fraction of major mergers result in formation of a seed SMBH. (abridged)

GEMS, Galaxy Evolution from Morphologies and SEDs, is a large-area (800 arcmin2) two-color (F606W and F850LP) imaging survey with the Advanced Camera for Surveys on HST. Centered on the Chandra Deep Field South, it covers an area of ~28'x28', or about 120 Hubble Deep Field areas, to a depth of m_AB(F606W)=28.3 (5sigma and m_AB(F850LP)=27.1 (5sigma) for compact sources. In its central ~1/4, GEMS incorporates ACS imaging from the GOODS project. Focusing on the redshift range 0.2<=z<=1.1, GEMS provides morphologies and structural parameters for nearly 10,000 galaxies where redshift estimates, luminosities and SEDs exist from COMBO-17. At the same time, GEMS contains detectable host galaxy images for several hundred faint AGN. This paper provides an overview of the science goals, the experiment design, the data reduction and the science analysis plan for GEMS.

We have used the 30'x30' Hubble Space Telescope image mosaic from the Galaxy Evolution from Morphology and SEDs (GEMS) project in conjunction with the COMBO-17 deep photometric redshift survey to define a sample of nearly 1500 galaxies with 0.65 < z < 0.75. With this sample, we can study the distribution of rest-frame V-band morphologies more than 6 Gyr ago, without differential bandpass shifting and surface brightness dimming across this narrow redshift slice. Focusing on red-sequence galaxies at z ~ 0.7, we find that 85% of their combined rest-frame V-band luminosity density comes from visually-classified E/S0/Sa galaxies (at V-band magnitudes brighter than -19.5). Similar results are obtained if automated classifiers are used. This fraction is identical to that found at the present day, and is biased by less than 10% by large scale structure and the morphology-density relation. Under the assumption that peculiar and edge-on disk galaxies are red by virtue of their dust content, we find that less than 13% of the total rest-frame V-band luminosity of the z ~ 0.7 red galaxy population is from dusty galaxies.

The integrated colors of distant galaxies provide a means for interpreting the properties of their stellar content. Here, we use rest-frame UV-to-optical colors to constrain the spectral-energy distributions and stellar populations of color-selected, B-dropout galaxies at z ~ 4 in the Great Observatories Origins Deep Survey. We combine the ACS data with ground-based near-infrared images, which extend the coverage of galaxies at z ~ 4 to the rest-frame B-band. We observe a color-magnitude trend in the rest-frame m(UV) - B versus B diagram for the z ~ 4 galaxies that has a fairly well-defined "blue-envelope", and is strikingly similar to that of color-selected, U-dropout galaxies at z ~ 3. We also find that although the co-moving luminosity density at rest-frame UV wavelengths (1600 Angstroms) is roughly comparable at z ~ 3 and z ~ 4, the luminosity density at rest-frame optical wavelengths increases by about one-third from z ~ 4 to z ~ 3. Although the star-formation histories of individual galaxies may involve complex and stochastic events, the evolution in the global luminosity density of the UV-bright galaxy population corresponds to an average star-formation history with a star-formation rate that is constant or increasing over these redshifts. This suggests that the evolution in the luminosity density corresponds to an increase in the stellar-mass density of more than 33%.

The redshift distribution of well-defined samples of distant early-type galaxies offers a means to test the predictions of monolithic and hierarchical galaxy formation scenarios. NICMOS maps of the entire Hubble Deep Field North in the F110W and F160W filters, when combined with the available WFPC2 data, allow us to calculate photometric redshifts and determine the morphological appearance of galaxies at rest-frame optical wavelengths out to z ~ 2.5. Here we report results for two subsamples of early-type galaxies, defined primarily by their morphologies in the F160W band, which were selected from the NICMOS data down to H160_{AB} < 24.0. The observed redshift distributions of our two early-type samples do not match that predicted by a monolithic collapse model, which shows an overabundance at z > 1.5. A hierarchical formation model better matches the redshift distribution of the HDF-N early-types at z > 1.5, but still does not adequately describe the observed early-types. The hierarchical model predicts significantly bluer colors on average than the observed early-type colors, and underpredicts the observed number of early-types at z < 1. [abridged]

Using U'- through Ks-band imaging data in the GOODS-South field, we construct a large, complete sample of 275 ``extremely red objects'' (EROs; K_s<22.0, R-K_s>3.35; AB), all with deep HST/ACS imaging in B_435, V_606, i_775, and z_850, and well-calibrated photometric redshifts. Quantitative concentration and asymmetry measurements fail to separate EROs into distinct morphological classes. We therefore visually classify the morphologies of all EROs into four broad types: ``Early'' (elliptical-like), ``Late'' (disk galaxies), ``Irregular'' and ``Other'' (chain galaxies and low surface brightness galaxies), and calculate their relative fractions and comoving space densities. For a broad range of limiting magnitudes and color thresholds, the relative number of early-type EROs is approximately constant at 33-44%, and the comoving space densities of Early- and Late-type EROs are comparable. Mean rest-frame spectral energy distributions (SEDs) at wavelengths between 0.1 and 1.2 um are constructed for all EROs. The SEDs are extremely similar in their range of shapes, independent of morphological type. The implication is that any differences between the broad-band SEDs of Early-type EROs and the other types are relatively subtle, and there is no robust way of photometrically distinguishing between different morphological types with usual optical/near-infrared photometry.

Cosmic variance is the uncertainty in observational estimates of the volume density of extragalactic objects such as galaxies or quasars arising from the underlying large-scale density fluctuations. This is often a significant source of uncertainty, especially in deep galaxy surveys, which tend to cover relatively small areas. We present estimates of the relative cosmic variance for one-point statistics (i.e. number densities) for typical scales and volumes sampled by the Great Observatories Origins Deep Survey (GOODS). We use two approaches: for objects with a known two-point correlation function that is well-approximated by a power law, one can use the standard analytic formalism to calculate the cosmic variance (in excess of shot noise). We use this approach to estimate the cosmic variance for several populations that are being studied in the GOODS program: Extremely Red Objects (ERO) at z~1, and Lyman Break Galaxies (LBG) at z~3 and z~4, using clustering information for similar populations in the literature. For populations with unknown clustering, one can use predictions from Cold Dark Matter theory to obtain a rough estimate of the variance as a function of number density. We present a convenient plot which allows one to use this approach to read off the cosmic variance for a population with a known mean redshift and estimated number density. We conclude that for the volumes sampled by GOODS, cosmic variance is a significant source of uncertainty for strongly clustered objects (~40-60% for EROs) and less serious for less clustered objects, ~10-20% for LBGs.

The redshift distribution of near-IR selected galaxies is often used to attempt to discriminate between the classical view of galaxy formation, in which present-day luminous galaxies were assembled at early times and evolve due to the passive aging of their stellar populations, and that of hierarchical structure formation, in which galaxies were assembled more recently via the merging of smaller objects. We carry out such a test here, by computing the distribution of photometric redshifts of K_AB<22 galaxies in the Great Observatories Origins Deep Field Survey GOODS Southern Field, and comparing the results with predictions from a semi-analytic model based on hierarchical structure formation, and a classical `passive evolution' model. We find that the redshift distributions at z <~ 1.5 of both the hierarchical and passive models are very similar to the observed one. At z>~1.5, the hierarchical model shows a deficit of galaxies, while the passive model predicts an excess. We investigate the nature of the observed galaxies in the redshift range where the models diverge, and find that the majority have highly disturbed morphologies, suggesting that they may be merger-induced starbursts. While the hierarchical model used here does not produce these objects in great enough numbers, the appearance of this population is clearly in better qualitative agreement with the hierarchical picture than with the classical passive evolution scenario. We conclude that the observations support the general framework of hierarchical formation, but suggest the need for new or modified physics in the models.

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later.

We investigate several fundamental properties of z ~ 4 Lyman-break galaxies by comparing observations with the predictions of a semi-analytic model based on the Cold Dark Matter theory of hierarchical structure formation. We use a sample of B_{435}-dropouts from the Great Observatories Origins Deep Survey, and complement the ACS optical B_{435}, V_{606}, i_{775}, and z_{850} data with the VLT ISAAC J, H, and K_{s} observations. We extract B_{435}-dropouts from our semi-analytic mock catalog using the same color criteria and magnitude limits that were applied to the observed sample. We find that the i_{775} - K_{s} colors of the model-derived and observed B_{435}-dropouts are in good agreement. However, we find that the i_{775}-z_{850} colors differ significantly, indicating perhaps that either too little dust or an incorrect extinction curve have been used. Motivated by the reasonably good agreement between the model and observed data we present predictions for the stellar masses, star formation rates, and ages for the z ~ 4 Lyman-break sample. We find that according to our model the color selection criteria used to select our z ~ 4 sample surveys 67% of all galaxies at this epoch down to z_{850} < 26.5. We find that our model predicts a roughly 40% mass build-up between the z \~ 4 and z ~ 3 epochs for the UV rest-frame L* galaxies. Furthermore, according to our model, at least 50% of the total stellar mass resides in relatively massive UV-faint objects that fall below our observational detection limit.

We investigate the physical significance of a new spectral parameter, eta. This parameter was defined using a Principal Component Analysis of the 2dF Galaxy Redshift Survey, to retain astrophysical information while minimising the effect of measurement uncertainties. We find that while eta is correlated with morphological type, there is a large scatter in this relationship. A tighter empirical relationship is found between eta and the equivalent width of the Halpha line, suggesting a connection with the relative star formation rate. We pursue this connection using spectral synthesis models. Using models in which the star formation history is parameterised in terms of an exponentially decreasing function of time, we find that there is a tight correlation between eta and the ratio of the present to the past-averaged rate of star-formation, often known as the ``birthrate'' parameter b. This correlation also holds in models with much more complicated star formation histories, generated by a semi- analytic model of galaxy formation based upon the hierarchical formation scenario. There are two possible causes for the tight correlations we find between eta and b in those galaxies with the most complex star formation histories; Firstly, the spectra themselves may be degenerate to the actual long-term star formation history of each galaxy in the optical wavelength range probed by the 2dFGRS. Secondly, the birthrate parameter b may represent a physically fundamental quality of galaxy halos -- their over-density relative to the background density. We conclude that the tight connection with b makes eta a physically meaningful, as well as convenient and robust, statistic for galaxy parameterisation and classification. (abridged)

Reducing the power on small scales relative to the `standard' LCDM model alleviates a number of possible discrepancies with observations, and is favored by the recent analysis of WMAP plus galaxy and Lyman-alpha forest data. Here, we investigate the epoch of reionization in several models normalized to WMAP on large scales, and with sufficiently reduced power on small scales to solve the halo concentration and substructure problems. These include a tilted model, the WMAP running-index model, and a Warm Dark Matter model. We assume that the Universe was reionized by stellar sources comprised of a combination of supermassive (~200 M_sun) Pop III stars and Pop II stars with a `normal' IMF. We find that in all of these models, structure formation and hence reionization occurs late, certainly at redshifts below ten, and more probably at z < 6. This is inconsistent (at 2 sigma) with the determination of z_reion ~ 17 from the WMAP temperature-polarization data and is only marginally consistent with SDSS quasar observations. The tension between the galactic-scale observations, which favor low-power models, and the early reionization favored by WMAP can only be resolved if the efficiency of Pop III star formation is dramatically higher than any current estimate, or if there is an exotic population of ionizing sources such as mini-quasars. Otherwise, we may have to live with the standard LCDM power spectrum, and solve the small-scale problems in some other way.

We calculate the global star formation rate density (SFRD) from z ~ 30-3 using a semi-analytic model incorporating the hierarchical assembly of dark matter halos, gas cooling via atomic hydrogen, star formation, supernova feedback, and suppression of gas collapse in small halos due to the presence of a photoionizing background. We compare the results with the predictions of simpler models based on the rate of dark matter halo growth and a fixed ratio of stellar-to-dark mass, and with observational constraints on the SFRD at 3 < z < 6. We also estimate the star formation rate due to very massive, metal-free Pop III stars using a simple model based on the halo formation rate, calibrated against detailed hydrodynamic simulations of Pop III star formation. We find that the total production rate of hydrogen-ionizing photons during the probable epoch of reionization (15 < z < 20) is approximately equally divided between Pop II and Pop III stars, and that if reionization is late (less than about 15, close to the lower limit of the range allowed by the WMAP results), then Pop II stars alone may be able to reionize the Universe.

We investigate a model for the high-ionization state gas associated with observed damped Lyman-alpha systems, based on a semi-analytic model of galaxy formation set within the paradigm of hierarchical structure formation. In our model, the hot gas in halos and sub-halos gives rise to CIV absorption, while the low-ionization state gas is associated with the cold gas in galaxies. The model matches the distribution of CIV column densities and leads naturally to kinematic properties that are in good agreement with the data. We examine the contribution of both hot and cold gas to sub-damped systems and suggest that the properties of these systems can be used as an important test of the model. We expect that sub-DLA systems will generally be composed of a single gas disk and thus predict that they should have markedly different kinematics than the damped systems. Finally, we find that hot halo gas produces less than one third of Lyman limit systems at redshift three. We model the contribution of mini-halos (halos with virial velocities < 35 km/s) to Lyman limit systems and find that they may contain as much gas as is observed in these systems. However, if we adopt realistic models of the gas density distribution we find that these systems are not a significant source of Lyman limit absorption. Instead we suggest that uncollapsed gas outside of virialized halos is responsible for most of the Lyman limit systems at high redshift.

A new approach to estimating photometric redshifts - using Artificial Neural Networks (ANNs) - is investigated. Unlike the standard template-fitting photometric redshift technique, a large spectroscopically-identified training set is required but, where one is available, ANNs produce photometric redshift accuracies at least as good as and often better than the template-fitting method. The Bayesian priors on the underlying redshift distribution are automatically taken into account. Furthermore, inputs other than galaxy colours - such as morphology, angular size and surface brightness - may be easily incorporated, and their utility assessed. Different ANN architectures are tested on a semi-analytic model galaxy catalogue and the results are compared with the template-fitting method. Finally the method is tested on a sample of ~ 20000 galaxies from the Sloan Digital Sky Survey. The r.m.s. redshift error in the range z < 0.35 is ~ 0.021.

In this paper, I discuss the capabilities and limitations of an 8-10 meter ultraviolet/optical telescope in space, the proposed successor to the Hubble Space Telescope, in the context of galaxy studies. The exquisite spatial resolution and excellent sensitivity of such a facility would open up new possibilities for the study of nearby dwarf galaxies ($z \la 0.5$), and for studying the internal structure and kinematics of more luminous galaxies at high redshift ($z \ga 2$). These applications are of particular importance because they would address areas in which the popular Cold Dark Matter theory is in potential conflict with observations.

We investigate the properties of a galaxy similar to the Milky Way within the context of standard disk formation theory in a LCDM universe. Using the standard assumption that baryons conserve specific angular momentum when they collapse, we conclude that the mean properties of the model galaxies are in good agreement with the Milky Way and other similar spiral galaxies, but the predicted scatter in disk scale lengths may be too large. A model in which half of the initial specific angular momentum is transfered to the dark matter may produce a smaller scatter, if very compact disks are unstable and evolve into spheroids or early type galaxies.

Wide field near-infrared (IR) surveys have revealed a population of galaxies with very red optical$-$IR colors, which have been termed ``Extremely Red Objects'' (EROs). Modeling suggests that such red colors (R-K > 5) could be produced by galaxies at z>~1 with either very old stellar populations or very high dust extinction. Recently it has been discovered that EROs are strongly clustered. Are these objects the high-redshift progenitors of present day giant ellipticals (gEs)? Are they already massive at this epoch? Are they the descendents of the $z\sim3$ Lyman Break Galaxies (LBG), which have also been identified as possible high redshift progenitors of giant ellipticals? We address these questions within the framework of the Cold Dark Matter paradigm using an analytic model that connects the number density and clustering or bias of an observed population with the halo occupation function (the number of observed galaxies per halo of a given mass). We find that EROs reside in massive dark matter halos, with average mass <M> > 1E13/h100 Msun. The occupation function that we derive for EROs is very similar to the one we derive for z=0, L>L* early type galaxies, whereas the occupation function for LBGs is skewed towards much smaller host halo masses (<M> ~ 1E11 to 1E12/h100 Msun. We then use the derived occupation function parameters to explore the possible evolutionary connections between these three populations.

The LCIR Survey, using the Cambridge IR Survey Instrument (CIRSI), reaches H \~ 20-21 over ~1 deg^2. We present results for 744 arcmin^2 in which public UBVRI data exist. We compare optical-IR colours with predictions of a semi-analytic hierarchical model (SAM) and find reasonable agreement. We also determine photometric redshifts, finding a median z of z_m ~ 0.55. We compare N(z) of different spectral types with models, showing that the observations are inconsistent with simple PLE models while the SAM provides a reasonable fit to the total N(z) but underestimates the number of z ~ 1 red spectral types. We also present N(z) for samples of red objects (EROs). We find that EROs with R - H > 4 and H < 20.5 have z_m ~ 1; redder EROs have higher z_m. For 19 < H < 20, EROs with R - H > 4 comprise ~18% of the observed galaxy population while in the SAM they contribute only ~4%. We also determine the angular correlation function w(theta) for magnitude, colour, spectral type and photo-z selected samples and use the estimated N(z) to derive the spatial clustering xi(r). Parametrizing xi(r) by xi(r_c,z)=(r_c/r_*(z))^(-1.8) (r_c comoving), we find that r_*(z) increases by ~1.5-2 times from z = 0 to z ~ 1.2. We interpret this as a selection effect - the galaxies selected at z ~ 1.2 are intrinsically very luminous. When galaxies are selected by absolute magnitude we find no evidence for evolution in r_* over this z range. Extrapolated to z = 0, we find r_*(0) ~ 6.5 h^-1Mpc for red galaxies and r_*(0) ~ 2-4 h^-1Mpc for blue galaxies. We also find that while the angular clustering amplitude of EROs with R - H > 4 or I - H > 3 is up to four times that of the whole galaxy population, the spatial clustering length r_*(z=1) is ~7.5-10.5 h^-1Mpc which is only a factor of ~1.7 times r_* for non-EROs lying in a similar z and luminosity range.

There is accumulating evidence that the faint end of the galaxy luminosity function might be very different in different locations. The luminosity function might be rising in rich clusters and flat or declining in regions of low density. If galaxies form according to the model of hierarchical clustering then there should be many small halos compared to the number of big halos. If this theory is valid then there must be a mechanism that eliminates at least the visible component of galaxies in low density regions. A plausible mechanism is photoionization of the intergalactic medium at a time before the epoch that most dwarf galaxies form in low density regions but after the epoch of formation for similar systems that ultimately end up in rich clusters. The dynamical timescales are found to accommodate this hypothesis in a flat universe with Omega_m < 0.4. If small halos exist but simply cannot be located because they have never become the sites of significant star formation, they still might have dynamical manifestations. These manifestations are hard to identify in normal groups of galaxies because small halos do not make a significant contribution to the global mass budget. However, it could be entertained that there are clusters of halos where there are only small systems, clusters that are at the low mass end of the hierarchical tree. There may be places where only a few small galaxies managed to form, enough for us to identify and use as test probes of the potential. It turns out that such environments might be common. Four probable groups of dwarfs are identified within 5 Mpc and the assumption they are gravitationally bound suggests M/L_B ~ 300 - 1200 M_sun/L_sun, 6 +/- factor 2 times higher than typical values for groups with luminous galaxies.

Data on galaxies at high redshift, identified by the Lyman-break photometric technique, can teach us about how galaxies form and evolve. The stellar masses and other properties of such Lyman break galaxies (LBGs) depend sensitively on the details of star formation. In this paper we consider three different star formation prescriptions, and use semi-analytic methods applied to the now-standard $\Lambda$CDM theory of hierarchical structure formation to show how these assumptions about star formation affect the predicted masses of the stars in these galaxies and the masses of the dark matter halos that host them. We find that, within the rather large uncertainties, recent estimates of the stellar masses of LBGs from multi-color photometry are consistent with the predictions of all three models. However, the estimated stellar masses are more consistent with the predictions of two of the models in which star formation is accelerated at high redshifts $z\gsim3$, and of these models the one in which many of the LBGs are merger-driven starbursts is also more consistent with indications that many high redshift galaxies are gas rich. The clustering properties of LBGs have put some constraints on the masses of their host halos, but due to similarities in the halo occupation of the three models we consider and degeneracies between model parameters, current constraints are not yet sufficient to distinguish between realistic models.

(Abridged) The Las Campanas Infrared Survey, based on broad-band optical and near-infrared photometry, is designed to robustly identify a statistically significant and representative sample of evolved galaxies at redshifts z>1. We have completed an H-band imaging survey over 1.1 square degrees of sky in six separate fields. The average 5 sigma detection limit in a four arcsecond diameter aperture is H ~ 20. Here we describe the design of the survey, the observation strategies, data reduction techniques, and object identification procedures. We present sample near-infrared and optical photometric catalogs for objects identified in two survey fields. We perform object detection in all bandpasses and identify ~ 54,000 galaxies over 1,408 square arcminutes of sky in the two fields. Of these galaxies, ~ 14,000 are detected in the H-band and ~ 2,000 have the colors of evolved galaxies, I - H >3, at z > 1. We find that (1) the differential number counts N(m) for the H-band detected objects has a slope of 0.44 at H < 19 and 0.27 for H > 19. In addition, we find that (2) the differential number counts for the H detected red objects has a slope of 0.85 at H < 20 and 0.32 at H > 20, with a mean surface density ~ 3,000 degree^{-2} mag^{-1} at H=20. Finally, we find that (3) galaxies with red optical to near-IR colors (I-H > 3) constitute ~ 20% of the H detected galaxies at H ~ 21, but only 2% at H = 19. We show that red galaxies are strongly clustered, which results in a strong field to field variation in their surface density. Comparisons of observations and predictions based on various formation scenarios indicate that these red galaxies are consistent with mildly evolving early-type galaxies at z ~ 1, although with a significant amount of on-going star formation as indicated by the large scatter in their V-I colors.

We present results from the Las Campanas Infrared Survey, designed to identify a statistically significant sample of z>=1 galaxies using photometric redshift techniques. Here we summarize the design and strategies of the survey and present the first estimate of the galaxy luminosity function at z>=1 based on H-band selected galaxies identified in our survey. Results of number count studies and luminosity function measurements indicate that most early-type galaxies were already in place by z~1.2 with a modest space density evolution and a mild luminosity evolution over that expected from passive evolution.

We apply standard disk formation theory with adiabatic contraction within cuspy halo models predicted by the standard LCDM cosmology. The resulting models score remarkably well when confronted with the broad range of observational data available for the Milky Way and M31 galaxies, giving a Milky Way virial mass of 1-2x10^12Msun and concentration C=12-17. We consider two types of models, in which: (A) baryons conserve angular momentum and (B) some of the angular momentum of the baryons is transferred to the dark matter. Type-A models produce good agreement with observed rotation curves and obey constraints in the solar neighborhood, but may have too much dark matter in the center to allow a fast rotating bar. The type-B models with angular momentum transport have a slightly more massive disk and less dark matter in the central part, allowing a fast rotating bar to persist. Both classes of models probably have sufficient baryonic mass in the central 3.5kpc to reproduce recent observational values of the optical depth to microlensing events towards the Galactic center. All models require that about 1/2 of all baryons expected inside the virial radius must not be in the disk or bulge. We investigate whether the range of virial masses allowed by our dynamical models is compatible with constraints from the galaxy luminosity function, and find a range of parameter space that is allowed by this constraint. We conclude that rotation curves and dynamical properties of ``normal'' high surface brightness spiral galaxies appear to be consistent with standard LCDM.

(Abridged) We have identified a population of faint red galaxies from a 0.62 square degree region of the Las Campanas Infrared Survey whose properties are consistent with their being the progenitors of early-type galaxies. The optical and IR colors, number-magnitude relation and angular clustering together indicate modest evolution and increased star formation rates among the early-type field population at redshifts between one and two. The counts of red galaxies with $H$ magnitudes between 17 and 20 rise with a slope that is much steeper than that of the total H sample. The surface density of red galaxies drops from roughly 3000 per square degree at H = 20.5, I-H > 3 to ~ 20 per square degree at H = 20, I-H > 5. The V-I colors are approximately 1.5 magnitudes bluer on average than a pure old population and span a range of more than three magnitudes. The colors, and photometric redshifts derived from them, indicate that the red galaxies have redshift distributions adequately described by Gaussians with sigma_z ~ 0.2$ centered near redshift one, with the exception that galaxies having $V-I<1.6$ and $I-H>3$ are primarily in the 1.5 < z < 2 range. We find co-moving correlation lengths of 9-10 Mpc at z ~ 1, comparable to, or larger than, those found for early-type galaxies at lower redshifts. A simple photometric evolution model reproduces the counts of the red galaxies, with only a ~ 30% decline in the underlying space density of early-type galaxies at z ~ 1.2. We suggest on the basis of the colors, counts, and clustering that these red galaxies are the bulk of the progenitors of present day early-type galaxies.

Cold Dark Matter theory predicts that the Local Group should contain many more dwarf-sized objects than the observed number of dwarf galaxies --- the so-called sub-structure problem. We investigate whether the suppression of star formation in these small objects due to the presence of a photoionizing background can resolve the problem. We make use of results from recent hydrodynamic simulations to build a recipe for the suppression of gas infall into semi-analytic galaxy formation models, and use these to predict the luminosity function of dwarf galaxies in the Local Group. In the models without photoionization ``squelching'', we predict a large excess of faint dwarf galaxies compared with the observed number in the Local Group --- thus, the usual recipe for supernovae feedback used in semi-analytic models does not solve the sub-structure problem on its own. When we include photoionization squelching, we find good agreement with the observations. We have neglected tidal destruction, which probably further reduces the number of dwarf galaxies. We conclude that photoionizing squelching easily solves the sub-structure problem. In fact, it is likely that once this effect is taken into account, models with reduced small-scale power (e.g. Warm Dark Matter) would underproduce dwarf galaxies.

We discuss how current and future data on the clustering and number density of z~3 Lyman-break galaxies (LBGs) can be used to constrain their relationship to dark matter haloes. We explore a three-parameter model in which the number of LBGs per dark halo scales like a power-law in the halo mass: N(M) = (M/M_1)^S for M>M_m. Here, M_m is the minimum mass halo that can host an LBG, M_1 is a normalization parameter, associated with the mass above which haloes host more than one observed LBG, and S determines the strength of the mass dependence. We show how these three parameters are constrained by three observable properties of LBGs: the number density, the large-scale bias, and the fraction of objects in close pairs. Given these three quantities, the three unknown model parameters may be estimated analytically, allowing a full exploration of parameter space. As an example, we assume an LCDM cosmology and consider the observed properties of a recent sample of spectroscopically confirmed LBGs. We find that the favored range for our model parameters is M_m = (0.4-8)x10^10 Msun/h, M_1 = (6-10)x10^12 Msun/h, and 0.9 < S < 1.1. The preferred region in M_m expands by an order of magnitude and slightly shallower slopes are acceptable if the allowed range of bias is permitted to span all recent observational estimates. We also discuss how the observed clustering of LBGs as a function of luminosity can be used to constrain halo occupation, although due to current observational uncertainties we are unable to reach any strong conclusions. Our methods and results can be used to constrain more realistic models that aim to derive the occupation function N(M) from first principles, and offer insight into how basic physical properties affect the observed properties of LBGs.

We model the acquisition of spin by dark-matter halos in semi-analytic merger trees. We explore two different algorithms; one in which halo spin is acquired from the orbital angular momentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells of material while still in the linear regime. We find that both scenarios produce the characteristic spin distribution of halos found in N-body simulations, namely, a log-normal distribution with mean ~0.04 and standard deviation ~0.5 in the log. A perfect match requires fine-tuning of two free parameters. Both algorithms also reproduce the general insensitivity of the spin distribution to halo mass, redshift and cosmology seen in N-body simulations. The spin distribution can be made strictly constant by physically motivated scalings of the free parameters. In addition, both schemes predict that halos which have had recent major mergers have systematically larger spin values. These algorithms can be implemented within semi-analytic models of galaxy formation based on merger trees. They yield detailed predictions of galaxy properties that strongly depend on angular momentum (such as size and surface brightness) as a function of merger history and environment.

We use N-body simulations combined with semi-analytic models to compute the clustering properties of modeled galaxies at z~3, and confront these predictions with the clustering properties of the observed population of Lyman-break galaxies (LBGs). Several scenarios for the nature of LBGs are explored, which may be broadly categorized into models in which high-redshift star formation is driven by collisional starbursts and those in which quiescent star formation dominates. For each model, we make predictions for the LBG overdensity distribution, the variance of counts-in-cells, the correlation length, and close pair statistics. Models which assume a one-to-one relationship between massive dark-matter halos and galaxies are disfavored by close pair statistics, as are models in which colliding halos are associated with galaxies in a simplified way. However, when modeling of gas consumption and star formation is included using a semi-analytic treatment, the quiescent and collisional starburst models predict similar clustering properties and none of these models can be ruled out based on the available clustering data. None of the ``realistic'' models predict a strong dependence of clustering amplitude on the luminosity threshold of the sample, in apparent conflict with some observational results.

We discuss how measurements of the absorption of $\gamma$-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. We use semi-analytic models (SAMs) of galaxy formation, set within the hierarchical structure formation scenario, to obtain predictions of the EBL for 0.1-1000$\mu$m. SAMs incorporate simplified physical treatments of the key processes of galaxy formation --- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production --- and have been shown to reproduce key observations at low and high redshift. Here we also introduce improved modelling of the spectral energy distributions in the mid-to-far-IR arising from emission by dust grains. Assuming a flat \lcdm cosmology with $\Omega_m=0.3$ and Hubble parameter $h=0.65$, we investigate the consequences of variations in input assumptions such as the stellar initial mass function (IMF) and the efficiency of converting cold gas into stars. We conclude that observational studies of the absorption of $\gamma$-rays with energies from 10s of Gev to 10s of TeV will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the IMF, the history of star formation, and the reprocessing of light by dust.

We assess the feasibility of detecting star formation in damped Lyman-alpha systems (DLAs) at z>1 through near-infrared spectroscopy using the forthcoming integral field units on 8m-class telescopes. Although their relation to galaxies is not well established, high-z DLAs contain most of the neutral gas in the Universe, and this reservoir is depleted with time - presumably through star formation. Line emission should be an indicator of star formation activity, but searches based on Lyman-alpha are unreliable because of the selective extinction of this resonant UV line. Using more robust lines such as H-alpha forces a move to the near-infrared at z>1. For line emission searches, spectroscopy is more sensitive than imaging, but previous long-slit spectroscopic searches have been hampered by the likelihood that any star forming region in the DLA galaxy disk would fall outside the narrow slit. The new integral field units such as CIRPASS on Gemini will cover sufficient solid angles to intercept these, even in the extreme case of large galactic disks at high redshift. On an 8m-class telescope, star formation rates of <1M_sun/yr will be reached at z~1.4 with H-alpha in the H-band. Such star formation rates are well below L* for the high-z Lyman-break population, and are comparable locally to the luminous giant HII complexes in M101. It appears that instruments such as CIRPASS on Gemini will have both the sensitivity and the survey area to measure star formation rates in z>1 DLAs. These observations will probe the nature of damped Lyman-alpha systems and address their relation to galaxies.

Do disk galaxies exist at redshifts much greater than unity, and how might they look different from local disks? How does the morphological mix of galaxies change with redshift? What can we learn from current observations about the properties of high redshift galactic disks? I present theoretical predictions based on semi-analytic hierarchical models, focussing on the role played by interactions, mergers, and starbursts in determining the observable properties of disk galaxies at high redshift, and discuss the interpretation of high redshift observations of possible proto-disks (damped Lyman-$\alpha$ systems and Lyman-break galaxies) in light of these predictions.

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

This paper has been withdrawn because the models on which it was based have undergone significant changes and improvements. A new paper with the same title, based on the improved models, is accepted for publication in MNRAS and is available at astro-ph/0006364. We have withdrawn this paper to avoid confusion.

Using semi-analytic models of galaxy formation set within the Cold Dark Matter (CDM) merging hierarchy, we investigate several scenarios for the nature of the high-redshift ($z \ga 2$) Lyman-break galaxies (LBGs). We consider a ``collisional starburst'' model in which bursts of star formation are triggered by galaxy-galaxy mergers, and find that a significant fraction of LBGs are predicted to be starbursts. This model reproduces the observed comoving number density of bright LBGs as a function of redshift and the observed luminosity function at $z\sim 3$ and $z\sim 4$, with a reasonable amount of dust extinction. Model galaxies at $z \sim 3$ have star formation rates, half-light radii, I-K colours, and internal velocity dispersions that are in good agreement with the data. Global quantities such as the star formation rate density and cold gas and metal content of the Universe as a function of redshift also agree well. Two ``quiescent'' models without starbursts are also investigated. In one, the star formation efficiency in galaxies remains constant with redshift, while in the other, it scales inversely with disc dynamical time, and thus increases rapidly with redshift. The first quiescent model is strongly ruled out as it does not produce enough high redshift galaxies once realistic dust extinction is accounted for. The second quiescent model fits marginally, but underproduces cold gas and very bright galaxies at high redshift. A general conclusion is that star formation at high redshift must be more efficient than locally. The collisional starburst model appears to accomplish this naturally without violating other observational constraints.

We present a method for constraining the evolution of the galaxy luminosity-velocity (LV) relation in hierarchical scenarios of structure formation. The comoving number density of dark-matter halos with circular velocity of 200 km/s is predicted in favored CDM cosmologies to be nearly constant over the redshift range 0<z<5. Any observed evolution in the density of bright galaxies implies in turn a corresponding evolution in the LV relation. We consider several possible forms of evolution for the zero-point of the LV relation and predict the corresponding evolution in galaxy number density. The Hubble Deep Field suggests a large deficit of bright (M_V < -19) galaxies at 1.4 < z < 2. If taken at face value, this implies a dimming of the LV zero-point by roughly 2 magnitudes. Deep, wide-field, near-IR selected surveys will provide more secure measurements to compare with our predictions.

A new method for classification of galaxy spectra is presented, based on a recently introduced information theoretical principle, the `Information Bottleneck'. For any desired number of classes, galaxies are classified such that the information content about the spectra is maximally preserved. The result is classes of galaxies with similar spectra, where the similarity is determined via a measure of information. We apply our method to approximately 6000 galaxy spectra from the ongoing 2dF redshift survey, and a mock-2dF catalogue produced by a Cold Dark Matter-based semi-analytic model of galaxy formation. We find a good match between the mean spectra of the classes found in the data and in the models. For the mock catalogue, we find that the classes produced by our algorithm form an intuitively sensible sequence in terms of physical properties such as colour, star formation activity, morphology, and internal velocity dispersion. We also show the correlation of the classes with the projections resulting from a Principal Component Analysis.

We investigate the properties of damped Lyman-alpha systems (DLAS) in semi-analytic models, focusing on whether the models can reproduce the kinematic properties of low-ionization metal lines described by Prochaska & Wolfe (1997b, 1998). We explore a variety of approaches for modelling the radial distribution of the cold neutral gas associated with the galaxies in our models, and find that our results are very sensitive to this ingredient. If we use an approach based on Fall & Efstathiou (1980), in which the sizes of the discs are determined by conservation of angular momentum, we find that the majority of the DLAS correspond to a single galactic disc. These models generically fail to reproduce the observed distribution of velocity widths. In alternative models in which the gas discs are considerably more extended, a significant fraction of DLAS arise from lines of sight intersecting multiple gas discs in a common halo. These models produce kinematics that fit the observational data, and also seem to agree well with the results of recent hydrodynamical simulations. Thus we conclude that Cold Dark Matter based models of galaxy formation can be reconciled with the kinematic data, but only at the expense of the standard assumption that DLAS are produced by rotationally supported gas discs whose sizes are determined by conservation of angular momentum. We suggest that the distribution of cold gas at high redshift may be dominated by another process, such as tidal streaming due to mergers.

Previously we have shown that it is possible to account for the kinematic properties of damped Lyman alpha systems (DLAS) in the context of semi-analytic models. In these models, hierarchical structure formation is approximated by constructing a merger tree for each dark matter halo. A natural consequence is that every virialized halo may contain not only a central galaxy, but also a number of satellite galaxies as determined by its merging history. Thus the kinematics of the DLAS arise from the combined effects of the internal rotation of gas disks and the motions between gas disks within a common halo. Here we investigate the sensitivity of this model to some of the assumptions made previously, including the modeling of satellite dynamics, the scale height of the gas, and the cosmology.

Previously we have shown (Maller et al, 1998) that the kinematics of Damped Lyman Alpha Systems (DLAS) as measured by Prochaska and Wolfe (1998) can be reproduced in a multiple disk model (MDM) if the gaseous disks are of sufficient radial extent. Here we discuss this model's predictions for the relationship between DLAS and Lyman break galaxies (LBGs), which we here take to be objects at z~3 brighter than R=25.5. We expect that future observations of the correlations between DLAS and LBGs will provide a new data set able to discriminate between different theoretical models of the DLAS. Djorgovski (1997) has already detected a few optical counterparts and more studies are underway.

We study the biasing relation between dark-matter halos or galaxies and the underlying mass distribution, using cosmological $N$-body simulations in which galaxies are modelled via semi-analytic recipes. The nonlinear, stochastic biasing is quantified in terms of the mean biasing function and the scatter about it as a function of time, scale and object properties. The biasing of galaxies and halos shows a general similarity and a characteristic shape, with no galaxies in deep voids and a steep slope in moderately underdense regions. At $\sim 8\hmpc$, the nonlinearity is typically $\lsim 10$ percent and the stochasticity is a few tens of percent, corresponding to $\sim 30$ percent variations in the cosmological parameter $\beta=\Omega^{0.6}/b$. Biasing depends weakly on halo mass, galaxy luminosity, and scale. The time evolution is rapid, with the mean biasing larger by a factor of a few at $z\sim 3$ compared to $z=0$, and with a minimum for the nonlinearity and stochasticity at an intermediate redshift. Biasing today is a weak function of the cosmological model, reflecting the weak dependence on the power-spectrum shape, but the time evolution is more cosmology-dependent, relecting the effect of the growth rate. We provide predictions for the relative biasing of galaxies of different type and color, to be compared with upcoming large redshift surveys. Analytic models in which the number of objects is conserved underestimate the evolution of biasing, while models that explicitly account for merging provide a good description of the biasing of halos and its evolution, suggesting that merging is a crucial element in the evolution of biasing.

We present V, V-I color-magnitude diagrams (CMDs) for three old star clusters in the Large Magellanic Cloud (LMC): NGC 1466, NGC 2257 and Hodge 11. Our data extend about 3 magnitudes below the main-sequence turnoff, allowing us to determine accurate relative ages and the blue straggler frequencies. Based on a differential comparison of the CMDs, any age difference between the three LMC clusters is less than 1.5 Gyr. Comparing their CMDs to those of M 92 and M 3, the LMC clusters, unless their published metallicities are significantly in error, are the same age as the old Galactic globulars. The similar ages to Galactic globulars are shown to be consistent with hierarchial clustering models of galaxy formation. The blue straggler frequencies are also similar to those of Galactic globular clusters. We derive a true distance modulus to the LMC of (m-M)=18.46 +/- 0.09 (assuming (m-M)=14.61 for M 92) using these three LMC clusters.

We investigate the clustering properties of high-redshift galaxies within three competing scenarios for assigning luminous galaxies to dark matter halos from N-body simulations: a one galaxy per massive halo model, a quiescent star formation model, and a collisional starburst model. We compare these models to observations of Lyman-Break galaxies at z~3$ With current data and the simple statistic used here, one cannot rule out any of these models, but we see potential for finding distinguishing features using statistics that are sensitive to the tails of the distribution, and statistics based on the number of multiple galaxies per halo, which we explore in an ongoing study.

Using high resolution N-body simulations with hydrodynamics and star formation, we investigate interactions and the resulting starbursts in galaxies with properties typical of $z\sim 3$. We apply spectral population models to produce mock-HST images, and discuss the observed magnitude, color, and morphological appearance of our simulated galaxies in both the rest-UV and rest-visual bands.

Using semi-analytic models of galaxy formation, we investigate galaxy properties such as the Tully-Fisher relation, the B and K-band luminosity functions, cold gas contents, sizes, metallicities, and colours, and compare our results with observations of local galaxies. We investigate several different recipes for star formation and supernova feedback, including choices that are similar to the treatment in Kauffmann, White & Guiderdoni (1993) and Cole et al. (1994) as well as some new recipes. We obtain good agreement with all of the key local observations mentioned above. In particular, in our best models, we simultaneously produce good agreement with both the observed B and K-band luminosity functions and the I-band Tully-Fisher relation. Improved cooling and supernova feedback modelling, inclusion of dust extinction, and an improved Press-Schechter model all contribute to this success. We present results for several variants of the CDM family of cosmologies, and find that models with values of $\Omega_0 \simeq 0.3$--0.5 give the best agreement with observations.

Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $\Lambda$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.

We present here the extragalactic background light (EBL) predicted by semi-analytic models of galaxy formation, and show how measurements of the absorption of gamma rays of $\sim$ TeV energies via pair production on the EBL can probe cosmology and the formation of galaxies. Semi-analytic models permit a physical treatment of the key processes of galaxy formation -- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production -- and have been shown to reproduce key observations at low and high redshift. Using this approach, we investigate the consequences of variations in input assumptions such as the stellar initial mass function and the underlying cosmology. We conclude that observational studies of the absorption of $\sim 10^{-2}-10^{2}$ TeV gamma rays will help to constrain the star formation history of the universe, and the nature and extent of the extinction of starlight due to dust and reradiation of the absorbed energy at infrared wavelengths.

The velocity profiles of weak metal absorption lines can be used to observationally probe the kinematic state of gas in damped Lyman-alpha systems. Prochaska and Wolfe have argued that the flat distribution of velocity widths combined with the asymmetric line profiles indicate that the DLAS are disks with large rotation velocities. An alternative explanation has been proposed by Haehnelt, Steinmetz, and Rauch, in which the observed large velocity widths and asymmetric profiles can be produced by lines of sight passing through two or more clumps each having relatively small internal velocity dispersions. We investigate the plausibility of this scenario in the context of semi-analytic models based on hierarchical merging trees and including simple treatments of gas dynamics, star formation, supernova feedback, and chemical evolution. We find that all the observed properties of the metal-line systems including the distribution of velocities and the asymmetric profiles, can be reproduced by lines of sight passing through sub-clumps that are bound within larger virialized dark matter halos. In order to produce enough multiple hits, we find that the cold gas must be considerably more extended than the optical radius of the proto-galaxies, perhaps even beyond the tidal radius of the sub-halo. This could occur due to tidal stripping or supernova-driven outflows.

We review our analysis of the clustering properties of ``Lyman-break'' galaxies (LBGs) at redshift z~3, previously discussed in Wechsler et al (1998). We examine the likelihood of spikes found by Steidel et al (1998) in the redshift distribution of LBGs, within a suite of models for the evolution of structure in the Universe. Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBGs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. The LBG correlation functions are less steep than galaxies today (gamma~1.4), but show similar or slightly longer correlation lengths. We have extened this analysis and include a preliminary comparison to the new data presented in Adelberger et al (1998). We also discuss work in progress, in which we use semi-analytic models to identify Lyman-break galaxies within dark-matter halos.

Gamma rays of ~TeV energies from distant sources suffer attenuation due to pair production off of ~1 micron EBL photons. We may exploit this process in order to indirectly measure the EBL and constrain models of galaxy formation. Here, using semi-analytic models of galaxy formation, we examine how gamma ray absorption may be used as an indirect probe of the stellar initial mass function (IMF), although there is a degeneracy with dust modeling. We point out that with the new generation of gamma ray telescopes including STACEE, MAGIC, HESS, VERITAS, and Milagro, we should soon possess a wealth of new data and a new method for probing the nature of the IMF.

Observations now probe the star formation history of the Universe back to a redshift of $z\sim5$. We investigate whether the predictions of semi-analytic models of galaxy formation based on hierarchical Cold Dark Matter (CDM) type models are in agreement with these direct observations and also with the ``fossil'' evidence contained in constraints on the ages of present day early-type galaxies. Previous models predicted that the star formation rate density falls off rather steeply at $z\ga 2$, and correspondingly that the majority of the stars in the Universe formed at relatively low redshift. We investigate the effect of including a bursting mode of star formation, assuming that galaxy-galaxy mergers trigger starbursts and using the merger rate that arises naturally in the CDM merging hierarchy. The resulting starbursts substantially increase the global star formation rate at high redshift, leading to predictions that are in good agreement with the star formation rate density at $z\sim3$ obtained from sub-millimeter observations (SCUBA) and optical/UV estimates after correction for dust extinction. The mass of stars formed at $z \ge 3$ is correspondingly in better agreement with the fossil evidence. We also investigate complementary global quantities such as the mass of cold gas and the average metallicity of cold gas as a function of redshift, and the integrated extra-galactic background light. We argue that these observations, taken together, provide strong constraints on the star formation history of the Universe, and that hierarchical models of the CDM type are in reasonable agreement with these observations when starbursts are included.

We evaluate the accuracy of semi-analytic merger-trees by comparing them with the merging histories of dark-matter halos in N-body simulations, focusing on the joint distribution of the number of progenitors and their masses. We first confirm that the halo mass function as predicted directly by the Press-Schechter (PS) model deviates from the simulations by up to 50% depending on the mass scale and redshift, while the means of the projected distributions of progenitor number and mass for a halo of a given mass are more accurately predicted by the Extended PS model. We then use the full merger trees to study the joint distribution as a function of redshift and parent-halo mass. We find that while the deviation of the mean quantities due to the inaccuracy of the Extended PS model partly propagates into the higher moments of the distribution, the merger-tree procedure does not introduce a significant additional source of error. In particular, certain properties of the merging history such as the mass ratio of the progenitors and the total accretion rate are reproduced quite accurately for galaxy sized halos ($\sim 10^{12}\msun$), and less so for larger masses. We conclude that although there could be $\sim 50%$ deviations in the absolute numbers and masses of progenitors and in the higher order moment of these distributions, the relative properties of progenitors for a given halo are reproduced fairly well by the merger trees. They can thus provide a useful framework for modelling galaxy formation once the above-mentioned limitations are taken into account.

Using semi-analytic models of galaxy formation, we investigate the properties of $z\sim3$ galaxies and compare them with the observed population of Lyman-break galaxies (LBGs). In addition to the usual quiescent mode of star formation, we introduce a physical model for starbursts triggered by galaxy-galaxy interactions. We find that with the merger rate that arises naturally in the CDM-based merging hierarchy, a significant fraction of bright galaxies identified at high redshift ($z \ga 2$) are likely to be low-mass, bursting satellite galaxies. The abundance of LBGs as a function of redshift and the luminosity function of LBGs both appear to be in better agreement with the data when the starburst mode is included, especially when effects of dust extinction are considered. The objects that we identify as LBGs have observable properties including low velocity dispersions that are in good agreement with the available data. In this ``Bursting Satellite'' scenario, quiescent star formation at $z\ga2$ is relatively inefficient and most of the observed LBGs are starbursts triggered by satellite mergers within massive halos. In high-resolution N-body simulations, we find that the most massive dark matter halos cluster at redshift $z\sim 3$ much as the LBGs are observed to do. This is true for both the $\Omega=1$ CHDM model and low-$\Omega$ \LCDM and OCDM models, all of which have fluctuation power spectra $P(k)$ consistent with the distribution of low-redshift galaxies. The Bursting Satellite scenario can resolve the apparent paradox of LBGs that cluster like massive dark matter halos but have narrow linewidths and small stellar masses.

We present two matched sets of five simulations each, covering five presently favored simple modifications to the standard cold dark matter (CDM) scenario. One simulation suite, with a linear box size of 75 Mpc/h, is designed for high resolution and good statistics on the group/poor cluster scale, and the other, with a box size of 300 Mpc/h, is designed for good rich cluster statistics. All runs had 57 million cold particles, and models with massive neutrinos had an additional 113 million hot particles. We consider separately models with massive neutrinos, tilt, curvature, and a nonzero cosmological constant in addition to the standard CDM model. We find that our tilted Omega+Omega_Lambda=1 (TLCDM) model produces too much small-scale power by a factor of ~3, and our open Lambda=0 (OCDM) model also exceeds observed small-scale power by a factor of 2. In addition, we take advantage of the large dynamic range in detectable halo masses our simulations allow to check the shape of the Press-Schechter approximation. We find good fits at cluster masses for delta_c=1.27--1.35 for a Gaussian filter and delta_c=1.57--1.73 for a tophat filter. However, Press-Schechter overpredicts the number density of halos compared to the simulations in the high resolution suite by a weakly cosmology-dependent factor of 1.5--2 at galaxy and group masses, which cannot be fixed by adjusting delta_c within reasonable bounds. An appendix generalizes the spherical collapse model to any isotropic cosmology.

We investigate several approaches for constructing Monte Carlo realizations of the merging history of virialized dark matter halos (``merger trees'') using the extended Press-Schechter formalism. We describe several unsuccessful methods in order to illustrate some of the difficult aspects of this problem. We develop a practical method that leads to the reconstruction of mean quantities such as the conditional and overall mass functions as given by the Press-Schechter model. This method is convenient, computationally efficient, and works for any power spectrum or background cosmology. In addition, we investigate statistics that describe the distribution of the number of progenitors and their masses as a function of redshift.

The combination of measurements of the local abundance of rich clusters of galaxies and its evolution to higher redshift offers the possibility of a direct measurement of $\Omega_0$ with little contribution from other cosmological parameters. We investigate the significance of recent claims that this evolution indicates that $\Omega_0$ must be small. The most recent cluster velocity dispersion function from a compilation including the ESO Northern Abell Cluster Survey (ENACS) results in a significantly higher normalization for models, corresponding to $\sigma_8\approx 0.6$ for $\Omega_0=1$, compared to the Eke, Cole, & Frenk result of $\sigma_8=0.52\pm 0.04$. Using the ENACS data for a $z=0$ calibration results in strong evolution in the abundance of clusters, and we find that the velocity dispersion function is consistent with $\Omega_0=1$. The results are dependent upon the choice and analysis of low-redshift and high-redshift data, so at present, the data is not good enough to determine $\Omega_0$ unambiguously.

We discuss variants of Cold Dark Matter (CDM) dominated cosmological models that give good agreement with a range of observations. We consider models with hot dark matter, tilt, $\Omega < 1$, or a cosmological constant. We also discuss the sensitivity of the results to other parameters, such as the Hubble parameter and the baryon fraction. We obtain constraints by combining the COBE data, cluster abundances, abundance of damped Lyman-$\alpha$ systems at $z\sim3$, the small-angle Cosmic Microwave Background anisotropy, and the small-scale non-linear power spectrum. We present non-linear power spectra from a new suite of N-body simulations for the ``best-bet'' models from each category.

The velocity dispersion of galaxies on small scales ($r\sim1h^{-1}$ Mpc), $\sigma_{12}(r)$, can be estimated from the anisotropy of the galaxy-galaxy correlation function in redshift space. We apply this technique to ``mock-catalogs'' extracted from N-body simulations of several different variants of Cold Dark Matter dominated cosmological models to obtain results which may be consistently compared to similar results from observations. We find a large variation in the value of $\sigma_{12}(1 h^{-1} Mpc)$ in different regions of the same simulation. We conclude that this statistic should not be considered to conclusively rule out any of the cosmological models we have studied. We attempt to make the statistic more robust by removing clusters from the simulations using an automated cluster-removing routine, but this appears to reduce the discriminatory power of the statistic. However, studying $\sigma_{12}$ as clusters with different internal velocity dispersions are removed leads to interesting information about the amount of power on cluster and subcluster scales. We also compute the pairwise velocity dispersion directly and compare this to the values obtained using the Davis-Peebles method, and find that the agreement is fairly good. We evaluate the models used for the mean streaming velocity and the pairwise peculiar velocity distribution in the original Davis-Peebles method by comparing the models with the results from the simulations.

The velocity dispersion of galaxies on scales of $r\sim1h^{-1}$ Mpc, $\sigma_{12}(r)$, may be estimated from the anisotropy of the galaxy-galaxy correlation function in redshift space. We present a reanalysis of the CfA1 survey, correct an error in the original analysis of Davis and Peebles (1983), and find that $\sigma_{12}(r)$ is extremely sensitive to the details of how corrections for infall into the Virgo cluster are applied. We conclude that a robust value of $\sigma_{12}$ cannot be obtained from this survey. We also discuss results from other redshift surveys, including the effect of removing clusters.