Source author record

Darren J. Croton

Darren J. Croton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2021arXiv

Exploring the relation between dust mass and galaxy properties using Dusty SAGE

We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two categories: ISM dust-poor and ISM dust-rich. We split the ISM dust-poor group into two subgroups: halo dust-rich and dust-poor (the latter contains galaxies that lack dust in both the ISM and halo). Halo dust-rich galaxies have high outflow rates of heated gas and dust and are more massive. We divide ISM dust-rich galaxies based on their specific star formation rate (sSFR) into star-forming and quenched subgroups. At redshift z=0, we find that ISM dust-rich galaxies have a relatively high sSFR, low bulge-to-total (BTT) mass ratio, and high gas metallicity. The high sSFR of ISM dust-rich galaxies allows them to produce dust in the stellar ejecta. Their metal-rich ISM enables dust growth via grain accretion. The opposite is seen in the ISM dust-poor group. Furthermore, ISM dust-rich galaxies are typically late-types, while ISM dust-poor galaxies resemble the early-type population, and we show how their ISM content evolves from being dust-rich to dust-poor. Finally, we investigate dust production from z=3 to z=0 and find that all groups evolve similarly, except for the quenched ISM dust-rich group.

preprint2020arXiv

The origin of dust in galaxies across cosmic time

We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage.

preprint2019arXiv

The 21cm bispectrum during reionization: a tracer of the ionization topology

We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.

preprint2016arXiv

Building disc structure and galaxy properties through angular momentum: The DARK SAGE semi-analytic model

We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find the discs naturally build a pseduobulge-like component. Our main results are focussed on predictions relating to the integrated mass--specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass--spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs.

preprint2016arXiv

Non-linear structure formation in the "Running FLRW" cosmological model

We present a suite of cosmological N-body simulations describing the "Running Friedmann-Lema{ï}tre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends ΛCDM with a time-evolving vacuum density, Λ(z), and time-evolving gravitational Newton's coupling, G(z). In this paper we review the model and introduce the necessary analytical treatment needed to adapt a reference N-body code. Our resulting simulations represent the first realisation of the full growth history of structure in the R-FLRW cosmology into the non-linear regime, and our normalisation choice makes them fully consistent with the latest cosmic microwave background data. The post-processing data products also allow, for the first time, an analysis of the properties of the halo and sub-halo populations. We explore the degeneracies of many statistical observables and discuss the steps needed to break them. Furthermore, we provide a quantitative description of the deviations of R-FLRW from ΛCDM, which could be readily exploited by future cosmological observations to test and further constrain the model.

preprint2016arXiv

Semi-Analytic Galaxy Evolution (SAGE): Model Calibration and Basic Results

This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a minimum set of basic halo properties. In this work we present the baryonic prescriptions implemented in SAGE to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling--radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.

preprint2016arXiv

The growth of disks and bulges during hierarchical galaxy formation. I: fast evolution vs secular processes

We present a theoretical model for the evolution of mass, angular momentum and size of galaxy disks and bulges, and we implement it into the semi-analytic galaxy formation code SAGE. The model follows both secular and violent evolutionary channels, including smooth accretion, disk instabilities, minor and major mergers. We find that the combination of our recipe with hierarchical clustering produces two distinct populations of bulges: merger-driven bulges, akin to classical bulges and ellipticals, and instability-driven bulges, akin to secular (or pseudo-)bulges. The model mostly reproduces the mass-size relation of gaseous and stellar disks, the evolution of the mass-size relation of ellipticals, the Faber-Jackson relation, and the magnitude-colour diagram of classical and secular bulges. The model predicts only a small overlap of merger-driven and instability-driven components in the same galaxy, and predicts different bulge types as a function of galaxy mass and disk fraction. Bulge type also affects the star formation rate and colour at a given luminosity. The model predicts a population of merger-driven red ellipticals that dominate both the low-mass and high-mass ends of the galaxy population, and span all dynamical ages; merger-driven bulges in disk galaxies are dynamically old and do not interfere with subsequent evolution of the star-forming component. Instability-driven bulges dominate the population at intermediate galaxy masses, especially thriving in massive disks. The model green valley is exclusively populated by instability-driven bulge hosts. Through the present implementation the mass accretion history is perceivable in the galaxy structure, morphology and colours.

