Source author record

Michael J. Hudson

Michael J. Hudson 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

47works
5topics
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

47 published item(s)

preprint2026arXiv

JWST Reveals a Candidate Jellyfish Galaxy at z=1.156

We report the discovery of COSMOS2020-635829 as a candidate jellyfish galaxy undergoing ram pressure stripping in a (proto)cluster at $z > 1$. High-resolution imaging from the James Webb Space Telescope reveals a symmetric stellar disk coupled to a unilateral tail of star-forming knots to the south. Using Gemini GMOS IFU observations, we show that these extra-planar continuum sources are embedded within an ionized gas tail that is kinematically connected to the disk of COSMOS2020-635829. If confirmed, this represents the highest-redshift discovery of a ram pressure stripped ionized gas tail. The tail sources are characterized by extremely young stellar populations ($\lesssim 100\,\mathrm{Myr}$), have stellar masses of ${\sim}10^8\,\mathrm{M_\odot}$, and star formation rates of $0.1\text{--}1\,\mathrm{M_\odot\,yr^{-1}}$. This work shows that ram pressure stripping can potentially perturb group and cluster galaxies at $z > 1$ and may contribute to environmental quenching even near Cosmic Noon.

preprint2022arXiv

Star formation characteristics of CNN-identified post-mergers in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

The importance of the post-merger epoch in galaxy evolution has been well-documented, but post-mergers are notoriously difficult to identify. While the features induced by mergers can sometimes be distinctive, they are frequently missed by visual inspection. In addition, visual classification efforts are highly inefficient because of the inherent rarity of post-mergers (~1% in the low-redshift Universe), and non-parametric statistical merger selection methods do not account for the diversity of post-mergers or the environments in which they appear. To address these issues, we deploy a convolutional neural network (CNN) which has been trained and evaluated on realistic mock observations of simulated galaxies from the IllustrisTNG simulations, to galaxy images from the Canada France Imaging Survey (CFIS), which is part of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS). We present the characteristics of the galaxies with the highest CNN-predicted post-merger certainties, as well as a visually confirmed subset of 699 post-mergers. We find that post-mergers with high CNN merger probabilities (p(x)>0.8) have an average star formation rate that is 0.1 dex higher than a mass- and redshift-matched control sample. The SFR enhancement is even greater in the visually confirmed post-merger sample, a factor of two higher than the control sample.

preprint2021arXiv

Assessing the accuracy of cosmological parameters estimated from velocity -- density comparisons via simulations

A promising method for measuring the cosmological parameter combination fsigma_8 is to compare observed peculiar velocities with peculiar velocities predicted from a galaxy density field using perturbation theory. We use N-body simulations and semi-analytic galaxy formation models to quantify the accuracy and precision of this method. Specifically, we examine a number of technical aspects, including the optimal smoothing length applied to the density field, the use of dark matter halos or galaxies as tracers of the density field, the effect of noise in the halo mass estimates or in the stellar-to-halo mass relation, and the effect of finite survey volumes. We find that for a Gaussian smoothing of 4 Mpc/h, the method has only small systematic biases at the level of 5%. Cosmic variance affects current measurements at the 5% level due to the volume of current redshift data sets.

preprint2021arXiv

Evolution of subhalo orbits in a smoothly-growing host halo potential

The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other relevant mechanisms, we study subhalo orbital evolution using numerical calculations in which subhaloes are modelled as massless particles orbiting in a time-varying spherical potential. We find that the radial action of the subhalo orbit decreases over the first few orbits, indicating that the response to the growth of the host halo is not adiabatic during this phase. The subhalo orbits can shrink by a factor of $\sim$1.5 in this phase. Subsequently, the radial action is well conserved and orbital contraction slows down. We propose a model accurately describing the orbital evolution. Given these results, we consider the spatial distribution of the population of subhaloes identified in high-resolution cosmological simulations. We find that it is consistent with this population having been accreted at z < 3, indicating that any subhaloes accreted earlier are unresolved in the simulations. We also discuss tidal stripping as a formation scenario for NGC1052-DF2, an ultra diffuse galaxy significantly lacking DM, and find that its expected DM mass could be consistent with observational constraints if its progenitor was accreted early enough, z > 1.5, although it should still be a relatively rare object.

preprint2020arXiv

How dark are filaments in the cosmic web?

The cold dark matter model predicts that dark matter haloes are connected by filaments. Direct measurements of the masses and structure of these filaments are difficult, but recently several studies have detected these dark-matter-dominated filaments using weak lensing. Here we study the efficiency of galaxy formation within the filaments by measuring their total mass-to-light ratios and stellar mass fractions. Specifically, we stack pairs of Luminous Red Galaxies (LRGs) with a typical separation on the sky of $8 h^{-1}$ Mpc. We stack background galaxy shapes around pairs to obtain mass maps through weak lensing, and we stack galaxies from the Sloan Digital Sky Survey (SDSS) to obtain maps of light and stellar mass. To isolate the signal from the filament, we construct two matched catalogues of physical and non-physical (projected) LRG pairs, with the same distributions of redshift and separation. We then subtract the two stacked maps. Using LRG pair samples from the BOSS survey at two different redshifts, we find that the evolution of the mass in filament is consistent with the predictions from perturbation theory. The filaments are not entirely dark: their mass-to-light ratios ($M/L = 351\pm137$ in solar units in the $r$-band) and stellar mass fractions ($M_{\rm stellar}/M = 0.0073\pm0.0030$) are consistent with the cosmic values (and with their redshift evolutions)

preprint2020arXiv

Joint analysis of 6dFGS and SDSS peculiar velocities for the growth rate of cosmic structure and tests of gravity

Measurement of peculiar velocities by combining redshifts and distance indicators is a powerful way to measure the growth rate of cosmic structure and test theories of gravity at low redshift. Here we constrain the growth rate of structure by comparing observed Fundamental Plane peculiar velocities for 15894 galaxies from the 6dF Galaxy Survey (6dFGS) and Sloan Digital Sky Survey (SDSS) with predicted velocities and densities from the 2M$++$ redshift survey. We measure the velocity scale parameter $β\equiv {Ω_m^γ}/b = 0.372^{+0.034}_{-0.050}$ and $0.314^{+0.031}_{-0.047}$ for 6dFGS and SDSS respectively, where $Ω_m$ is the mass density parameter, $γ$ is the growth index, and $b$ is the bias parameter normalized to the characteristic luminosity of galaxies, $L^*$. Combining 6dFGS and SDSS we obtain $β= 0.341\pm0.024$, implying that the amplitude of the product of the growth rate and the mass fluctuation amplitude is $fσ_8 = 0.338\pm0.027$ at an effective redshift $z=0.035$. Adopting $Ω_m = 0.315\pm0.007$ as favoured by Planck and using $γ=6/11$ for General Relativity and $γ=11/16$ for DGP gravity, we get $S_8(z=0) = σ_8 \sqrt{Ω_m/0.3} =0.637 \pm 0.054$ and $0.741\pm0.062$ for GR and DGP respectively. This measurement agrees with other low-redshift probes of large scale structure but deviates by more than $3σ$ from the latest Planck CMB measurement. Our results favour values of the growth index $γ> 6/11$ or a Hubble constant $H_0 > 70$\,km\,s$^{-1}$\,Mpc$^{-1}$ or a fluctuation amplitude $σ_8 < 0.8$ or some combination of these. Imminent redshift surveys such as Taipan, DESI, WALLABY, and SKA1-MID will help to resolve this tension by measuring the growth rate of cosmic structure to 1\% in the redshift range $0 < z < 1$.

