Source author record

August E. Evrard

August E. Evrard 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

33works
7topics
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

33 published item(s)

preprint2022arXiv

HSC-XXL : Baryon budget of the 136 XXL Groups and Clusters

We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass--total mass ($M_{500}$) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass--total mass relation is shallower than unity, $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of $M_{500}$ agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution.

preprint2021arXiv

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available.

preprint2021arXiv

The Role of Machine Learning in the Next Decade of Cosmology

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

preprint2020arXiv

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $γ= -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $γ\simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$.

preprint2016arXiv

Galaxy Cluster Mass Estimation from Stacked Spectroscopic Analysis

We use simulated galaxy surveys to study: i) how galaxy membership in redMaPPer clusters maps to the underlying halo population, and ii) the accuracy of a mean dynamical cluster mass, $M_σ(λ)$, derived from stacked pairwise spectroscopy of clusters with richness $λ$. Using $\sim\! 130,000$ galaxy pairs patterned after the SDSS redMaPPer cluster sample study of Rozo et al. (2015 RMIV), we show that the pairwise velocity PDF of central--satellite pairs with $m_i < 19$ in the simulation matches the form seen in RMIV. Through joint membership matching, we deconstruct the main Gaussian velocity component into its halo contributions, finding that the top-ranked halo contributes $\sim 60\%$ of the stacked signal. The halo mass scale inferred by applying the virial scaling of Evrard et al. (2008) to the velocity normalization matches, to within a few percent, the log-mean halo mass derived through galaxy membership matching. We apply this approach, along with mis-centering and galaxy velocity bias corrections, to estimate the log-mean matched halo mass at $z=0.2$ of SDSS redMaPPer clusters. Employing the velocity bias constraints of Guo et al. (2015), we find $\langle \ln(M_{200c})|λ\rangle = \ln(M_{30}) + α_m \ln(λ/30)$ with $M_{30} = 1.56 \pm 0.35 \times 10^{14} M_\odot$ and $α_m = 1.31 \pm 0.06_{stat} \pm 0.13_{sys}$. Systematic uncertainty in the velocity bias of satellite galaxies overwhelmingly dominates the error budget.

preprint2016arXiv

Host Galaxy Identification for Supernova Surveys

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.

preprint2016arXiv

RHAPSODY-G simulations II - Baryonic growth and metal enrichment in massive galaxy clusters

We study the evolution of the stellar component and the metallicity of both the intracluster medium and of stars in massive ($M_{\rm vir}\approx 6\times 10^{14}$ M$_{\odot}/h$) simulated galaxy clusters from the Rhapsody-G suite in detail and compare them to observational results. The simulations were performed with the AMR code RAMSES and include the effect of AGN feedback at the sub-grid level. AGN feedback is required to produce realistic galaxy and cluster properties and plays a role in mixing material in the central regions and regulating star formation in the central galaxy. In both our low and high resolution runs with fiducial stellar yields, we find that stellar and ICM metallicities are a factor of two lower than in observations. We find that cool core clusters exhibit steeper metallicity gradients than non-cool core clusters, in qualitative agreement with observations. We verify that the ICM metallicities measured in the simulation can be explained by a simple "regulator" model in which the metallicity is set by a balance of stellar yield and gas accretion. It is plausible that a combination of higher resolution and higher metal yield in AMR simulation would allow the metallicity of simulated clusters to match observed values; however this hypothesis needs to be tested with future simulations. Comparison to recent literature highlights that results concerning the metallicity of clusters and cluster galaxies might depend sensitively on the scheme chosen to solve the hydrodynamics.

preprint2015arXiv

Persistent, Global Identity for Scientists via ORCID

Scientists have an inherent interest in claiming their contributions to the scholarly record, but the fragmented state of identity management across the landscape of astronomy, physics, and other fields makes highlighting the contributions of any single individual a formidable and often frustratingly complex task. The problem is exacerbated by the expanding variety of academic research products and the growing footprints of large collaborations and interdisciplinary teams. In this essay, we outline the benefits of a unique scholarly identifier with persistent value on a global scale and we review astronomy and physics engagement with the Open Researcher and Contributor iD (ORCID) service as a solution.

preprint2015arXiv

Rhapsody-G simulations: galaxy clusters as baryonic closed boxes and the covariance between hot gas and galaxies