preprint2016arXiv

The Theoretical Astrophysical Observatory: Cloud-Based Mock Galaxy Catalogues

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.

preprint2015arXiv

Beyond the Halo: Redefining environment with unbound matter in N-body simulations

Approximately half of the matter in the Universe is "unbound" at z = 0, according to N-body simulations such as the Millennium Run. Here, we use the milli-Millennium simulation to examine the distribution of unbound matter in relation to the dark matter halos which host galaxies. We measure the unbound matter within two types of windows, using a halo dependent radius and a fixed radius at several different scales. We also consider the timescales over which a halo can accrete the local unbound matter at z = 2 and z = 0. Finally, we compare the unbound matter to observable properties of galaxies, such as local galaxy count environment and stellar mass. We find that halos at z = 2 can accrete far more of the nearby unbound matter over a Hubble time than halos at z = 0 and that 78% of particles within 5 $R_{vir}$ of a halo at z = 2 will be accreted by z = 0, compared to 36% of particles within 5 $h^{-1}$ Mpc of the halo. We also find that galaxy count environment is closely related to the amount of nearby unbound matter when measured on the same scale.

preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

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

preprint2015arXiv

The Gigaparsec WiggleZ Simulations: Characterising scale dependant bias and associated systematics in growth of structure measurements

We use the Gigaparsec WiggleZ (GiggleZ) simulations to characterise galaxy bias and its scale dependence for a range of redshifts and halo masses in a standard $Λ$LCDM cosmology. Assuming bias converges to a scale independent form at large scales, we develop a phenomenological model which fully expresses the mass and redshift dependence of bias and its scale dependence in real or redshift space. We then use this to illustrate how scale-dependent bias can systematically skew measurements of the growth-rate of cosmic structure obtained from redshift-space distortion measurements. When data is fit only to scales $k_{\rm max}{\le}0.1$ $[h^{-1} \rm{Mpc}]^{-1}$, we find that these effects are significant only for large biases ($b{\gtrsim}3$) at large redshifts ($z{\gtrsim}1$). However, when smaller scales are incorporated ($k_{\rm max}{\le}0.2$ $[h^{-1} \rm{Mpc}]^{-1}$) to increase measurement precision, the combination of reduced statistical uncertainty and increased scale dependent bias can result in highly significant systematics for most large halos across all redshifts. We identify several new interesting aspects of bias, including a significant large-scale bias boost for small halos at low-redshifts due to substructure effects ($\sim$20\% for Milky Way-like systems) and a nearly redshift-independent halo mass (corresponding to a redshift-space bias of ${\sim}1.5$) for which halo bias has little-or-no scale dependence on scales greater than $3$ $[h^{-1} {\rm Mpc}]$. This suggests an optimal strategy of targeting bias ${\sim}{1.5}$ systems for clustering studies which are dominated by systematic uncertainties in how observed halo (or galaxy) distributions map to their underlying mass distribution, such as cosmological measurements of neutrino masses. Code for generating our fitting formula is publicly available at http://gbpoole.github.io/Poole_2014a_code/ (Abridged)

preprint2015arXiv

The WiggleZ Dark Energy Survey: Probing the Epoch of Radiation Domination using Large Scale Structure