preprint2019arXiv

Neural physical engines for inferring the halo mass distribution function

An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques.

preprint2016arXiv

Satellite quenching timescales in clusters from projected phase space measurements matched to simulated orbits

We measure the star formation quenching efficiency and timescale in cluster environments. Our method uses N-body simulations to estimate the probability distribution of possible orbits for a sample of observed SDSS galaxies in and around clusters based on their position and velocity offsets from their host cluster. We study the relationship between their star formation rates and their likely orbital histories via a simple model in which star formation is quenched once a delay time after infall has elapsed. Our orbit library method is designed to isolate the environmental effect on the star formation rate due to a galaxy's present-day host cluster from `pre-processing' in previous group hosts. We find that quenching of satellite galaxies of all stellar masses in our sample ($10^{9}-10^{11.5}\,{\rm M}_\odot$) by massive ($> 10^{13}\,{\rm M}_\odot$) clusters is essentially $100$ per cent efficient. Our fits show that all galaxies quench on their first infall, approximately at or within a Gyr of their first pericentric passage. There is little variation in the onset of quenching from galaxy-to-galaxy: the spread in this time is at most $\sim 2$ Gyr at fixed $M_*$. Higher mass satellites quench earlier, with very little dependence on host cluster mass in the range probed by our sample.

preprint2015arXiv

Beyond Sérsic + exponential disc morphologies in the Coma Cluster