Within a sufficiently large cosmic volume, conservation of baryons implies a simple `closed box' view in which the sum of the baryonic components must equal a constant fraction of the total enclosed mass. We present evidence from Rhapsody-G hydrodynamic simulations of massive galaxy clusters that the closed-box expectation may hold to a surprising degree within the interior, non-linear regions of haloes. At a fixed halo mass, we find a significant anti-correlation between hot gas mass fraction and galaxy mass fraction (cold gas + stars), with a rank correlation coefficient of -0.69 within $R_{500c}$. Because of this anti-correlation, the total baryon mass serves as a low-scatter proxy for total cluster mass. The fractional scatter of total baryon fraction scales approximately as $0.02 (Δ_c/100)^{0.6}$, while the scatter of either gas mass or stellar mass is larger in magnitude and declines more slowly with increasing radius. We discuss potential observational tests using cluster samples selected by optical and hot gas properties; the simulations suggest that joint selection on stellar and hot gas has potential to achieve 5% scatter in total halo mass.

preprint2014arXiv

A Model for Multi-property Galaxy Cluster Statistics

The massive dark matter halos that host groups and clusters of galaxies have observable properties that appear to be log-normally distributed about power-law mean scaling relations in halo mass. Coupling this assumption with either quadratic or cubic approximations to the mass function in log space, we derive closed-form expressions for the space density of halos as a function of multiple observables as well as forms for the low-order moments of properties of observable-selected samples. Using a Tinker mass function in a ΛCDM cosmology, we show that the cubic analytic model reproduces results obtained from direct, numerical convolution at the 10 percent level or better over nearly the full range of observables covered by current observations and for redshifts extending to z = 1.5. The model provides an efficient framework for estimating effects arising from selection and covariance among observable properties in survey samples.

preprint2013arXiv

Cluster Cosmology at a Crossroads: Neutrino Masses

Galaxy clusters --- in combination with CMB and BAO data --- can provide precise constraints on the sum of neutrino masses. However, these constraints depend on the calibration of the mass--observable relation. For instance, the mass calibration employed in Planck Collaboration (2011a,b) rules out the minimal 6-parameter ΛCDM model at 3.7σ, and implies a sum of neutrino masses \sum m_ν= 0.39 \pm 0.10. By contrast, the mass calibration favored by Rozo et al. (2012b) from a self-consistent analysis of X-ray, SZ, and optical scaling relations is consistent with a minimal flat ΛCDM model with no massive neutrinos (1.7σ), and is a better fit to additional data (e.g. H0). We discuss these results in light of the recent SPT and ACT analysis, noting that the Rozo et al. (2012b) model suggests the current mild "tension" (<2σ) between CMB and BAO+H_0 data will decrease as the uncertainties in these measurements decrease.

preprint2013arXiv

Problem Roulette: Studying Introductory Physics in the Cloud

Problem Roulette (PR) is a web-based study service at the University of Michigan that offers random-within-topic access to a large library of past exam problems in introductory physics courses. Built on public-private cloud computing infrastructure, PR served nearly 1000 students during Fall 2012 term, delivering more than 60,000 problem pages. The service complements that of commercial publishing houses by offering problems authored by local professors. We describe the service architecture, including reporting and analytical capabilities, and present an initial evaluation of the impact of its use. Among roughly 500 students studying electromagnetism, we find that the 229 students who worked fifty or more problems over the term outperformed their complement by 0.40 grade points (on a 4.0 scale). This improvement partly reflects a selection bias that better students used the service more frequently. Adjusting for this selection bias by using overall grade point average (GPA) as a classifier, we find a grade point improvement of 0.22 for regular PR users, significantly above the random noise level of 0.04. Simply put, students who worked one or more additional problem per weekday earned nearly a quarter-letter average grade improvement irrespective of GPA. Student comments emphasize the importance of randomness in helping them synthesize concepts. The PR source code is publicly available on the bitbucket repository.

preprint2013arXiv

Virial scaling of galaxies in clusters: bright to faint is cool to hot

By combining galaxy tracers from high-resolution N-body and hydrodynamical simulations, we present a consistent picture of the behaviour of galaxy velocities in massive clusters. In haloes above ~ 10^14 Msun, the brightest satellite galaxies are slightly cooler compared to the dark matter, while fainter satellites are hotter. Within the virial radius of a cluster, the mean velocity dispersion based on the 100 brightest galaxies is a factor of 1.065 +/- 0.005 (stat) +/- 0.027 (sys) higher than that of the dark matter (corresponding to a ~10-15 per cent bias in the dynamical mass estimate) while that based on only the five brightest galaxies is 0.868 +/- 0.039 (stat) +/- 0.035 (sys). These trends are approximately independent of redshift. The velocity structure is sensitive to the modelling of galaxies in clusters, indicative of the complex interplay of tidal stripping, dynamical friction, and merging. Velocity dispersions derived from instantaneous subhalo properties are larger than those employing either peak subhalo properties or hydrodynamical galaxy tracers. The latter two methods are consistent, implying that stacked spectroscopic analysis of cluster samples should, after correction for projection, show a trend towards slightly higher velocities when fainter galaxies are included, with an unbiased measure of dark matter velocity dispersion coming from approximately 30 galaxies per cluster. We show evidence that the velocity distribution function of bright galaxies near the cluster centre has a low-velocity tail due to strong dynamical friction.

preprint2012arXiv

A Comparative Study of Local Galaxy Clusters: I. Derived X-ray Observables

We examine systematic differences in the derived X-ray properties of galaxy clusters as reported by three different groups: Vikhlinin et al. (2009a), Mantz et al. (2010b), and Planck Collaboration (2011b). The sample overlap between any two pairs of works ranges between 16 to 28 galaxy clusters in common. We find systematic differences in most reported properties, including the total cluster mass, M500. The most extreme case is an average 45% \pm 5% difference in cluster mass between the Planck Collaboration (2011b) and Mantz et al. (2010b), for clusters at z > 0.13 (averaged over 16 clusters). These mass differences induce differences in cluster observables defined within an R500 aperture. After accounting for aperture differences, we find very good agreement in gas mass estimates between the different groups. However, the soft-band X-ray luminosity, LX, core-excised spectroscopic temperature, TX, and gas thermal energy, YX = MgasTX display mean differences at the 5%-15% level. We also find that the low (z \leq 0.13) and high (z \geq 0.13) galaxy cluster samples in Planck Collaboration (2011b) appear to be systematically different: the YSZ/YX ratio for these two sub- samples is ln(YSZ/YX) = -0.06 \pm 0.04 and ln(YSZ/YX) = 0.08 \pm 0.04 for the low and high redshift sub-samples respectively.

preprint2012arXiv

A Comparative Study of Local Galaxy Clusters: II: X-ray and SZ Scaling Relations

We compare cluster scaling relations published for three different samples selected via X-ray and Sunyaev-Zel'dovich (SZ) signatures. We find tensions driven mainly by two factors: i) systematic differences in the X-ray cluster observables used to derive the scaling relations, and ii) uncertainty in the modeling of how the gas mass of galaxy clusters scales with total mass. All scaling relations are in agreement after accounting for these two effects. We describe a multivariate scaling model that enables a fully self-consistent treatment of multiple observational catalogs in the presence of property covariance, and apply this formalism when interpreting published results. The corrections due to scatter and observable covariance can be significant. For instance, our predicted Ysz-Lx scaling relation differs from that derived using the naive "plug in" method by \approx 25%. Finally, we test the mass normalization for each of the X-ray data sets we consider by applying a space density consistency test: we compare the observed REFLEX luminosity function to expectations from published Lx-M relations convolved with the mass function for a WMAP7 flat ΛCDM model.

preprint2012arXiv

A High Throughput Workflow Environment for Cosmological Simulations

The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at the percent-level in order to fully leverage the power of next-generation catalogs. Simulations of large-scale structure provide the means to characterize these uncertainties. We are using XSEDE resources to produce multiple synthetic sky surveys of galaxies and large-scale structure in support of science analysis for the Dark Energy Survey. In order to scale up our production to the level of fifty 10^10-particle simulations, we are working to embed production control within the Apache Airavata workflow environment. We explain our methods and report how the workflow has reduced production time by 40% compared to manual management.

preprint2012arXiv

Closing the Loop: A Self-Consistent Model of Optical, X-ray, and SZ Scaling Relations for Clusters of Galaxies

We demonstrate that optical data from SDSS, X-ray data from ROSAT and Chandra, and SZ data from Planck, can be modeled in a fully self-consistent manner. After accounting for systematic errors and allowing for property covariance, we find that scaling relations derived from optical and X-ray selected cluster samples are consistent with one another. Moreover, these clusters scaling relations satisfy several non-trivial spatial abundance constraints and closure relations. Given the good agreement between optical and X-ray samples, we combine the two and derive a joint set of LX-M and YSZ-M relations. Our best fit YSZ-M relation is in good agreement with the observed amplitude of the thermal SZ power spectrum for a WMAP7 cosmology, and is consistent with the masses for the two CLASH galaxy clusters published thus far. We predict the halo masses of the remaining z \leq 0.4 CLASH clusters, and use our scaling relations to compare our results with a variety of X-ray and weak lensing cluster masses from the literature.

preprint2011arXiv

Cosmological Parameters from Observations of Galaxy Clusters

Studies of galaxy clusters have proved crucial in helping to establish the standard model of cosmology, with a universe dominated by dark matter and dark energy. A theoretical basis that describes clusters as massive, multi-component, quasi-equilibrium systems is growing in its capability to interpret multi-wavelength observations of expanding scope and sensitivity. We review current cosmological results, including contributions to fundamental physics, obtained from observations of galaxy clusters. These results are consistent with and complementary to those from other methods. We highlight several areas of opportunity for the next few years, and emphasize the need for accurate modeling of survey selection and sources of systematic error. Capitalizing on these opportunities will require a multi-wavelength approach and the application of rigorous statistical frameworks, utilizing the combined strengths of observers, simulators and theorists.

preprint2011arXiv

Influence of Projection in Cluster Cosmology Studies

Projection tends to skew the mass-observable relation of galaxy clusters by creating a small fraction of severely blended systems, those for which the measured observable property of a cluster is strongly boosted relative to the value of its primary host halo. We examine the bias in cosmological parameter estimates caused by incorrectly assuming a Gaussian (projection-free) mass-observable relation when the true relation is non-Gaussian due to projection. We introduce a mixture model for projection and explore Fisher forecasts for a survey of 5000 sq. deg. to z=1.1 and an equivalent mass threshold of 10^13.7 h^-1 solar masses. Using a blended fraction motivated by optical cluster finding applied to the Millennium Simulation and applying Planck and otherwise weak priors, we find that the biases in Omega_DE and w are significant, being factors of 2.8 and 2.4, respectively, times previous forecast uncertainties. Incorporating eight new degrees of freedom to describe cluster selection with projection increases the forecast uncertainty in Omega_DE and w by similar factors. Knowledge of these additional parameters at the 5% level limits degradation in dark energy constraints to <10% relative to projection-free forecasts. We discuss strategies for using simulations and complementary observations to characterize the fraction of blended clusters and their mass selection properties.

preprint2011arXiv

The Evolution of Radio Galaxies and X-ray Point Sources in Coma Cluster Progenitors Since z~1.2

Using Chandra imaging spectroscopy and Very Large Array (VLA) L-band radio maps, we have identified radio sources at P_{1.4GHz} >=5x10^{23} W Hz^{-1} and X-ray point sources (XPSs) at L_{0.3-8keV}>=5x10^{42} erg s^{-1} in L>L* galaxies in 12 high-redshift (0.4<z<1.2) clusters of galaxies. The radio galaxies and XPSs in this cluster sample, chosen to be consistent with Coma Cluster progenitors at these redshifts, are compared to those found at low-z analyzed in Hart et al. (2009). Within a projected radius of 1 Mpc of the cluster cores, we find 17 cluster radio galaxies (11 with secure redshifts, including one luminous FR II radio source at z=0.826, and 6 more with host galaxy colors similar to cluster ellipticals). The radio luminosity function (RLF) of the cluster radio galaxies as a fraction of the cluster red sequence (CRS) galaxies reveals significant evolution of this population from high-z to low-z, with higher power radio galaxies situated in lower temperature clusters at earlier epochs. Additionally, there is some evidence that cluster radio galaxies become more centrally concentrated than CRS galaxies with cosmic time. Within this same projected radius, we identify 7 spectroscopically-confirmed cluster XPSs, all with CRS host galaxy colors. Consistent with the results from Martini et al. (2009), we estimate a minimum X-ray active fraction of 1.4+/-0.8% for CRS galaxies in high-z clusters, corresponding to an approximate 10-fold increase from 0.15+/-0.15% at low-z. Although complete redshift information is lacking for several XPSs in z>0.4 cluster fields, the increased numbers and luminosities of the CRS radio galaxies and XPSs suggest a substantial (9-10 fold) increase in the heat injected into high redshift clusters by AGN compared to the present epoch.

preprint2010arXiv

LoCuSS: A Comparison of Cluster Mass Measurements from XMM-Newton and Subaru - Testing Deviation from Hydrostatic Equilibrium and Non-Thermal Pressure Support

We compare X-ray hydrostatic and weak-lensing mass estimates for a sample of 12 clusters that have been observed with both XMM-Newton and Subaru. At an over-density of Δ=500, we obtain 1-M^{X}/M^{WL}=0.01+/-0.07 for the whole sample. We also divided the sample into undisturbed and disturbed sub-samples based on quantitative X-ray morphologies using asymmetry and fluctuation parameters, obtaining 1-M^{X}/M^{WL}=0.09+/-0.06 and -0.06+/-0.12 for the undisturbed and disturbed clusters, respectively. In addition to non-thermal pressure support, there may be a competing effect associated with adiabatic compression and/or shock heating which leads to overestimate of X-ray hydrostatic masses for disturbed clusters, for example, in the famous merging cluster A1914. Despite the modest statistical significance of the mass discrepancy, on average, in the undisturbed clusters, we detect a clear trend of improving agreement between M^{X} and M^{WL} as a function of increasing over-density, M^{X}/M^{WL}=(0.908+/-0.004)+(0.187+/-0.010) \log_{10}(Δ/500). We also examine the gas mass fractions, f_{gas}=M^{gas}/M^{WL}, finding that they are an increasing function of cluster radius, with no dependence on dynamical state, in agreement with predictions from numerical simulations. Overall, our results demonstrate that XMM-Newton and Subaru are a powerful combination for calibrating systematic uncertainties in cluster mass measurements.

preprint2010arXiv

Sensitivity of galaxy cluster dark energy constraints to halo modeling uncertainties

We perform a sensitivity study of dark energy constraints from galaxy cluster surveys to uncertainties in the halo mass function, bias and the mass-observable relation. For a set of idealized surveys, we evaluate cosmological constraints as priors on sixteen nuisance parameters in the halo modeling are varied. We find that surveys with a higher mass limit are more sensitive to mass-observable uncertainties while surveys with low mass limits that probe more of the mass function shape and evolution are more sensitive to mass function errors. We examine the correlations among nuisance and cosmological parameters. Mass function parameters are strongly positively (negatively) correlated with Omega_DE (w). For the mass-observable parameters, Omega_DE is most sensitive to the normalization and its redshift evolution while w is more sensitive to redshift evolution in the variance. While survey performance is limited mainly by mass-observable uncertainties, the current level of mass function error is responsible for up to a factor of two degradation in ideal cosmological constraints. For surveys that probe to low masses (10^13.5 h^-1 M_sun), even percent-level constraints on model nuisance parameters result in a degradation of ~ sqrt{2} (2) on Omega_DE (w) relative to perfect knowledge.

preprint2009arXiv

Accretion Shocks in Clusters of Galaxies and their SZ Signature from Cosmological Simulations

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

preprint2009arXiv

Cosmological Constraints from the SDSS maxBCG Cluster Catalog

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find σ_8(Ω_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find σ_8=0.807\pm 0.020 and Ω_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

preprint2008arXiv

The X-ray Cluster Normalization of the Matter Power Spectrum

The number density of galaxy clusters provides tight statistical constraints on the matter fluctuation power spectrum normalization, traditionally phrased in terms of sigma_8, the root mean square mass fluctuation in spheres with radius 8 h^-1 Mpc. We present constraints on sigma_8 and the total matter density Omega_m0 from local cluster counts as a function of X-ray temperature, taking care to incorporate and minimize systematic errors that plagued previous work with this method. In particular, we present new determinations of the cluster luminosity - temperature and mass - temperature relations, including their intrinsic scatter, and a determination of the Jenkins mass function parameters for the same mass definition as the mass - temperature calibration. Marginalizing over the 12 uninteresting parameters associated with this method, we find that the local cluster temperature function implies sigma_8 (Omega_m0/0.32)^alpha = 0.86+/-0.04 with alpha = 0.30 (0.41) for Omega_m0 < 0.32 (Omega_mo > 0.32) (68% confidence for two parameters). This result agrees with a wide range of recent independent determinations, and we find no evidence of any additional sources of systematic error for the X-ray cluster temperature function determination of the matter power spectrum normalization. The joint WMAP5 + cluster constraints are: Omega_m0 = 0.30+0.03/-0.02 and sigma_8 = 0.85+0.04/-0.02 (68% confidence for two parameters).

preprint2007arXiv

MaxBCG: A Red Sequence Galaxy Cluster Finder

Measurements of galaxy cluster abundances, clustering properties, and mass to- light ratios in current and future surveys can provide important cosmological constraints. Digital wide-field imaging surveys, the recently-demonstrated fidelity of red-sequence cluster detection techniques, and a new generation of realistic mock galaxy surveys provide the means for construction of large, cosmologicallyinteresting cluster samples, whose selection and properties can be understood in unprecedented depth. We present the details of the "maxBCG" algorithm, a cluster-detection technique tailored to multi-band CCD-imaging data. MaxBCG primarily relies on an observational cornerstone of massive galaxy clusters: they are marked by an overdensity of bright, uniformly red galaxies. This detection scheme also exploits classical brightest cluster galaxies (BCGs), which are often found at the center of these same massive clusters. (ABRIDGED)

preprint2006arXiv

Cluster-Supercluster Alignments

We study correlations in spatial orientation between galaxy clusters and their host superclusters using a Hubble Volume N-body realization of a concordance cosmology and an analytic model for tidally-induced alignments. We derive an analytic form for distributions of the alignment angle as functions of halo mass (M), ellipticity (epsilon), distance (r) and velocity (v) and show that the model, after tuning of three parameters, provides a good fit to the numerical results. The parameters indicate a high degree of alignment along anisotropic, collapsed filaments. The degree of alignment increases with M and epsilon while it decreases with r and is independent of v. We note the possibility of using the cluster-supercluster alignment effect as a cosmological probe to constrain the slope of the initial power spectrum.

preprint2006arXiv

The Asymptotic Form of Cosmic Structure: Small Scale Power and Accretion History

We explore the effects of small scale structure on the formation and equilibrium of dark matter halos in a universe dominated by vacuum energy. We present the results of a suite of four N-body simulations, two with a LCDM initial power spectrum and two with WDM-like spectra that suppress the early formation of small structures. All simulations are run into to far future when the universe is 64Gyr/h old, long enough for halos to essentially reach dynamical equilibrium. We quantify the importance of hierarchical merging on the halo mass accretion history, the substructure population, and the equilibrium density profile. We modify the mass accretion history function of Wechsler et al. (2002) by introducing a parameter, γ, that controls the rate of mass accretion, dln(M) / dln(a) ~ a^(-γ), and find that this form characterizes both hierarchical and monolithic formation. Subhalo decay rates are exponential in time with a much shorter time scale for WDM halos. At the end of the simulations, we find truncated Hernquist density profiles for halos in both the CDM and WDM cosmologies. There is a systematic shift to lower concentration for WDM halos, but both cosmologies lie on the same locus relating concentration and formation epoch. Because the form of the density profile remains unchanged, our results indicate that the equilibrium halo density profile is set independently of the halo formation process.

preprint2004arXiv

The Ultimate Halo Mass in a LCDM Universe

In the far future of an accelerating LCDM cosmology, the cosmic web of large-scale structure consists of a set of increasingly isolated halos in dynamical equilibrium. We examine the approach of collisionless dark matter to hydrostatic equilibrium using a large N-body simulation evolved to scale factor a = 100, well beyond the vacuum--matter equality epoch, a_eq ~ 0.75, and 53/h Gyr into the future for a concordance model universe (Omega_m ~ 0.3, Omega_Lambda ~ 0.7). The radial phase-space structure of halos -- characterized at a < a_eq by a pair of zero-velocity surfaces that bracket a dynamically active accretion region -- simplifies at a > 10 a_eq when these surfaces merge to create a single zero-velocity surface, clearly defining the halo outer boundary, rhalo, and its enclosed mass, mhalo. This boundary approaches a fixed physical size encompassing a mean interior density ~ 5 times the critical density, similar to the turnaround value in a classical Einstein-deSitter model. We relate mhalo to other scales currently used to define halo mass (m200, mvir, m180b) and find that m200 is approximately half of the total asymptotic cluster mass, while m180b follows the evolution of the inner zero velocity surface for a < 2 but becomes much larger than the total bound mass for a > 3. The radial density profile of all bound halo material is well fit by a truncated Hernquist profile. An NFW profile provides a somewhat better fit interior to r200 but is much too shallow in the range r200 < r < rhalo.

preprint2003arXiv

Future Evolution of Structure in an Accelerating Universe

Current cosmological data indicate that our universe contains a substantial component of dark vacuum energy that is driving the cosmos to accelerate. We examine the immediate and longer term consequences of this dark energy (assumed here to have a constant density). Using analytic calculations and supporting numerical simulations, we present criteria for test bodies to remain bound to existing structures. We show that collapsed halos become spatially isolated and dynamically relax to a particular density profile with logarithmic slope steeper than -3 at radii beyond r_200. The asymptotic form of the space-time metric is then specified. We develop this scenario further by determining the effects of the accelerating expansion on the background radiation fields and individual particles. In an appendix, we generalize these results to include quintessence.

preprint2001arXiv

Clustering of dark matter halos on the light-cone: scale-, time- and mass-dependence of the halo biasing in the Hubble volume simulations

We develop a phenomenological model to predict the clustering of dark matter halos on the light-cone by combining several existing theoretical models. Assuming that the velocity field of halos on large scales is approximated by linear theory, we propose an empirical prescription of a scale-, mass-, and time-dependence of halo biasing. We test our model against the Hubble Volume $N$-body simulation and examine its validity and limitations. We find a good agreement in two-point correlation functions of dark matter halos between the phenomenological model predictions and measurements from the simulation for $R>5h^{-1}$Mpc both in the real and redshift spaces. Although calibrated on the mass scale of groups and clusters and for redshifts up to $z\sim2$, the model is quite general and can be applied to a wider range of astrophysical objects, such as galaxies and quasars, if the relation between dark halos and visible objects is specified.

preprint1997arXiv

The Intracluster Gas Fraction in X-ray Clusters : Constraints on the Clustered Mass Density

The mean intracluster gas fraction of X-ray clusters within their hydrostatic regions is derived from recent observational compilations of David, Jones & Forman and White & Fabian. At radii encompassing a mean density 500 times the critical value, the individual sample bi-weight means are moderately ($2.4 σ$) discrepant; revising binding masses with a virial relation calibrated by numerical simulations removes the discrepancy and results in a combined sample mean and standard error $\fbargas(r_{500}) = (0.060 \pm 0.003) \h^{-3/2}$. For hierarchical clustering models with an extreme physical assumption to maximize cluster gas content, this value constrains the universal ratio of total, clustered to baryonic mass $Ω_m/Ω_b \le 23.1 \h^{3/2}$; combining with the primordial nucleosynthesis upper limit on $Ω_b$ results in $Ω_m \h^{1/2} < 0.60$. A less conservative, physically plausible approach based on low D/H inferences from quasar absorption spectra and accounting for baryons within cluster galaxies yields an estimate $Ω_m \h^{2/3} = 0.28 \pm 0.07$ with sources of systematic error involved in the derivation providing approximately $35%$ uncertainty. Additional effects which could provide consistency with the Einstein-deSitter case $Ω_m \se 1$ are presented, and their observable signatures discussed.

preprint1995arXiv

ROSAT Observations of Compact Groups of Galaxies

We have systematically analyzed a sample of 13 new and archival ROSAT PSPC observations of compact groups of galaxies: 12 Hickson Compact Groups plus the NCG 2300 group. We find that approximately two-thirds of the groups have extended X-ray emission and, in four of these, the emission is resolved into diffuse emission from gas at a temperature of $k$T $\sim 1$ keV in the group potential. All but one of the groups with extended emission have a spiral fraction of less than 50\%. The baryon fraction of groups with diffuse emission is 5--19\%, similar to the values in clusters of galaxies. However, with a single exception (HCG 62), the gas-to-stellar-mass ratio in our groups has a median value near 5\%, somewhat greater than the values for individual early-type galaxies and two orders of magnitude lower than in clusters of galaxies. The X-ray luminosities of individual group galaxies are comparable to those of similar field galaxies, although the L$_X$-L$_B$ relation for early-type galaxies may be flatter in compact groups than in the field.