Darren J. Croton

Simulations of the formation, evolution and clustering of galaxies and quasars

The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.Numerical simulations are a primary theoretical tool to study the nonlinear gravitational growth of structure in the Universe, and to link the initial conditions of cold dark matter (CDM) cosmogonies to observations of galaxies at the present day. Without direct numerical simulation, the hierarchical build-up of structure with its threedimensional dynamics would be largely inaccessible. Since the dominant mass component, the dark matter, is assumed to consist of weakly interacting elementary particles that interact only gravitationally, such simulations use a set of discrete point particles to represent the collisionless dark matter fluid. This representation as an N-body system is obviously only a coarse approximation, and im-

The WiggleZ Dark Energy Survey: mapping the distance–redshift relation with baryon acoustic oscillations

We present measurements of the baryon acoustic peak at redshifts z = 0.44, 0.6 and 0.73 in the galaxy correlation function of the final dataset of the WiggleZ Dark Energy Survey. We combine our correlation function with lower-redshift measurements from the 6-degree Field Galaxy Survey and Sloan Digital Sky Survey, producing a stacked survey correlation function in which the statistical significance of the detection of the baryon acoustic peak is 4.9-σ relative to a zero-baryon model with no peak. We fit cosmological models to this combined baryon acoustic oscillation (BAO) dataset comprising six distance-redshift data points, and compare the results to similar fits to the latest compilation of supernovae (SNe) and Cosmic Microwave Background (CMB) data. The BAO and SNe datasets produce consistent measurements of the equation-ofstate w of dark energy, when separately combined with the CMB, providing a powerful check for systematic errors in either of these distance probes. Combining all datasets we determine w = 1.03 ± 0.08 for a flat Universe, consistent with a cosmological constant model. Assuming dark energy is a cosmological constant and varying the spatial curvature, we find k = 0.004± 0.006.

The All-wavelength Extended Groth Strip International Survey (AEGIS) Data Sets

In this the first of a series of Letters, we present a panchromatic data set in the Extended Groth Strip region of the sky. Our survey, the All-Wavelength Extended Groth Strip International Survey (AEGIS), aims to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray, GALEX ultraviolet, CFHT/MegaCam Legacy Survey optical, CFHT/CFH12K optical, Hubble Space Telescope/ACS optical and NICMOS near-infrared, Palomar/WIRC near-infrared, Spitzer/IRAC mid-infrared, Spitzer/MIPS far-infrared, and VLA radio continuum. In addition, this region of the sky has been targeted for extensive spectroscopy using the Deep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

The WiggleZ Dark Energy Survey: Joint measurements of the expansion and growth history at z < 1

We perform a joint determination of the distance–redshift relation and cosmic expansion rate at redshifts z = 0.44, 0.6 and 0.73 by combining measurements of the baryon acoustic peak and Alcock–Paczynski distortion from galaxy clustering in the WiggleZ Dark Energy Survey, using a large ensemble of mock catalogues to calculate the covariance between the measurements. We find that D_A(z) = (1205 ± 114, 1380 ± 95, 1534 ± 107) Mpc and H(z) = (82.6 ± 7.8, 87.9 ± 6.1, 97.3 ± 7.0) km s^(−1) Mpc^(−1) at these three redshifts. Further combining our results with other baryon acoustic oscillation and distant supernovae data sets, we use a Monte Carlo Markov Chain technique to determine the evolution of the Hubble parameter H(z) as a stepwise function in nine redshift bins of width Δz = 0.1, also marginalizing over the spatial curvature. Our measurements of H(z), which have precision better than 7 per cent in most redshift bins, are consistent with the expansion history predicted by a cosmological constant dark energy model, in which the expansion rate accelerates at redshift z < 0.7.

The Evolving Luminosity Function of Red Galaxies