[abridged] We explore the diversity of internal galaxy structures in the Coma cluster across a wide range of luminosities ($-17$\,$>$\,$M_g$\,$>$\,$-22$) and cluster-centric radii ($0$\,$<$\,$r_{\rm{cluster}}$\,$<$\,1.3 $r_{200}$) through analysis of deep Canada-France-Hawaii Telescope $i$ band imaging. We present 2D multi-component decomposition via GALFIT, encompassing a wide range of candidate model morphologies with up to three photometric components. Particular focus is placed on early-type galaxies with outer discs (i.e. S0s), and deviations from simple (`unbroken') exponential discs. Rigorous filtering ensures that each model component provides a statistically significant improvement to the goodness-of-fit. The majority of Coma cluster members in our sample (478 of 631) are reliably fitted by symmetric structural models. Of these, 134 ($28\%$) are single Sérsic objects, 143 ($30\%$) are well-described by 2 component structures, while 201 ($42\%$) require more complex models. Multi-component Sérsic galaxies resemble compact psuedobulges ($n\sim$\,2, $R_e \sim$\, 4 kpc) surrounded by extended Gaussian-like outer structures ($R_e > 10$ kpc). 11\% of galaxies ($N=52$) feature a break in their outer profiles, indicating `truncated' or `anti-truncated' discs. Beyond the break radius, truncated galaxies are structurally consistent with exponential discs, disfavouring physical truncation as their formation mechanism. Bulge luminosity in anti-truncated galaxies correlates strongly with galaxy luminosity, indicating a bulge-enhancing origin for these systems. Both types of broken disc are found overwhelmingly ($>70\%$) in `barred' galaxies, despite a low measured bar fraction for Coma ($20\pm2\%$). Thus, galaxy bars play an important role in formation of broken disc structures. No strong variation in galaxy structure is detected with projected cluster-centric radius.

preprint2015arXiv

CFHTLenS: Testing the Laws of Gravity with Tomographic Weak Lensing and Redshift Space Distortions

Dark energy may be the first sign of new fundamental physics in the Universe, taking either a physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and galaxy peculiar velocities provide complementary probes of General Relativity, and in combination allow us to test modified theories of gravity in a unique way. We perform such an analysis by combining measurements of cosmic shear tomography from the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the strongest existing joint constraints on the metric potentials that describe general theories of gravity. For scale-independent modifications to the metric potentials which evolve linearly with the effective dark energy density, we find present-day cosmological deviations in the Newtonian potential and curvature potential from the prediction of General Relativity to be (Delta Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = -0.05 \pm 0.3 respectively (68 per cent CL).

preprint2015arXiv

CFHTLenS: Weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.

preprint2015arXiv

Cosmological parameters from the comparison of peculiar velocities with predictions from the 2M++ density field

Peculiar velocity measurements are the only tool available in the low-redshift Universe for mapping the large-scale distribution of matter and can thus be used to constrain cosmology. Using redshifts from the 2M++ redshift compilation, we reconstruct the density of galaxies within 200 Mpc/h, allowing for the first time good sampling of important superclusters such as the Shapley Concentration. We compare the predicted peculiar velocities from 2M++ to Tully-Fisher and SNe peculiar velocities. We find a value of $β^* \equiv Ω_{\rm{m}}^{0.55}/b^* = 0.431 \pm 0.021$, suggesting $Ω_{\rm{m}}^{0.55}σ_{\rm{8,lin}} = 0.401 \pm 0.024$, in good agreement with other probes. The predicted peculiar velocity of the Local Group arising from the 2M++ volume alone is $540 \pm 40$ km/s, towards $l = 268 \pm 4, b= 38 \pm 6$, only $10^\circ$ out of alignment with the Cosmic Microwave Background dipole. To account for velocity contributions arising from sources outside the 2M++ volume, we fit simultaneously for $β^*$ and an external bulk flow in our analysis. We find that an external bulk flow is preferred at the 5.1$σ$ level, and the best fit has a velocity of $159\pm23$ km/s towards $l=304 \pm 11, b = 6 \pm 13$. Finally, the predicted bulk flow of a 50 Mpc/h Gaussian-weighted volume centred on the Local Group is $230 \pm 30$ km/s, in the direction $l=293\pm 8, b = 14 \pm 10$, in agreement with predictions from $Λ$CDM.

preprint2015arXiv

Kinematic Bias in Cosmological Distance Measurement

Recent calculations using non-linear relativistic cosmological perturbation theory show biases in the mean luminosity distance and distance modulus at low redshift. We show that these effects may be understood very simply as a non-relativistic, and purely kinematic, Malmquist-like bias, and we describe how the effect changes if one averages over sources that are limited by apparent magnitude. This effect is essentially identical to the distance bias from small-scale random velocities that has previously been considered by astronomers, though we find that the standard formula overestimates the homogeneous bias by a factor 2.

preprint2015arXiv

On the perturbation of the luminosity distance by peculiar motions

We consider some aspects of the perturbation to the luminosity distance $d(z)$ that are of relevance for SN1a cosmology and for future peculiar velocity surveys at non-negligible redshifts. 1) Previous work has shown that the correction to the lowest order perturbation $δd / d = -δv / c z$ has the peculiar characteristic that it appears to depend on the absolute state of motion of sources, rather than on their motion relative to that of the observer. The resolution of this apparent violation of the equivalence principle is that it is necessary to allow for evolution of the velocities with time, and also, when considering perturbations on the scale of the observer-source separation, to include the gravitational redshift effect. We provide an expression for $δd / d$ that provides a physically consistent way to measure peculiar velocities and determine their impact for SN1a cosmology. 2) We then calculate the perturbation to the redshift as a function of source flux density, which has been proposed as an alternative probe of large-scale motions. We show how the inclusion of surface brightness modulation modifies the relation between $δz(m)$ and the peculiar velocity, and that, while the noise properties of this method might appear promising, the velocity signal is swamped by the effect of galaxy clustering for most scales of interest. 3) We show how, in linear theory, peculiar velocity measurements are biased downwards by the effect of smaller scale motions or by measurement errors (such as in photometric redshifts). Our results nicely explain the effects seen in simulations by Koda et al.\ 2013. We critically examine the prospects for extending peculiar velocity studies to larger scales with near-term future surveys.

preprint2014arXiv

CFHTLenS: A Weak Lensing Shear Analysis of the 3D-Matched-Filter Galaxy Clusters

We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$σ$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$σ$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^β$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $β\sim 1.4 \pm 0.1$, both of which are in 1$σ$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org.

preprint2014arXiv

CFHTLenS: Co-evolution of galaxies and their dark matter haloes

Galaxy-galaxy weak lensing is a direct probe of the mean matter distribution around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield statistically significant galaxy halo mass measurements over a much wider range of stellar masses ($10^{8.75}$ to $10^{11.3} M_{\odot}$) and redshifts ($0.2 < z < 0.8$) than previous weak lensing studies. At redshift $z \sim 0.5$, the stellar-to-halo mass ratio (SHMR) reaches a maximum of $4.0\pm0.2$ percent as a function of halo mass at $\sim 10^{12.25} M_{\odot}$. We find, for the first time from weak lensing alone, evidence for significant evolution in the SHMR: the peak ratio falls as a function of cosmic time from $4.5 \pm 0.3$ percent at $z \sim 0.7$ to $3.4 \pm 0.2$ percent at $z \sim 0.3$, and shifts to lower stellar mass haloes. These evolutionary trends are dominated by red galaxies, and are consistent with a model in which the stellar mass above which star formation is quenched "downsizes" with cosmic time. In contrast, the SHMR of blue, star-forming galaxies is well-fit by a power law that does not evolve with time. This suggests that blue galaxies form stars at a rate that is balanced with their dark matter accretion in such a way that they evolve along the SHMR locus. The redshift dependence of the SHMR can be used to constrain the evolution of the galaxy population over cosmic time.

preprint2014arXiv

CFHTLenS: Cosmological constraints from a combination of cosmic shear two-point and three-point correlations

Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2σ$. The combined constraint on $Σ_8 = σ_8 (Ω_{\rm m}/0.27)^α$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $Ω_{\rm m}$ and $σ_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $Ω_{\rm m}$ and $σ_8$ by 26% for flat $Λ$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $Ω_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function.

preprint2014arXiv

Dark Matter Halos in Galaxies and Globular Cluster Populations

We combine a new, comprehensive database for globular cluster populations in all types of galaxies with a new calibration of galaxy halo masses based entirely on weak lensing. Correlating these two sets of data, we find that the mass ratio $η\equiv M_{GCS}/M_{h}$ (total mass in globular clusters, divided by halo mass) is essentially constant at $\langle η\rangle \sim 4 \times 10^{-5}$, strongly confirming earlier suggestions in the literature. Globular clusters are the only known stellar population that formed in essentially direct proportion to host galaxy halo mass. The intrinsic scatter in $η$ appears to be at most 0.2 dex; we argue that some of this scatter is due to differing degrees of tidal stripping of the globular cluster systems between central and satellite galaxies. We suggest that this correlation can be understood if most globular clusters form at very early stages in galaxy evolution, largely avoiding the feedback processes that inhibited the bulk of field-star formation in their host galaxies. The actual mean value of $η$ also suggests that about $1/4$ of the \emph{initial} gas mass present in protogalaxies collected into GMCs large enough to form massive, dense star clusters. Finally, our calibration of $\langle η\rangle$ indicates that the halo masses of the Milky Way and M31 are $(1.2\pm0.5)\times 10^{12} M_{\odot}$ and $(3.9\pm1.8)\times 10^{12} M_{\odot}$ respectively.

preprint2014arXiv

Dissecting the Red Sequence: The Bulge and Disc Colours of Early-Type Galaxies in the Coma Cluster

We explore the internal structure of red sequence galaxies in the Coma cluster across a wide range of luminosities ($-17>M_g>-22$) and cluster-centric radii ($0<r_{\rm{cluster}}<1.3 r_{200}$). We present the 2D bulge-disc decomposition of galaxies in deep Canada-France-Hawaii Telescope $u,g,i$ imaging using GALFIT. Rigorous filtering is applied to identify an analysis sample of 200 galaxies which are well described by an `archetypal' S0 structure (central bulge + outer disc). We consider internal bulge and/or disc colour gradients by allowing component sizes to vary between bands. Gradients are required for $30\%$ of analysis sample galaxies. Bulge half-light radii are found to be uncorrelated with galaxy luminosity ($R_e \sim 1$ kpc, $n\sim2$) for all but the brightest galaxies ($M_g<-20.5$). The S0 discs are brighter (at fixed size, or smaller at fixed luminosity) than those of star-forming spirals. A similar colour-magnitude relation is found for both bulges and discs. The global red sequence for S0s in Coma hence results from a combination of both component trends. We measure an average bulge $-$ disc colour difference of $0.09\pm0.01$ mag in $g-i$, and $0.16\pm0.01$ mag in $u-g$. Using simple stellar population models, bulges are either $\sim2$-$3\times$ older, or $\sim2\times$ more metal-rich than discs. The trend towards bluer global S0 colours observed further from Coma's core is driven by a significant correlation in disc colour with cluster-centric radius. An equivalent trend is detected in bulge colours at a marginal significance level. Our results therefore favour environment-mediated mechanisms of disc fading as the dominant factor in S0 formation.

preprint2014arXiv

The Phase Space and Stellar Populations of Cluster Galaxies at z ~ 1: Simultaneous Constraints on the Location and Timescale of Satellite Quenching

We investigate the velocity vs. position phase space of z ~ 1 cluster galaxies using a set of 424 spectroscopic redshifts in 9 clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories: quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commonly found at small clustercentric radii with high clustercentric velocities, and appear to trace a coherent ``ring" in phase space. Using several zoom simulations of clusters we show that the coherent distribution of the poststarbursts can be reasonably well-reproduced using a simple quenching scenario. Specifically, the phase space is best reproduced if satellite quenching occurs on a rapid timescale (0.1 < tau_{Q} < 0.5 Gyr) after galaxies make their first passage of R ~ 0.5R_{200}, a process that takes a total time of ~ 1 Gyr after first infall. We compare this quenching timescale to the timescale implied by the stellar populations of the poststarburst galaxies and find that the poststarburst spectra are well-fit by a rapid quenching (tau_{Q} = 0.4^{+0.3}_{-0.4} Gyr) of a typical star-forming galaxy. The similarity between the quenching timescales derived from these independent indicators is a strong consistency check of the quenching model. Given that the model implies satellite quenching is rapid, and occurs well within R_{200}, this would suggest that ram-pressure stripping of either the hot or cold gas component of galaxies are the most plausible candidates for the physical mechanism. The high cold gas consumption rates at z ~ 1 make it difficult to determine if hot or cold gas stripping is dominant; however, measurements of the redshift evolution of the satellite quenching timescale and location may be capable of distinguishing between the two.

preprint2013arXiv

CFHTLenS tomographic weak lensing cosmological parameter constraints: Mitigating the impact of intrinsic galaxy alignments

We present a finely-binned tomographic weak lensing analysis of the Canada-France-Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, with accurate shear and photometric redshifts for a galaxy sample with a median redshift of zm =0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six redshift bins, each spanning the angular range of 1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological probes: the cosmic microwave background with data from WMAP7, baryon acoustic oscillations with data from BOSS, and a prior on the Hubble constant from the HST distance ladder. This leads to constraints on the normalisation of the matter power spectrum sigma_8 = 0.799 +/- 0.015 and the matter density parameter Omega_m = 0.271 +/- 0.010 for a flat Lambda CDM cosmology. For a flat wCDM cosmology we constrain the dark energy equation of state parameter w = -1.02 +/- 0.09. We also provide constraints for curved Lambda CDM and wCDM cosmologies. We find the intrinsic alignment contamination to be galaxy-type dependent with a significant intrinsic alignment signal found for early-type galaxies, in contrast to the late-type galaxy sample for which the intrinsic alignment signal is found to be consistent with zero.

preprint2013arXiv

CFHTLenS: Combined probe cosmological model comparison using 2D weak gravitational lensing

We present cosmological constraints from 2D weak gravitational lensing by the large-scale structure in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate photometric redshifts and measured shapes for 4.2 million galaxies between redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function over angular scales ranging between 0.8 and 350 arcmin. Using non-linear models of the dark-matter power spectrum, we constrain cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. The best constraints from lensing alone are obtained for the small-scale density-fluctuations amplitude sigma_8 scaled with the total matter density Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 = 0.79+-0.03. We combine the CFHTLenS data with WMAP7, BOSS and an HST distance-ladder prior on the Hubble constant to get joint constraints. For a flat LambdaCDM model, we find Omega_m = 0.283+-0.010 and sigma_8 = 0.813+-0.014. In the case of a curved wCDM universe, we obtain Omega_m = 0.27+-0.03, sigma_8 = 0.83+-0.04, w_0 = -1.10+-0.15 and Omega_K = 0.006+0.006-0.004. We calculate the Bayesian evidence to compare flat and curved LambdaCDM and dark-energy CDM models. From the combination of all four probes, we find models with curvature to be at moderately disfavoured with respect to the flat case. A simple dark-energy model is indistinguishable from LambdaCDM. Our results therefore do not necessitate any deviations from the standard cosmological model.

preprint2013arXiv

CFHTLenS: Mapping the Large Scale Structure with Gravitational Lensing

We present a quantitative analysis of the largest contiguous maps of projected mass density obtained from gravitational lensing shear. We use data from the 154 deg2 covered by the Canada-France-Hawaii Telescope Lensing Survey. Our study is the first attempt to quantitatively characterize the scientific value of lensing maps, which could serve in the future as a complementary approach to the study of the dark universe with gravitational lensing. We show that mass maps contain unique cosmological information beyond that of traditional two-points statistical analysis techniques. Using a series of numerical simulations, we first show how, reproducing the CFHTLenS observing conditions, gravitational lensing inversion provides a reliable estimate of the projected matter distribution of large scale structure. We validate our analysis by quantifying the robustness of the maps with various statistical estimators. We then apply the same process to the CFHTLenS data. We find that the 2-points correlation function of the projected mass is consistent with the cosmological analysis performed on the shear correlation function discussed in the CFHTLenS companion papers. The maps also lead to a significant measurement of the third order moment of the projected mass, which is in agreement with analytic predictions, and to a marginal detection of the fourth order moment. Tests for residual systematics are found to be consistent with zero for the statistical estimators we used. A new approach for the comparison of the reconstructed mass map to that predicted from the galaxy distribution reveals the existence of giant voids in the dark matter maps as large as 3 degrees on the sky. Our analysis shows that lensing mass maps can be used for new techniques such as peak statistics and the morphological analysis of the projected dark matter distribution.

preprint2013arXiv

CFHTLenS: The Environmental Dependence of Galaxy Halo Masses from Weak Lensing

We use weak gravitational lensing to analyse the dark matter halos around satellite galaxies in galaxy groups in the CFHTLenS dataset. This dataset is derived from the CFHTLS-Wide survey, and encompasses 154 sq. deg of high-quality shape data. Using the photometric redshifts, we divide the sample of lens galaxies with stellar masses in the range 10^9 Msun to 10^10.5 Msun into those likely to lie in high-density environments (HDE) and those likely to lie in low-density environments (LDE). Through comparison with galaxy catalogues extracted from the Millennium Simulation, we show that the sample of HDE galaxies should primarily (~61%) consist of satellite galaxies in groups, while the sample of LDE galaxies should consist of mostly (~87%) non-satellite (field and central) galaxies. Comparing the lensing signals around samples of HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE galaxies clearly shows a positive contribution from their host groups on their lensing signals at radii of ~500--1000 kpc, the typical separation between satellites and group centres. More importantly, the subhalos of HDE galaxies are less massive than those around LDE galaxies by a factor 0.65 +/- 0.12, significant at the 2.9 sigma level. A natural explanation is that the halos of satellite galaxies are stripped through tidal effects in the group environment. Our results are consistent with a typical tidal truncation radius of ~40 kpc.

preprint2013arXiv

CFHTLenS: The relation between galaxy dark matter haloes and baryons from weak gravitational lensing

We present a study of the relation between dark matter halo mass and the baryonic content of host galaxies, quantified via luminosity and stellar mass. Our investigation uses 154 deg2 of Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) lensing and photometric data, obtained from the CFHT Legacy Survey. We employ a galaxy-galaxy lensing halo model which allows us to constrain the halo mass and the satellite fraction. Our analysis is limited to lenses at redshifts between 0.2 and 0.4. We express the relationship between halo mass and baryonic observable as a power law. For the luminosity-halo mass relation we find a slope of 1.32+/-0.06 and a normalisation of 1.19+0.06-0.07x10^13 h70^-1 Msun for red galaxies, while for blue galaxies the best-fit slope is 1.09+0.20-0.13 and the normalisation is 0.18+0.04-0.05x10^13 h70^-1 Msun. Similarly, we find a best-fit slope of 1.36+0.06-0.07 and a normalisation of 1.43+0.11-0.08x10^13 h70^-1 Msun for the stellar mass-halo mass relation of red galaxies, while for blue galaxies the corresponding values are 0.98+0.08-0.07 and 0.84+0.20-0.16x10^13 h70^-1 Msun. For red lenses, the fraction which are satellites tends to decrease with luminosity and stellar mass, with the sample being nearly all satellites for a stellar mass of 2x10^9 h70^-2 Msun. The satellite fractions are generally close to zero for blue lenses, irrespective of luminosity or stellar mass. This, together with the shallower relation between halo mass and baryonic tracer, is a direct confirmation from galaxy-galaxy lensing that blue galaxies reside in less clustered environments than red galaxies. We also find that the halo model, while matching the lensing signal around red lenses well, is prone to over-predicting the large-scale signal for faint and less massive blue lenses. This could be a further indication that these galaxies tend to be more isolated than assumed. [abridged]