We place the most robust constraint to date on the scale of the turnover in the cosmological matter power spectrum using data from the WiggleZ Dark Energy Survey. We find this feature to lie at a scale of $k_0=0.0160^{+0.0041}_{-0.0035}$ [h/Mpc] (68% confidence) for an effective redshift of 0.62 and obtain from this the first-ever turnover-derived distance and cosmology constraints: a measure of the cosmic distance-redshift relation in units of the horizon scale at the redshift of radiation-matter equality (r_H) of D_V(z=0.62)/r_H=18.3 (+6.3/-3.3) and, assuming a prior on the number of extra relativistic degrees of freedom $N_{eff}=3$, constraints on the matter density parameter $Ω_Mh^2=0.136^{+0.026}_{-0.052}$ and on the redshift of matter-radiation equality $z_{eq}=3274^{+631}_{-1260}$. All results are in excellent agreement with the predictions of standard LCDM models. Our constraints on the logarithmic slope of the power spectrum on scales larger than the turnover is bounded in the lower limit with values only as low as -1 allowed, with the prediction of standard LCDM models easily accommodated by our results. Lastly, we generate forecasts for the achievable precision of future surveys at constraining $k_0$, $Ω_Mh^2$, $z_{eq}$ and $N_{eff}$. We find that BOSS should substantially improve upon the WiggleZ turnover constraint, reaching a precision on $k_0$ of $\pm$9% (68% confidence), translating to precisions on $Ω_Mh^2$ and $z_{eq}$ of $\pm$10% (assuming a prior $N_{eff}=3$) and on $N_{eff}$ of (+78/-56)% (assuming a prior $Ω_Mh^2=0.135$). This is sufficient precision to sharpen the constraints on $N_{eff}$ from WMAP, particularly in its upper limit. For Euclid, we find corresponding attainable precisions on $(k_0, Ω_Mh^2, N_eff)$ of (3,4,+17/-21)%. This represents a precision approaching our forecasts for the Planck Surveyor.

preprint2015arXiv

Thinking outside the halo: Tracing the large-scale distribution of diffuse cosmic metals with semi-analytic models

With the installation of the Cosmic Origins Spectrograph on the Hubble Space Telescope, measurements of the metal content of the low redshift intergalactic medium (IGM) are now available. Using a new grid-based model for diffuse gas coupled to the SAGE semi-analytic model of galaxy formation, we examine the impact of supernova feedback on the pollution of the IGM. We consider different assumptions for the reheating and ejection of gas by supernovae and their dependence on galaxy circular velocity and gas surface density. Where metals are present, we find the most likely metallicity to be $-1.5 < $log$_{10}$(Z/Z$_{\odot}$)$< -1.0$ at $z = 0$, consistent with both observations and more sophisticated hydrodynamic simulations. Our model predicts that the regions of the IGM with the highest metallicities will be near galaxies with M$_{\star} \sim 10^{10.5}h^{-1}$M$_{\odot}$ and in environments of densities $\sim 10 \times$ the mean. We also find that 90% of IGM metals at $z = 0$ are ejected by galaxies with stellar masses less than $10^{10.33}h^{-1}$M$_{\odot}$.

preprint2014arXiv

The WiggleZ Dark Energy Survey: Improved Distance Measurements to z = 1 with Reconstruction of the Baryonic Acoustic Feature

We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in the redshift range 0.2<$z$<1, and in all three find improvement after reconstruction in the detection of the baryonic acoustic feature and its usage as a standard ruler. We measure model independent distance measures $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ of 1716 $\pm$ 83 Mpc, 2221 $\pm$ 101 Mpc, 2516 $\pm$ 86 Mpc (68% CL) at effective redshifts z = 0.44, 0.6, 0.73, respectively, where $D_{\mathrm V}$ is the volume-average-distance, and $r_{\mathrm s}$ is the sound horizon at the end of the baryon drag epoch. These significantly improved 4.8, 4.5 and 3.4 percent accuracy measurements are equivalent to those expected from surveys with up to 2.5 times the volume of WiggleZ. These measurements are fully consistent with cosmologies allowed by the analyses of the Planck Collaboration and the Sloan Digital Sky Survey.We provide the $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ posterior probability distributions and their covariances. When combining these measurements with temperature fluctuations measurements of Planck, the polarization of WMAP9, and the 6dF Galaxy Survey baryonic acoustic feature, we do not detect deviations from a flat LCDM model. Assuming this model we constrain the current expansion rate to $H_0$ = 67.15 $\pm$ 0.98 kms$^{-1}$Mpc$^{-1}$. Allowing the equation of state of dark energy to vary we obtain $w_\mathrm{DE}$ = -1.080 $\pm$ 0.135. When assuming a curved LCDM model we obtain a curvature value of $Ω_{\mathrm K}$ = -0.0043 $\pm$ 0.0047.

preprint2014arXiv

Where do galaxies end? Comparing measurement techniques of hydrodynamic-simulation galaxies' integrated properties

