We study the incidence of nuclear activity in a large sample of massive post-starburst galaxies at z~0.7 selected from the Sloan Digital Sky Survey, and identify active galactic nuclei based on radio continuum and optical emission lines. Over our mass range of 10^10.6-10^11.5 Msun, the incidence of radio activity is weakly dependent on stellar mass and independent of stellar age, while radio luminosity depends strongly on stellar mass. Optical nuclear activity incidence depends most strongly on the Dn4000 line index, a proxy for stellar age, with an active fraction that is ~ten times higher in the youngest versus oldest post-starburst galaxies. Since a similar trend is seen between age and molecular gas fractions, we argue that, like in local galaxies, the age trend reflects a peak in available fueling rather than feedback from the central black hole on the surrounding galaxy.

We use the FIRE-2 cosmological simulations to study the formation of a virial temperature, quasi-static gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that virialization of the inner CGM coincides with abrupt changes in the properties of the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results hence suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed.

High-redshift Lyman-alpha blobs (LABs) are an enigmatic class of objects that have been the subject of numerous observational and theoretical investigations. It is of particular interest to determine the dominant power sources for the copious luminosity, as direct emission from HII regions, cooling gas, and fluorescence due to the presence of active galactic nuclei (AGN) can all contribute significantly. In this paper, we present the first theoretical model to consider all of these physical processes in an attempt to develop an evolutionary model for the origin of high-z LABs. This is achieved by combining a series of high-resolution cosmological zoom-in simulations with ionization and Lyman-alpha (Lya) radiative transfer models. We find that massive galaxies display a range of Lya luminosities and spatial extents (which strongly depend on the limiting surface brightness used) over the course of their lives, though regularly exhibit luminosities and sizes consistent with observed LABs. The model LABs are typically powered from a combination of recombination in star-forming galaxies, as well as cooling emission from gas associated with accretion. When AGN are included in the model, the fluorescence caused by AGN-driven ionization can be a significant contributor to the total Lya luminosity as well. We propose that the presence of an AGN may be predicted from the Gini coefficient of the blob's surface brightness. Within our modeled mass range, there are no obvious threshold physical properties that predict appearance of LABs, and only weak correlations of the luminosity with the physical properties of the host galaxy. This is because the emergent Lya luminosity from a system is a complex function of the gas temperature, ionization state, and Lya escape fraction.

We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds.

The flat star formation (SF) history of the Milky Way (MW) requires gas in the Galactic disc to be replenished, most likely from a reservoir outside the Galaxy. Such a replenishment may be achieved by a form of `positive' feedback, whereby SF feedback creates a Galactic fountain cycle that collects and cools additional gas from the hot halo surrounding the Galaxy. In this paper we present a model of this process for the MW. A section of the Galactic disc is allowed to form stars which subsequently explode as supernovae (SNe) and send gas out into the hot halo. The gas that is sent out is colder than the hot halo gas and, as it mixes, the halo gas is cooled, providing fuel for further SF as the mixture falls back onto the Galactic disc. We find that this process can be sufficient to maintain a roughly-constant cold gas mass in the MW over at least 3 Gyr. Our results further suggest that there is a positive feedback trend whereby increasing SF leads to an increase in the cold gas budget at average SF rates below 0.5 solar masses per year but above which becomes negative, where further increasing the SFR causes the cold gas budget to decrease. We have constructed an analytical model for this that reproduces the data well and could have profound implications for galaxy evolution in feedback-dominated regimes.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

Dark matter haloes play a fundamental role in cosmological structure formation. The most common approach to model their assembly mechanisms is through N-body simulations. In this work we present an innovative pathway to predict dark matter halo formation from the initial density field using a Deep Learning algorithm. We implement and train a Deep Convolutional Neural Network (DCNN) to solve the task of retrieving Lagrangian patches from which dark matter halos will condense. The volumetric multi-label classification task is turned into a regression problem by means of the euclidean distance transformation. The network is complemented by an adaptive version of the watershed algorithm to form the entire protohalo identification pipeline. We show that splitting the segmentation problem into two distinct sub-tasks allows for training smaller and faster networks, while the predictive power of the pipeline remains the same. The model is trained on synthetic data derived from a single full N-body simulation and achieves deviations of ~10% when reconstructing the dark matter halo mass function at z=0. This approach represents a promising framework for learning highly non-linear relations in the primordial density field. As a practical application, our method can be used to produce mock dark matter halo catalogues directly from the initial conditions of N-body simulations.

Data with uncertain, missing, censored, and correlated values are commonplace in many research fields including astronomy. Unfortunately, such data are often treated in an ad hoc way in the astronomical literature potentially resulting in inconsistent parameter estimates. Furthermore, in a realistic setting, the variables of interest or their errors may have non-normal distributions which complicates the modeling. I present a novel approach to compute the likelihood function for such data sets. This approach employs Gaussian copulas to decouple the correlation structure of variables and their marginal distributions resulting in a flexible method to compute likelihood functions of data in the presence of measurement uncertainty, censoring, and missing data. I demonstrate its use by determining the slope and intrinsic scatter of the star forming sequence of nearby galaxies from observational data. The outlined algorithm is implemented as the flexible, easy-to-use, open-source Python package LEO-Py.

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a wide range of galaxy properties, from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a well-ordered disk is present, we calculate the rotational profile <v_phi>(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM_enc/r) assuming spherical symmetry. In the simulated massive high-redshift galaxies, the gas rotation traces the circular velocity reasonably well at intermediate radii r~1-3 kpc, but the two quantities diverge significantly outside that range. At larger radii, gradients in the turbulent pressure can bias dynamical mass measurements low by ~10-40%. In the interior, the assumption of a spherically-symmetric gravitational potential becomes increasingly poor owing to a massive disk component, reducing the gas rotational velocities by >~10%. Finally, in the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies.

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimeter broadband flux to total infrared luminosity (L_IR), and hence star formation rate, in high-z galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g., redshift and star formation rate. In this paper, we analyse the dust temperatures of massive (M* > 10^10 Msun) z=2-6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. At z~2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature does not strongly evolve with redshift over z=2-6 at fixed IR luminosity but is tightly correlated with L_IR at fixed z. The mass-weighted temperature is important for accurately estimating the total dust mass. We also analyse an 'equivalent' dust temperature for converting (sub)millimeter flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature ('warmer dust') do not necessarily have higher mass-weighted temperatures. A 'two-phase' picture for interstellar dust can explain the different scaling relations of the various dust temperatures.

We present the first detailed study of the spatially-resolved dust continuum emission of simulated galaxies at 1<z<5. We run the radiative transfer code SKIRT on a sample of submillimeter-bright galaxies drawn from the Feedback in Realistic Environments (FIRE) project. These simulated galaxies reach Milky Way masses by z=2. Our modelling provides predictions for the full rest-frame far-ultraviolet-to-far-infrared spectral energy distributions of these simulated galaxies, as well as 25-pc-resolution maps of their emission across the wavelength spectrum. The derived morphologies are notably different in different wavebands, with the same galaxy often appearing clumpy and extended in the far-ultraviolet yet an ordered spiral at far-infrared wavelengths. The observed-frame 870-$\mu$m half-light radii of our FIRE-2 galaxies are ~0.5-4kpc, consistent with existing ALMA observations of galaxies with similarly high redshifts and stellar masses. In both simulated and observed galaxies, the dust continuum emission is generally more compact than the cold gas and the dust mass, but more extended than the stellar component. The most extreme cases of compact dust emission seem to be driven by particularly compact recent star-formation, which generates steep dust temperature gradients. Our results confirm that the spatial extent of the dust continuum emission is sensitive to both the dust mass and SFR distributions.

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) at z>=5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust extinction, dust emission, and dust temperature within these simulated z>=5 galaxies. Our sample forms a tight correlation between infrared excess (IRX=F_IR/F_UV) and ultraviolet (UV)-continuum slope (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L_IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L_IR, the peak wavelength of the dust spectral energy distributions for z>=5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z=0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (LFs) at z=5-10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IRLFs and UV luminosity densities at z=5-10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z>=5 galaxies.

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.

Properties of galaxies vary systematically with the mass of their parent dark matter halos. This basic galaxy - halo connection shows a fair amount of scatter whose origin is not fully understood. Here, we study how differences in the halo assembly history affect central galaxies in low mass (M_halo < 10^12 M_sun) halos at z=2-6 with the help of the MassiveFIRE suite of cosmological simulations. In contrast to previous works that tie galaxy properties to halo concentration and halo formation redshift, we focus on halo growth rate as a measure of assembly history. We find that, at fixed halo mass, faster growing halos tend to have lower stellar masses and higher SFRs per unit stellar mass but similar overall SFRs. We provide a simple explanation for these findings with the help of an analytic model that captures approximately the behavior of our hydrodynamical simulations. Specifically, among halos of a given current mass, quickly growing halos have lower stellar masses (and thus higher sSFRs) because they were less massive and had comparably lower cold gas masses and SFRs in the past than slowly growing halos. By combining these findings with estimates for the scatter of the halo growth rate, we show that variations in growth rate at fixed halo mass may largely explain the scatter of the stellar mass - halo mass relation. In contrast, halo growth variations likely play only a minor role in the scatter of the star forming sequence in low mass galaxies.

We explore the radial variation of star formation histories in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 9 low-mass field dwarfs with M_ star = 10^5 - 10^7 M_sun from previous FIRE results, and a new suite of 17 higher mass field dwarfs with M_star = 10^7 - 10^9 M_sun introduced here. We find that age gradients are common in our dwarfs, with older stars dominant at large radii. The strength of the gradient correlates with overall galaxy age such that earlier star formation produces a more pronounced gradient. The relation between formation time and strength of the gradient is driven by both mergers and star-formation feedback. Mergers can both steepen and flatten the age gradient depending on the timing of the merger and star formation history of the merging galaxy. In galaxies without significant mergers, early feedback pushes stars to the outskirts at early times. Interestingly, among galaxies without mergers, those with large dark matter cores have flatter age gradients because these galaxies have more late-time feedback. If real galaxies have age gradients as we predict, stellar population studies that rely on sampling a limited fraction of a galaxy can give a biased view of its global star formation history. We show that central fields can be biased young by a few Gyrs while outer fields are biased old. Fields positioned near the 2D half-light radius will provide the least biased measure of a dwarf galaxy's global star formation history.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

We present an analysis of the global and spatially-resolved Kennicutt-Schmidt (KS) star formation relation in the FIRE (Feedback In Realistic Environments) suite of cosmological simulations, including halos with $z = 0$ masses ranging from $10^{10}$ -- $10^{13}$ M$_{\odot}$. We show that the KS relation emerges and is robustly maintained due to the effects of feedback on local scales regulating star-forming gas, independent of the particular small-scale star formation prescriptions employed. We demonstrate that the time-averaged KS relation is relatively independent of redshift and spatial averaging scale, and that the star formation rate surface density is weakly dependent on metallicity and inversely dependent on orbital dynamical time. At constant star formation rate surface density, the `Cold \& Dense' gas surface density (gas with $T < 300$~K and $n > 10$~cm$^{-3}$, used as a proxy for the molecular gas surface density) of the simulated galaxies is $\sim$0.5~dex less than observed at $\sim$kpc scales. This discrepancy may arise from underestimates of the local column density at the particle-scale for the purposes of shielding in the simulations. Finally, we show that on scales larger than individual giant molecular clouds, the primary condition that determines whether star formation occurs is whether a patch of the galactic disk is thermally Toomre-unstable (not whether it is self-shielding): once a patch can no longer be thermally stabilized against fragmentation, it collapses, becomes self-shielding, cools, and forms stars, regardless of epoch or environment.