preprint2013arXiv

Disentangling satellite galaxy populations using orbit tracking in simulations

Physical processes regulating star formation in satellite galaxies represent an area of ongoing research, but the projected nature of observed coordinates makes separating different populations of satellites (with different processes at work) difficult. The orbital history of a satellite galaxy leads to its present-day phase space coordinates; we can also work backwards and use these coordinates to statistically infer information about the orbital history. We use merger trees from the MultiDark Run 1 N-body simulation to compile a catalog of the orbits of satellite haloes in cluster environments. We parameterize the orbital history by the time since crossing within 2.5 rvir of the cluster centre and use our catalog to estimate the probability density over a range of this parameter given a set of present-day projected (i.e. observable) phase space coordinates. We show that different populations of satellite haloes, e.g. infalling, backsplash and virialized, occupy distinct regions of phase space, and semi-distinct regions of projected phase space. This will allow us to probabilistically determine the time since infall of a large sample of observed satellite galaxies, and ultimately to study the effect of orbital history on star formation history (the topic of a future paper). We test the accuracy of our method and find that we can reliably recover this time within +/-2.58 Gyr in 68 per cent of cases by using all available phase space coordinate information, compared to +/-2.64 Gyr using only position coordinates and +/-3.10 Gyr guessing 'blindly', i.e. using no coordinate information, but with knowledge of the overall distribution of infall times. In some regions of phase space, the accuracy of the infall time estimate improves to +/-1.85 Gyr. Although we focus on time since infall, our method is easily generalizable to other orbital parameters (e.g. pericentric distance and time).