Using the suite of high-resolution zoom re-simulations of individual haloes by Martig et al., and the large-scale simulation \emph{MassiveBlack-II}, we examine the differences in measured galaxy properties from techniques with various aperture definitions of where galaxies end. We perform techniques popular in the literature and present a new technique of our own, where the aperture radius is based on the baryonic mass profiles of simulated (sub)haloes. For the average Milky-Way-mass system, we find the two most popular techniques in the literature return differences of order 30 per cent for stellar mass, a factor of 3 for gas mass, 40 per cent for star formation rate, and factors of several for gas accretion and ejection rates. Individual cases can show variations greater than this, with the severity dependent on the concentration of a given system. The average difference in integrated properties for a more general galaxy population are not as striking, but are still significant for stellar and gas mass, especially for optical-limit apertures. The large differences that can occur are problematic for comparing results from various publications. We stress the importance of both defining and justifying a technique choice and discourage using popular apertures that use an exact fraction of the virial radius, due to the unignorable variation in galaxy-to-(sub)halo size. Finally, we note that technique choice does not greatly affect simulated galaxies from lying within the scatter of observed scaling relations, but it can alter the derived best-fit slope for the Kennicutt-Schmidt relation.

preprint2013arXiv

Measures of Galaxy Environment - III. Difficulties in identifying proto-clusters at z ~ 2

Galaxy environment is frequently discussed, but inconsistently defined. It is especially difficult to measure at high redshift where only photometric redshifts are available. With a focus on early forming proto-clusters, we use a semi-analytical model of galaxy formation to show how the environment measurement around high redshift galaxies is sensitive to both scale and metric, as well as to cluster viewing angle, evolutionary state, and the availability of either spectroscopic or photometric data. We use two types of environment metrics (nearest neighbour and fixed aperture) at a range of scales on simulated high-z clusters to see how "observed" overdensities compare to "real" overdensities. We also "observationally" identify z = 2 proto-cluster candidates in our model and track the growth histories of their parent halos through time, considering in particular their final state at z = 0. Although the measured environment of early forming clusters is critically dependent on all of the above effects (and in particular the viewing angle), we show that such clusters are very likely (< 90%) to remain overdense at z = 0, although many will no longer be among the most massive. Object to object comparisons using different methodologies and different data, however, require much more caution.

preprint2013arXiv

Measures of galaxy environment -- II. Rank-ordered mark correlations

We analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+δ$, while such a small dependence would be difficult to detect by traditional methods. We then show that rank ordering the marks and using the rank as a weight is a simple way of comparing the correlation signals for different marks. With this we quantify to what extent fixed-aperture overdensities are sensitive to large-scale halo environments, nearest-neighbor overdensities are sensitive to small-scale environments within haloes, and colour is a better tracer of overdensity than is luminosity.

preprint2013arXiv

The simplest model of galaxy formation I: A formation history model of galaxy stellar mass growth

We introduce a simple model to self-consistently connect the growth of galaxies to the formation history of their host dark matter haloes. Our model is defined by two simple functions: the "baryonic growth function" which controls the rate at which new baryonic material is made available for star formation, and the "physics function" which controls the efficiency with which this material is converted into stars. Using simple, phenomenologically motivated forms for both functions that depend only on a single halo property, we demonstrate the model's ability to reproduce the z=0 red and blue stellar mass functions. Furthermore, by adding redshift as a second input variable to the physics function we show that the reproduction of the global stellar mass function out to z=3 is improved. We conclude by discussing the general utility of our new model, highlighting its usefulness for creating mock galaxy samples which have a number of key advantages over those generated by other techniques.

preprint2013arXiv

The Stellar Masses of ~40,000 UV Selected Galaxies from the WiggleZ Survey at 0.3<z<1.0: Analogues of Lyman Break Galaxies?