The process by which massive galaxies transition from blue, star-forming disks into red, quiescent galaxies remains one of the most poorly-understood aspects of galaxy evolution. In this investigation, we attempt to gain a better understanding of how star formation is quenched by focusing on a massive post-starburst galaxy at z = 0.747. The target has a high stellar mass and a molecular gas fraction of ~30% -- unusually high for its low star formation rate. We look for indicators of star formation suppression mechanisms in the stellar kinematics and age distribution of the galaxy obtained from spatially resolved Gemini Integral-Field spectra and in the gas kinematics obtained from ALMA. We find evidence of significant rotation in the stars, but we do not detect a stellar age gradient within 5 kpc. The molecular gas is aligned with the stellar component, and we see no evidence of strong gas outflows. Our target may represent the product of a merger-induced starburst or of morphological quenching; however, our results are not completely consistent with any of the prominent quenching models.

Recent long wavelength observations on the thermal dust continuum suggest that the Rayleigh-Jeans (RJ) tail can be used as a time-efficient quantitative probe of the dust and ISM mass in high-$z$ galaxies. We use high-resolution cosmological simulations from the Feedback in Realistic Environment (FIRE) project to analyze the dust emission of $M_*>10^{10}\;M_{\odot}$ galaxies at $z=2-4$. Our simulations (MassiveFIRE) explicitly include various forms of stellar feedback, and they produce the stellar masses and star formation rates of high-$z$ galaxies in agreement with observations. Using radiative transfer modelling, we show that sub-millimeter (sub-mm) luminosity and molecular ISM mass are tightly correlated and that the overall normalization is in quantitative agreement with observations. Notably, sub-mm luminosity traces molecular ISM mass even during starburst episodes as dust mass and mass-weighted temperature evolve only moderately between $z=4$ and $z=2$, including during starbursts. Our finding supports the empirical approach of using broadband sub-mm flux as a proxy for molecular gas content in high-$z$ galaxies. We thus expect single-band sub-mm observations with ALMA to dramatically increase the sample size of high-$z$ galaxies with reliable ISM masses in the near future.

We present a suite of cosmological zoom-in simulations at z>5 from the Feedback In Realistic Environments project, spanning a halo mass range M_halo~10^8-10^12 M_sun at z=5. We predict the stellar mass-halo mass relation, stellar mass function, and luminosity function in several bands from z=5-12. The median stellar mass-halo mass relation does not evolve strongly at z=5-12. The faint-end slope of the luminosity function steepens with increasing redshift, as inherited from the halo mass function at these redshifts. Below z~6, the stellar mass function and ultraviolet (UV) luminosity function slightly flatten below M_star~10^4.5 M_sun (fainter than M_1500~-12), owing to the fact that star formation in low-mass halos is suppressed by the ionizing background by the end of reionization. Such flattening does not appear at higher redshifts. We provide redshift-dependent fitting functions for the SFR-M_halo, SFR-M_star, and broad-band magnitude-stellar mass relations. We derive the star formation rate density and stellar mass density at z=5-12 and show that the contribution from very faint galaxies becomes more important at z>8. Furthermore, we find that the decline in the z~6 UV luminosity function brighter than M_1500~-20 is largely due to dust attenuation. Approximately 37% (54%) of the UV luminosity from galaxies brighter than M_1500=-13 (-17) is obscured by dust at z~6. Our results broadly agree with current data and can be tested by future observations.

We study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of $M_{\rm halo} \approx10^{10}\,M_{\odot}$ dark matter halos as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the $z=0$ stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, $z_{\rm f}$, in both CDM and WDM models. WDM halos follow the same $M_{\star}(z=0)-V_{\rm max}(z_{\rm f})$ relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest $M_{\star}$ WDM galaxies form all of their stars at late times. The discovery of young ultra-faint dwarf galaxies with no ancient star formation -- which do not exist in our CDM simulations -- would therefore provide evidence in support of WDM.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We quantify the gas-phase abundance of deuterium and fractional contribution of stellar mass loss to the gas in cosmological zoom-in simulations from the Feedback In Realistic Environments project. At low metallicity, our simulations confirm that the deuterium abundance is very close to the primordial value. The chemical evolution of the deuterium abundance that we derive here agrees quantitatively with analytical chemical evolution models. We furthermore find that the relation between the deuterium and oxygen abundance exhibits very little scatter. We compare our simulations to existing high-redshift observations in order to determine a primordial deuterium fraction of 2.549 +/- 0.033 x 10^-5 and stress that future observations at higher metallicity can also be used to constrain this value. At fixed metallicity, the deuterium fraction decreases slightly with decreasing redshift, due to the increased importance of mass loss from intermediate-mass stars. We find that the evolution of the average deuterium fraction in a galaxy correlates with its star formation history. Our simulations are consistent with observations of the Milky Way's interstellar medium: the deuterium fraction at the solar circle is 85-92 per cent of the primordial deuterium fraction. We use our simulations to make predictions for future observations. In particular, the deuterium abundance is lower at smaller galactocentric radii and in higher mass galaxies, showing that stellar mass loss is more important for fuelling star formation in these regimes (and can even dominate). Gas accreting onto galaxies has a deuterium fraction above that of the galaxies' interstellar medium, but below the primordial fraction, because it is a mix of gas accreting from the intergalactic medium and gas previously ejected or stripped from galaxies.