We trace the assembly history of red galaxies since z = 1 by measuring their evolving space density with the B-band luminosity function. Our sample of 39,599 red galaxies, selected from 6.96 deg2 of imaging from the NOAO Deep Wide-Field Survey and the Spitzer IRAC Shallow Survey, is an order of magnitude larger, in size and volume, than comparable samples in the literature. We measure a higher space density of z ~ 0.9 red galaxies than some of the recent literature, in part because we account for the faint yet significant galaxy flux that falls outside of our photometric aperture. The B-band luminosity density of red galaxies, which effectively measures the evolution of ~L* galaxies, increases by only 36% ± 13% from z = 0 to z = 1. If red galaxy stellar populations have faded by 1.24 B-band magnitudes since z = 1, the stellar mass contained within the red galaxy population has roughly doubled over the past 8 Gyr. This is consistent with star-forming galaxies being transformed into L* red galaxies after a decline in their star formation rates. In contrast, the evolution of 4L* red galaxies differs only slightly from a model with negligible z < 1 star formation and no galaxy mergers. If this model approximates the luminosity evolution of red galaxy stellar populations, then 80% of the stellar mass contained within todays 4L* red galaxies was already in place at z = 0.7. While red galaxy mergers have been observed, such mergers do not produce rapid growth of 4L* red galaxy stellar masses between z = 1 and the present day.

MERGERS AND BULGE FORMATION IN ΛCDM: WHICH MERGERS MATTER?

We use a suite of semi-empirical models to predict the galaxy-galaxy merger rate and relative contributions to bulge growth as a function of mass (both halo and stellar), redshift, and mass ratio. The models use empirical constraints on the halo occupation distribution, evolved forward in time, to robustly identify where and when galaxy mergers occur. Together with the results of high-resolution merger simulations, this allows us to quantify the relative contributions of mergers with different properties (e.g., mass ratios, gas fractions, redshifts) to the bulge population. We compare with observational constraints, and find good agreement. We also provide useful fitting functions and make public a code to reproduce the predicted merger rates and contributions to bulge mass growth. We identify several robust conclusions. (1) Major mergers dominate the formation and assembly of ~L * bulges and the total spheroid mass density, but minor mergers contribute a non-negligible ~30%. (2) This is mass dependent: bulge formation and assembly is dominated by more minor mergers in lower-mass systems. In higher-mass systems, most bulges originally form in major mergers near ~L *, but assemble in increasingly minor mergers. (3) The minor/major contribution is also morphology dependent: higher B/T systems preferentially form in more major mergers, with B/T roughly tracing the mass ratio of the largest recent merger; lower B/T systems preferentially form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more mergers to reach the same B/T as high-mass systems. Gas-richness dramatically suppresses the absolute efficiency of bulge formation, but does not strongly influence the relative contribution of major versus minor mergers. (5) Absolute merger rates at fixed mass ratio increase with galaxy mass. (6) Predicted merger rates agree well with those observed in pair and morphology-selected samples, but there is evidence that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates also agree with the integrated growth in bulge mass density with cosmic time, but with a factor ~2 uncertainty in both—up to half the bulge mass density could come from non-merger processes. We systematically vary the model assumptions, totaling ~103 model permutations, and quantify the resulting uncertainties. Our conclusions regarding the importance of different mergers for bulge formation are very robust to these changes. The absolute predicted merger rates are systematically uncertain at the factor ~2 level; uncertainties grow at the lowest masses and high redshifts.

The DEEP2 Galaxy Redshift Survey: Color and Luminosity Dependence of Galaxy Clustering at z ∼ 1