We characterise the stellar masses and star formation rates in a sample of almost 40000 spectroscopically confirmed UV luminous galaxies at 0.3<z<1.0 selected from within the WiggleZ Dark Energy Survey. In particular, we match this UV bright population to wide-field infrared surveys such as the near infrared UKIDSS LAS and the mid infrared WISE All-Sky Survey. We find that ~30% of the UV luminous WiggleZ galaxies are detected at >5sigma in the UKIDSS-LAS at all redshifts. An even more luminous subset of 15% are also detected in the WISE 3.4 and 4.6um bands. We compute stellar masses for this very large sample of extremely blue galaxies and quantify the sensitivity of the stellar mass estimates to various assumptions made during the SED fitting. The median stellar masses are log10(M*/M0)=9.6\pm0.7, 10.2\pm0.5 and 10.4\pm0.4 for the IR-undetected, UKIDSS detected and UKIDSS+WISE detected galaxies respectively. We demonstrate that the inclusion of NIR photometry can lead to tighter constraints on the stellar masses. The mass estimates are found to be most sensitive to the inclusion of secondary bursts of star formation as well as changes in the stellar population synthesis models, both of which can lead to median discrepancies of the order of 0.3dex in the stellar masses. We find that the best-fit M/LK is significantly lower (by ~0.4 dex) than that predicted by simple optical colour based estimators, in particular for the bluer galaxies with younger best-fit ages. The WiggleZ galaxies have star formation rates of 3-10 M0/yr and mostly lie at the upper end of the main sequence of star-forming galaxies at these redshifts. Their rest-frame UV luminosities and stellar masses are comparable to both local compact UV-luminous galaxies as well as Lyman break galaxies at z~2-3.(abridged)

preprint2012arXiv

A diversity of progenitors and histories for isolated spiral galaxies

We analyze a suite of 33 cosmological simulations of the evolution of Milky Way-mass galaxies in low-density environments. Our sample spans a broad range of Hubble types at z=0, from nearly bulgeless disks to bulge-dominated galaxies. Despite the fact that a large fraction of the bulge is typically in place by z=1, we find no significant correlation between the morphology at z=1 and at z=0. The z=1 progenitors of disk galaxies span a range of morphologies, including smooth disks, unstable disks, interacting galaxies and bulge-dominated systems. By z=0.5, spiral arms and bars are largely in place and the progenitor morphology is correlated with the final morphology. We next focus on late-type galaxies with a bulge-to-total ratio B/T<0.3 at z=0. These show a correlation between B/T at z=0 and the mass ratio of the largest merger at z<2, as well as with the gas accretion rate at z>1. We find that the galaxies with the lowest B/T tend to have a quiet baryon input history, with no major mergers at z<2, and with a low and constant gas accretion rate that keeps a stable angular-momentum direction. More violent merger or gas accretion histories lead to galaxies with more prominent bulges. Most disk galaxies have a bulge Sersic index n<2. The galaxies with the highest bulge Sersic index tend to have histories of intense gas accretion and disk instability rather than active mergers.

preprint2012arXiv

CANDELS: The progenitors of compact quiescent galaxies at z~2

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

preprint2012arXiv

Constraining the last 7 billion years of galaxy evolution in semi-analytic models

We investigate the ability of the Croton et al. (2006) semi-analytic model to reproduce the evolution of observed galaxies across the final 7 billion years of cosmic history. Using Monte-Carlo Markov Chain techniques we explore the available parameter space to produce a model which attempts to achieve a statistically accurate fit to the observed stellar mass function at z=0 and z~0.8, as well as the local black hole-bulge relation. We find that in order to be successful we are required to push supernova feedback efficiencies to extreme limits which are, in some cases, unjustified by current observations. This leads us to the conclusion that the current model may be incomplete. Using the posterior probability distributions provided by our fitting, as well as the qualitative details of our produced stellar mass functions, we suggest that any future model improvements must act to preferentially bolster star formation efficiency in the most massive halos at high redshift.

preprint2012arXiv

Measures of Galaxy Environment - I. What is "Environment"?

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

preprint2012arXiv

Predictions for ASKAP Neutral Hydrogen Surveys