We study the morphologies and sizes of galaxies at z>5 using high-resolution cosmological zoom-in simulations from the Feedback In Realistic Environments project. The galaxies show a variety of morphologies, from compact to clumpy to irregular. The simulated galaxies have more extended morphologies and larger sizes when measured using rest-frame optical B-band light than rest-frame UV light; sizes measured from stellar mass surface density are even larger. The UV morphologies are usually dominated by several small, bright young stellar clumps that are not always associated with significant stellar mass. The B-band light traces stellar mass better than the UV, but it can also be biased by the bright clumps. At all redshifts, galaxy size correlates with stellar mass/luminosity with large scatter. The half-light radii range from 0.01 to 0.2 arcsec (0.05-1 kpc physical) at fixed magnitude. At z>5, the size of galaxies at fixed stellar mass/luminosity evolves as (1+z)^{-m}, with m~1-2. For galaxies less massive than M_star~10^8 M_sun, the ratio of the half-mass radius to the halo virial radius is ~10% and does not evolve significantly at z=5-10; this ratio is typically 1-5% for more massive galaxies. A galaxy's "observed" size decreases dramatically at shallower surface brightness limits. This effect may account for the extremely small sizes of z>5 galaxies measured in the Hubble Frontier Fields. We provide predictions for the cumulative light distribution as a function of surface brightness for typical galaxies at z=6.

The physical mechanisms that quench star formation, turning blue star-forming galaxies into red quiescent galaxies, remain unclear. In this Letter, we investigate the role of gas supply in suppressing star formation by studying the molecular gas content of post-starburst galaxies. Leveraging the wide area of the SDSS, we identify a sample of massive intermediate-redshift galaxies that have just ended their primary epoch of star formation. We present ALMA CO(2-1) observations of two of these post-starburst galaxies at z~0.7 with M* ~ 2x10^11 Msun. Their molecular gas reservoirs of (6.4 +/- 0.8) x 10^9 Msun and (34.0 +/- 1.6) x 10^9 Msun are an order of magnitude larger than comparable-mass galaxies in the local universe. Our observations suggest that quenching does not require the total removal or depletion of molecular gas, as many quenching models suggest. However, further observations are required both to determine if these apparently quiescent objects host highly obscured star formation and to investigate the intrinsic variation in the molecular gas properties of post-starburst galaxies.

Accurate measurements of galaxy masses and sizes are key to tracing galaxy evolution over time. Cosmological zoom-in simulations provide an ideal test bed for assessing the recovery of galaxy properties from observations. Here, we utilize galaxies with $M_*\sim10^{10}-10^{11.5}M_{\odot}$ at z~1.7-2 from the MassiveFIRE cosmological simulation suite, part of the Feedback in Realistic Environments (FIRE) project. Using mock multi-band images, we compare intrinsic galaxy masses and sizes to observational estimates. We find that observations accurately recover stellar masses, with a slight average underestimate of ~0.06 dex and a ~0.15 dex scatter. Recovered half-light radii agree well with intrinsic half-mass radii when averaged over all viewing angles, with a systematic offset of ~0.1 dex (with the half-light radii being larger) and a scatter of ~0.2 dex. When using color gradients to account for mass-to-light variations, recovered half-mass radii also exceed the intrinsic half-mass radii by ~0.1 dex. However, if not properly accounted for, aperture effects can bias size estimates by ~0.1 dex. No differences are found between the mass and size offsets for star-forming and quiescent galaxies. Variations in viewing angle are responsible for ~25% of the scatter in the recovered masses and sizes. Our results thus suggest that the intrinsic scatter in the mass-size relation may have previously been overestimated by ~25%. Moreover, orientation-driven scatter causes the number density of very massive galaxies to be overestimated by ~0.5 dex at $M_*\sim10^{11.5}M_{\odot}$.

We introduce massive black holes (BHs) in the Feedback In Realistic Environments project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass halos ($M_{\rm halo}(z=2) \approx 10^{12.5}\,\rm{M}_{\odot}$) down to $z=1$. These simulations model stellar feedback by supernovae, stellar winds, and radiation, and BH growth using a gravitational torque-based prescription tied to resolved properties of galactic nuclei. We do not include BH feedback. We show that early BH growth occurs through short ($\lesssim 1\,$Myr) accretion episodes that can reach or even exceed the Eddington rate. In this regime, BH growth is limited by bursty stellar feedback continuously evacuating gas from galactic nuclei, and BHs remain under-massive relative to the local $M_{\rm BH}$-$M_{\rm bulge}$ relation. BH growth is more efficient at later times, when the nuclear stellar potential retains a significant gas reservoir, star formation becomes less bursty, and galaxies settle into a more ordered state, with BHs rapidly converging onto the scaling relation when the host reaches $M_{\rm bulge} \sim 10^{10}\,\rm{M}_{\odot}$. Our results are not sensitive to the details of the accretion model so long as BH growth is tied to the gas content within $\sim 100\,$pc of the BH. Our simulations imply that bursty stellar feedback has strong implications for BH and AGN demographics, especially in the early Universe and for low-mass galaxies.

Star formation rate (SFR) distributions of galaxies are often assumed to be bimodal with modes corresponding to star-forming and quiescent galaxies, respectively. Both classes of galaxies are typically studied separately and SFR distributions of star-forming galaxies are commonly modelled as lognormals. Using both observational data and results from numerical simulations, I argue that this division into star-forming and quiescent galaxies is unnecessary from a theoretical point of view and that the SFR distributions of the whole population can be well fit by zero-inflated negative binomial distributions. This family of distributions has 3 parameters that determine the average SFR of the galaxies in the sample, the scatter relative to the star-forming sequence, and the fraction of galaxies with zero SFRs, respectively. The proposed distributions naturally account for (i) the discrete nature of star formation, (ii) the presence of 'dead' galaxies with zero SFRs, and (iii) asymmetric scatter. Excluding 'dead' galaxies, the distribution of log SFR is unimodal with a peak at the star forming sequence and an extended tail towards low SFRs. However, uncertainties and biases in the SFR measurements can create the appearance of a bimodal distribution.

We use cosmological hydrodynamic simulations with stellar feedback from the FIRE project to study the physical nature of Lyman limit systems (LLSs) at z<1. At these low redshifts, LLSs are closely associated with dense gas structures surrounding galaxies, such as galactic winds, dwarf satellites, and cool inflows from the intergalactic medium. Our analysis is based on 14 zoom-in simulations covering the halo mass range M_h~10^9-10^13 Msun at z=0, which we convolve with the dark matter halo mass function to produce cosmological statistics. We find that the majority of cosmologically-selected LLSs are associated with halos in the mass range 10^10 < M_h < 10^12 Msun. The incidence and HI column density distribution of simulated absorbers with columns 10^16.2 < N_HI < 2x10^20 cm^-2 are consistent with observations. High-velocity outflows (with radial velocity exceeding the halo circular velocity by a factor >~2) tend to have higher metallicities ([X/H] ~ -0.5) while very low metallicity ([X/H] < -2) LLSs are typically associated with gas infalling from the intergalactic medium. However, most LLSs occupy an intermediate region in metallicity-radial velocity space, for which there is no clear trend between metallicity and radial kinematics. Metal-enriched inflows arise in the FIRE simulations as a result of galactic winds that fall back onto galaxies at low redshift. The overall simulated LLS metallicity distribution has a mean (standard deviation) [X/H] = -0.9 (0.4) and does not show significant evidence for bimodality, in contrast to recent observational studies but consistent with LLSs arising from halos with a broad range of masses and metallicities.

