Galaxy evolution in groups and clusters: satellite star formation histories and quenching timescales in a hierarchical Universe

Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFR) and lie on the red-sequence than central (field) galaxies. Using galaxy group/cluster catalogs from SDSS DR7, together with a cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching timescales of satellites of M_star > 5 x 10^9 M_sun at z=0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z~0.5, or ~5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs and quiescent fractions. With this, we constrain the importance and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at M_star < 10^10 M_sun. We then constrain satellite star formation histories, finding a 'delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2-4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of < 0.8 Gyr. These quenching timescales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in satellite quiescent fraction with halo mass arises simply because of satellites quenching in a lower mass group prior to infall (group preprocessing), which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact provides key physical insight into the subhalo abundance matching method.Comment: 25 pages, 13 figures. Accepted for publication in MNRAS, matches published versio

Galaxy evolution near groups and clusters: ejected satellites and the spatial extent of environmental quenching

Galaxies that are several virial radii beyond groups/clusters show preferentially quiescent star formation rates. Using a galaxy group/cluster catalog from the Sloan Digital Sky Survey, together with a cosmological N-body simulation, we examine the origin of this environmental quenching beyond the virial radius. Accounting for the clustering of groups/clusters, we show that central galaxies show enhanced SFR quenching out to 2.5 virial radii beyond groups/clusters, and we demonstrate that this extended environmental enhancement can be explained simply by 'ejected' satellite galaxies that orbit beyond their host halo's virial radius. We show that ejected satellites typically orbit for several Gyr beyond the virial radius before falling back in, and thus they compose up to 40% of all central galaxies near groups/clusters. We show that a model in which ejected satellites experience the same SFR quenching as satellites within a host halo can explain essentially all environmental dependence of galaxy quenching. Furthermore, ejected satellites (continue to) lose significant halo mass, an effect that is potentially observable via gravitational lensing. The SFRs/colors and stellar-to-halo masses of ejected satellites highlight the importance of environmental history and present challenges to models of galaxy occupation that ignore such history.Comment: 15 pages, 9 figures. Accepted for publication in MNRAS, matches published versio

A gravitational lensing explanation for the excess of strong Mg-II absorbers in GRB afterglow spectra

GRB afterglows offer a probe of the intergalactic medium out to high redshift which complements observations along more abundant quasar lines-of-sight. Although both quasars and GRB afterglows should provide a-priori random sight-lines through the intervening IGM, it has been observed that strong Mg-II absorbers are twice as likely to be found along sight-lines toward GRBs. Several proposals to reconcile this discrepancy have been put forward, but none has been found sufficient to explain the magnitude of the effect. In this paper we estimate the effect of gravitational lensing by galaxies and their surrounding mass distributions on the statistics of Mg-II absorption. We find that the multi-band magnification bias could be very strong in the spectroscopic GRB afterglow population and that gravitational lensing can explain the discrepancy in density of absorbers, for plausibly steep luminosity functions. The model makes the prediction that approximately 20%-60% of the spectroscopic afterglow sample (i.e. ~ 5-15 of 26 sources) would have been multiply imaged, and hence result in repeating bursts. We show that despite this large lensing fraction it is likely that none would yet have been identified by chance owing to the finite sky coverage of GRB searches. We predict that continued optical monitoring of the bright GRB afterglow locations in the months and years following the initial decay would lead to identification of lensed GRB afterglows. A confirmation of the lensing hypothesis would allow us to constrain the GRB luminosity function down to otherwise inaccessibly faint levels, with potential consequences for GRB models.Comment: 8 pages, 3 figures. Submitted to MNRAS

The Stellar Mass Components of Galaxies: Comparing Semi-Analytical Models with Observation

We compare the stellar masses of central and satellite galaxies predicted by three independent semianalytical models with observational results obtained from a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. In particular, we compare the stellar mass functions of centrals and satellites, the relation between total stellar mass and halo mass, and the conditional stellar mass functions, which specify the average number of galaxies of stellar mass M_* that reside in a halo of mass M_h. The semi-analytical models only predict the correct stellar masses of central galaxies within a limited mass range and all models fail to reproduce the sharp decline of stellar mass with decreasing halo mass observed at the low mass end. In addition, all models over-predict the number of satellite galaxies by roughly a factor of two. The predicted stellar mass in satellite galaxies can be made to match the data by assuming that a significant fraction of satellite galaxies are tidally stripped and disrupted, giving rise to a population of intra-cluster stars in their host halos. However, the amount of intra-cluster stars thus predicted is too large compared to observation. This suggests that current galaxy formation models still have serious problems in modeling star formation in low-mass halos.Comment: 12 pages, 6 figures, accepted for publication in Ap

Cosmological Constraints from a Combination of Galaxy Clustering and Lensing -- I. Theoretical Framework