The Australian Square Kilometer Array Pathfinder (ASKAP) will revolutionise our knowledge of gas-rich galaxies in the Universe. Here we present predictions for two proposed extragalactic ASKAP neutral hydrogen (HI) emission-line surveys, based on semi-analytic models applied to cosmological N-body simulations. The ASKAP HI All-Sky Survey, known as WALLABY, is a shallow 3 Pi survey (z = 0 - 0.26) which will probe the mass and dynamics of over 600,000 galaxies. A much deeper small-area HI survey, called DINGO, aims to trace the evolution of HI from z = 0 - 0.43, a cosmological volume of 40 million Mpc^3, detecting potentially 100,000 galaxies. The high-sensitivity 30 antenna ASKAP core (diameter ~2 km) will provide an angular resolution of 30 arcsec (at z=0). Our simulations show that the majority of galaxies detected in WALLABY (87.5%) will be resolved. About 5000 galaxies will be well resolved, i.e. more than five beams (2.5 arcmin) across the major axis, enabling kinematic studies of their gaseous disks. This number would rise to 160,000 galaxies if all 36 ASKAP antennas could be used; the additional six antennas provide baselines up to 6 km, resulting in an angular resolution of 10 arcsec. For DINGO this increased resolution is highly desirable to minimise source confusion; reducing confusion rates from a maximum of 10% of sources at the survey edge to 3%. We estimate that the sources detected by WALLABY and DINGO will span four orders of magnitude in total halo mass (from 10^{11} to 10^{15} Msol) and nearly seven orders of magnitude in stellar mass (from 10^{5} to 10^{12} Msol), allowing us to investigate the process of galaxy formation across the last four billion years.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

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

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

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

preprint2011arXiv

Statistical Analysis of Galaxy Surveys-IV: An objective way to quantify the impact of superstructures on galaxy clustering statistics

For galaxy clustering to provide robust constraints on cosmological parameters and galaxy formation models, it is essential to make reliable estimates of the errors on clustering measurements. We present a new technique, based on a spatial Jackknife (JK) resampling, which provides an objective way to estimate errors on clustering statistics. Our approach allows us to set the appropriate size for the Jackknife subsamples. The method also provides a means to assess the impact of individual regions on the measured clustering, and thereby to establish whether or not a given galaxy catalogue is dominated by one or several large structures, preventing it to be considered as a "fair sample". We apply this methodology to the two- and three-point correlation functions measured from a volume limited sample of M* galaxies drawn from data release seven of the Sloan Digital Sky Survey (SDSS). The frequency of jackknife subsample outliers in the data is shown to be consistent with that seen in large N-body simulations of clustering in the cosmological constant plus cold dark matter cosmology. We also present a comparison of the three-point correlation function in SDSS and 2dFGRS using this approach and find consistent measurements between the two samples.

preprint2011arXiv

The mid-life crisis of the Milky Way and M31