preprint2013arXiv

First measurement of the bulk flow of nearby galaxies using the cosmic microwave background

Peculiar velocities in the nearby Universe can be measured via the kinetic Sunyaev-Zel'dovich (kSZ) effect. Using a statistical method based on an optimised cross-correlation with nearby galaxies, we extract the kSZ signal generated by plasma halo of galaxies from the Cosmic Microwave Background (CMB) temperature anisotropies observed by the Wilkinson Microwave Anisotropy Probe (WMAP). Marginalising over the thermal Sunyaev-Zel'dovich contribution from clusters of galaxies, possible unresolved point source contamination, and Galactic foregrounds (dust, synchrotron and free-free emission), we report a kSZ bulk flow signal present at the 90% confidence level in the seven-year WMAP data. When only galaxies within 50 Mpc/h are included in the kSZ template we find a bulk flow in the CMB frame of |V|=533 +/- 263 km/s, in the direction l=324 +/- 27, b=-7 +/- 17, consistent with bulk flow measurements on a similar scale using classical distance indicators. We show how this comparison constrains the (ionised) baryonic budget in the local universe. On very large (~ 500 Mpc/h) scales, we find a 95% upper limit of 470 km/s, inconsistent with some analyses of bulk flow of clusters from the kSZ. We estimate that the significance of the bulk flow signal may increase to 3-5 sigma using data from the PLANCK probe.

preprint2013arXiv

Quenching star formation in cluster galaxies

In order to understand the processes that quench star formation within rich clusters, we construct a library of subhalo orbits drawn from $Λ$CDM cosmological N-body simulations of four rich clusters. The orbits are combined with models of star formation followed by quenching in the cluster environment. These are compared with observed bulge and disc colours and stellar absorption linestrength indices of satellite galaxies. Models in which the bulge stellar populations depend only on the galaxy subhalo mass while the disc quenching depends on the cluster environment are acceptable fits to the data. An exponential disc quenching timescale of 3 - 3.5 Gyr is preferred. Models with short ($\lesssim 1$ Gyr) quenching timescales yield cluster-centric gradients in disc colours and Balmer line indices that are too steep compared to observations. We also examine models in which there is quenching in lower mass groups prior to cluster infall ("pre-processing"), finding that such models are a better fit to the data than models without pre-processing and require similar quenching times. The data slightly prefer models where quenching occurs only for galaxies falling within about 0.5 $r_{200}$. Finally, we have examined models with short quenching timescales of 1 Gyr, but a long delay time of 3 Gyr prior to quenching. All models with short quenching timescales, even such "delayed-then-rapid" quenching models, produce excessively red galaxies near the cluster core and are strongly disfavoured by the data. These results imply that the environments of rich clusters must impact star formation rates of infalling galaxies on relatively long timescales -- several times longer than a typical halo spends within the virial radius of a cluster. This scenario favours gentler quenching mechanisms such as slow "strangulation" over more rapid ram-pressure stripping.

preprint2012arXiv

CFHTLenS tomographic weak lensing: Quantifying accurate redshift distributions

The Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) comprises deep multi-colour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate photometric redshifts and shape measurements. We demonstrate that the redshift probability distribution function summed over galaxies provides an accurate representation of the galaxy redshift distribution accounting for random and catastrophic errors for galaxies with best fitting photometric redshifts z_p < 1.3. We present cosmological constraints using tomographic weak gravitational lensing by large-scale structure. We use two broad redshift bins 0.5 < z_p <= 0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and measure the shear correlation function on angular scales in the range ~1-40 arcmin. We show that the problematic redshift scaling of the shear signal, found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For a flat Lambda-CDM model and a fixed matter density Omega_m=0.27, we find the normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When combined with cosmic microwave background data (WMAP7), baryon acoustic oscillation data (BOSS), and a prior on the Hubble constant from the HST distance ladder, we find that CFHTLenS improves the precision of the fully marginalised parameter estimates by an average factor of 1.5-2. Combining our results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074 and sigma_8=0.802 \pm 0.013.

preprint2012arXiv

CFHTLenS: The Canada-France-Hawaii Telescope Lensing Survey

We present the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) that accurately determines a weak gravitational lensing signal from the full 154 square degrees of deep multi-colour data obtained by the CFHT Legacy Survey. Weak gravitational lensing by large-scale structure is widely recognised as one of the most powerful but technically challenging probes of cosmology. We outline the CFHTLenS analysis pipeline, describing how and why every step of the chain from the raw pixel data to the lensing shear and photometric redshift measurement has been revised and improved compared to previous analyses of a subset of the same data. We present a novel method to identify data which contributes a non-negligible contamination to our sample and quantify the required level of calibration for the survey. Through a series of cosmology-insensitive tests we demonstrate the robustness of the resulting cosmic shear signal, presenting a science-ready shear and photometric redshift catalogue for future exploitation.

preprint2012arXiv

Power Spectrum Estimation from Peculiar Velocity Catalogues