We analyse the star formation rates (SFRs), colours, and dust extinctions of galaxies in massive (10^12.5-10^13.5 Msun) haloes at z~2 in high-resolution, cosmological zoom-in simulations as part of the Feedback in Realistic Environments (FIRE) project. The simulations do not model feedback from active galactic nuclei (AGN) but reproduce well the observed relations between stellar and halo mass and between stellar mass and SFR. About half (a third) of the simulated massive galaxies (massive central galaxies) at z~2 have broad-band colours classifying them as 'quiescent', and the fraction of quiescent centrals is steeply decreasing towards higher redshift, in agreement with observations. The progenitors of z~2 quiescent central galaxies are, on average, more massive, have lower specific SFRs, and reside in more massive haloes than the progenitors of similarly massive star forming centrals. The simulations further predict a morphological mix of galaxies that includes disk-dominated, irregular, and early-type galaxies. However, our simulations do not reproduce the reddest of the quiescent galaxies observed at z~2. We also do not find evidence for a colour bimodality, but are limited by our modest sample size. In our simulations, the star formation activity of central galaxies of moderate mass (Mstar~10^10-10^11 Msun) is affected by a combination of two distinct physical processes. Outflows powered by stellar feedback result in a short-lived (<100 Myr), but almost complete, suppression of star formation activity after which many galaxies quickly recover and continue to form stars at normal rates. In addition, galaxies residing in slowly growing haloes tend to experience a moderate reduction of their SFRs ('cosmological starvation'). The relative importance of these processes and AGN feedback is uncertain and will be explored in future work.

Recent spatially resolved observations of galaxies at z=0.6-3 reveal that high-redshift galaxies show complex kinematics and a broad distribution of gas-phase metallicity gradients. To understand these results, we use a suite of high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, which include physically motivated models of the multi-phase ISM, star formation, and stellar feedback. Our simulations reproduce the observed diversity of kinematic properties and metallicity gradients, broadly consistent with observations at z=0-3. Strong negative metallicity gradients only appear in galaxies with a rotating disk, but not all rotationally supported galaxies have significant gradients. Strongly perturbed galaxies with little rotation always have flat gradients. The kinematic properties and metallicity gradient of a high-redshift galaxy can vary significantly on short time-scales, associated with starburst episodes. Feedback from a starburst can destroy the gas disk, drive strong outflows, and flatten a pre-existing negative metallicity gradient. The time variability of a single galaxy is statistically similar to the entire simulated sample, indicating that the observed metallicity gradients in high-redshift galaxies reflect the instantaneous state of the galaxy rather than the accretion and growth history on cosmological time-scales. We find weak dependence of metallicity gradient on stellar mass and specific star formation rate (sSFR). Low-mass galaxies and galaxies with high sSFR tend to have flat gradients, likely due to the fact that feedback is more efficient in these galaxies. We argue that it is important to resolve feedback on small scales in order to produce the diverse metallicity gradients observed.

Galaxy formation models are now able to reproduce observed relations such as the relation between galaxies' star formation rates (SFRs) and stellar masses ($M_*$) and the stellar mass--halo mass relation. We demonstrate that comparisons of the short-timescale variability in galaxy SFRs with observational data provide an additional useful constraint on the physics of galaxy formation feedback. We apply SFR indicators with different sensitivity timescales to galaxies from the Feedback in Realistic Environments (FIRE) simulations. We find that the SFR--$M_*$ relation has a significantly greater scatter when the H$\alpha$-derived SFR is considered compared with when the far-ultraviolet (FUV)-based SFR is used. This difference is a direct consequence of bursty star formation because the FIRE galaxies exhibit order-of-magnitude SFR variations over timescales of a few Myr. We show that the difference in the scatter between the simulated H$\alpha$- and FUV-derived SFR--$M_*$ relations at $z=2$ is consistent with observational constraints. We also find that the H$\alpha$/FUV ratios predicted by the simulations at $z=0$ are similar to those observed for local galaxies except for a population of low-mass ($M_* \lesssim 10^{9.5} {\rm M}_\odot$) simulated galaxies with lower H$\alpha$/FUV ratios than observed. We suggest that future cosmological simulations should compare the H$\alpha$/FUV ratios of their galaxies with observations to constrain the feedback models employed.

The morphology of massive star-forming galaxies at high redshift is often dominated by giant clumps of mass ~10^8-10^9 Msun and size ~100-1000 pc. Previous studies have proposed that giant clumps might have an important role in the evolution of their host galaxy, particularly in building the central bulge. However, this depends on whether clumps live long enough to migrate from their original location in the disc or whether they get disrupted by their own stellar feedback before reaching the centre of the galaxy. We use cosmological hydrodynamical simulations from the FIRE (Feedback in Realistic Environments) project that implement explicit treatments of stellar feedback and ISM physics to study the properties of these clumps. We follow the evolution of giant clumps in a massive (stellar mass ~10^10.8 Msun at z=1), discy, gas-rich galaxy from redshift z>2 to z=1. Even though the clumpy phase of this galaxy lasts over a gigayear, individual gas clumps are short-lived, with mean lifetime of massive clumps of ~20 Myr. During that time, they turn between 0.1% and 20% of their gas into stars before being disrupted, similar to local GMCs. Clumps with M>10^7 Msun account for ~20% of the total star formation in the galaxy during the clumpy phase, producing ~10^10 Msun of stars. We do not find evidence for net inward migration of clumps within the galaxy. The number of giant clumps and their mass decrease at lower redshifts, following the decrease in the overall gas fraction and star-formation rate.