Upcoming next generation galactic surveys, such as GAIA and HERMES, will deliver unprecedented detail about the structure and make-up of our Galaxy, the Milky Way, and promise to radically improve our understanding of it. However, to benefit our broader knowledge of galaxy formation and evolution we first need to quantify how typical the Galaxy is with respect to other galaxies of its type. Through modeling and comparison with a large sample of galaxies drawn from the Sloan Digital Sky Survey and Galaxy Zoo, we provide tentative yet tantalizing evidence to show that both the Milky Way and nearby M31 are undergoing a critical transformation of their global properties. Both appear to possess attributes that are consistent with galaxies midway between the distinct blue and red bimodal color populations. In extragalactic surveys, such `green valley' galaxies are transition objects whose star formation typically will have all but extinguished in less than ~5 Gyrs. This finding reveals the possible future of our own galactic home, and opens a new window of opportunity to study such galactic transformations up close.

preprint2010arXiv

Absence of Evidence Is Not Evidence of Absence: The Color-Density Relation at Fixed Stellar Mass Persists to z ~ 1

We use data drawn from the DEEP2 Galaxy Redshift Survey to investigate the relationship between local galaxy density, stellar mass, and rest-frame galaxy color. At z ~ 0.9, we find that the shape of the stellar mass function at the high-mass (log (M*/Msun) > 10.1) end depends on the local environment, with high-density regions favoring more massive systems. Accounting for this stellar mass-environment relation (i.e., working at fixed stellar mass), we find a significant color-density relation for galaxies with 10.6 < log(M*/Msun) < 11.1 and 0.75 < z < 0.95. This result is shown to be robust to variations in the sample selection and to extend to even lower masses (down to log(M*/Msun) ~ 10.4). We conclude by discussing our results in comparison to recent works in the literature, which report no significant correlation between galaxy properties and environment at fixed stellar mass for the same redshift and stellar mass domain. The non-detection of environmental dependence found in other data sets is largely attributable to their smaller samples size and lower sampling density, as well as systematic effects such as inaccurate redshifts and biased analysis techniques. Ultimately, our results based on DEEP2 data illustrate that the evolutionary state of a galaxy at z ~ 1 is not exclusively determined by the stellar mass of the galaxy. Instead, we show that local environment appears to play a distinct role in the transformation of galaxy properties at z > 1.

preprint2010arXiv

Where do Wet, Dry, and Mixed Galaxy Mergers Occur? A Study of the Environments of Close Galaxy Pairs in the DEEP2 Galaxy Redshift Survey

We study the environment of wet, dry, and mixed galaxy mergers at 0.75 < z < 1.2 using close pairs in the DEEP2 Galaxy Redshift Survey. We find that the typical environment of mixed and dry merger candidates is denser than that of wet mergers, mostly due to the color-density relation. While the galaxy companion rate (Nc) is observed to increase with overdensity, using N-body simulations we find that the fraction of pairs that will eventually merge decreases with the local density, predominantly because interlopers are more common in dense environments. After taking into account the merger probability of pairs as a function of local density, we find only marginal environment dependence of the fractional merger rate for wet mergers over the redshift range we have probed. On the other hand, the fractional dry merger rate increases rapidly with local density due to the increased population of red galaxies in dense environments. We also find that the environment distribution of K+A galaxies is similar to that of wet mergers alone and of wet+mixed mergers, suggesting a possible connection between K+A galaxies and wet and/or wet+mixed mergers. We conclude that, as early as z ~ 1, high-density regions are the preferred environment in which dry mergers occur, and that present-day red-sequence galaxies in overdense environments have, on average, undergone 1.2+-0.3 dry mergers since this time, accounting for (38+-10)% of their mass accretion in the last 8 billion years. Our findings suggest that dry mergers are crucial in the mass-assembly of massive red galaxies in dense environments, such as Brightest Cluster Galaxies (BCGs) in galaxy groups and clusters.

preprint2009arXiv

A simple model to link the properties of quasars to the properties of dark matter halos out to high redshift

We present a simple model of how quasars occupy dark matter halos from z=0 to z=5 using the observed mBH-sigma relation and quasar luminosity functions. This provides a way for observers to statistically infer host halo masses for quasar observations using luminosity and redshift alone. Our model is deliberately simple and sidesteps any need to explicitly describe the physics. In spite of its simplicity, the model reproduces many key observations and has predictive power: 1) model quasars have the correct luminosity function (by construction) and spatial clustering (by consequence); 2) we predict high redshift quasars of a given luminosity live in less massive dark matter halos than the same luminosity quasars at low redshifts; 3) we predict a factor of ~5 more 10^8.5Msun black holes at z~2 than is currently observed; 4) we predict a factor of ~20 evolution in the amplitude of the mBH-Mhalo relation between z=5 and the present day; 5) we expect luminosity dependent quasar lifetimes of between tQ~10^(7-8)yr, but which may become as short as 10^(5-6)yr for quasars brighter than L*; 6) while little luminosity dependent clustering evolution is expected at z<1, increasingly strong evolution is predicted for L>L* quasars at higher redshifts. These last two results arise from the narrowing distribution of halo masses that quasars occupy as the Universe ages. We also deconstruct both "downsizing" and "upsizing" trends predicted by the model at different redshifts and space densities. Importantly, this work illustrates how current observations cannot distinguish between more complicated physically motivated quasar models and our simple phenomenological approach. It highlights the opportunities such methodologies provide.

preprint2008arXiv

Red Galaxy Growth and the Halo Occupation Distribution

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, using the observed luminosity function and clustering of 40,696 0.2<z<1.0 red galaxies in Bootes. Half of 10^{11.9} Msun/h halos host a red central galaxy, and this fraction increases with increasing halo mass. We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that the stellar mass contained within the red population has doubled since z=1, with the stellar mass within red satellite galaxies tripling over this redshift range. In cluster mass halos most of the stellar mass resides within satellite galaxies and the intra-cluster light, with a minority of the stellar mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are proportional to halo mass to the power of a third. We thus conclude that halo mergers do not always lead to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, the stellar masses of their central galaxies typically grow by only 30%.

preprint2008arXiv

Statistical Analysis of Galaxy Surveys - I. Robust error estimation for 2-point clustering statistics

We present a test of different error estimators for 2-point clustering statistics, appropriate for present and future large galaxy redshift surveys. Using an ensemble of very large dark matter LambdaCDM N-body simulations, we compare internal error estimators (jackknife and bootstrap) to external ones (Monte-Carlo realizations). For 3-dimensional clustering statistics, we find that none of the internal error methods investigated are able to reproduce neither accurately nor robustly the errors of external estimators on 1 to 25 Mpc/h scales. The standard bootstrap overestimates the variance of xi(s) by ~40% 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 the use of galaxy clustering in placing constraints on cosmological parameters. For example, in a 2-parameter fit to the projected correlation function, we find that the standard bootstrap systematically overestimates the 95% confidence interval, while the jackknife method remains biased, but to a lesser extent. The scatter we find between realizations, for Gaussian statistics, implies that a 2-sigma confidence interval, as inferred from an internal estimator, could correspond in practice to anything from 1-sigma to 3-sigma. Finally, by an oversampling of sub-volumes, it is possible to obtain bootstrap variances and confidence intervals that agree with external error estimates, but it is not clear if this prescription will work for a general case.

preprint2007arXiv

The Stellar Mass Tully-Fisher Relation to z=1.2 from AEGIS

We combine newly measured rotation velocities, velocity dispersions, and stellar masses to construct stellar mass Tully-Fisher relations (M*TFRs) for 544 galaxies with strong emission lines at 0.1<z<1.2 from the All Wavelength Extended Groth Strip International Survey (AEGIS) and the Deep Extragalactic Evolutionary Probe 2 Survey (DEEP2). The conventional M*TFR using only rotation velocity (Vrot) shows large scatter (~1.5 dex in velocity). The scatter and residuals are correlated with morphology in the sense that disturbed, compact, and major merger galaxies have lower velocities for their masses. We construct an M*TFR using the kinematic estimator S_0.5 which is defined as sqrt(0.5Vrot^2 + sigma_g^2) and accounts for disordered or non-circular motions through the gas velocity dispersion (sigma_g). The new M*TFR, termed S_0.5/M*TFR, is remarkably tight over 0.1<z<1.2 with no detectable evolution of its intercept or slope with redshift. The average best fit relation has 0.47 dex scatter in stellar mass, corresponding to ~1.2 'magnitudes,' assuming a constant mass-to-light ratio. Interestingly, the S_0.5/M*TFR is consistent with the absorption-line based stellar mass Faber-Jackson relation for nearby elliptical galaxies in terms of slope and intercept, which might suggest a physical connection between the two relations.

preprint2006arXiv

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 39599 red galaxies, selected from 6.96 square degrees of imaging from the NOAO Deep Wide-Field and Spitzer IRAC Shallow surveys, 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 which 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 percent 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 star formation and no galaxy mergers since z=1. If this model approximates the luminosity evolution of red galaxy stellar populations, then 80 percent of the stellar mass contained within today's 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.

preprint2005arXiv

The 2dF Galaxy Redshift Survey: luminosity functions by density environment and galaxy type

We use the 2dF Galaxy Redshift Survey to measure the dependence of the bJ-band galaxy luminosity function on large-scale environment, defined by density contrast in spheres of radius 8h-1Mpc, and on spectral type, determined from principal component analysis. We find that the galaxy populations at both extremes of density differ significantly from that at the mean density. The population in voids is dominated by late types and shows, relative to the mean, a deficit of galaxies that becomes increasingly pronounced at magnitudes brighter than M_bJ-5log10h <-18.5. In contrast, cluster regions have a relative excess of very bright early-type galaxies with M_bJ-5log10h < -21. Differences in the mid to faint-end population between environments are significant: at M_bJ-5log10h=-18 early and late-type cluster galaxies show comparable abundances, whereas in voids the late types dominate by almost an order of magnitude. We find that the luminosity functions measured in all density environments, from voids to clusters, can be approximated by Schechter functions with parameters that vary smoothly with local density, but in a fashion which differs strikingly for early and late-type galaxies. These observed variations, combined with our finding that the faint-end slope of the overall luminosity function depends at most weakly on density environment, may prove to be a significant challenge for models of galaxy formation.