The peculiar velocities of galaxies are an inherently valuable cosmological probe, providing an unbiased estimate of the distribution of matter on scales much larger than the depth of the survey. Much research interest has been motivated by the high dipole moment of our local peculiar velocity field, which suggests a large scale excess in the matter power spectrum, and can appear to be in some tension with the LCDM model. We use a composite catalogue of 4,537 peculiar velocity measurements with a characteristic depth of 33 h-1 Mpc to estimate the matter power spectrum. We compare the constraints with this method, directly studying the full peculiar velocity catalogue, to results from Macaulay et al. (2011), studying minimum variance moments of the velocity field, as calculated by Watkins, Feldman & Hudson (2009) and Feldman, Watkins & Hudson (2010). We find good agreement with the LCDM model on scales of k > 0.01 h Mpc-1. We find an excess of power on scales of k < 0.01 h Mpc-1, although with a 1 sigma uncertainty which includes the LCDM model. We find that the uncertainty in the excess at these scales is larger than an alternative result studying only moments of the velocity field, which is due to the minimum variance weights used to calculate the moments. At small scales, we are able to clearly discriminate between linear and nonlinear clustering in simulated peculiar velocity catalogues, and find some evidence (although less clear) for linear clustering in the real peculiar velocity data.

preprint2012arXiv

Probing Satellite Halos with Weak Gravitational Lensing

We demonstrate the possibility of detecting tidal stripping of dark matter subhalos within galaxy groups using weak gravitational lensing. We have run ray-tracing simulations on galaxy catalogues from the Millennium Simulation to generate mock shape catalogues. The ray-tracing catalogues assume a halo model for galaxies and groups, using various models with different distributions of mass between galaxy and group halos to simulate different stages of group evolution. Using these mock catalogues, we forecast the lensing signals that will be detected around galaxy groups and satellite galaxies, as well as test two different methods for isolating the satellites' lensing signals. A key challenge is to determine the accuracy to which group centres can be identified. We show that with current and ongoing surveys, it will possible to detect stripping in groups of mass 10^12--10^15 Msun.

preprint2012arXiv

The growth rate of cosmic structure from peculiar velocities at low and high redshifts

Peculiar velocities are an important probe of the growth rate of mass density fluctuations in the Universe. Most previous studies have focussed exclusively on measuring peculiar velocities at intermediate ($0.2 < z < 1$) redshifts using statistical redshift-space distortions. Here we emphasize the power of peculiar velocities obtained directly from distance measurements at low redshift ($z \lesssim 0.05$), and show that these data break the usual degeneracies in the Ω_{m,0} -- $σ_{8,0}$ parameter space. Using only peculiar velocity data, we find $Ω_{m,0} = 0.259\pm0.045$ and $σ_{8,0} = 0.748\pm0.035$. Fixing the amplitude of fluctuations at very high redshift using observations of the Cosmic Microwave Background (CMB), the same data can be used to constrain the growth index $γ$, with the strongest constraints coming from peculiar velocity measurements in the nearby Universe. We find $γ= 0.619\pm 0.054$, consistent with LCDM. Current peculiar velocity data already strongly constrain modified gravity models, and will be a powerful test as data accumulate.

preprint2011arXiv

A Slight Excess of Large Scale Power from Moments of the Peculiar Velocity Field

The peculiar motions of galaxies can be used to infer the distribution of matter in the Universe. It has recently been shown that measurements of the peculiar velocity field indicates an anomalously high bulk flow of galaxies in our local volume. In this paper we find the implications of the high bulk flow for the power spectrum of density fluctuations. We find that analyzing only the dipole moment of the velocity field yields an average power spectrum amplitude which is indeed much higher than the LCDM value. However, by also including shear and octupole moments of the velocity field, and marginalizing over possible values for the growth rate, an average power spectrum amplitude which is consistent with LCDM is recovered. We attempt to infer the shape of the matter power spectrum from moments of the velocity field, and find a slight excess of power on scales ~ h-1 Gpc.

preprint2011arXiv

Cosmic flows in the nearby universe from Type Ia Supernovae

Peculiar velocities are one of the only probes of very large-scale mass density fluctuations in the nearby Universe. We present new "minimal variance" bulk flow measurements based upon the "First Amendment" compilation of 245 Type Ia supernovae (SNe) peculiar velocities and find a bulk flow of 249 +/- 76 km/s in the direction l= 319 +/- 18 deg, b = 7 +/- 14 deg. The SNe bulk flow is consistent with the expectations of ΛCDM. However, it is also marginally consistent with the bulk flow of a larger compilation of non-SNe peculiar velocities (Watkins, Feldman, & Hudson 2009). By comparing the SNe peculiar velocities to predictions of the IRAS Point Source Catalog Redshift survey (PSCz) galaxy density field, we find Ω_{m}^{0.55} σ_{8,lin} = 0.40 +/- 0.07, which is in agreement with ΛCDM. However, we also show that the PSCz density field fails to account for 150 +/- 43 km/s of the SNe bulk motion.

preprint2011arXiv

Environmental quenching and hierarchical cluster assembly: Evidence from spectroscopic ages of red-sequence galaxies in Coma

We explore the variation in stellar population ages for Coma cluster galaxies as a function of projected cluster-centric distance, using a sample of 362 red-sequence galaxies with high signal-to-noise spectroscopy. The sample spans a wide range in luminosity (0.02-4 L*) and extends from the cluster core to near the virial radius. We find a clear distinction in the observed trends of the giant and dwarf galaxies. The ages of red-sequence giants are primarily determined by galaxy mass, with only weak modulation by environment, in the sense that galaxies at larger cluster-centric distance are slightly younger. For red-sequence dwarfs (with mass <10^10 Msun), the roles of mass and environment as predictors of age are reversed: there is little dependence on mass, but strong trends with projected cluster-centric radius are observed. The average age of dwarfs at the 2.5 Mpc limit of our sample is approximately half that of dwarfs near the cluster centre. The gradient in dwarf galaxy ages is a global cluster-centric trend, and is not driven by the ongoing merger of the NGC 4839 group to the south west of Coma. We interpret these results using environmental histories extracted from the Millennium Simulation for members of massive clusters. Hierarchical cluster assembly naturally leads to trends in the accretion times of galaxies as a function of projected cluster-centric radius. On average, simulated galaxies now located in cluster cores joined halos above any given mass threshold earlier than those now located in the outskirts of clusters. We test environmental quenching models, in which star formation is halted in galaxies when they enter halos of a given mass, or become satellite galaxies. The models broadly reproduce the gradients observed in Coma, but for dwarf galaxies the efficiency of environmental quenching must be very high to match the strong trends observed.

preprint2011arXiv

The 2M++ galaxy redshift catalogue

Peculiar velocities arise from gravitational instability, and thus are linked to the surrounding distribution of matter. In order to understand the motion of the Local Group with respect to the Cosmic Microwave Background, a deep all-sky map of the galaxy distribution is required. Here we present a new redshift compilation of 69~160 galaxies, dubbed 2M++, to map large-scale structures of the Local Universe over nearly the whole sky, and reaching depths of K <= 12.5, or 200 Mpc/h. The target catalogue is based on the Two-Micron-All-Sky Extended Source Catalog (2MASS-XSC). The primary sources of redshifts are the 2MASS Redshift Survey, the 6dF galaxy redshift survey and the Sloan Digital Sky Survey (DR7). We assess redshift completeness in each region and compute the weights required to correct for redshift incompleteness and apparent magnitude limits, and discuss corrections for incompleteness in the Zone of Avoidance. We present the density field for this survey, and discuss the importance of large-scale structures such as the Shapley Concentration.

preprint2011arXiv