We present a new method that simultaneously solves for cosmology and galaxy bias on non-linear scales. The method uses the halo model to analytically describe the (non-linear) matter distribution, and the conditional luminosity function (CLF) to specify the halo occupation statistics. For a given choice of cosmological parameters, this model can be used to predict the galaxy luminosity function, as well as the two-point correlation functions of galaxies, and the galaxy-galaxy lensing signal, both as function of scale and luminosity. In this paper, the first in a series, we present the detailed, analytical model, which we test against mock galaxy redshift surveys constructed from high-resolution numerical $N$-body simulations. We demonstrate that our model, which includes scale-dependence of the halo bias and a proper treatment of halo exclusion, reproduces the 3-dimensional galaxy-galaxy correlation and the galaxy-matter cross-correlation (which can be projected to predict the observables) with an accuracy better than 10 (in most cases 5) percent. Ignoring either of these effects, as is often done, results in systematic errors that easily exceed 40 percent on scales of \sim 1 h^{-1}\Mpc, where the data is typically most accurate. Finally, since the projected correlation functions of galaxies are never obtained by integrating the redshift space correlation function along the line-of-sight out to infinity, simply because the data only cover a finite volume, they are still affected by residual redshift space distortions (RRSDs). Ignoring these, as done in numerous studies in the past, results in systematic errors that easily exceed 20 perent on large scales (r_\rmp \gta 10 h^{-1}\Mpc). We show that it is fairly straightforward to correct for these RRSDs, to an accuracy better than $\sim 2$ percent, using a mildly modified version of the linear Kaiser formalism

Cosmological Constraints from a Combination of Galaxy Clustering and Lensing -- III. Application to SDSS Data

We simultaneously constrain cosmology and galaxy bias using measurements of galaxy abundances, galaxy clustering and galaxy-galaxy lensing taken from the Sloan Digital Sky Survey. We use the conditional luminosity function (which describes the halo occupation statistics as function of galaxy luminosity) combined with the halo model (which describes the non-linear matter field in terms of its halo building blocks) to describe the galaxy-dark matter connection. We explicitly account for residual redshift space distortions in the projected galaxy-galaxy correlation functions, and marginalize over uncertainties in the scale dependence of the halo bias and the detailed structure of dark matter haloes. Under the assumption of a spatially flat, vanilla {\Lambda}CDM cosmology, we focus on constraining the matter density, {\Omega}m, and the normalization of the matter power spectrum, {\sigma}8, and we adopt WMAP7 priors for the spectral index, the Hubble parameter, and the baryon density. We obtain that \Omegam = 0.278_{-0.026}^{+0.023} and {\sigma}8 = 0.763_{-0.049}^{+0.064} (95% CL). These results are robust to uncertainties in the radial number density distribution of satellite galaxies, while allowing for non-Poisson satellite occupation distributions results in a slightly lower value for {\sigma}8 (0.744_{-0.047}^{+0.056}). These constraints are in excellent agreement (at the 1{\sigma} level) with the cosmic microwave background constraints from WMAP. This demonstrates that the use of a realistic and accurate model for galaxy bias, down to the smallest non-linear scales currently observed in galaxy surveys, leads to results perfectly consistent with the vanilla {\Lambda}CDM cosmology.Comment: 21 pages, 9 figures, 5 tables, submitted to MNRA

Satellite Kinematics I: A New Method to Constrain the Halo Mass-Luminosity Relation of Central Galaxies

Satellite kinematics can be used to probe the masses of dark matter haloes of central galaxies. In order to measure the kinematics with sufficient signal-to-noise, one uses the satellite galaxies of a large number of central galaxies stacked according to similar properties (e.g., luminosity). However, in general the relation between the luminosity of a central galaxy and the mass of its host halo will have non-zero scatter. Consequently, this stacking results in combining the kinematics of satellite galaxies in haloes of different masses, which complicates the interpretation of the data. In this paper we present an analytical framework to model satellite kinematics, properly accounting for this scatter and for various selection effects. We show that in the presence of scatter in the halo mass-luminosity relation, the commonly used velocity dispersion of satellite galaxies can not be used to infer a unique halo mass-luminosity relation. In particular, we demonstrate that there is a degeneracy between the mean and the scatter of the halo mass-luminosity relation. We present a new technique that can break this degeneracy, and which involves measuring the velocity dispersions using two different weighting schemes: host-weighting (each central galaxy gets the same weight) and satellite-weighting (each central galaxy gets a weight proportional to its number of satellites). The ratio between the velocity dispersions obtained using these two weighting schemes is a strong function of the scatter in the halo mass-luminosity relation, and can thus be used to infer a unique relation between light and mass from the kinematics of satellite galaxies.Comment: 8 pages, 3 figures, MNRAS submitte

Constraints on the relationship between stellar mass and halo mass at low and high redshift

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

An analytical model for the accretion of dark matter subhalos

An analytical model is developed for the mass function of cold dark matter subhalos at the time of accretion and for the distribution of their accretion times. Our model is based on the model of Zhao et al. (2009) for the median assembly histories of dark matter halos, combined with a simple log-normal distribution to describe the scatter in the main-branch mass at a given time for halos of the same final mass. Our model is simple, and can be used to predict the un-evolved subhalo mass function, the mass function of subhalos accreted at a given time, the accretion-time distribution of subhalos of a given initial mass, and the frequency of major mergers as a function of time. We test our model using high-resolution cosmological $N$-body simulations, and find that our model predictions match the simulation results remarkably well. Finally, we discuss the implications of our model for the evolution of subhalos in their hosts and for the construction of a self-consistent model to link galaxies and dark matter halos at different cosmic times.Comment: 14 pages, 10 figures (caption for figure 10 fixed). Accepted for publication in Ap

Assessment of horse owners’ ability to recognise equine laminitis: A cross-sectional study of 93 veterinary diagnosed cases in Great Britain

The primary objective was to establish whether cases of owner‐suspected laminitis would be confirmed as laminitis by the attending veterinary surgeon. Secondary objectives were to compare owner‐ and veterinary‐reported information from veterinary‐confirmed cases of equine laminitis