The thermal Sunyaev-Zel'dovich (SZ) effect and soft X-ray emission are routinely observed around massive galaxies and in galaxy groups and clusters. We study these observational diagnostics of galaxy haloes for a suite of cosmological `zoom-in' simulations from the `Feedback In Realistic Environments' project, which spans a large range in halo mass 10^10-10^13 Msun). We explore the effect of stellar feedback on the hot gas observables. The properties of our simulated groups, such as baryon fractions, SZ flux, and X-ray luminosities (L_X), are broadly consistent with existing observations, even though feedback from active galactic nuclei is not included. We make predictions for future observations of lower-mass objects for both SZ and diffuse X-ray measurements, finding that they are not just scaled-down versions of massive galaxies, but more strongly affected by galactic winds driven by star formation. Low-mass haloes (<~10^11 Msun) retain a low fraction of their baryons, which results in a strong suppression of the SZ signal. Our simulations therefore predict a scaling with halo mass that is steeper than self-similar for haloes less massive than 10^13 Msun. For halo masses <~10^12 Msun, L_X is time-variable and correlated primarily with the star formation rate (SFR). For these objects, the diffuse X-ray emission is powered mostly by galactic winds and the gas dominating the X-ray emission is flowing out with radial velocities close to the halo's circular velocity. For halo masses >~10^13 Msun, on the other hand, L_X is much less variable and not correlated with the SFR, because the emission originates from the quasi-hydrostatic, virialized halo gas.

We use cosmological simulations from the Feedback In Realistic Environments (FIRE) project, which implement a comprehensive set of stellar feedback processes, to study ultra-violet (UV) metal line emission from the circum-galactic medium of high-redshift (z=2-4) galaxies. Our simulations cover the halo mass range Mh ~ 2x10^11 - 8.5x10^12 Msun at z=2, representative of Lyman break galaxies. Of the transitions we analyze, the low-ionization C III (977 A) and Si III (1207 A) emission lines are the most luminous, with C IV (1548 A) and Si IV (1394 A) also showing interesting spatially-extended structures. The more massive halos are on average more UV-luminous. The UV metal line emission from galactic halos in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ~2 dex. The peaks of UV metal line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic halos to power the metal line emission. The strong time variability implies that even some relatively low-mass halos may be detectable. Conversely, flux-limited samples will be biased toward halos whose central galaxy has recently experienced a strong burst of star formation. Spatially-extended UV metal line emission around high-redshift galaxies should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope and Keck Cosmic Web Imager (KCWI).

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

Observations reveal that quasar host halos at z~2 have large covering fractions of cool dense gas (>~60% for Lyman limit systems within a projected virial radius). Most simulations have so far have failed to explain these large observed covering fractions. We analyze a new set of 15 simulated massive halos with explicit stellar feedback from the FIRE project, covering the halo mass range M_h~2x10^12-10^13 Msun at z=2. This extends our previous analysis of the circum-galactic medium of high-redshift galaxies to more massive halos. AGN feedback is not included in these simulations. We find Lyman limit system covering fractions consistent with those observed around quasars. The large HI covering fractions arise from star formation-driven galactic winds, including winds from low-mass satellite galaxies that interact with cosmological filaments. We show that it is necessary to resolve these satellite galaxies and their winds to reproduce the large Lyman limit system covering fractions observed in quasar-mass halos. Our simulations predict that galaxies occupying dark matter halos of mass similar to quasars but without a luminous AGN should have Lyman limit system covering fractions comparable to quasars.

The cosmic noon (z~1.5-3) marked a period of vigorous star formation for most galaxies. However, about a third of the more massive galaxies at those times were quiescent in the sense that their observed stellar populations are inconsistent with rapid star formation. The reduced star formation activity is often attributed to gaseous outflows driven by feedback from supermassive black holes, but the impact of black hole feedback on galaxies in the young Universe is not yet definitively established. We analyze the origin of quiescent galaxies with the help of ultra-high resolution, cosmological simulations that include feedback from stars but do not model the uncertain consequences of black hole feedback. We show that dark matter halos with specific accretion rates below ~0.25-0.4 per Gyr preferentially host galaxies with reduced star formation rates and red broad-band colors. The fraction of such halos in large dark matter only simulations matches the observed fraction of massive quiescent galaxies (~10^10-10^11 Msun). This strongly suggests that halo accretion rate is the key parameter determining which massive galaxies at z~1.5-3 become quiescent. Empirical models that connect galaxy and halo evolution, such as halo occupation distribution or abundance matching models, assume a tight link between galaxy properties and the masses of their parent halos. These models will benefit from adding the specific accretion rate of halos as a second model parameter.

Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).

Most galaxy evolution simulations as well as a variety of observational methods assume a linear scaling between the (galaxy-averaged) dust-to-gas ratio D and metallicity Z of the interstellar medium (ISM). Indeed, nearby galaxies with solar or moderately sub-solar metallicities clearly follow this trend albeit with significant scatter. However, a growing number of observations show that the linear scaling breaks down for metal-poor galaxies (Z<0.2 Z_sun), highlighting the need for a more sophisticated modeling of the dust-to-metal ratio of galaxies. Here we study the co-evolution of dust and metal abundances in galaxies with the help of a dynamical, one-zone model that incorporates dust formation and destruction processes in addition to gas inflows, outflows, and metal enrichment. The dynamical model is consistent with various observational constraints, including the stellar mass -- metallicity relation, the stellar mass -- halo mass relation, and the observed Z -- D relation for both metal-poor and metal-rich galaxies. The functional form of the Z -- D relation follows from a basic equilibrium ansatz, similar to the ideas used previously to model the stellar mass -- metallicity relation. Galactic outflows regulate the inflow rate of gas from the cosmic web for galaxies of a given star formation rate. The mass loading factor of outflows thus dictates the rate at which the dust and metal content of the ISM is diluted. The stellar mass dependence of the mass loading factor drives the evolution of metallicities, dust-to-gas ratios, and dust-to-metal ratios in galaxies.

The Hubble sequence is a common classification scheme for the structure of galaxies. Despite the tremendous usefulness of this diagnostic, we still do not fully understand when, where, and how this morphological ordering was put in place. Here, we investigate the morphological evolution of a sample of 22 high redshift ($z\geq3$) galaxies extracted from the Argo simulation. Argo is a cosmological zoom-in simulation of a group-sized halo and its environment. It adopts the same high resolution ($\sim10^4$ M$_\odot$, $\sim100$ pc) and sub-grid physical model that was used in the Eris simulation but probes a sub-volume almost ten times bigger with as many as 45 million gas and star particles in the zoom-in region. Argo follows the early assembly of galaxies with a broad range of stellar masses ($\log M_{\star}/{\rm M}_{\odot}\sim8-11$ at $z\simeq3$), while resolving properly their structural properties. We recover a diversity of morphologies, including late-type/irregular disc galaxies with flat rotation curves, spheroid dominated early-type discs, and a massive elliptical galaxy, already established at $z\sim3$. We identify major mergers as the main trigger for the formation of bulges and the steepening of the circular velocity curves. Minor mergers and non-axisymmetric perturbations (stellar bars) drive the bulge growth in some cases. The specific angular momenta of the simulated disc components fairly match the values inferred from nearby galaxies of similar $M_{\star}$ once the expected redshift evolution of disc sizes is accounted for. We conclude that morphological transformations of high redshift galaxies of intermediate mass are likely triggered by processes similar to those at low redshift and result in an early build-up of the Hubble sequence.

Observations show a prevalence of high redshift galaxies with large stellar masses and predominantly passive stellar populations. A variety of processes have been suggested that could reduce the star formation in such galaxies to observed levels, including quasar mode feedback, virial shock heating, or galactic winds driven by stellar feedback. However, the main quenching mechanisms have yet to be identified. Here we study the origin of star formation quenching using Argo, a cosmological, hydrodynamical zoom-in simulation that follows the evolution of a massive galaxy at $z\geq{}2$. This simulation adopts the same sub-grid recipes of the Eris simulations, which have been shown to form realistic disk galaxies, and, in one version, adopts also a mass and spatial resolution identical to Eris. The resulting galaxy has properties consistent with those of observed, massive (M_* ~ 1e11 M_sun) galaxies at z~2 and with abundance matching predictions. Our models do not include AGN feedback indicating that supermassive black holes likely play a subordinate role in determining masses and sizes of massive galaxies at high z. The specific star formation rate (sSFR) of the simulated galaxy matches the observed M_* - sSFR relation at early times. This period of smooth stellar mass growth comes to a sudden halt at z=3.5 when the sSFR drops by almost an order of magnitude within a few hundred Myr. The suppression is initiated by a leveling off and a subsequent reduction of the cool gas accretion rate onto the galaxy, and not by feedback processes. This "cosmological starvation" occurs as the parent dark matter halo switches from a fast collapsing mode to a slow accretion mode. Additional mechanisms, such as perhaps radio mode feedback from an AGN, are needed to quench any residual star formation of the galaxy and to maintain a low sSFR until the present time.

Cold Dark Matter (CDM) theory, a pillar of modern cosmology and astrophysics, predicts the existence of a large number of starless dark matter halos surrounding the Milky Way (MW). However, clear observational evidence of these "dark" substructures remains elusive. Here, we present a detection method based on the small, but detectable, velocity changes that an orbiting substructure imposes on the stars in the MW disk. Using high-resolution numerical simulations we estimate that the new space telescope Gaia should detect the kinematic signatures of a few starless substructures provided the CDM paradigm holds. Such a measurement will provide unprecedented constraints on the primordial matter power spectrum at low-mass scales and offer a new handle onto the particle physics properties of dark matter.

Observations show that star formation in galaxies is closely correlated with the abundance of molecular hydrogen. Modeling this empirical relation from first principles proves challenging, however, and many questions regarding its properties remain open. For instance, the exact functional form of the relation is still debated and it is also unknown whether it applies at z>4, where CO observations are sparse. Here, we analyze how the shape of the star formation -- gas relation affects the cosmic star formation history and global galaxy properties using an analytic model that follows the average evolution of galaxies in dark matter halos across cosmic time. We show that a linear relation with an H2 depletion time of ~2.5 Gyr, as found in studies of nearby galaxies, results in good agreement with current observations of galaxies at both low and high redshift. These observations include the evolution of the cosmic star formation rate density, the z~4-9 UV luminosity function, the evolution of the mass -- metallicity relation, the relation between stellar and halo mass, and the gas-to-stellar mass ratios of galaxies. In contrast, the short depletion times that result from adopting a highly super-linear star formation -- gas relation lead to large star formation rates, substantial metal enrichment (~0.1 solar), and low gas-to-stellar mass ratios already at z~10, in disagreement with observations. These results can be understood in terms of an equilibrium picture of galaxy evolution in which gas inflows, outflows, and star formation drive the metallicities and gas fractions toward equilibrium values that are determined by the ratio of the accretion time to the gas depletion time. In this picture, the cosmic modulation of the accretion rate is the primary process that drives the evolution of stellar masses, gas masses, and metallicities of galaxies from high redshift until today.

Molecular hydrogen (H2) is the primary component of the reservoirs of cold, dense gas that fuel star formation in our galaxy. While the H2 abundance is ultimately regulated by physical processes operating on small scales in the interstellar medium (ISM), observations have revealed a tight correlation between the ratio of molecular to atomic hydrogen in nearby spiral galaxies and the pressure in the mid-plane of their disks. This empirical relation has been used to predict H2 abundances in galaxies with potentially very different ISM conditions, such as metal-deficient galaxies at high redshifts. Here, we test the validity of this approach by studying the dependence of the pressure -- H2 relation on environmental parameters of the ISM. To this end, we follow the formation and destruction of H2 explicitly in a suite of hydrodynamical simulations of galaxies with different ISM parameters. We find that a pressure -- H2 relation arises naturally in our simulations for a variety of dust-to-gas ratios or strengths of the interstellar radiation field in the ISM. Fixing the dust-to-gas ratio and the UV radiation field to values measured in the solar neighborhood results in fair agreement with the relation observed in nearby galaxies with roughly solar metallicity. However, the parameters (slope and normalization) of the pressure -- H2 relation vary in a systematical way with ISM properties. A particularly strong trend is the decrease of the normalization of the relation with a lowering of the dust-to-gas ratio of the ISM. We show that this trend and other properties of the pressure -- H2 relation are natural consequences of the transition from atomic to molecular hydrogen with gas surface density.

Interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way provide the majority of the gamma rays observed by the Fermi Gamma Ray Space Telescope. In addition to the gas which is densely concentrated along the Galactic Disk, hydrodynamical simulations and observational evidence favor the presence of a halo of hot (T~10^6 K) ionized hydrogen (H_II), extending with non-negligible densities out to the virial radius of the Milky Way. We show that cosmic ray collisions with this circum-galactic gas should be expected to provide a significant flux of gamma rays, on the order of 10% of the observed isotopic gamma ray background at energies above 1 GeV. In addition, gamma rays originating from the extended H_II halos of other galaxies along a given line-of-sight should contribute to this background at a similar level.

There is ample observational evidence that the star formation rate (SFR) surface density, Sigma_SFR, is closely correlated with the surface density of molecular hydrogen, Sigma_H2. This empirical relation holds both for galaxy-wide averages and for individual >=kpc sized patches of the interstellar medium (ISM), but appears to degrade substantially at a sub-kpc scale. Identifying the physical mechanisms that determine the scale-dependent properties of the observed Sigma_H2-Sigma_SFR relation remains a challenge from a theoretical perspective. To address this question, we analyze the slope and scatter of the Sigma_H2-Sigma_SFR relation using a set of cosmological, galaxy formation simulations with a peak resolution of ~100 pc. These simulations include a chemical network for molecular hydrogen, a model for the CO emission, and a simple, stochastic prescription for star formation that operates on ~100 pc scales. Specifically, star formation is modeled as a Poisson process in which the average SFR is directly proportional to the present mass of H2. The predictions of our numerical model are in good agreement with the observed Kennicutt-Schmidt and Sigma_H2-Sigma_SFR relations. We show that observations based on CO emission are ill suited to reliably measure the slope of the latter relation at low (<20 M_sun pc^-2) H2 surface densities on sub-kpc scales. Our models also predict that the inferred Sigma_H2-Sigma_SFR relation steepens at high H2 surface densities as a result of the surface density dependence of the CO/H2 conversion factor. Finally, we show that on sub-kpc scales most of the scatter in the relation is a consequence of discreteness effects in the star formation process. In contrast, variations of the CO/H2 conversion factor are responsible for most of the scatter measured on super-kpc scales.

Massive compact galaxies seem to be more common at high redshift than in the local universe, especially in denser environments. To investigate the fate of such massive galaxies identified at z~2 we analyse the evolution of their properties in three cosmological hydrodynamical simulations that form virialised galaxy groups of mass ~10^13 Msun hosting a central massive elliptical/S0 galaxy by redshift zero. We find that at redshift ~2 the population of galaxies with M_*> 2 10^10 Msun is diverse in terms of mass, velocity dispersion, star formation and effective radius, containing both very compact and relatively extended objects. In each simulation all the compact satellite galaxies have merged into the central galaxy by redshift 0 (with the exception of one simulation where one of such satellite galaxy survives). Satellites of similar mass at z = 0 are all less compact than their high redshift counterparts. They form later than the galaxies in the z = 2 sample and enter the group potential at z < 1, when dynamical friction times are longer than the Hubble time. Also, by z = 0 the central galaxies have increased substantially their characteristic radius via a combination of in situ star formation and mergers. Hence in a group environment descendants of compact galaxies either evolve towards larger sizes or they disappear before the present time as a result of the environment in which they evolve. Since the group-sized halos that we consider are representative of dense environments in the LambdaCDM cosmology, we conclude that the majority of high redshift compact massive galaxies do not survive until today as a result of the environment.

Characterizing the conversion factor between CO emission and column density of molecular hydrogen, X_CO, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds (GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there are physical grounds to expect that it should also depend on the gas metallicity, surface density, and strength of the interstellar radiation field. The XCO factor may also depend on the scale on which CO emission is averaged due to effects of limited resolution. We study the dependence of X_CO on gas properties and averaging scale using a model that is based on a combination of results of sub-pc scale magneto-hydrodynamic simulations and on the gas distribution from self-consistent cosmological simulations of galaxy formation. Our model predicts a value of X_CO that is consistent with the Galactic value for interstellar medium conditions typical for the Milky Way. For such conditions the predicted X_CO varies by only a factor of two for gas surfaced densities in the range \sim 50 - 500 M_sun / pc^2. However, the model also predicts that more generally on the scale of GMCs, X_CO is a strong function of metallicity, and depends on the column density and the interstellar UV flux. We show explicitly that neglecting these dependencies in observational estimates can strongly bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset range compared to the true distribution. We find that when averaged on \sim kpc scales the X-factor depends only weakly on radiation field and column density, but is still a strong function of metallicity. The predicted metallicity dependence can be approximated as X_CO \sim Z^{-{\gamma}} with {\gamma} ~ 0.5 - 0.8.

We present measurements of two types of cluster galaxy alignments based on a volume limited and highly pure ($\ge$ 90%) sample of clusters from the GMBCG catalog derived from SDSS DR7. We detect a clear BCG alignment (the alignment of major axis of the BCG toward the distribution of cluster satellite galaxies). We find that the BCG alignment signal becomes stronger as the redshift and BCG absolute magnitude decrease, and becomes weaker as BCG stellar mass decreases. No dependence of the BCG alignment on cluster richness is found. We can detect a statistically significant ($\ge$ 3 sigma) satellite alignment (the alignment of the major axes of the cluster satellite galaxies toward the BCG) only when we use the isophotal fit position angles (PAs, hereafter), and the satellite alignment depends on the apparent magnitudes rather than the absolute magnitudes of the BCGs. This suggests the detected satellite alignment based on isophotoal PAs from the SDSS pipeline is possibly due to the contamination from the diffuse light of nearby BCGs. We caution that this should not be simply interpreted as non-existence of the satellite alignment, but rather that we cannot detect them with our current photometric SDSS data. We perform our measurements on both SDSS $r$ band and $i$ band data, but did not observe a passband dependence of the alignments.

The physical mechanisms and timescales that determine the morphological signatures and the quenching of star formation of typical (~L*) elliptical galaxies are not well understood. To address this issue, we have simulated the formation of a group of galaxies with sufficient resolution to track the evolution of gas and stars inside about a dozen galaxy group members over cosmic history. Galaxy groups, which harbor many elliptical galaxies in the universe, are a particularly promising environment to investigate morphological transformation and star formation quenching, due to their high galaxy density, their relatively low velocity dispersion, and the presence of a hot intragroup medium. Our simulation reproduces galaxies with different Hubble morphologies and, consequently, enables us to study when and where the morphological transformation of galaxies takes place. The simulation does not include feedback from active galactic nuclei showing that it is not an essential ingredient for producing quiescent, red elliptical galaxies in galaxy groups. Ellipticals form, as suspected, through galaxy mergers. In contrast with what has often been speculated, however, these mergers occur at z>1, before the merging progenitors enter the virial radius of the group and before the group is fully assembled. The simulation also shows that quenching of star formation in the still star-forming elliptical galaxies lags behind their morphological transformation, but, once started, is taking less than a billion years to complete. As long envisaged the star formation quenching happens as the galaxies approach and enter the finally assembled group, due to quenching of gas accretion and (to a lesser degree) stripping. A similar sort is followed by unmerged, disk galaxies, which, as they join the group, are turned into the red-and-dead disks that abound in these environments.

It is a well established empirical fact that the surface density of the star formation rate, Sigma_SFR, strongly correlates with the surface density of molecular hydrogen, Sigma_H2, at least when averaged over large (~kpc) scales. Much less is known, however, if (and how) the Sigma_SFR-Sigma_H2 relation depends on environmental parameters, such as the metallicity or the UV radiation field in the interstellar medium (ISM). Furthermore, observations indicate that the scatter in the Sigma_SFR-Sigma_H2 relation increases rapidly with decreasing averaging scale. How the scale-dependent scatter is generated and how one recovers a tight ~ kpc scale Sigma_SFR-Sigma_H2 relation in the first place is still largely debated. Here, these questions are explored with hydrodynamical simulations that follow the formation and destruction of H2, include radiative transfer of UV radiation, and resolve the ISM on ~60 pc scales. We find that within the considered range of H2 surface densities (10-100 Msun/pc^2) the Sigma_SFR-Sigma_H2 relation is steeper in environments of low metallicity and/or high radiation fields (compared to the Galaxy), that the star formation rate at a given H2 surface density is larger, and the scatter is increased. Deviations from a "universal" Sigma_SFR-Sigma_H2 relation should be particularly relevant for high redshift galaxies or for low-metallicity dwarfs at z~0. We also find that the use of time-averaged SFRs produces a large, scale dependent scatter in the Sigma_SFR-Sigma_H2 relation. Given the plethora of observational data expected from upcoming surveys such as ALMA the scale-scatter relation may indeed become a valuable tool for determining the physical mechanisms connecting star formation and H2 formation.

A star formation efficiency per free fall time that evolves over the life time of giant molecular clouds (GMCs) may have important implications for models of supersonic turbulence in molecular clouds or for the relation between star formation rate and H2 surface density. We discuss observational data that could be interpreted as evidence of such a time variability. In particular, we investigate a recent claim based on measurements of H2 and stellar masses in individual GMCs. We show that this claim depends crucially on the assumption that H2 masses do not evolve over the life times of GMCs. We exemplify our findings with a simple toy model that uses a constant star formation efficiency and, yet, is able to explain the observational data.

Here we introduce PHAT, the PHoto-z Accuracy Testing programme, an international initiative to test and compare different methods of photo-z estimation. Two different test environments are set up, one (PHAT0) based on simulations to test the basic functionality of the different photo-z codes, and another one (PHAT1) based on data from the GOODS survey. The accuracy of the different methods is expressed and ranked by the global photo-z bias, scatter, and outlier rates. Most methods agree well on PHAT0 but produce photo-z scatters that can differ by up to a factor of two even in this idealised case. A larger spread in accuracy is found for PHAT1. Few methods benefit from the addition of mid-IR photometry. Remaining biases and systematic effects can be explained by shortcomings in the different template sets and the use of priors on the one hand and an insufficient training set on the other hand. Scatters of 4-8% in Delta_z/(1+z) were obtained, consistent with other studies. However, somewhat larger outlier rates (>7.5% with Delta_z/(1+z)>0.15; >4.5% after cleaning) are found for all codes. There is a general trend that empirical codes produce smaller biases than template-based codes. The systematic, quantitative comparison of different photo-z codes presented here is a snapshot of the current state-of-the-art of photo-z estimation and sets a standard for the assessment of photo-z accuracy in the future. The rather large outlier rates reported here for PHAT1 on real data should be investigated further since they are most probably also present (and possibly hidden) in many other studies. The test data sets are publicly available and can be used to compare new methods to established ones and help in guiding future photo-z method development. (abridged)

We trace the evolution of central galaxies in three ~10^13 M_sun galaxy groups simulated at high resolution in cosmological hydrodynamical simulations. The evolution in the group potential leads, at z=0, to central galaxies that are massive, gas-poor early-type systems supported by stellar velocity dispersion resembling either elliptical or S0 galaxies. Their z~2-2.5 main progenitors are massive M* ~ 3-10 x 10^10 M_sun, star forming (20-60 M_sun/yr) galaxies which host substantial reservoirs of cold gas (~5 x 10^9 M_sun) in extended gas disks. Our simulations thus show that star forming galaxies observed at z~2 are likely the main progenitors of central galaxies in galaxy groups at z=0. Their central stellar densities stay approximately constant from z~1.5 down to z=0. Instead, the galaxies grow inside-out, by acquiring a stellar envelope outside the innermost ~2 kpc. Consequently the density within the effective radius decreases by up to two orders of magnitude. Both major and minor mergers contribute to most of the mass accreted outside the effective radius and thus drive the evolution of the half-mass radii. In one of the three simulated groups the short central cooling time leads to a dramatic rejuvenation of the central group galaxy at z<1, affecting its morphology, kinematics and colors. This episode is eventually terminated by a group-group merger. Our simulations demonstrate that, in galaxy groups, the interplay between halo mass assembly, galaxy merging and gas accretion has a substantial influence on the star formation histories and z=0 morphologies of central galaxies.[Abridged]

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

We use a set of high-resolution N-body simulations of binary galaxy mergers to show that the morphologies of the tidal features that are seen around a large fraction of nearby, massive ellipticals in the field, cannot be reproduced by equal-mass dissipationless mergers; rather, they are well explained by the accretion of disk-dominated galaxies. In particular, the arm- and looplike morphologies of the observed tidal debris can only be produced by the kinematically cold material of the disk components of the accreted galaxies. The tidal features that arise from such "cold-accretion" events onto a massive elliptical are visible for significantly longer timescales than the features produced by elliptical-elliptical mergers (about 1-2 Gyr vs. a few hundred million years). Mass ratios of the order of 1:10 between the accreting elliptical and the accreted disk are sufficient to match the brightness of the observed debris. Furthermore, stellar population synthesis models and simple order-of-magnitude calculations indicate that the colors of the tidal features generated in such minor cold-accretion events are relatively red, in agreement with the observations. The minor cold-accretion events that explain the presence, brightness, and structural and color properties of the tidal debris cause only a modest mass and luminosity increase in the accreting massive elliptical. These results, coupled with the relative statistical frequencies of disk- and bulge-dominated galaxies in the field, suggest that massive ellipticals assemble most of their mass well before their tidal debris forms through the accretion of relatively little, kinematically cold material rather than in very recent, dissipationless major mergers.

ZEBRA, the Zurich Extragalactic Bayesian Redshift Analyzer, is a tool for estimating redshifts and template types of galaxies using medium- and broad-band photometric data. ZEBRA employs novel techniques within the template-fitting approach to produce high-quality Maximum-Likelihood and Bayesian redshift estimates. This manuscript serves as a user guide to ZEBRA. It explains how to use ZEBRA, specifies input and output formats, and gives a short account of the available options. ZEBRA is a free and open-source software distributed under the GNU Public License 3 and available at http://www.exp-astro.phys.ethz.ch/ZEBRA . Several upgrades are currently being implemented in ZEBRA. An updated documentation will be provided at each new release. Any problems, comments and suggestions on the code and the manual should be sent via e-mail to zebra@phys.ethz.ch.

ABRIDGED: We study the evolution since z~1 of the rest-frame B luminosity function of the early-type galaxies (ETGs) in ~0.7 deg^2 in the COSMOS field. In order to identify ALL progenitors of local ETGs we construct the sample of high-z galaxies using two complementary criteria: (i) A morphological selection based on the Zurich Estimator of Structural Types, and (ii) A photometric selection based on the galaxy properties in the (U-V)-M_V color-magnitude diagram. We furthermore constrain both samples so as to ensure that the selected progenitors of ETGs are compatible with evolving into systems which obey the mu_B-r_{hl} Kormendy relation. Assuming the luminosity evolution derived from studies of the fundamental plane for high-z ETGs, our analysis shows no evidence for a decrease in the number density of the most massive ETGs out to z~ 0.7: Both the morphologically- and the photometrically-selected sub-samples show no evolution in the number density of bright (~L>2.5L*) ETGs. Allowing for different star formation histories, and cosmic variance, we estimate a maximum decrease in the number density of massive galaxies at that redshift of ~30%. We observe, however, in both the photometrical and morphological samples, a deficit of up to ~2-3 of fainter ETGs over the same cosmic period. Our results argue against a significant contribution of recent dissipationless ``dry'' mergers to the formation of the most massive ETGs. We suggest that the mass growth in low luminosity ETGs can be explained with a conversion from z~0.7 to z=0 of blue, irregular and disk galaxies into low- and intermediate-mass ``red'' ETGs, possibly also through gas rich mergers.

We measure photometric redshifts and spectral types for galaxies in the COSMOS survey. We use template fitting technique combined with luminosity function priors and with the option to simultaneously estimate dust extinction (i.e. E(B-V)) for each galaxy.Our estimated redshifts are accurate to i<25 and z~1.2. Using simulations with sampling and noise characteristics similar to those in COSMOS, the accuracy and reliability is estimated for the photometric redshifts as a function of the magnitude limits of the sample, S/N ratios and the number of bands used. From the simulations we find that the ratio of derived 95% confidence interval in the redshift probability distribution to the estimated photometric redshift (D95) can be used to identify and exclude the catastrophic failures in the photometric redshift estimates. We compare the derived redshifts with high-reliability spectroscopic redshifts for a sample of 868 normal galaxies with z < 1.2 from zCOSMOS. Considering different scenarios, depending on using prior, no prior and/or extinction, we compare the photometric and spectroscopic redshifts for this sample. This corresponds to an rms scatter of 0.031, with a small number of outliers (<2.5%). We also find good agreement (rms=0.10) between photometric and spectroscopic redshifts for Type II AGNs. We compare results from our photometric redshift procedure with three other independent codes and find them in excellent agreement. We show preliminary results, based on photometric redshifts for the entire COSMOS sample (to i < 25 mag.).

zCOSMOS is a large redshift survey that is being undertaken in the COSMOS field using 600 hours of observation with the VIMOS spectrograph on the 8-m VLT. The survey is designed to characterise the environments of COSMOS galaxies from the 100 kpc scales of galaxy groups up to the 100 Mpc scale of the cosmic web and to produce diagnostic information on galaxies and active galactic nuclei. The zCOSMOS survey consists of two parts: (a) zCOSMOS-bright, a magnitude-limited I-band IAB < 22.5 sample of about 20,000 galaxies with 0.1 < z < 1.2 covering the whole 1.7 deg2 COSMOS ACS field and designed to mimic the parameters of the 2dfGRS; and (b) zCOSMOS-deep, a survey of approximately 10,000 galaxies selected through colour-selection criteria to have 1.4 < z < 3.0, within the central 1 deg2. This paper describes the survey design and the construction of the target catalogues, and briefly outlines the observational program and the data pipeline. In the first observing season, spectra of 1303 zCOSMOS-bright targets and of 977 zCOSMOS-deep targets have been obtained. These are briefly analysed to demonstrate the characteristics that may be expected from zCOSMOS, and particularly zCOSMOS-bright, when it is finally completed between 2008-2009. The power of combining spectroscopic and photometric redshifts is demonstrated, especially in correctly identifying the emission line in single-line spectra and in determining which of the less reliable spectroscopic redshifts are correct and which are incorrect. Our zCOSMOS-deep spectra demonstrate the effectiveness of our selection techniques to isolate high redshift galaxies at 1.4 < z < 3.0 and of VIMOS to measure their redshifts using ultraviolet absorption lines.

(ABRIDGED) Motivated by the desire to reliably and automatically classify structure of thousands of COSMOS galaxies, we present ZEST, the Zurich Estimator of Structural Types. To classify galaxy structure, ZEST uses: (i) Five non-parametric diagnostics: asymmetry, concentration, Gini coefficient, 2nd-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (ii) The exponent n of single--Sersic fits to the 2D surface brightness distributions. To fully exploit the wealth of information while reducing the redundancy present in these diagnostics, ZEST performs a principal component (PC) Analysis. We use a sample of ~56,000 I<24 COSMOS galaxies to show that the first three PCs fully describe the key aspects of the galaxy structure, i.e., to calibrate a three-dimensional classification grid of axis PC_1, PC_2, and PC_3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. (1996). The ZEST classification breaks most of the degeneracy between different galaxy populations that affects morphological classifications based on only some of the diagnostics included in ZEST. As a first application, we present the evolution since z~1 of the Luminosity Functions of COSMOS galaxies of early, disk and irregular galaxies and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find that the LF up to a redshift z=1 is consistent with a pure-luminosity evolution (of about 0.95 magnitudes at z \~0.7). We highlight however two trends, that are in general agreement with a down-sizing scenario for galaxy formation: (1.) A deficit of a factor of about two at z~0.7 of MB>-20.5 structurally--classified early--type galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of irregular galaxies.

Simulations of nearby (0.015 < z < 0.025) SDSS galaxies have been used to reproduce as accurately as possible the appearance that they would have on COSMOS ACS images if they had been observed at z ~ 0.7 and z ~ 1.2. By adding the SDSS galaxies to random locations in the COSMOS images, we simulate the effects of chance superpositions of high redshift galaxies with unrelated foreground or background objects. We have used these simulated images, together with those of real COSMOS galaxies at these same redshifts, to undertake a "blind" morphological classification of galaxies to identify those that appear to be undergoing mergers and thus to estimate the change in merger fraction with redshift. We find that real mergers are harder to recognize at high redshift, and also that the chance superposition of unrelated galaxies often produces the appearance of mergers where in reality none exists. In particular, we estimate that 1.5 - 2.0% of objects randomly added to ACS images are misclassified as mergers due to projection with unrelated objects, and as a result, that 40% of the apparent mergers in COSMOS at z=0.7 are likely to be spurious. We find that the fraction of galaxies undergoing mergers increases as (1+z)^3.8+/-1.2 to z ~ 0.7 and that this trend appears to continue to z = 1.2. Merger candidates at z ~ 0.7 are bluer than the parent population, especially when the statistical effects of the chance projections are accounted for. Merger candidates are more asymmetric than the population as a whole, and are often associated with irregular morphology. Nevertheless, the majority (~60%) of the merger candidates appear to be associated with spiral galaxies although in this case we cannot correct for the effects of chance projections.

We present ZEBRA, the Zurich Extragalactic Bayesian Redshift Analyzer. The current version of ZEBRA combines and extends several of the classical approaches to produce accurate photometric redshifts down to faint magnitudes. In particular, ZEBRA uses the template-fitting approach to produce Maximum Likelihood and Bayesian redshift estimates based on: (1.) An automatic iterative technique to correct the original set of galaxy templates to best represent the SEDs of real galaxies at different redshifts; (2.) A training set of spectroscopic redshifts for a small fraction of the photometric sample; and (3.) An iterative technique for Bayesian redshift estimates, which extracts the full two-dimensional redshift and template probability function for each galaxy. We demonstrate the performance of ZEBRA by applying it to a sample of 866 I_AB <= 22.5 COSMOS galaxies with available u*, B, V, g', r', i', z' and K_s photometry and zCOSMOS spectroscopic redshifts in the range 0 < z < 1.3. Adopting a 5-sigma-clipping that excludes less than 10 galaxies, both the Maximum Likelihood and Bayesian estimates for this sample have an accuracy dz/1+z better than 0.03. Similar accuracies are recovered using mock galaxies.

We study a sample of approximately 16,500 galaxies with I_AB <= 22.5 in the COSMOS field. Structural information on the galaxies is derived by fitting single Sersic models to their two-dimensional surface brightness distributions. We investigate the evolution of the number density of disk galaxies larger than 5 kpc between redshift z~1 and the present epoch. To this end, we use the measurements of the half-light radii to construct, as a function of redshift, the size function of both the total disk galaxy population and of disk galaxies split in four bins of bulge-to-disk ratio. Furthermore, we use a selected sample of roughly 1800 SDSS galaxies to calibrate our results with respect to the local universe. We find that: (i) The number density of disk galaxies with intermediate sizes (r_{1/2}~5-7 kpc) remains nearly constant from z~1 to today. (ii) The number density of the largest disks (r_{1/2}>7 kpc) decreases by a factor of about two out to z~1. (iii) There is a constancy in the number density of large bulgeless disks out to z~1; the deficit of large disks at early epochs seems to arise from a smaller number of bulged disks. Our results indicate that the bulk of the large disk galaxy population has completed its growth by z~1, and support the hypothesis that secular evolution processes produce - or at least add stellar mass to - the bulge components of disk galaxies.