The HST/ACS Coma Cluster Survey IV. Intergalactic Globular Clusters and the Massive Globular Cluster System at the Core of the Coma Galaxy Cluster

Intracluster stellar populations are a natural result of tidal interactions in galaxy clusters. Measuring these populations is difficult, but important for understanding the assembly of the most massive galaxies. The Coma cluster is one of the nearest truly massive galaxy clusters, and is host to a correspondingly large system of globular clusters (GCs). We use imaging from the HST/ACS Coma Cluster Survey to present the first definitive detection of a large population of intracluster GCs (IGCs) that fills the Coma cluster core and is not associated with individual galaxies. The GC surface density profile around the central massive elliptical galaxy, NGC 4874, is dominated at large radii by a population of IGCs that extend to the limit of our data (R<520 kpc). We estimate that there are 47000+/-1600 (random) +4000/-5000 (systematic) IGCs out to this radius, and that they make up ~70% of the central GC system, making this the largest GC system in the nearby Universe. Even including the GC systems of other cluster galaxies, IGCs still make up ~30-45% of the GCs in the cluster core. Observational limits from previous studies of the intracluster light (ICL) suggest that the IGC population has a high specific frequency. If the IGC population has a specific frequency similar to high-S_N dwarf galaxies, then the ICL has a total stellar mass of ~10^12 M_sun within the cluster core. The ICL makes up approximately half of the stellar luminosity and one-third of the stellar mass of the central (NGC4874+ICL) system. The color distribution of the IGC population is bimodal, with blue, metal-poor GCs outnumbering red, metal-rich GCs by a ratio of 4:1. The fraction of red IGCs (20%), and the red color of those GCs, implies that IGCs can originate from the halos of relatively massive, L* galaxies, and not solely from the disruption of dwarf galaxies. (Abridged)

preprint2010arXiv

Group-finding with photometric redshifts: The Photo-z Probability Peaks algorithm

We present a galaxy group-finding algorithm, the Photo-z Probability Peaks (P3) algorithm, optimized for locating small galaxy groups using photometric redshift data by searching for peaks in the signal-to-noise of the local overdensity of galaxies in a three-dimensional grid. This method is an improvement over similar two-dimensional matched-filter methods in reducing background contamination through the use of redshift information, allowing it to accurately detect groups at lower richness. We present the results of tests of our algorithm on galaxy catalogues from the Millennium Simulation. Using a minimum S/N of 3 for detected groups, a group aperture size of 0.25 Mpc/h, and assuming photometric redshift accuracy of sigma_z = 0.05 it attains a purity of 84% and detects ~295 groups/deg.^2 with an average group richness of 8.6 members. Assuming photometric redshift accuracy of sigma_z = 0.02, it attains a purity of 97% and detects ~143 groups/deg.^2 with an average group richness of 12.5 members. We also test our algorithm on data available for the COSMOS field and the presently-available fields from the CFHTLS-Wide survey, presenting preliminary results of this analysis.

preprint2010arXiv

The HST/ACS Coma Cluster Survey. II. Data Description and Source Catalogs

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

preprint2010arXiv

Ultraviolet tails and trails in cluster galaxies: A sample of candidate gaseous stripping events in Coma