We present measurements of the color and luminosity dependence of galaxy clustering at -->z ~ 1 in the DEEP2 Galaxy Redshift Survey. Using volume-limited subsamples in bins of both color and luminosity, we find the following: (1) The clustering dependence is much stronger with color than with luminosity and is as strong with color at -->z ~ 1 as is found locally. We find no dependence of the clustering amplitude on color for galaxies on the red sequence, but a significant dependence on color for galaxies within the blue cloud. (2) For galaxies in the range -->L/L* ~ 0.7–2, a stronger large-scale luminosity dependence is seen for all galaxies than is seen for red and blue galaxies separately. The small-scale clustering amplitude depends significantly on luminosity for blue galaxies, with brighter samples having a stronger rise on scales -->rp h−1 Mpc. (3) Redder galaxies exhibit stronger small-scale redshift-space distortions (fingers of god), and both red and blue populations show large-scale distortions in -->ξ (rp,π) due to coherent infall. (4) While the clustering length, -->r0, increases smoothly with galaxy color (in narrow bins), its power-law exponent, γ, exhibits a sharp jump from the blue cloud to the red sequence. The intermediate-color green galaxy population likely includes transitional galaxies moving from the blue cloud to the red sequence; on large scales green galaxies are as clustered as red galaxies but show infall kinematics and a small-scale correlation slope akin to the blue galaxy population. (5) We compare our results to a semianalytic galaxy formation model applied to the Millennium Run simulation. Differences between the data and the model suggest that in the model star formation is shut down too efficiently in satellite galaxies.

CANDELS: The Progenitors of Compact Quiescent Galaxies at z 2

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

THE DEEP2 GALAXY REDSHIFT SURVEY: THE RELATIONSHIP BETWEEN GALAXY PROPERTIES AND ENVIRONMENT AT z 1

We study the mean environment of galaxies in the DEEP2 Galaxy Redshift Survey as a function of rest-frame color, luminosity, and [OII] 3727u equivalent width. The local galaxy overdensity for > 14,000 galaxies at 0.75 < z < 1.35 is estimated using the projected 3 rd -nearest-neighbor surface density. Of the galaxy properties studied, mean environment is found to depend most strongly on galaxy color; all major features of the correlation between mean overdensity and rest-frame color observed in the local universe were already in place at z ∼ 1. In contrast to local results, we find a substantial slope in the mean dependence of environment on luminosity for blue, star-forming galaxies at z ∼ 1, with brighter blue galaxies being found on average in regions of greater overdensity. We discuss the roles of galaxy clusters and groups in establishing the observed correlations between environment and galaxy properties at high redshift, and we also explore the evidence for a “downsizing of quenching” from z ∼ 1 to z ∼ 0. Our results add weight to existing evidence that the mechanism(s) that result in star-formation quenching are efficient in group environments as well as clusters. This work is the first of its kind at high redshift and represents the first in a series of papers addressing the role of environment in galaxy formation at 0 < z < 1. Subject headings: galaxies:high-redshift, galaxies:evolution, galaxies:statistics, galaxies:fundamental parameters, large-scale structure of universe

Statistical analysis of galaxy surveys ― I. Robust error estimation for two-point clustering statistics

We present a test of different error estimators for two-point clustering statistics, appropriate for present and future large galaxy redshift surveys. Using an ensemble of very large dark matter ACDM N-body simulations, we compare internal error estimators (jackknife and bootstrap) to external ones (Monte Carlo realizations). For three-dimensional clustering statistics, we find that none of the internal error methods investigated is able to reproduce either accurately or robustly the errors of external estimators on 1 to 25 h ―1 Mpc scales. The standard bootstrap overestimates the variance of ξ (s) by ∼40 per cent on all scales probed, but recovers, in a robust fashion, the principal eigenvectors of the underlying covariance matrix. The jackknife returns the correct variance on large scales, but significantly overestimates it on smaller scales. This scale dependence in the jackknife affects the recovered eigenvectors, which tend to disagree on small scales with the external estimates. Our results have important implications for fitting models to galaxy clustering measurements. For example, in a two-parameter fit to the projected correlation function, we find that the standard bootstrap systematically overestimates the 95 per cent confidence interval, while the jackknife method remains biased, but to a lesser extent. Ignoring the systematic bias, the scatter between realizations, for Gaussian statistics, implies that a 2σ confidence interval, as inferred from an internal estimator, corresponds in practice to anything from 1σ to 3σ. By oversampling the subvolumes, we find that it is possible, at least for the cases we consider, to obtain robust bootstrap variances and confidence intervals that agree with external error estimates. Our results are applicable to two-point statistics, like ξ(s) and w p (r p ), measured in large redshift surveys, and show that the interpretation of clustering measurements with internally estimated errors should be treated with caution.