We have used new deep observations of the Coma cluster from GALEX to identify 13 star-forming galaxies with asymmetric morphologies in the ultraviolet. Aided by optical broad-band and H-alpha imaging, we interpret the asymmetric features as being due to star formation within gas stripped from the galaxies by interaction with the cluster environment. The selected objects display a range of structures from broad fan-shaped systems of filaments and knots (`jellyfish') to narrower and smoother tails extending up to 100 kpc in length. Some of the features have been discussed previously in the literature, while others are newly identified here. As an ensemble, the candidate stripping events are located closer to the cluster centre than other star-forming galaxies; their radial distribution is similar to that of all cluster members, dominated by passive galaxies. The fraction of blue galaxies which are undergoing stripping falls from 40% in the central 500 kpc, to less than 5% beyond 1 Mpc. We find that tails pointing away from (i.e. galaxies moving towards) the cluster centre are strongly favoured (11/13 cases). From the small number of `outgoing' galaxies with stripping signatures we conclude that the stripping events occur primarily on first passage towards the cluster centre, and are short-lived compared to the cluster crossing time. Using infall trajectories from simulations, the observed fraction of blue galaxies undergoing stripping can be reproduced if the events are triggered at a threshold radius of ~1 Mpc and detectable for ~500 Myr. HST images are available for two galaxies from our sample and reveal compact blue knots coincident with UV and H-alpha emission, apparently forming stars within the stripped material. Our results confirm that stripping of gas from infalling galaxies, and associated star formation in the stripped material, is a widespread phenomenon in rich clusters.

preprint2009arXiv

Abundance ratios in red-sequence galaxies over a wide mass range: The "X-planes" for magnesium, calcium, carbon and nitrogen

We analyse the abundance ratios of the light elements Mg, Ca, C and N, relative to Fe, for 147 red-sequence galaxies in the Coma cluster and the Shapley Supercluster. The sample covers a six-magnitude range in luminosity, from giant ellipticals to dwarfs at M^*+4. We exploit the wide mass range to investigate systematic trends in the abundance ratios Mg/Fe, Ca/Fe, C/Fe and N/Fe. We find that each of these ratios can be well modelled using two-parameter relations of the form [X/Fe] = a0 + a1 log sigma + a2 [Fe/H], where sigma is the velocity dispersion. Analysing these "X-planes" reveals new structure in the abundance patterns, beyond the traditional one-parameter (e.g. Mg/Fe-sigma) correlations. The X-planes for the alpha elements, Mg and Ca, indicate a positive correlation with velocity dispersion, and simultaneously an anti-correlation with Fe/H (i.e. a1>0 and a2<0). Taking both effects into account dramatically reduces the scatter, compared to the traditional X/Fe-sigma relations. For C and N, a similar correlation with velocity dispersion is recovered, but there is no additional dependence on Fe/H (i.e. a1>0 and a2~0). The explicit dependence of X/Fe on two parameters is evidence that at least two physical processes are at work in setting the abundance patterns. The Fe/H dependence of Mg/Fe and Ca/Fe, at fixed sigma, may result from different durations of star formation, from galaxy to galaxy. The absence of corresponding Fe/H dependence for C and N is consistent with these elements being generated in lower-mass stars. The increase with sigma, at fixed Fe/H, is similar for elements Mg, C and N, and slightly shallower for Ca. This pattern of trends cannot be explained solely by a systematic variation of star-formation time-scale with sigma.

preprint2009arXiv

Ages and metallicities for quiescent galaxies in the Shapley Supercluster: Driving parameters of the stellar populations

We use high signal-to-noise spectroscopy for a sample of 232 quiescent galaxies in the Shapley Supercluster, to investigate how their stellar populations depend on velocity dispersion, luminosity and stellar mass. The sample spans a large range in velocity dispersion (sigma from 30-300 km/s) and in luminosity (M_R from -18.7 to -23.2). Estimates of age, total metallicity (Z/H) and alpha-element abundance ratio (a/Fe) were derived from absorption-line analysis, using single-burst models. Age, Z/H and a/Fe are all correlated positively with velocity dispersion, but we also find significant residual trends with luminosity: at given sigma, the brighter galaxies are younger, less alpha-enriched, and have higher Z/H. At face value, these results might suggest that the stellar populations depend on stellar mass as well as on velocity dispersion. However, we show that the observed trends can be reproduced by models in which the stellar populations depend systematically only on sigma, and are independent of stellar mass Mstel. For age, the observed luminosity correlation arises because young galaxies are brighter, at fixed Mstel. For metallicity, the observed luminosity dependence arises because metal-rich galaxies, at fixed mass, tend also to be younger, and hence brighter. We find a good match to the observed luminosity correlations with Age ~ sigma^0.40, Z/H ~ sigma^0.35, a/Fe ~ sigma^0.20, where the slopes are close to those found when fitting traditional scaling relations. We conclude that the star formation and enrichment histories of galaxies are determined primarily by the depth of their gravitational potential wells. The observed residual correlations with luminosity do not imply a corresponding dependence on stellar mass.

preprint2008arXiv

Consistently Large Cosmic Flows on Scales of 100 Mpc/h: a Challenge for the Standard LCDM Cosmology

Peculiar velocity surveys have non-uniform spatial distributions of tracers, so that the bulk flow estimated from them does not correspond to that of a simple volume such as a sphere. Thus bulk flow estimates are generally not strictly comparable between surveys, even those whose effective depths are similar. In addition, the sparseness of typical surveys can lead to aliasing of small scale power into what is meant to be a probe of the largest scales. Here we introduce a new method of calculating bulk flow moments where velocities are weighted to give an optimal estimate of the bulk flow of an idealized survey, with the variance of the difference between the estimate and the actual flow being minimized. These "minimum variance" estimates can be designed to estimate the bulk flow on a particular scale with minimal sensitivity to small scale power, and are comparable between surveys. We compile all major peculiar velocity surveys and apply this new method to them. We find that most surveys we studied are highly consistent with each other. Taken together the data suggest that the bulk flow within a Gaussian window of radius 50 Mpc/h is 407 km/s toward l=287 and b=8. The large-scale bulk motion is consistent with predictions from the local density field. This indicates that there are significant density fluctuations on very large scales. A flow of this amplitude on such a large scale is not expected in the WMAP5-normalized LCDM cosmology, for which the predicted one-dimensional r.m.s. velocity is ~110 km/s. The large amplitude of the observed bulk flow favors the upper values of the WMAP5 error-ellipse, but even the point at the top of the WMAP5 95% confidence ellipse predicts a bulk flow which is too low compared to that observed at >98% confidence level.

preprint2007arXiv

The Peculiar Velocities of Local Type Ia Supernovae and their Impact on Cosmology

We quantify the effect of supernova Type Ia peculiar velocities on the derivation of cosmological parameters. The published distant and local Ia SNe used for the Supernova Legacy Survey first-year cosmology report form the sample for this study. While previous work has assumed that the local SNe are at rest in the CMB frame (the No Flow assumption), we test this assumption by applying peculiar velocity corrections to the local SNe using three different flow models. The models are based on the IRAS PSCz galaxy redshift survey, have varying beta = Omega_m^0.6/b, and reproduce the Local Group motion in the CMB frame. These datasets are then fit for w, Omega_m, and Omega_Lambda using flatness or LambdaCDM and a BAO prior. The chi^2 statistic is used to examine the effect of the velocity corrections on the quality of the fits. The most favored model is the beta=0.5 model, which produces a fit significantly better than the No Flow assumption, consistent with previous peculiar velocity studies. By comparing the No Flow assumption with the favored models we derive the largest potential systematic error in w caused by ignoring peculiar velocities to be Delta w = +0.04. For Omega_Lambda, the potential error is Delta Omega_Lambda = -0.04 and for Omega_m, the potential error is Delta Omega_m < +0.01. The favored flow model (beta=0.5) produces the following cosmological parameters: w = -1.08 (+0.09,-0.08), Omega_m = 0.27 (+0.02,-0.02) assuming a flat cosmology, and Omega_Lambda = 0.80 (+0.08,-0.07) and Omega_m = 0.27 (+0.02,-0.02) for a w = -1 (LambdaCDM) cosmology.

preprint2005arXiv

NOAO Fundamental Plane Survey -- II. Age and Metallicity along the Red Sequence

We present spectroscopic linestrength data for 4097 red-sequence galaxies in 93 low-redshift galaxy clusters, and use these to investigate variations in average stellar populations as a function of galaxy mass. Our analysis includes an improved treatment of nebular emission contamination, which affects ~10% of the sample galaxies. Using the stellar population models of D. Thomas and collaborators, we simultaneously fit twelve observed linestrength-sigma relations in terms of common underlying trends of age, [Z/H] (total metallicity) and a/Fe (alpha-element enhancement). We find that the observed linestrength-sigma relations can be explained only if higher-mass red-sequence galaxies are, on average, older, more metal rich, and more alpha-enhanced than lower-mass galaxies. Quantitatively, the scaling relations are age=sigma^(0.59+/-0.13), Z/H=sigma^(0.53+/-0.08) and a/Fe=sigma^(0.31+/-0.06), where the errors reflect the range obtained using different subsets of indices. We conclude that although the stars in giant red galaxies in clusters formed early, most of the galaxies at the faint end joined the red sequence only at recent epochs. This "down-sizing" trend is in good qualitative agreement with observations of the red sequence at higher redshifts, but is not predicted by semi-analytic models of galaxy formation.

preprint2000arXiv

On the dependence of spectroscopic indices of early-type galaxies on age, metallicity and velocity dispersion

We investigate the Mg-sigma and <Fe>-sigma relations in a sample of 72 early-type galaxies drawn mostly from cluster and group environments using a homogeneous data-set which is well-calibrated onto the Lick/IDS system. The small intrinsic scatter in Mg at a given sigma gives upper limits on the spread in age and metallicity of 49% and 32% respectively, if the spread is attributed to one quantity only and if the variations in age and metallicity are uncorrelated. The age/metallicity distribution as inferred from the Hbeta vs <Fe> diagnostic diagram reinforces this conclusion, as we find mostly galaxies with large luminosity weighted ages spanning a range in metallicity. In our sample we do not find significant evidence for an anti-correlation of ages and metallicities which would keep the index-sigma relations tight while hiding a large spread in age and metallicity. As a result of correlated errors in the age-metallicity plane, a mild age-metallicity anti-correlation cannot be completely ruled out given the current data. Correcting the line-strengths indices for non-solar abundance ratios following the recent paper by Trager et al., leads to higher mean metallicity and slightly younger age estimates while preserving the metallicity sequence. The [Mg/Fe] ratio is mildly correlated with the central velocity dispersion and ranges from [Mg/Fe]=0.05 to 0.3 for galaxies with sigma > 100 km/s. Under the assumption that there is no age gradient along the index-sigma relations, the abundance-ratio corrected Mg-sigma, <Fe>-sigma and Hbeta-sigma relations give consistent estimates of Delta [M/H]/ Delta log sigma = 0.9 (+- 0.1). The slope of the Hbeta-sigma relation limits a potential age trend as a function of sigma to 2-3 Gyrs along the sequence.(abriged)