Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the statistical character of our sky s population of clusters of galaxies , viewed from radio to gamma - ray wavelengths , is sensitive to models of cosmology , astrophysics , and large - scale gravity . galaxy clusters are cosmographic buoys that signal locations of peaks in the large - scale matter density . the population is shallow and finite . surveys in the coming decades will definitively map our universe s terrain as defined by the highest @xmath0 peaks . current maps have advanced to the stage where abell 2163 , a cluster at redshift @xmath1 with a plasma virial temperature @xmath2 @xcite and galaxy velocity dispersion @xmath3 @xcite , has been nominated a candidate for _ the _ most massive cluster in the universe @xcite , the cosmic equivalent of mount everest . physically , galaxy clusters are manifested in the most massive of the bound structures  termed _ halos _ ( or _ haloes _ )  that emerge in the cosmic web of large - scale structure ( lss ) . the lss web is a gravitationally amplified descendant of a weak noise field seeded by quantum fluctuations during an early , inflationary epoch @xcite . its evolutionary dynamics have been well studied into the non - linear regime by n - body simulations @xcite . locally bound regions ( the halos ) emerge , initially via coherent infall within a narrow mass range and , subsequently , via a combination of infall and hierarchical merging that widen the dynamic range , pushing to increasingly larger halo masses . the merging process is of considerable interest for cluster studies , driving astrophysical signatures that can test physical models from the nature of dark matter @xcite to the magnetohydrodynamics of hot , dilute plasmas ( e.g. * ? ? ? merging also potentially confuses cosmological studies , by creating close halo pairs that may appear as one cluster in projection and by introducing variance into observable signals . halos are multi - component systems consisting of dark matter and baryons in several phases : black holes ; stars ; cold , molecular gas ; warm / hot gas ; and non - thermal plasma . after decades of study via n - body and hydrodynamic simulation and related methods ( see recent review by * ? ? ? * ) , models for the detailed evolution of the baryons in clusters are growing in capability to describe an increasingly large and rich volume of observations . what is clear empirically is that the galaxy formation process is globally inefficient : a recent study by @xcite finds that stellar mass accounts for only @xmath4 percent of the total baryon budget in the most massive halos . radiative cooling of gas is overcome by feedback from various sources , including mechanical and radiative input from supernova winds and black hole jets , thermal conduction and other plasma processes , and ablation and harassment during gravitational encounters . while the hierarchical nature of structure formation implies that galaxy and cluster formation are deeply intertwined and , therefore , that detailed understanding of cluster structure and evolution _ requires _ that we understand galaxy formation , the scales separating the most massive clusters from the largest galaxies  roughly a factor of 100 in length and 1000 in mass  allow progress to be made by approximate physical treatments . the dark matter kinematic structure , including remnant , fine - scale _ sub - halos _ @xcite , as well as the morphology and scaling behaviors of the hot , intracluster medium ( icm ) that dominates the baryonic component @xcite , are examples of areas where direct simulations made good , early progress . a key aspect of their multi - component nature is the fact that clusters offer multiple , observable signals across the electromagnetic spectrum ( e.g. * ? ? ? * ) . at x - ray wavelengths , the hot icm emits thermal bremsstrahlung and line emission from ionized metals injected into the plasma by stripping and feedback processes . stellar emission from galaxies and intracluster light dominates the optical and near - infrared . at millimeter wavelengths , inverse compton scattering within clusters distorts the spectrum of the cosmic microwave background ( cmb ) . gravitational lensing offers a unique probe into the total matter distributions in clusters . synchrotron emission from relativistic electrons is visible at radio frequencies . these and other signatures discussed below provide physically coupled , and often observationally independent , lines of evidence with which to test astrophysical models of cluster evolution . a challenge to cluster cosmology is the construction of accurate statistical models that address survey observables explicitly while incorporating intrinsic property covariance . the use of clusters to study cosmology has a history dating to zwicky s discovery of dark matter in the coma cluster @xcite . brightest cluster galaxies were later employed as standard candles to study the local expansion history of the universe ; @xcite actually derived ( with low significance ) a negative deceleration parameter using this approach , implying accelerated expansion consistent with present findings . in the 1980 s , measurement of the enhanced spatial clustering of clusters relative to galaxies supported the model of gaussian random initial conditions expected from inflation @xcite . in the early 1990 s , an apparent discrepancy between local baryon fraction measurements of clusters @xcite with primordial nucleosynthesis expectations helped rule out a model with critical matter density @xcite . the revelation of hot clusters at high redshift later that decade @xcite presaged the ultimate discovery of dark energy from type ia supernova ( snia ) surveys . the turn of the millennium witnessed a flurry of activity aimed at measuring the amplitude of the matter power spectrum from cluster counts . x - ray studies in particular showed that the amplitude was lower than had been accepted previously ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , a result later confirmed by cmb and cosmic shear measurements . these studies also exposed the importance of understanding systematic effects associated with the use of directly observable quantities as proxies for mass @xcite . recent studies have used cluster counts or the icm mass fraction in very massive systems ( both methods described in more detail below ) to constrain cosmological parameters . these studies are consistent with other observations that find a universe dominated by dark energy ( @xmath5 ) , with sub - dominant dark matter ( @xmath6 ) , and a small minority of baryonic material ( @xmath7 ) @xcite . a detailed pedagogical treatment of how cluster studies helped establish this reference cosmology is given in the review of @xcite . @xcite review x - ray studies of clusters from the rosat satellite era .    explaining the nature of the dark energy and dark matter are core problems of physics . the consensus ` concordance ' cosmological model , @xmath8cdm , postulates that dark energy ( de ) is associated with a small , non - zero vacuum energy , equivalent to a cosmological constant term in einstein s equations . another possibility is that de arises from a light scalar field ( or fields ) that evolves over cosmic time . a third option is that de is essentially an apparition , not a source term of einstein s general relativistic equations but a reflection of their breakdown at length and time scales of cosmic dimensions ( e.g. * ? ? ? sky surveys of cosmic systems , from supernovae to galaxies to clusters of galaxies , provide the means to discriminate among these alternatives . forthcoming cluster surveys at mm , optical / near - infrared , and x - ray wavelengths , discussed in section  [ sec : new_surveys ] , have the potential to find hundreds of thousands of groups and clusters . figure  [ fig : surveyhistory ] puts these efforts into historical perspective , by plotting size against year of publication for cluster samples that generated cosmological constraints discussed in this review . symbol size is proportional to median sample redshift , and symbol types encode the selection method . the stars at far right show theoretical estimates of the all - sky number and median redshift of halos with masses above @xmath9 and @xmath10 . the former mass limit roughly marks the transition from galaxy groups to galaxy clusters , while the latter marks the deepest potential wells with icm temperatures @xmath11 . current surveys have made good progress , but the full population of clusters remains largely undiscovered . optical and x - ray surveys have the longest histories , but these traditional methods are being complemented by new approaches . space - based surveys in the near - infrared extend optical methods to @xmath12 @xcite , and the first few clusters identified by their gravitational lensing signature have been published @xcite . ongoing mm surveys have released the first sets of clusters discovered through the sunyaev - zeldovich ( sz ) effect @xcite , with the promise of much more to come . panoramic , multi - wavelength surveys of common sky areas offer profound improvements to our understanding of clusters as astrophysical systems , which in turn further empowers their use for cosmological studies . and while considerable challenges to interpretation and modeling of survey data certainly exist , a halo model framework , discussed in section  [ sec : theory ] , is rising to meet this task .      a feature common to many techniques that study de is the nature of the input data , which consists of catalogs of properties , @xmath13 , of discrete objects that lie along our past light - cone . upcoming wide - field surveys will generate @xmath13-catalogs of large dimension that will be distilled to constrain perhaps tens of cosmological and astrophysical parameters . such catalogs may contain internal support through the use of complementary methods : besides galaxy clusters ( cl ) , the same data set can be analyzed for baryon acoustic oscillations ( bao ) and , for optical surveys , weak lensing ( wl ) . ( in the case of repeat observations , optical surveys can also be analyzed for type ia supernovae and gravitational time delay signatures . ) science processing leads to a compressed set of statistical signals , @xmath14 , where @xmath15 indicates an aforementioned method . for large cluster surveys , @xmath16 might consist of counts of clusters binned by sky area , redshift and detected signal . extracting accurate constraints on a set of cosmological model parameters , @xmath17 , from these surveys requires sophisticated likelihood analyses . the critical ingredient is @xmath18 , the underlying likelihood that the cl ( and other ) statistics of the observed sky would be realized within a particular universe . key capabilities that enable such likelihood analysis are :    1 .   to predict statistical expectations , @xmath18 , for _ many _ universes , @xmath17 ; + 2 .   to extract _ unbiased _ statistical signals from the raw catalog , @xmath19 ; + 3 .   to understand the expected signal covariance , cov(@xmath20    genuine understanding of cosmological models from observed cluster data is dependent on the degree to which theory and simulation can provide robust predictions for the observed signals . while numerical simulations of lss can predict catalog - level yields for a given cosmology ( e.g. * ? ? * ) , such predictions necessarily entail additional astrophysical assumptions , meaning @xmath18 is actually @xmath21 , where @xmath22 represents degrees of freedom introduced by an assumed astrophysical model . recovering cosmological information from survey data therefore necessitates marginalization over a reasonable range of astrophysical assumptions . on the other hand , as cosmological constraints from all methods improve , the cluster community can potentially invert the problem , recovering constraints on astrophysical models after marginalizing over cosmology . we begin this review by describing the theoretical basis for cluster cosmology ( section  [ sec : theory ] ) , and include there an opening discussion of important sources of systematic error . key observational windows are described in section  [ sec : observations ] , and recent cosmological constraints are reviewed in section  [ sec : cosmo_constraints ] . in section [ sec : physics ] , cluster contributions to particle physics and gravity are examined . in section  [ sec : opportunities ] , we highlight opportunities for important , near - term progress . in closing , we emphasize some essential considerations in survey modeling and analysis ( section  [ sec : modeling ] ) before presenting our conclusions ( section  [ sec : conclusions ] ) . this section sketches a theoretical description of the halo framework that supports cluster cosmology . there is considerable richness to the galaxy formation problem that we omit here ; the recent review of @xcite provides substantial detail . simulation studies of halo evolution into the strongly non - linear regime are becoming increasingly powerful , but finite resolution and uncertainty in astrophysical treatments limit predictive power . while not yet robust enough to offer sharp prior characterization of the astrophysics required for cosmological studies , simulations offer key insights into the structure and physics sensitivity of the functions that relate observable signals to halo mass and epoch . ample evidence now supports the picture that lss formed via gravitational amplification of initially small density fluctuations , @xmath23 . cosmic microwave background anisotropy measurements are consistent with expectations from a large class of basic inflationary models ( e.g. * ? ? ? such models are characterized by an instantaneous primordial power spectrum , @xmath24 , with spectral index , @xmath25 , expected to be close to unity . here , @xmath26 is the fourier transform of the density fluctuations , @xmath27 . after inflation ceases , fluctuations in the coupled photon - baryon - dark matter fluid evolve in ways that are now well understood from linearized boltzmann treatments @xcite . for the standard case of adiabatic fluctuations , and on scales above the baryon jeans mass , the post - recombination matter ( dark matter and baryons ) power spectrum exhibits a growing mode that scales with the cosmic expansion parameter , @xmath28 , as @xmath29 here , @xmath30 is a transfer function that encapsulates evolution before recombination at @xmath31 , @xmath32 is the density perturbation growth factor from linear theory , and @xmath33 is the controlling parameter set of the background cosmological model . dark energy models that involve modifications to general relativity may introduce @xmath34-dependence into the growth function above . the power spectrum in the reference @xmath8cdm  model is set by present - epoch energy densities , @xmath35 , where @xmath36 is the dimensionless hubble constant and @xmath37 is the density of component @xmath38 relative to the critical density , @xmath39 . the curvature density , @xmath40 , is zero to within @xmath41 @xcite , consistent with a flat spatial metric on cosmic scales . our notation uses ` b ' for baryons , ` c ' for cold , dark matter ( cdm ) , and ` m ' for all matter : @xmath42 . in the minimal model , the dark energy is a vacuum energy with equation of state , @xmath43 . we use @xmath44 for this case and employ @xmath45 when referring to models wherein the dark energy equation of state , @xmath46 , differs from @xmath47 . [ tab : wmapcosm ]    cc + parameter & value +   + @xmath44 & @xmath48 + @xmath49 & @xmath50 + @xmath51 & @xmath52 + @xmath53 & @xmath54 + @xmath25 & @xmath55 + @xmath56 & @xmath57 + @xmath58 & @xmath59 + @xmath60 from @xcite . current constraints for a flat @xmath8cdm  model from cmb measurements , combined with angular clustering of red galaxies and local measurements of @xmath61 , are shown in table  [ tab : wmapcosm ] @xcite . the parameter @xmath62 is the variance in density fluctuations evaluated at horizon crossing , which is independent of @xmath34 for @xmath63 , and the wavenumber @xmath64 corresponds to a large comoving length scale , @xmath65 .    in the minimal model , the matter density fluctuations filtered within a sphere of comoving radius @xmath66 are gaussian distributed with zero mean . the comoving radius defines a mass , @xmath67 , of matter within that radius in the young universe , when @xmath68 . early observations that the variance in galaxy counts is near unity on a scale of @xmath69 led to this as a conventional choice of scale at which to quote the fluctuation amplitude ( see table  [ tab : wmapcosm ] ) . the corresponding mass , @xmath70 , is characteristic of rich clusters of galaxies . the variance of linearly - evolved , cdm fluctuations , filtered on mass scale @xmath71 , has the form @xmath72 where the filter function is @xmath73 $ ] for the typical case of sharp ( or top - hat ) spatial filtering within radius @xmath66 . evaluating equation  [ eq : sigmasq ] at @xmath74 and @xmath75 produces the oft - quoted matter power spectrum normalization parameter , @xmath58 . we will see below that @xmath76 serves as a similarity variable for expressing model - independent forms of the halo space density and clustering . the evolution of the fluctuation spectrum , equation  [ eq : powspect ] , is valid at early times or at scales sufficiently large so that @xmath77 at all times . on small scales , where cdm power spectra are generically maximum , fluctuation growth produces @xmath78 , and linear theory breaks down . mode - mode coupling terms become important to the dynamics , and solutions in fourier space become difficult . while higher - order perturbation theory solutions can extend analytic evolution to later times than linear theory ( e.g. * ? ? ? * ; * ? ? ? * ) , the full problem is typically treated using n - body simulations , discussed below . a recent analytical advance considers lss as an effective fluid . @xcite show that integrating out small - scale , non - linear structures renormalizes the cosmological background and introduces dissipative terms , of order @xmath79 , into the dynamics of large - scale modes , with @xmath80 the typical velocity dispersion of collapsed halos . since even the most massive halos have @xmath81 , the magnitude of these effects is very small . furthermore , @xcite show that virialized halos decouple completely from large - scale dynamics , at all orders in the post - newtonian expansion . astrophysical structures , from the first stars at high redshift to galaxy clusters at low redshift , tend to emerge from local maxima of the filtered density field . while density peaks are generally non - spherical @xcite , a first - order description considers them spherical and isolated from their surroundings . birkhoff s theorem then implies that the expansion histories of radial mass shells within a peak follow trajectories perturbed from the overall background , with sufficiently dense shells expanding to a maximum size and then contracting . the traditional _ ansatz _ assumes collapse by a radial factor of two @xcite , after which a quasi - virialized and quasi - hydrostatic structure  a perfectly spherical _ halo _  is born . the collapse criterion is that the _ linearly - evolved _ perturbation amplitude reach a critical value , @xmath82 , with @xmath83 the conventional choice . applying this idea to the cdm spectrum , equation  [ eq : sigmasq ] , leads to a characteristic mass scale , @xmath84 , defined by @xmath85 = { { \ensuremath{\delta_{\mathrm{c}}}}}$ ] . at a given epoch , a spectrum of halo masses exist , with masses above ( below ) @xmath84 forming from perturbations with amplitudes above ( below ) the _ rms _ level of the filtered gaussian spectrum . considerable literature ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * and many others ) has established this picture as the _ halo model _ of large - scale structure . we review here only aspects relevant for cluster cosmology ; a more thorough review can be found in @xcite . the basic element of the halo model is the population mean space density , @xmath86 , in units of number per unit comoving volume , commonly referred to as the _ mass function_. expressed as a differential function of mass , it takes the form @xmath87 where @xmath88 is the comoving mean matter density and @xmath89 is a model - dependent function of the filtered perturbation spectrum , equation  [ eq : sigmasq ] . analytic forms for @xmath89 capture much , but not all , of the behavior seen in n - body simulations , as discussed below . the spatial clustering of halos is described by a modified version of the matter power spectrum . on large spatial scales , or low wavenumbers , the halo autocorrelation power spectrum is modified , @xmath90 where @xmath91 , the halo _ bias function _ , is independent of @xmath34 , for the case of gaussian fluctuations , but dependent on mass and epoch . while this expression applies to the spatial autocorrelation of systems with fixed mass @xmath71 , it generalizes to the cross - correlation between sets of halos at different masses , @xmath92 . the theory of peaks in gaussian random fields expresses the bias as a function of the normalized peak height , @xmath93 @xcite . below , we show that n - body simulations support the forms of equations  ( [ eq : massftnbasic ] ) and ( [ eq : corrftnbasic ] ) , but a precise fit to the mass function requires that @xmath89 be adjusted to include explicit redshift dependence , @xmath94 . there are subtleties to the definition of mass in simulations that must also be taken into account . various astrophysical processes play out within the photon - baryon components of the evolving cosmic web , including hydrodynamic , magnetohydrodynamic and radiative transfer effects ; star and black hole formation with associated feedback of momentum , energy and entropy ; and so on . except for the immediate vicinity of black holes , these processes involve classical physics that is largely known . but the fully three - dimensional and non - linear nature of the problem , the wide dynamic range in length and time scales , and the non - trivial couplings among the constituent physical processes introduce tremendous complexity into baryon evolution . galaxy formation is truly a grand challenge computational problem . we touch on select issues relevant to the observable features of galaxy clusters . _ shocks and turbulent mhd heating . _ during halo formation , gravitational potential energy in the baryonic component is thermalized via shocks . the highest mach numbers , of tens or more , should occur in the accretion shocks at the edges of clusters ( e.g. * ? ? ? * ) . while these strong shocks are expected to be efficient particle accelerators , recent observations place tight limits on the volume - averaged pressure contributions from relativistic particles @xcite shocks with mach numbers of a few are also associated with major mergers : a spectacular example is the narrow radio relic in the cluster ciza j2242.8 + 5301 , for which @xcite use multi - frequency radio and polarization observations to infer a mach number @xmath95 in a shock located @xmath96 from the cluster center . most of the energy thermalized during cluster formation , however , is dissipated in weak shocks that are persistently driven by dissipating sub - structures and ongoing minor mergers . shocks are also driven by jets from agns and , at earlier times , by winds from star forming galaxies . details of the nano - parsec scale physics that drive thermalization remain under active study , especially the roles of magnetic fields , turbulence and plasma instabilities ( e.g. * ? ? ? observations and simulations discussed below indicate that thermalization is efficient ; thermal pressure supplies the bulk of support against gravity within the halo potential except during brief periods near periapsis of major mergers . _ radiative cooling . _ since the intracluster plasma ( and , to a lesser extent , its interstellar counterpart ) is optically thin at most wavelengths , radiation loss is the primary cooling mechanism for the baryonic component of halos . indeed , the classic criterion for setting an upper bound on galaxy size comes from balancing the gas cooling time against the halo dynamical time @xcite . the first generation of stars form at @xmath97 , aided by molecular hydrogen line emission , within halos of mass @xmath98 @xcite . by @xmath99 , atomic line cooling in halos with virial temperatures above @xmath100  k produces the first generation of galaxies , which grow hierarchically for a time determined by the large - scale environment . proto - cluster regions have more efficient cooling at high redshift than do proto - voids , but the feedback from vigorous , early production of compact sources helps to quench star formation before a large fraction of baryons are converted to stars . the cooling timescale of the gas in massive halos is typically longer than a hubble time , except for a subset of systems that exhibit cool cores . the central @xmath101 region of such systems tends to be x - ray bright and typically contains a dominant elliptical galaxy . we discuss aspects of cool core phenomenology in section  [ sec : future : coreevol ] .    _ star and black hole formation . _ cold , molecular gas fuels star formation . the star formation rate can roughly be considered as proportional to the local rate of gas cooling below @xmath100  k , but there are other considerations . different venues for star formation exist , ranging from quiescent disks to the bulges of tidally - triggered starburst galaxies , and it is not yet clear whether a single model based on local gas conditions captures the full range of observed behavior . supermassive black hole ( smbh ) growth occurs through mergers and accretion in galactic cores , and these central engines drive quasar and radio jet activity ( e.g. * ? ? ? sloan digital sky survey ( sdss ) quasar studies indicate that smbhs of mass @xmath102 exist at @xmath103 @xcite , and processes for forming such large black holes in the first few hundred million years of the universe have been proposed @xcite .    in gaussian random fields , small - scale peaks are more abundant when embedded within large - scale peaks , so the largest galaxies and quasars at high redshift represent the progenitors of massive galaxies observed in low redshift clusters .    _ feedback from compact sources . _ feedback of mass , momentum , and entropy from stellar / smbh sources is important at all stages of the lss hierarchy . photoionization and supernova - driven winds serve to limit cooling and star formation in low - mass halos @xcite . jets driven by accretion onto the central smbh appear to be required to limit the maximum size of galaxies ( e.g. * ? ? ? * ; * ? ? ? formulations for this feedback typically tie the energy input to the mass accretion rate which , in turn , is governed by the local rate of cooling and/or cold accretion . the end result of this competition between cooling and heating is that heating largely wins . the overall efficiency of star formation is small , and peaks in halos of roughly galactic scale ( e.g. * ? ? ? figure  [ fig : dai2010fstarvcirc ] shows a recent compilation of stellar mass fraction ( @xmath104 ) measurements as a function of halo circular velocity , @xmath105 , with @xmath71 the total halo mass and @xmath106 its radius @xcite . the horizontal lines show the cosmic baryon fraction , @xmath107 , derived from wilkinson microwave anisotropy probe ( _ wmap _ ) data analysis @xcite . the stellar mass fraction is maximized at a few tens of percent of the cosmic mean in halos with @xmath108 , equivalent to a mass of @xmath109 at @xmath110 . in cluster - sized halos , the stellar fraction declines with mass , taking on values @xmath111 of the global baryon fraction at the highest masses . yet , the largest galaxies are found in the cores of massive clusters , and their very old stellar populations produce a characteristically narrow ` red sequence ' in a color - magnitude diagram of cluster members .    _ dynamical and thermodynamical equilibrium . _ in the context of the evolving cosmic web , the processes above must contend with conditions imposed by halo merging . at any given time , major mergers , such as those involving progenitor pair mass ratios larger than @xmath112 , occur in @xmath111 of the population , concentrated toward the highest masses . these rare events can drive the mass contents of a halo considerably out of equilibrium . minor mergers , while much more frequent , are also less damaging . current simulations and observations indicate that the dynamical and thermodynamical response of halos is quite fast . hydrostatic and virial equilibrium assumptions are typically valid to within roughly ten percent for the majority of the cluster population ( e.g. * ? ? ? * ; * ? ? ? all this astrophysical evolution offers a treasure trove of observational possibilities . uniquely in massive clusters , all of the matter is readily observed , allowing a complete census to be taken . stars make up 13% @xcite ; @xmath113 resides in the hot , diffuse , intergalactic gas @xcite ; and the rest is in the form of non - baryonic cdm ( section  [ sec : physics_dm ] ) .      as tracers of massive halos , galaxy clusters provide a number of signatures that are sensitive to the underlying cosmology . we review here the principles underlying key methods . a typical set of cosmological parameters for such studies might consist of the primordial spectrum amplitude and slope , the present - epoch densities of the three energy components dominant at late times , the dimensionless hubble constant , and the de equation of state parameters , @xmath114 where the last two parameters define a linearly - evolving de equation of state , @xmath115 this particular set is meant to be illustrative . there is considerable variation in the literature , and many works restrict analysis to a flat cosmology , which removes one degree of freedom from the above through the condition @xmath116 .      the yield of upcoming cluster surveys will be sufficiently large to enable disaggregation by angular position , redshift , and the observed signal , @xmath117 . ( note the latter is also referred to in the literature as the _ mass proxy _ , or sometimes the _ observable mass _ , @xmath118 . ) complications associated with the signal  mass likelihood and with redshift estimation are discussed below . as a starting point , consider a perfect tracer of mass , @xmath119 , with error - free redshifts , @xmath120 . within a given survey , the expected number of halos , @xmath121 , in a cell described by mass bin @xmath28 and redshift bin @xmath15 with solid angle @xmath122 is @xmath123 cosmology enters this expression through the mass function and the volume element , @xmath124 . the counts in each large spatial bin will deviate from the mean by an excess number , @xmath125 , determined by the local large - scale density field , @xmath27 . following @xcite , the spatial covariance of the counts is @xmath126 where @xmath127 describes the spatial correlation between mass - redshift bins , @xmath128 here , @xmath129 is the window function for cell @xmath15 ( that , when present , can include the effects of redshift estimate uncertainties ) and @xmath130 is a geometric term that depends on the comoving separation , @xmath131 , between cells @xmath15 and @xmath132 . when cells @xmath15 and @xmath132 sample different redshifts , an accurate approximation uses their geometric mean to evaluate @xmath133 @xcite .    combining the spatial clustering with a diagonal shot noise term forms the full covariance for a survey sample . derivatives of the mean counts and covariance with respect to model parameters form the fisher information matrix used in survey forecasts . expressions for the fisher matrix can be found in @xcite . equations  ( [ eq : mbar ] ) through ( [ eq : xipixel ] ) serve as the foundation of likelihood analysis of large cluster surveys . to be useful in practice , these expressions must undergo a number of modifications , including : transformation from mass to the signal used for cluster detection , @xmath134 ; inclusion of counting errors arising from incompleteness ( missed sources ) and impurities ( false sources ) ; and inclusion of photometric uncertainties , @xmath135 . we discuss these issues in section  [ sec : theorytopractice ] and summarize current results in section  [ sec : constraints : growth ] . the mass fraction of hot gas , @xmath136 , measured within a characteristic radius of a halo at redshift @xmath137 can be written as @xmath138 where @xmath139 accounts for star formation and other baryon effects within that radius . at large radii in the most massive halos , where the hot icm dominates the baryon budget and the impacts of feedback processes are modest , baryon losses are small and @xmath140 is a reasonable expectation .    motivated by the growing body of measurements of @xmath136 from the _ rosat _ x - ray satellite , @xcite and @xcite recognized that a mismatch in the dependence on metric distance , @xmath141 , between gas mass ( @xmath142 ) and total mass ( @xmath143 ) measured from x - ray observations implied that gas fraction measurements in massive clusters could be exploited as a distance estimator , with @xmath144 . like type ia supernovae , massive clusters serve as standard calibration sources that test the expansion history of the universe . key benefits , relative to survey counts , are the ability to perform this test with a relatively small number of clusters and the relative insensitivity to cluster selection . we summarize results from this exercise in section  [ sec : constraints : fgas ] .      in a similar vein , @xcite noted that x - ray and sz measurements could be combined to determine distances to clusters . the cmb spectral shift is governed by the compton y - parameter , a measure of the electron pressure along the line of sight , @xmath145 . given an observed sz signal , @xmath146 , and a predicted signal based on x - ray measurements of the icm density and temperature , @xmath147 , the angular diameter distance scales as @xmath148 the cosmological constraint originates from the distance dependence of the x - ray measurements , @xmath149 , and the requirement that @xmath150 . accurate sz and x - ray flux and temperature calibration are particularly important to this method , referred to below as xsz .      the thermal and kinetic sz signals from clusters ( section  [ sec : obs_sz ] ) cause distortions in the cmb at small angular scales ( @xmath151 ) . if the distortion pattern from a single halo of mass @xmath71 at redshift @xmath137 is described by an angular fourier transform , @xmath152 , then the full halo population will generate a fluctuation spectrum @xcite @xmath153 adding halo spatial correlations gives a small correction to this estimate @xcite . this approach to testing cosmology is limited by degeneracy with astrophysical assumptions , as the interplay between @xmath152 and @xmath154 makes clear . measurements of the cosmic peculiar velocity field contain additional cosmological information ( e.g. @xcite and references therein ) . the kinetic sz effect ( section  [ sec : obs_sz ] ) in principle offers a way to measure the peculiar velocities of galaxy clusters . although some initial results based on such measurements have been reported ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , the technique has not yet reached the maturity of those discussed above and is not discussed further in this review . n - body simulations of a single , collisionless dark matter fluid offer the means to investigate non - linear evolution of lss under an implicit ` light - traces - mass ' assumption . the technology supporting such simulations has advanced to the state where @xmath155 is available ( but not yet realized ) on peta - scale computational platforms @xcite . employing larger-@xmath156 simply to model bigger volumes is a natural mode of growth , since parallelization is relatively simple ( large - volume domain decomposition minimizes the particle transfer among computational nodes ) , the number of timesteps is independent of @xmath156 , and the light - traces - mass assumption is easier to justify under modest mass and force resolution . large - volume simulations produce generous halo population realizations with which to calibrate the mass function and clustering of halos , and current state - of - the - art studies employ ensembles of @xmath157-particle simulations . coupled n - body and gas dynamic simulation methods enable multi - fluid studies that break free of the light - traces - mass assumption . indeed , the first application of this class of codes tested the possible separation of baryons and neutrinos within clusters formed in a universe dominated by massive neutrinos @xcite . the field has advanced considerably since then , and we refer the reader to @xcite for a recent review . we discuss primarily dark matter simulations here , with some relevant multi - fluid simulations results presented in the next section .      through the mass function , halo mass provides the critical measure that connects observables to the underlying cosmology . but halos are complex , dynamic structures that confound attempts at a unique definition of mass .    in the model of spherical collapse applied to initial density peaks , the halo edge and interior mass are readily defined by the outermost caustic in dark matter or by the location of the shock in cold baryonic accretion @xcite . in both cases , this radius marks an abrupt transition in the mean radial velocity , separating a nearly hydrostatic interior from an infall - dominated exterior . halos forming in 3-d simulations deviate from this ideal case in important ways , some of which can be described by higher - order analytic approaches to peak evolution @xcite . the collapse process is more ellipsoidal than spherical , and merging competes with smooth accretion as the dominant mode of halo growth ( e.g. , * ? ? ? defining centers and boundaries in this complex environment has become a matter of convention . two common algorithmic conventions have emerged : ( i ) percolation , also known as friends - of - friends ( fof ) , and ( ii ) spherical overdensity ( so ) . fof first links all pairs of particles within a given distance , @xmath158 , then merges them into groups based on a shared link condition ( ` a friend of a friend is a friend ' ) . the so approach first filters the particle field to identify peaks , then grows spheres around peaks with sizes determined by an interior density threshold , @xmath159 . the threshold density @xmath160 is typically chosen to be either the background matter density , @xmath161 , or the critical density , @xmath162 . unless otherwise specified , we adopt the latter convention in this article . = 4.75 cm = 4.75 cm = 4.75 cm    several studies discuss the relative merits of these approaches and argue values for the parameters @xmath158 and @xmath163 ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? figure  [ fig : lukic09 ] provides a visualization of three simulated halos spanning a range of dynamical and morphological behaviors . in each panel , white particles are members of the fof halo with @xmath164 while green are so members using @xmath165 against @xmath166 . these are typical of parameter values used in the literature . the left panel shows a relatively isolated system where the two methods give fairly consistent results . the other panels show two discrepant cases ; in the middle is a highly - structured , active merger while , at right , percolation across a filamentary bridge links two similarly sized systems that are just beginning to merge .    the discrepant cases do not dominate in number , but neither are they uncommon . for cosmological studies , what is important is to establish an accurate accounting process to enumerate observable halo features . roughly speaking , observers viewing the systems in figure  [ fig : lukic09 ] would be likely to identify one dominant cluster in the left and middle panels , and two in the right . an fof accounting system would need to admit a non - unitary condition ( one halo maps to two clusters ) when converting mass to observable signals . in contrast , so masses map to integrated aperture observations more directly . for this reason , so masses see more frequent use for survey data analysis . the original multiplicity function paper of @xcite used the clustering of particles in n - body experiments with @xmath167 to support their analytic form for @xmath89 in equation  [ eq : massftnbasic ] . later , @xcite used @xmath168 simulations to set free parameters of their @xmath89 model derived using an ellipsoidal , rather than spherical , collapse approximation . using a suite of simulations of open and flat cosmologies with @xmath169 ranging from @xmath112 to 1 , @xcite found a unique , three - parameter form for @xmath89 that produced a mass function accurate to @xmath170 across the suite of models . a recent study by @xcite employs 22 large ( @xmath171 ) simulations produced with three independent n - body codes to calibrate a functional form motivated by @xcite , @xmath172 \ , e^{-c/\sigma^2 }   . \label{eq : massftntinker}\ ] ] this study was the first to open the density threshold degree of freedom ; their fitting parameters are published as functions of @xmath163 ( against @xmath173 ) for @xmath174 $ ] . with the high statistical power of their simulation ensemble , @xcite achieve a fit with @xmath175 statistical precision in halo number at @xmath176 for a @xmath8cdm  cosmology . maintaining this precision for redshifts @xmath177 requires the introduction of mild redshift dependence into the fit parameters , @xmath178 , @xmath179 and @xmath180 . the theoretically expected halo counts above masses @xmath181 and @xmath10 in the reference @xmath8cdm  cosmology , shown in figure  [ fig : surveyhistory ] , are based on the tinker form for threshold @xmath165 against the mean mass density ( see fitting formulae in * ? ? ? on the other hand , the bias function measured in the same simulation ensemble shows no need for such redshift - dependent corrections . framed in terms of the normalized linear perturbation amplitude , @xmath182 , @xcite find a robust fit of the form @xmath183 with a single set of parameters @xmath184 that are written only as functions of @xmath163 . for the case @xmath185 ( i.e. , @xmath186 peaks ) , the value of the bias is large , @xmath187 , for @xmath165 . the cluster power spectrum , equation  [ eq : corrftnbasic ] , can be enhanced by factors of several tens over the mass power spectrum . the very massive end of the fof mass function was recently revised by @xcite using @xmath188-particle simulations in @xmath8cdm  cubic volumes up to @xmath189 in scale . above @xmath190 , their fit lies up to @xmath191 above prior calibrations @xcite . gravitational relaxation drives the phase - space structure of halos to a common structure that applies from small galactic satellites to the most massive galaxy clusters . the form of the radial density , @xmath192 is known as the navarro - frenk - white ( nfw ) profile @xcite . here , @xmath193 is the scale radius , @xmath194 is the concentration parameter ( with @xmath195 ) and @xmath196}$ ] . simulations show that concentration and mass are weakly correlated . in the mass range of galaxies to clusters , @xmath197 , with @xmath198 at @xmath110 and @xmath199 at @xmath200 ( e.g. , * ? ? ? that study finds that a fixed concentration , @xmath201 , applies in the mean to high mass halos , independent of redshift . tracking the mass accretion histories of halos in simulations , @xcite find a common functional form , and show that the formation epoch correlates strongly with concentration . the concentration  mass relation can be understood as a result of adiabatic contraction of differently - shaped peaks in the linear density field @xcite . cluster cosmology originates from phenomena observed on the sky , in the 2 + 1 space of angular coordinates and redshift . the observables employed for a likelihood analysis must be predicted under a set of combined cosmological and astrophysical parameters , @xmath202 . for constraints based on cluster counts , the mass function , @xmath86 , written in terms of spherical overdensity or percolation measures from simulations needs to be translated into a signal function , @xmath203 , for one or more signals , @xmath204 . we use the terms _ signal _ and _ observable _ interchangeably , and generically they refer to bulk measures at mm ( sz decrement @xmath205 ) , optical ( richness , @xmath206 , or velocity dispersion , @xmath207 ) , or x - ray ( luminosity , @xmath208 ; temperature , @xmath209 ; gas mass , @xmath210 ; and/or gas thermal energy , @xmath211 ) wavelengths ( see section  [ sec : observations ] ) . an ideal experiment would measure all of these observables within apertures optimally matched to the underlying halo sizes , @xmath212 . this ideal is often frustrated by signal - to - noise constraints and confused by projection effects and foreground / background contamination .      scaling relations for cluster signals , based on assumptions of virial equilibrium and self - similar internal structure , were first published by @xcite . in this model , halos at fixed mass and redshift are identical , and scalings with mass and redshift follow calculated power - law behaviors . observations generally support power - law behavior , but not always with the self - similar slope ( section  [ sec : constraints : scaling ] ) . we describe here a non - self - similar model that incorporates arbitrary mass scaling and allows for variations at fixed mass and redshift .    for compactness of notation , let @xmath213 , for each of the @xmath156 observables , @xmath204 , and let @xmath214 . the power - law assumption transforms to log - linear scaling @xmath215 where the average is over a very large cosmic volume . the elements of @xmath216 are the slopes of the individual mass - observable relations , and the intercepts @xmath217 reflect the evolution at fixed mass . at a fixed epoch , we can always choose units such that @xmath218 ( as we do below ) . for cosmological studies , a measure of merit is the equivalent mass scatter in each signal , @xmath219 . various processes , including different formation histories and the stochastic nature of mergers , generate deviations from the mean . taking these as gaussian in the log leads to a form for the _ conditional signal likelihood _ , @xmath220^\dagger \mathbf{\psi}^{-1 } \left [ { \mathbf{s}}-{\mathbf{\bar{s}}}(\mu , z)\right ] \right\ } . \label{eq : pofs}\ ] ] the elements of the covariance matrix , @xmath221 , could have mass or redshift dependence , but a first - order approach considers them as constants .    when the mass variance of signals is small , @xmath222 , then the above expressions can be convolved with a locally power - law approximation to the mass function , @xmath223 , to obtain the local _ signal space density _ function , @xmath224 , \label{eq : nofs}\ ] ] where @xmath225 is the variance about the mean log - mass selected by the set of signals @xmath226 , @xmath227 the first term above is the mean mass for the case of a flat mass function , @xmath228 . the second term , represents the ( eddington ) mass bias induced by asymmetry in the mass function convolution . upscattering of low - mass systems dominates when @xmath229 , and the high - mass end of the @xmath8cdm   mass function is steep , @xmath230 @xcite . these equations make explicit the degeneracy between cosmology ( e.g.  @xmath231 and @xmath28 ) and astrophysics ( e.g.  @xmath216 and @xmath232 ) inherent in cluster counts . they provide the means to compute biases , relative to a mass complete sample , associated with signal - limited cluster samples ( discussed further in section  [ sec : theory : selection ] ) . figure  [ fig : stanek10 ] provides support for this model from millennium gas simulation analysis @xcite . the left panel shows deviations about the mean behavior of four intrinsic ( 3-dimensional ) properties measured within @xmath233 for @xmath234 halos with mass @xmath235 at @xmath110 . the lower diagonal and red histograms show results from a cooling and preheating ( ph ) treatment of the baryons , where the entropy is instantaneously raised to @xmath236 at @xmath237 . only a small fraction of baryons cool into stars in this model @xcite . the upper diagonal and blue histograms are from a gravity - only ( go ) treatment , where the gas is heated only by shocks and does not cool . the internal properties generally have modest variance , and pairs tend to be positively correlated with typical correlation coefficient @xmath238 . halos identified by a pair of properties will have mass variance @xmath239 shown by the areas of the off - diagonal circles in the right panel of figure  [ fig : stanek10 ] . individual properties lie along the diagonal . the intrinsic gas thermal energy , @xmath205 , selects mass with @xmath240 dispersion , the best individual measure for both physics cases . this level is also seen in the simulations of @xcite which include cooling , star formation and feedback . pairs of intrinsic measurements always improve mass selection , and the strong correlation between @xmath241 and @xmath205 combines with the large mass variance of @xmath241 to achieve mass selection with @xmath242 scatter in the ph model . applying bayes theorem to this model allows one to write the likelihood of mass and an observable , @xmath243 , for a sample selected on observable @xmath244 . when the two signals are correlated , one can show that the scaling with mass of the non - selection signal will be @xmath245   , \label{eq : s2bargs1}\ ] ] which is biased relative to the naive expectation of @xmath246 . the intrinsic correlation between signals at fixed mass is relatively challenging to constrain from current data , but first measurements have been made for samples selected using optical @xcite and x - ray @xcite observations . clusters on the sky relate to halos through selection on one or more observables . matching cluster detections ( which originally reside in a 2 + 1 space of angular position and signal - to - noise ) to halos can sometimes be complex ; two halos along nearly the same line of sight may be blended into a single cluster , or a single halo may be fragmented into more than one cluster . the frequency of these occurrences is typically not large , @xmath247 , but the exact values are sensitive to a number of factors , particularly detection method and mass , and so are best modeled via direct sky realizations ( e.g. * ? ? ? * ) . the selection observable can be distorted from its intrinsic value ( equation  [ eq : pofs ] ) by triaxiality , by additional sources along the line - of - sight , by mis - centering and/or mis - estimation of the radial scale , and by other effects . telescope / instrument calibration and data processing methods also contribute to the error budget . for upcoming studies using cluster counts , photometric redshift errors have an important , but not dominant , effect ( section  [ sec : new_surveys ] ) . testing cosmology with halo counts and clustering requires that the theoretical mass function be transformed , via the scaling relations and a model of the selection process , into a prediction for the distribution of clusters in the space of survey observables ( e.g. redshift and x - ray flux ) . the scaling relation parameters set the space density portion of the survey yield ( equation  [ eq : nofs ] ) in terms of the ( cosmologically dependent ) local amplitude , @xmath248 , and logarithmic slope , @xmath249 , of the mass function . sample selection must be well understood to avoid perturbing @xmath248 and @xmath249 from their true values , biasing cosmological results . fortunately , such effects can be mitigated by survey self - calibration @xcite or by calibration using follow - up observations , as discussed below . the task of empirically constraining the scaling relations is complicated by the fact that the clusters targeted for follow - up observations are themselves subject to selection effects related to their original discovery . in an x - ray flux - limited sample , for example , higher x - ray luminosity at a given mass leads to a larger probability of detection ( commonly known as malmquist bias ) . the effects of selection bias must therefore be accounted for in the calibration of scaling relations , much as in the cosmological analysis ( e.g. @xcite ) . figure  [ fig : constraints : selection ] illustrates the influence of selection on observed scaling relation data in a cartoon case . the full population ( black crosses and green points ) obeys a scaling law ( red line ) with non - trivial intrinsic scatter . in the simple case where detection requires a particular threshold luminosity ( the dashed , blue line ) , it can be seen that , even if every detected cluster is followed up to obtain precise measurements of the mass and luminosity , the resulting data set will be a biased representation of the full population . while complete at the highest masses , the sample is increasingly incomplete at low masses , with the low - luminosity systems absent . a closely related consideration is the effect of the underlying mass function on the observed scaling relation data . the distribution of the relation s independent variable(s ) ( in this case cluster masses ) within the full population generically influences constraints on scaling laws ( e.g. @xcite ) . neglecting to account for this influence corresponds to the assumption of uniformly distributed independent variables ; often this approximation is sufficient , but the exponentially steep slope of the cluster mass function suggests that we should take the issue seriously in the context of cluster cosmology @xcite . figure  [ fig : constraints : selection ] illustrates how the steepness of the mass function influences the fraction of the observed data which are strongly biased relative to the underlying scaling relation . given the need to solve for both the slope and scatter of the scaling relation , accounting for the disparity in the number of high - mass and low - mass systems is critical . we note that simply conditioning the sampling distribution on cluster detection , as some authors have done , is not sufficient to rigorously recover all the scaling information . note that this effect has a floor set by non - zero intrinsic scatter in the scaling relations , but the effect can in principle be enhanced by measurement error . however , measurement errors in current x - ray and optical cluster surveys are typically smaller than the intrinsic dispersion , even at the survey limit . thus , re - measurement of the survey observables through deeper , follow - up observations ( e.g. to improve the signal - to - noise of x - ray or sz flux ) does not circumvent the issue of selection bias in the scaling relation analysis . while selection bias clearly influences scaling relations involving the selection observable , it also influences relations of other signals with which the selection observable has non - zero intrinsic correlation ( equation  [ eq : s2bargs1 ] ) . this is illustrated in figure  [ fig : constraints : selection2 ] , for a signal which is correlated with the selection observable with coefficient 0.1 ( left panel ) and 0.9 ( right panel ) . the red line shows the true scaling law and the points shown correspond to the detected clusters from figure  [ fig : constraints : selection ] ( right panel ) . with relatively mild intrinsic correlation , as has been found for temperatures and soft x - ray flux detection @xcite , the distribution of data points closely follows the underlying relation ; for more extreme values of the correlation coefficient , as might be expected , e.g. , between temperature and sz signal , deviations due to selection bias become evident . note that the severity of the effect also depends on the covariance of the signals rather than only on the correlation coefficient ( i.e. the size of the marginal scatter in each signal is also important ) . cluster samples are often characterized in terms of _ completeness _ and _ purity _ @xcite . completeness is used in many ways , but its simplest form for cluster cosmology refers to the fraction of halos above mass @xmath71 at redshift @xmath137 that are identified in a survey with some observable limit , @xmath250 . completeness of unity is achievable at high masses when the survey limit , @xmath250 , lies sufficiently far in the signal likelihood s negative tail . impurity is a measure of false positive sources in the sample . fewer conventions for its definition exist in the literature . generically , one can write the observed counts above some signal limit @xmath117 as a sum , @xmath251 , where the first term represents genuine cluster systems  manifestations of a single massive halo along the line of sight  and the second expresses detections of other origin . zero impurity , @xmath252 , is a desired goal . telescopes aimed at a distant halo necessarily collect photons that originate elsewhere along the multi - gigaparsec sightline than within the target system . due to their softer angular profiles , sz , lensing and optical cluster signals can be blended more readily than x - ray . chance orientations of two or more halos within local supercluster regions create an asymmetric tail to high signal values . considered in terms of mass selection , the effect produces a tail to low masses in the distribution of halo mass selected at a given signal ( e.g. * ? ? ? * ) .    since the matter components of halos are generally ellipsoidal rather than spherical , orientation variations also produce scatter in signals observed in halos of fixed mass . signals are generally maximized when viewed along the long axis and minimized along the short axis . orientation can affect cluster selection , with prolate systems oriented along the line - of - sight being preferentially included . since its collisional nature drives the x - ray emitting gas toward equipotential surfaces , it tends to be rounder than the dark matter and so less susceptible to orientation bias .    as discussed in section  [ sec : observations ] , the density squared dependence of the x - ray emissivity means that x - ray selection is less prone to projected confusion . optical richness measurements roughly trace mass density and are therefore more easily confused by projection and orientation effects . sz measurements are intermediate , since the sz effect depends on electron pressure , the product of density and temperature . it is important to keep in mind that theory offers many potential deviations from the reference @xmath8cdm  cosmology sketched above . key model assumptions  that the dark matter is a weakly interacting massive particle , that inflation produced a gaussian spectrum of initial density fluctuations with a power - law initial spectrum , that small - amplitude metric perturbations are well described by newtonian , weak field expansions in general relativity , and so on  need to be rigorously tested . in section  [ sec : physics ] , we discuss ways in which clusters can be used to test a number of proposed modifications to the reference model . in this section we review briefly the physics underlying multiwavelength observations of galaxy clusters . we summarize efforts to construct cluster catalogs , with an emphasis on surveys that have led to cosmological constraints . we discuss techniques used to measure the masses of clusters , and observable proxies that correlate tightly with mass . most of the baryons in the universe are in diffuse gas . typically , this gas is very difficult to observe . within galaxy clusters , however , gravity squeezes the gas , heating it to virial temperatures of @xmath253@xmath254 , which causes it to shine brightly in x - rays . galaxy clusters therefore ` light up ' at x - ray wavelengths as luminous , continuous , spatially - extended sources ( figure  [ fig : multiwavelength ] ) . = 4.5 cm = 4.5 cm = 4.5 cm    the primary x - ray emission mechanisms from the diffuse icm are collisional : free - free emission ( bremsstrahlung ) ; free - bound ( recombination ) emission ; and bound - bound emission ( mostly line radiation ) . the emissivities of these processes are proportional to the square of the electron density , which ranges from @xmath255 in the centers of bright ` cool core ' clusters to @xmath256 in cluster outskirts . at these low densities , the x - ray emitting plasma is optically thin and in the coronal limit , which makes modeling straightforward . for survey observations , the primary x - ray observables are flux , spectral hardness and spatial extent . using deeper , follow - up observations of individual clusters , modern x - ray satellites allow the spatially - resolved spectra of clusters to be determined precisely , permitting measurements of the density , temperature and metallicity profiles of the icm , and a host of derived thermodynamic quantities . for reviews of the principles underlying x - ray observations of clusters see , e.g. , @xcite and @xcite . the optical and near - ir emission from galaxy clusters is predominantly starlight . the galaxy populations of clusters are dominated by ellipticals and lenticulars ( i.e. early - type galaxies ) . this is particularly true in the central regions , where the largest and most luminous galaxies are found ( figure  [ fig : multiwavelength ] ) . the old and relatively homogeneous nature of their stellar populations leads to the majority of the galaxies in clusters occupying relatively tight loci in color - magnitude diagrams ( e.g. * ? ? ? this characteristic has proved important to modern cluster finding algorithms .    for optical surveys of clusters , the main observables are the richness ( i.e. the number of galaxies within the detection aperture ) , luminosity and color . for follow - up observations of individual clusters , aimed in particular at measuring their masses , the primary observables are the galaxy number density , luminosity , and velocity dispersion profiles . typical velocity dispersions for large clusters are of order 1000@xmath257 .    for reviews of optical studies of galaxy clusters including discussions of the development of the field , see @xcite and @xcite . as cmb photons pass through a galaxy cluster they have a non - negligible chance to inverse compton scatter off the hot icm electrons . this scattering boosts the photon energy and gives rise to a small but significant frequency - dependent shift in the cmb spectrum observed through the cluster known as the thermal sunyaev - zeldovich ( hereafter sz or tsz ) effect @xcite . the magnitude of the effect is proportional to the line of sight integral of the product of the gas density and temperature . the kinetic sz ( ksz ) effect is an additional , smaller distortion of the cmb spectrum due to the peculiar motion of a cluster with respect to the hubble flow ( i.e. the cmb rest frame ) . the magnitude of the ksz effect is proportional to the peculiar velocity . for a review see @xcite . according to general relativity , the gravity associated with a mass concentration will bend light rays passing near to it in a phenomenon known as gravitational lensing . this can both magnify and distort the images of background galaxies . with modern data , gravitational lensing can be detected clearly in the statistical appearance of background galaxies observed through clusters ( weak lensing ) , and in the field ( often termed cosmic shear ) . occasionally , lensing can also lead to strong distortions and multiple images of individual sources ( strong lensing ) . for a galaxy cluster and background galaxies of known redshifts , the measured gravitational shear can be used to infer the cluster mass . for a recent review of gravitational lensing , see @xcite . a well - designed cluster survey should meet requirements in terms of angular scale , flux sensitivity and redshift coverage . the survey should be as complete ( i.e. not have missed clusters that it should have detected ) and pure ( i.e. not have detected spurious clusters ) as possible and the selection function describing the completeness and purity as a function of signal , position and redshift should be known precisely . the survey observables should correlate as tightly as possible with mass . in tension with these requirements , surveys must also be constructed within the context of limited resources . x - ray observations currently offer the most mature and powerful technique for constructing cluster catalogs . the primary advantages of x - ray surveys are their exquisite purity and completeness , and the tight correlations between x - ray observables and mass . galaxy clusters are simple to identify at x - ray wavelengths , being the only x - ray luminous , continuous , spatially extended , extragalactic x - ray sources . clusters have typical soft x - ray band luminosities of @xmath258 erg / s or more , and spatial extents of several arcmin or larger , even at high redshifts . given modest angular resolution , e.g. @xmath259  arcmin ( full width half maximum , fwhm ; easily achievable ) and tens of detected counts , the x - ray emission from galaxy clusters can be detected against a background populated otherwise only sparsely with point - like active galactic nuclei . the first x - ray cluster catalogs constructed for cosmological work ( ( * ? ? ? * often called the  brightest 50  or b50 catalog ) ; @xcite ) were based on the ariel v and heao-1 all - sky surveys , and pointed observations made with the _ einstein observatory _ and _ exosat _ ( see also @xcite ) . these catalogs provided early evidence for evolution in the x - ray luminosity function of clusters @xcite and were used subsequently in a series of pioneering cosmological works ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? these catalogs were eventually superseded by surveys carried out with the _ rosat _ satellite . this mission , launched in june 1990 , had two main parts : the _ rosat _ all - sky survey ( rass ; * ? ? ? * ) , spanning the first 6 months ; and pointed observations , which took place over the next 8 years . the main instrument aboard _ rosat _ , the position sensitive proportional counter , had a modest point spread function ( psf ; @xmath260  arcmin fwhm in survey mode ) , but low background and a wide field of view ( @xmath261 degree diameter ) . the main cluster catalogs constructed from the rass and used in cosmological studies include the rosat brightest cluster sample ( bcs ; * ? ? ? * ) , which covered the northern hemisphere at high galactic latitudes and low redshifts ( @xmath262 ) to a flux limit of @xmath263 @xmath264(0.12.4 kev ) ; the rosat - eso flux - limited x - ray galaxy cluster survey ( reflex ; * ? ? ? * ) , which covered the southern sky at low redshifts to a flux limit of @xmath265 @xmath264   in the same band ; the hiflugcs sample @xcite of the x - ray brightest clusters at high galactic latitudes , with @xmath266 @xmath264(0.12.4kev ) ; and the massive cluster survey ( macs ; * ? ? ? * ) , which extended this work to higher redshifts @xmath267 and slightly fainter fluxes @xmath268 @xmath264 . other cluster surveys have been constructed from the rass , or are in the process of being constructed , but have not yet been used to derive rigorous cosmological constraints .    a number of x - ray cluster catalogs have also been constructed based on serendipitous discoveries in the pointed phase of the rosat mission . notable among these are the rosat deep cluster survey and the 400 square degree rosat pspc galaxy cluster survey ( 400d ; * ? ? ? * ) , which have been used to derive cosmological constraints . these catalogs cover much smaller areas than the rass , but reach an order of magnitude or more fainter in flux ( figure  [ fig : xray_surveys ] ) .    a second major advantage of x - ray surveys is the observed strong correlation between x - ray luminosity and mass across the entire flux and redshift range of interest . these quantities follow a simple power law relation ( section  [ sec : constraints : scaling ] ) , with a dispersion in luminosity at a given mass of @xmath269 and no significant outliers @xcite . the density - squared dependence also makes the x - ray survey signal from clusters relatively insensitive to projection effects . thus , an x - ray survey of sufficient depth can be translated straightforwardly into statistical knowledge of the distribution of massive halos .    in principle , x - ray surveys could be constructed using even lower - scatter mass proxies ( section  [ sec : obs_masses_proxies ] ) such as temperature or center - excised luminosity as the survey observable . however , given the ease and depth to which total x - ray luminosity can be measured , these lower - scatter mass proxies are typically used as auxiliary data ( e.g. * ? ? ? * ; * ? ? ? the primary disadvantage of x - ray cluster surveys is that they can only be carried out from space , which makes their construction relatively expensive . the first extensive cluster catalog was constructed at optical wavelengths by george abell @xcite based on visual inspection of photographic plates from the palomar observatory sky survey . abell identified clusters as concentrations of 50 or more galaxies in a magnitude range m3 to m3 + 2 ( where m3 is the magnitude of the third brightest cluster member ) and radius @xmath270 mpc ( with distance estimated based on the magnitude of the tenth brightest galaxy ) . clusters were further characterized into richness and distance classes . abell s catalog was updated and extended to the southern sky by @xcite ( hereafter aco ) . the final aco sample has more than 4000 clusters . an additional , early optical cluster catalog extending to poorer systems was compiled by zwicky and collaborators ( see e.g. * ? ? ? * ) , although the search criteria were less strict than abell s . @xcite applied a percolation algorithm to an early cfa redshift catalog to identify a set of 92 nearby groups and clusters . using 4-m class telescopes and a mix of photographic plate and ccd observations , @xcite opened high redshift cluster studies by identifying 418 systems over @xmath271 deg@xmath272 extending out to @xmath273 . spatial and photometric matched filter methods ( e.g. * ? ? ? * ) as well as the introduction of n - body simulations to calibrate projection effects ( e.g. * ? ? ? * ) marked the beginning of the modern era of optical cluster cosmology . because the cores of galaxy clusters are dominated by red , early - type galaxies , an effective way to reduce the impact of projection effects is to use color information to select for overdensities of red galaxies ( e.g. * ? ? ? * and references therein ) . the red - sequence cluster survey ( rcs ) , a sample of 956 clusters identified with a single ( @xmath274 ) color , provided the first modern cosmological constraints using optical selection @xcite .    to cover a broad range of redshifts , multi - color photometry is needed to track the intrinsic 4000 angstrom break feature of old stellar populations as it reddens . the five - band photometry of the sloan digital sky survey ( sdss ) enabled such selection . the maxbcg catalog @xcite of 13,823 clusters with optical richness @xmath275 was produced using @xmath276 colors and spans the redshift range @xmath277 . cosmological constraints from this sample @xcite are discussed below . recently , larger sdss clusters samples have become available , identified using photo - z clustering @xcite , a gaussian mixture modeling extension of the maxbcg method @xcite , and an adaptive matched filtering approach @xcite . these catalogs contain between 40000 and 69000 clusters spanning @xmath278 , and cover roughly 8000 deg@xmath272 of sky . a primary challenge to cosmological analysis using such catalogs is the definition of robust mass proxies that possess minimal and well - understood scatter across the full mass and redshift ranges of interest . projection of filamentary structures and small groups along the line of sight has a greater impact on optical cluster catalogs than x - ray , and these effects introduce a degree of skewness into the mass - observable relations @xcite . uncertainty in modeling this and other selection effects currently limits the constraining power offered by the large sample sizes of optical cluster catalogs . the first large catalogs of galaxy clusters selected from observations of the sz effect are currently under construction , using measurements made with the south pole telescope ( spt ; * ? ? ? * ; * ? ? ? * ) , the atacama cosmology telescope ( act ; * ? ? ? * ; * ? ? ? * ) and the _ planck _ satellite @xcite . the primary advantage of sz surveys is that , in contrast to x - ray and optical measurements , the sz signal of a cluster does not undergo surface brightness dimming . sz surveys are therefore well - suited , in principle , to searches for massive clusters at high redshifts . the surveys mentioned above are each expected to produce catalogs of hundreds of massive systems at intermediate - to - high redshifts . challenges for these projects include determining the optimal observables ( i.e. the best mass proxies ) to measure from the survey data in the current , low signal - to - noise ratio regime ; calibrating the mass scaling for these observables ; and understanding in detail the impact of contamination by radio and infrared sources @xcite . projection effects are also expected to be more significant for sz surveys than for x - rays @xcite . accurate measurements of cluster masses provide a cornerstone of cosmological work . x - ray mass measurements are based on the assumption of hydrostatic equilibrium ( hse ) in the icm . for a spherically symmetric system in hse , the measured gas density and temperature profiles can be related to the total mass ( e.g. * ? ? ? * ) , @xmath279,\ ] ] where @xmath280 is the mass within radius @xmath106 , @xmath281 is the icm temperature , @xmath282 is the gas particle density , @xmath283 is newton s constant , @xmath34 is the boltzmann constant , and @xmath284 is the mean molecular weight . note that the mass within radius @xmath106 depends more strongly on the temperature than the density at that radius . hydrostatic equilibrium requires that the gravitational potential remain stationary on a sound crossing time ; that all motions in the gas be subsonic ; and that forces other than gas pressure and gravity are unimportant . the hydrostatic method can therefore not be applied robustly to systems undergoing major merger events , nor to regions of otherwise relaxed clusters where these assumptions break down , e.g. in their central regions where strong agn feedback effects are commonly observed @xcite .    out to intermediate radii , measurements of the gas temperature and density profiles with _ chandra _ or _ xmm - newton _ are straightforward . at large radii ( @xmath285 ) , however , where the x - ray emission is faint , such measurements become challenging . recent advances in this regard have been made with the _ suzaku _ satellite , and opportunities for additional progress remain ( section  [ sec : future_outskirts ] ) . potentially increased levels of non - thermal pressure support ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) and gas clumping @xcite can also complicate measurements at large radii . a number of approaches have been used in implementing the hydrostatic method . the most common , which employs relatively strong priors , uses parameterized fits to the observed , projected surface brightness and temperature profiles ; these are then used to calculate the appropriate partial derivatives at each radius to determine the mass profile ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? a second , arguably preferable , approach employs a non - parametric deprojection of the brightness and temperature data , but assumes that the mass distribution follows a well - motivated parameterized form ( e.g. equation  [ eq : nfw ] ; @xcite , @xcite ) ; this approach simultaneously provides a framework for testing the validity of various mass models . in the case of very high quality x - ray data , a fully non - parametric deprojection of the surface brightness and temperature data can be employed , without additional , regularizing assumptions ( e.g. * ? ? ? x - ray mass measurements are relatively insensitive to triaxiality @xcite . for dynamically relaxed clusters , and for measurements out to intermediate radii , simulations indicate that hydrostatic x - ray masses should exhibit modest scatter ( @xmath286% ) and be biased low by @xmath287% ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , due primarily to kinetic pressure arising from residual gas motions .      like the x - ray method , optical - dynamical mass measurements are based on the assumption of dynamical equilibrium , with the galaxies used as test particles in the cluster . the mass enclosed within radius @xmath106 is given by the jeans equation ( e.g. * ? ? ? * ; * ? ? ? * )    @xmath288 , \label{eq : jeans}\ ] ]    where @xmath289 is the galaxy number density , @xmath290 the 3-dimensional velocity dispersion , and @xmath291 the velocity anisotropy parameter . these quantities can be determined under model assumptions from the projected galaxy number density and velocity dispersion profiles . an advantage of the optical dynamical method over the x - ray method is that it is insensitive to several forms of non - thermal pressure support that affect x - ray mass measurements ( e.g. magnetic fields , turbulence , and cosmic ray pressure ) . the galaxy population can also be observed at high contrast out to large radii . however , where the x - ray gas is a collisional fluid that returns rapidly to equilibrium following a disruption , the galaxies are collisionless and relax on a longer timescale @xcite . whereas x - ray mass measurements are relatively insensitive to triaxiality , the galaxy velocity anisotropy must be accounted for . the precision of optical dynamical measurements is also limited by the finite number of galaxies . while identifying the center of a cluster is straightforward at x - ray wavelengths , at optical wavelengths this can be a source of significant uncertainty . the infall regions of clusters form a characteristic trumpet - shaped pattern in radius - redshift phase - space diagrams , the edges of which are termed caustics . the identification of these caustics enables mass measurements out to larger radii ( up to @xmath292 mpc ) , to an accuracy of @xmath293% ( e.g. * ? ? ? in contrast to x - ray and optical dynamical methods , gravitational lensing offers a way to measure the masses of clusters that is free of assumptions regarding the dynamical state of the gravitating matter @xcite . weak lensing methods have an important role in cosmological work : while triaxiality is expected to introduce scatter in individual ( deprojected ) mass measurements at the level of tens of per cent @xcite , for statistical samples of clusters and using suitable mass estimators , working over optimized radial ranges and with good knowledge of the redshift distribution of the background population , weak lensing measurements are expected to provide almost unbiased results on the mean mass @xcite . the most common technique employed in weak lensing mass measurements is fitting the observed , azimuthally - averaged gravitational shear profile with a simple parameterized mass model ( e.g. * ? ? ? stacking analyses of clusters detected in survey fields have also proved successful in calibrating the mean mass  observable scaling relations down to relatively low masses ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? strong lensing enables precise measurements of the projected masses through regions enclosed by gravitational arcs . in combination with weak lensing , strong lensing constraints can improve significantly the absolute calibration of projected mass maps ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? . deprojected strong lensing measurements are particularly sensitive to triaxiality ( e.g. * ? ? ? * ; * ? ? ? nonetheless , detailed strong lensing studies of small samples of highly relaxed clusters have reported mass measurements in good agreement with x - ray results ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? some recent works have used clusters with multiple systems of strongly - lensed arcs to constrain the geometry of the universe ( * ? ? ? * and references therein ) . this approach is related to the shear ratio test also discussed for weak lensing measurements with clusters @xcite . a good mass proxy should be straightforward to measure and correlate tightly with mass , exhibiting minimal dispersion across the mass and redshift range of interest . robust , low - scatter mass proxy information for just a small fraction ( appropriately selected ) of the clusters in a survey can boost significantly its constraining power with respect to self - calibration alone ( e.g. * ? ? ? * ; * ? ? ? the total x - ray luminosity has an intrinsic dispersion at fixed mass of @xmath269 @xcite . the scatter in optical richness at fixed mass for the maxbcg catalog is also @xmath294 , with a modest inferred non - gaussian tail toward low masses @xcite . the scatter in the projected , integrated sz flux at fixed mass is expected to be somewhat smaller ( @xmath295% ; @xcite ) , although this is yet to be measured robustly from data . ( note that the scatter in the observed , projected sz flux is significantly larger than for the predicted , intrinsic signal ; figure  [ fig : stanek10 ] . ) since total x - ray luminosity , optical richness and integrated sz flux are all survey observables , their scatter versus mass tends to play a significant role in tests based on cluster counts even when more precise mass measurements for a fraction of the clusters are available from follow - up data . although not the basis for cluster surveys , other observables that correlate tightly with mass can provide powerful mass proxies . for the most massive clusters , the x - ray emitting gas mass is strongly correlated with total mass , with an observed scatter @xmath296 at fixed mass @xcite . the x - ray temperature , and the product of gas mass and temperature ( @xmath297 ) , have observed scatters @xmath298 ( e.g. @xcite ) ; it is thought that the tightness of the @xmath299@xmath71 relation may extend to lower masses than for either the @xmath209@xmath71 or @xmath210@xmath71 relations @xcite . center - excised x - ray luminosity also traces mass extremely well , with an observed scatter @xmath296 ( @xcite ; see section  [ sec : future : coreevol ] ) . weak lensing mass measurements exhibit larger scatter ( tens of per cent ) . however , the minimal predicted bias in the mean mass for statistical samples of lensing measurements ( section  [ sec : obs_masses_lensing ] ) makes the combination with x - ray measurements particularly promising .    while the relation between the three - dimensional dark matter velocity dispersion and mass is predicted to be tight @xcite , velocity anisotropy and projection effects cause the one - dimensional velocity dispersion to exhibit larger scatter ( @xmath300% , implying a scatter in mass at fixed velocity dispersion of @xmath269 ; @xcite ) . the relationship between galaxy and dark matter velocity dispersion , the velocity bias , also remains uncertain . the past decade has seen marked improvements in the data and analysis techniques employed in cluster cosmology . _ chandra _ observations of relaxed clusters provided important measurements of the distance scale , confirming the recent acceleration of the cosmic expansion . improvements in mass measurements led to a convergence in estimates of @xmath58 from cluster counts , and x - ray surveys extending to @xmath301 provided the first constraints on dark energy from the growth of structure . we begin this section with a review of the latest results from measurements of the number and growth of clusters , as well as the closely related problem of constraining scaling relations . we then review the state of other cluster - based probes of cosmology . we conclude by summarizing the latest constraints on dark energy from the combination of galaxy cluster data with other , independent probes . the application of cluster measurements to other areas of fundamental physics are discussed in section  [ sec : physics ] .        because the cosmological model influences predictions for scaling relation observations , and vice versa ( section  [ sec : theory : selection ] ) , the two must be constrained simultaneously . the statistical framework needed for this approach , in which a subset of detected clusters are targeted for more detailed observations including mass estimation , is described by @xcite . without reproducing the details here , we enumerate the components of the statistical model :    1 .   the mass function and expansion history , which together predict the number of clusters as a function of mass and redshift , @xmath302 . 2 .   stochastic scaling relations which describe the ( multivariate ) distribution of observables as a function of mass and redshift . measurements for each cluster and associated sampling distributions , describing the probability of obtaining particular measured values given true cluster properties ( i.e. a model for the measurement errors ) , accounting for any covariances in measured quantities . we note that all quantities entering the scaling relations need not be measured for every cluster , apart from ( necessarily ) those determining cluster selection ( e.g. , mass follow - up need not be complete ) . 4 .   the selection function , quantifying the probability of a cluster being detected and included in the sample .    using the likelihood function corresponding to this formalism ( see @xcite ) , constraints on cosmological and scaling relation parameters are obtained by marginalizing over the true , unknowable properties of each cluster , constrained by the measured values . straightforward adaptations of this approach , for example incorporating a cluster bias model ( section  [ sec : theory : countsclustering ] ) or stacked gravitational lensing data for poor clusters and groups , are possible . we concentrate our review on three recent , independent results based on measurements of the local abundance and evolution of clusters from x - ray and optical surveys , as well as preliminary results from new sz surveys . @xcite studied the bcs @xcite , reflex @xcite and bright macs @xcite cluster samples compiled from the rass . this survey strategy , covering a large fraction of the sky to relatively shallow depth , is optimized to the task of finding the largest clusters at the expense of depth in redshift . the mantz et  al . data set consists of 238 clusters with masses @xmath303 distributed over the redshift range @xmath304 . pointed _ chandra _ or _ rosat _ follow - up observations were used to measure @xmath210 for 94 clusters , while the remaining 144 had measurements only of the cluster redshift and survey flux . the @xmath305@xmath306 relation was marginalized over using hydrostatic mass estimates for relaxed clusters from @xcite . the data set was used to simultaneously constrain cosmological parameters and the @xmath208@xmath71 and @xmath209@xmath71 relations using the method described above and detailed in @xcite . @xcite also used x - ray selected clusters , but pursued a different survey strategy . their data set consists of disjoint low - redshift and high - redshift samples , with 49 clusters ( originally culled from several rass samples ) at redshifts @xmath307 ( nearly all @xmath308 ) and an additional 36 clusters serendipitously detected in the 400 square degree survey @xcite at @xmath309 . this survey strategy , covering a smaller area of the sky to greater depth , naturally finds fewer of the most massive clusters , but extends to lower masses ( @xmath310 ) and higher redshifts . all 85 clusters in the full data set were followed up with _ chandra _ , and their masses were estimated using either @xmath210 or @xmath299 as a proxy ; the proxy  mass relations were calibrated using hydrostatic mass estimates for a sample of well observed , low - redshift clusters . the analysis includes empirical constraints on the scaling relations and corrections for selection effects , though not using exactly the approach described in section  [ sec : constraints : procedure ] . the top panels of figure  [ fig : constraints : mfcn_fgas ] provide a simplified visualization of how constraints on dark energy arise from these data , comparing the observed mass function in two redshift ranges to model predictions for different cosmological models . in particular , an open universe with no dark energy clearly under - predicts the evolution of the mass function over the redshift range of the data . = 5.63 cm = 5.6 cm    the optically selected maxbcg sample @xcite employed by @xcite probes a different part of the cluster population ; it is restricted to lower redshifts than the x - ray samples described above ( @xmath277 ) , but extends to lower masses ( @xmath311 ) . this lower effective mass limit , which changes less strongly with redshift compared to x - ray surveys , makes the maxbcg sample significantly larger than the others , with @xmath312 clusters divided into 9 bins based on optical richness . mean masses for 5 richness ranges were estimated through a weak gravitational lensing analysis of stacked clusters , providing information from which to constrain the richness  mass relation . the cosmological analysis accounts for the covariance between cluster counts in each richness bin and the mean lensing mass estimates . the results obtained by these three groups on flat @xmath8cdm and constant @xmath313 models are summarized in table  [ tab : constraints : basic ] . note that , for the two works which fit @xmath313 models , the results on @xmath169 and @xmath58 are dominated by the low - redshift data and so are not degraded noticeably by the introduction of @xmath313 as a free parameter ; thus all three sets of constraints are directly comparable . the agreement between the different works , as well as others listed in table  [ tab : constraints : basic ] , is encouraging ; in particular , the close agreement in the constraints on @xmath58 reflects the relatively recent convergence in cluster mass estimates using different techniques , and our improved understanding of the relevant systematics ( section  [ sec : obs_masses ] ; see also , e.g. , @xcite ) . importantly , the concordance @xmath8cdm model provides an acceptable fit to the data in each case . [ tab : constraints : basic ]    ccccccc + reference@xmath314 & & @xmath58 & @xmath315 & @xmath316 & @xmath313 & @xmath53 +   +   + m10 & x - ray & @xmath317 & @xmath318 & @xmath319 & @xmath320 & + v09 & x - ray & @xmath321 & @xmath322 & @xmath319 & @xmath323 & +   + r10 & optical & @xmath324 & @xmath325 & @xmath319 & @xmath47 & + h09 & x - ray & @xmath326 & @xmath112 & @xmath319 & @xmath47 & +   + s03 & x - ray & @xmath327 & @xmath328 & @xmath319 & @xmath47 & +   + a08 & x - ray & & @xmath329 & @xmath330 & @xmath47 & + a08 & x - ray & & @xmath331 & @xmath319 & @xmath332 & + e09 & x - ray & & @xmath333 & @xmath319 & @xmath334 & + l06 & x - ray+sz & & @xmath335 & @xmath319 & @xmath47 & +   + b06 & x - ray+sz & & @xmath112 & @xmath319 & @xmath47 & @xmath336 + s04 & x - ray+sz & & @xmath112 & @xmath319 & @xmath47 & @xmath337 +     + @xmath60 entries @xmath338 indicate the assumption of global spatial flatness ( or use of _ wmap _ priors requiring flatness at the few percent level ) ; other entries without error bars indicate parameters that were fixed in the corresponding analysis . @xmath339 error bars are marginalized ( single - parameter ) 68.3% confidence intervals , and include each author s estimate of the systematic uncertainties ( with the exception of s04 ) . @xmath314 a08  =  @xcite ; b06  =  @xcite ; e09  =  @xcite ; h09  =  @xcite ; l06  =  @xcite ; m10  =  @xcite ; r10  =  @xcite ; s03  =  @xcite ; s04  =  @xcite ; v09  =  @xcite . @xmath340 cluster surveys extending to redshifts @xmath341 are required to constrain @xmath313 from the evolution of the mass function . note that the @xmath315 and @xmath58 constraints from these works are essentially unchanged , whether or not @xmath313 is allowed to differ from @xmath47 . figure  [ fig : constraints : rozo_wm_s8 ] ( left panel ) shows the joint constraints on @xmath169 and @xmath58 obtained by ( * ? ? ? * solid lines ) , which display the typical degeneracy between those parameters from cluster survey data . ( the degeneracy can be broken , for example , by including cluster @xmath342 data ; see section  [ sec : constraints : fgas ] . ) also shown are results from 5 years of _ wmap _ data ( * ? ? ? * dashed lines ) , which are tight for the assumed flat @xmath8cdm model . nevertheless , it is evident that the combination of the two types of data ( shaded regions ) is significantly improved over either one individually . = 5.5 cm    the constraints on constant @xmath313 models from @xcite and @xcite are respectively shown in figures  [ fig : constraints : rozo_wm_s8 ] ( right panel ) and [ fig : constraints : mantz_w ] , along with results from various other cosmological data sets . the cluster results are in good agreement with one another , as well as with the other , independent data . the @xmath343 precision constraint on @xmath313 from cluster growth alone ( including systematic uncertainties ) is clearly competitive , and further constraining power may be available from theoretical advances ( see section  [ sec : physics : lcdm ] ) . the strong degeneracy in cmb constraints evident in figure  [ fig : constraints : mantz_w ] also arises in many other models that are more complex than flat @xmath8cdm ; the ability of cluster data to break this degeneracy by providing precise and independent constraints on @xmath58 ( right panel of figure  [ fig : constraints : mantz_w ] ) makes the combination particularly powerful ( e.g. section  [ sec : physics : nu ] ) . although their respective surveys are not yet complete , preliminary results have been reported from sz - detected clusters from the south pole telescope ( spt ; * ? ? ? * ) and atacama cosmology telescope ( act ; * ? ? ? as the number of cluster detections used in these works is small , respectively 21 and 9 at redshifts @xmath344 , the data are not capable of producing interesting constraints on their own , but both groups demonstrate consistency with results from the _ wmap _ satellite @xcite . neither group has yet obtained simultaneous constraints on cosmology and the relevant sz scaling relation , making their final error budgets strongly dependent on the priors chosen to constrain the scaling relation .      to obtain scaling relations appropriate for constraining cosmology , the analysis of the two must be simultaneous and must account properly for the influence of the survey selection function ( e.g. @xcite and @xcite ) and the cluster mass function , as in ( * ? ? ? * see also sections  [ sec : constraints : procedure ] and section  [ sec : theory : selection ] ) . some authors have included corrections for the expected malmquist bias given a flux limit or selection function @xcite , but most analyses of scaling relations in the literature employs least - squares regression without detailed consideration of sample selection . since various cluster samples , with different selection functions , have been used , it is not surprising that results on the slopes , scatters and evolution of the scaling relations have varied widely compared with the formal uncertainties . another statistical issue is that the covariance that comes about when the scaling quantities are measured from the same observations , for example temperature and hydrostatic masses from x - ray data , is commonly ignored ( see section  [ sec : modeling ] ) . the method used to estimate masses has also varied . most authors have used the assumption of hydrostatic equilibrium applied to x - ray data , regardless of the dynamical state of the cluster , which must introduce spurious scatter due to departures from equilibrium and non - thermal support ( e.g. @xcite ) . more recently , mass proxies such as gas mass ( e.g. @xcite ) and x - ray thermal energy ( @xmath299 ; e.g. @xcite ) , or gravitational lensing signal ( e.g. @xcite ) have been employed . these variations in mass estimation and analysis methods , in addition to changes to instrument calibration over the years , makes a comprehensive and fair census of scaling relation results problematic . here we focus on the cosmological importance of scaling relation measurements , citing examples from recent work where the issues mentioned above are at least partially mitigated .    for x - ray and sz observables , under the assumption of strict self - similarity ( no additional heating or cooling ) , @xcite derived specific slopes and redshift dependences for the power - law form of equation  [ eq : sbarmu ] : @xmath345^{4/3 } , \nonumber \\    { \ensuremath{kt_{\mathrm{mw } } } } & \propto & \left [ e(z)m \right]^{2/3 } , \nonumber \\ e(z)\,y & \propto & \left [ e(z)m \right]^{5/3},\end{aligned}\ ] ] where the factors of @xmath346 are appropriate for measurements made at a fixed critical - overdensity radius . the subscripts ` bol ' and ` mw ' reflect the fact that these predictions apply to the bolometric luminosity and mass - weighted temperature . optical richness is more complex to predict , but empirical studies that map galaxies to sub - halos in simulations support a power law richness  mass relation for groups and clusters @xcite . figures  [ fig : constraints : scaling ] and [ fig : constraints : other_scaling ] show a few examples of recent scaling relation measurements . @xcite present a @xmath71@xmath208 relation for x - ray selected clusters in the cosmos field by measuring stacked weak lensing masses ( left panel of figure  [ fig : constraints : scaling ] ) . under the common assumption of symmetric scatter in the log , stacking on @xmath208 allows the mean @xmath347 relation to be recovered , at the cost of losing information about the intrinsic scatter . as the authors note , transforming these results to an @xmath348 relation introduces a dependence on the mass function . stacked lensing was also used to determine the mass  richness relation of optically selected clusters in sdss by @xcite . their results are shown in the right panel of figure  [ fig : constraints : scaling ] ; red points in the figure show that compatible masses were derived from galaxy velocity dispersion measurements . = 5.8 cm = 5.8 cm    for more massive systems , lensing can provide mass estimates on a cluster - by - cluster basis . the left panel of figure  [ fig : constraints : other_scaling ] shows such a mass  temperature relation from @xcite for x - ray selected clusters . since temperature only weakly influences x - ray detectability , selection bias should be relatively unimportant here , although intrinsic correlation with luminosity can still produces subtle biases relative to a mass - limited sample ( section  [ sec : theory : selection ] ) . = 5.9 cm    an sz scaling relation from @xcite , using clusters detected by spt , appears in the right panel of figure  [ fig : constraints : other_scaling ] ; in this case , masses are estimated from measurements of the x - ray thermal energy , @xmath299 , and an assumed @xmath299@xmath71 relation . this case is instructive in that the plotted @xmath205 values are derived from the survey data , with suggestions of malmquist bias in the flattening observed at the lowest masses ( compare with figure  [ fig : constraints : selection ] ; this potential bias is accounted for in the fit in figure  [ fig : constraints : other_scaling ] ) . recently , the _ planck _ collaboration released an early sz - selected cluster catalog together with an analysis of sz scalings for existing x - ray and optical cluster samples ( * ? ? ? * and references therein ) . measurements of the scaled sz signal from _ planck _ as a function of x - ray luminosity for @xmath349 x - ray selected galaxy clusters are shown in the left panel of figure  [ fig : constraints : planck_scaling ] , along with the stacked , mean sz signal in luminosity bins and model expectations derived from the observed x - ray properties @xcite . there is good agreement between the model expectations and the _ planck _ observations , validating the simple model description of the hot icm over the mass and redshift range probed . = 7.5 cm = 7.5 cm    the right panel of figure  [ fig : constraints : planck_scaling ] shows the result of a similar exercise , stacking the _ planck _ sz signal as a function of richness for the optically selected maxbcg catalog . also shown is the expected signal , based on the observed mean mass  richness relation for the maxbcg sample @xcite and a standard x - ray @xmath205@xmath71 scaling relation @xcite . in this case , the observed _ planck _ sz signal lies well below model predictions . for an x - ray detected sub - sample of the maxbcg catalog , however , relatively good agreement between the _ planck _ sz measurements and model predictions is observed @xcite . an important difference between the @xmath208 and richness comparisons is that the centers and radii used to extract the clusters sz signals from the ( relatively low resolution ) _ planck _ data were respectively based on x - ray and optical estimates . besides optical center and scale mis - estimates , other potential sources for the discrepancy in figure  [ fig : constraints : planck_scaling]b include lower than expected purity in the optical sample ( with lower mass halos being enhanced in richness by projection ) ; unmodeled biases in the scaling relations used ( e.g. the impact of non - thermal pressure support on the x - ray relations ) ; and the presence of an intrinsic anti - correlation between the galactic and hot gas mass fractions in the clusters . these and other possibilities are under investigation .    in general , most current scaling results are consistent with power - law relations over a wide range in mass ( e.g. @xcite ) . however , the exponents often differ from the values in equation  [ eq : constraints : virial_scaling ] , due at least in part to the impact of astrophysical processes , mainly star / smbh formation and associated feedback , discussed in section  [ sec : theory_astrophysics ] . for example , the @xmath208@xmath71 slope is commonly found to exceed @xmath350 , with the value depending on the mass range explored . there is less consensus on the @xmath209@xmath71 slope , with most estimates consistent with self - similarity ( e.g. @xcite ; @xmath351 ) , but some finding a shallower slope ( e.g. @xcite ; @xmath352 ) ; systematics related to correlated observables , mass estimation ( see section  [ sec : modeling ] ) , instrument calibration , selection , and differences in the mass range studied may play a role . there are fewer empirical estimates of the @xmath205@xmath71 or @xmath299@xmath71 slopes , but most are broadly compatible with the predicted value . we note that departure from self - similar scaling with mass does not necessarily imply departure from self - similar evolution with redshift ; for example , the high - mass halos in the preheated simulations of @xcite have x - ray luminosities that evolve within @xmath353 of the self - similar expectation at @xmath354 , despite having a slope with mass that is @xmath355 steeper than self - similar . the marginal intrinsic scatter in each relation determines the degree of selection bias in surveys and each observable s usefulness as a mass proxy ; current results are reviewed in section  [ sec : obs_masses_proxies ] . the evolution of the scaling relations is of great importance , since it is potentially degenerate with the cosmological signal of cluster growth . because the impact of survey limits can vary strongly with redshift for x - ray surveys , it is particularly crucial to account for them in this context . both @xcite and @xcite investigate evolution in the @xmath208@xmath71 relation , finding only marginal @xmath356 departures from the self - similar redshift dependence in equation  [ eq : constraints : virial_scaling ] . @xcite additionally found no evidence for departures from self - similar evolution in the @xmath209@xmath71 or center - excised @xmath208@xmath71 relations ( see section  [ sec : future : coreevol ] ) . as current and upcoming surveys expand the reach of cluster samples to @xmath357 , obtaining precise constraints on scaling relation evolution will be imperative ( see also section  [ sec : modeling ] ) .    the effects of survey bias on scaling relations can be _ partially _ mitigated by selecting clusters based on an observable other than the observable of immediate interest ( e.g. @xcite ) . however , as sections  [ sec : theory : mvscaling ] and [ sec : theory : selection ] make clear , the amount of residual bias depends on the intrinsic correlation between the two observables at fixed mass . to date , the only empirical estimates of such correlation are for mass and x - ray luminosity at fixed optical richness , @xmath358 @xcite ; and for soft - band x - ray luminosity and temperature at fixed mass , @xmath359 @xcite . large , overlapping surveys , and/or surveys coupled with multi - wavelength follow - up campaigns , are required to better constrain the property covariance of clusters ; ultimately , such constraints will provide a new level of robustness to cosmological work .      in principle , the confirmed existence of even a single galaxy cluster of implausibly high mass would challenge the standard @xmath8cdm model with gaussian initial conditions . recently , the discovery of high redshift , massive systems such as xmmuj2235.3 - 2557 at @xmath360 @xcite and spt - clj0546 - 5345 at @xmath361 @xcite have led to reports of possible tension with gaussian @xmath8cdm @xcite . such a test is in some sense an attractively simple alternative to the more involved work discussed above . however , despite the fact that the test involves only one or a few very massive clusters detected in a survey rather than the complete sample , a robust assessment of the likelihood still requires a detailed understanding of the selection function and survey biases , as well as a full accounting for the effects of scatter in the mass  observable scaling relations ( sections  [ sec : theory : selection ] and [ sec : constraints : procedure ] ) . the accuracy and precision of the cluster mass measurements are also critical , due to the steepness of the high mass tail of the cluster mass function . errors in mass measurements at the tens of per cent level , for example ( as might be expected for weak lensing measurements of an individual cluster ) , can modify the likelihood of such a cluster being observed by up to an order of magnitude . @xcite estimate confidence limits for the exclusion of the gaussian @xmath8cdm model based on the properties of the most massive galaxy cluster , or @xmath156 most massive galaxy clusters , detected in a given survey , employing constraints on the expansion history from current data . these authors conclude that none of the presently known high mass , high redshift clusters are in significant tension with the standard gaussian @xmath8cdm paradigm . the framework described in section  [ sec : constraints : procedure ] for simultaneously constraining cosmology and scaling relations has yet to be applied to the spatial clustering of clusters . however , @xcite have obtained cosmological constraints using a method based on the karhunen - love eigenvectors of 428 clusters from the x - ray selected reflex sample above a luminosity of @xmath362 ( see @xcite for the theory underlying this method of estimating the power spectrum ) . the constraints are relatively weak compared with more recent results from other methods : @xmath363 and @xmath364 at 95.4% confidence , fixing @xmath365 and without including the systematic uncertainty due to cosmic variance . in @xcite , the analysis was extended to include cluster abundance using an empirical x - ray luminosity  mass scaling relation from @xcite , demonstrating that the combination of the two methods breaks the degeneracy between @xmath169 and @xmath58 , in particular reducing the uncertainty on @xmath58 to the point where systematics related to mass estimation dominate ( table  [ tab : constraints : basic ] ) .    in principle , the spatial distribution of clusters can be employed in a simpler way , as has been done with individual galaxies , by using the baryon acoustic oscillation signature in the power spectrum as a probe of cosmic distance . @xcite and @xcite respectively analyzed the correlation function and power spectrum of clusters in the optically selected maxbcg catalog , finding a weak ( @xmath366 ) detection of the bao peak . more recently , @xcite found no significant evidence for a bao feature in the smaller , x - ray selected reflex  ii catalog . cluster gas mass fractions can be measured robustly using x - ray , or the combination of x - ray and sz , data for dynamically relaxed clusters ( section  [ sec : observations ] ) . when measured from x - ray data , @xmath342 values within a given angular aperture depend on cosmology as @xmath367 , while the predicted @xmath342 for a given cosmology is given by equation  [ eq : constraints : fgas_model ] . ( a more detailed expression for the comparison of measured and predicted @xmath342 values is given by @xcite , who include terms accounting for instrument calibration , non - thermal pressure , and the relationship between the characteristic radius of the model , e.g. @xmath368 , and the aperture of the measurement . ) the apparent redshift dependence of @xmath342 measurements on the cosmological background is illustrated in the bottom panels of figure  [ fig : constraints : mfcn_fgas ] .    in principle , @xmath342 measurements can be made at radii corresponding to any overdensity . in practice , this overdensity should be low enough ( i.e. the radius large enough ) that non - gravitational feedback effects do not introduce prohibitive scatter ; but not so small ( i.e. the radius not so large ) that the measurements become dominated by the systematic limitations of the instruments . a variety of simulations indicate that radii @xmath369 are sufficiently large to benefit from low @xmath342 scatter ( e.g. @xcite ) , while at radii @xmath370 uncertainties in the x - ray background and the impact of gas clumping can become a concern ( @xcite ; see also section  [ sec : future_outskirts ] ) .    using this method , ( * ? ? ? * see also @xcite ) obtained cosmological constraints using @xmath342 measurements at @xmath368 from _ chandra _ observations of 42 hot ( @xmath371 ) , relaxed clusters covering the redshift range @xmath372 . their analysis incorporated weak priors on @xmath53 and @xmath373 from hubble key project @xcite and big bang nucleosynthesis data , priors on the baryonic depletion and stellar mass fraction of clusters and their evolution ( the term @xmath139 in equation  [ eq : constraints : fgas_model ] ) , and marginalized over systematic allowances accounting for instrument calibration and non - thermal pressure . from equation  [ eq : constraints : fgas_model ] , we see that the normalization of the @xmath374 curve , combined with the priors on @xmath53 and @xmath373 , provides a constraint on @xmath169 ; while the shape of @xmath374 allows dark energy parameters to be constrained via the apparent dependence of @xmath342 on distance . the results for spatially flat , constant @xmath313 models and non - flat @xmath8cdm models are respectively shown in figures  [ fig : constraints : mantz_w ] and [ fig : constraints : allen_wm_wl ] . for the non - flat models , the presence of dark energy is detected at high ( @xmath375 ) confidence , comparable to current snia results ; the constraints , including systematic uncertainties , are @xmath376 and @xmath377 , with the flat @xmath8cdm model yielding an acceptable goodness of fit . for constant @xmath313 models , the @xmath342 data are again competitive with other cosmological results , obtaining @xmath378 . the systematic , cluster - to - cluster scatter in @xmath342 is small , @xmath379 , corresponding to only 5% in distance ; this high precision results from the restriction to hot , dynamically relaxed systems for which total masses can be accurately estimated ( section  [ sec : obs_masses_xray ] ) . interestingly , for clusters with @xmath380 , the measured @xmath342 values show no dependence on temperature ( right panel of figure  [ fig : constraints : allen_wm_wl ] ) , indicating that , for the most massive clusters , the self - similar expectation of constant @xmath342 with mass is realized . at lower temperature and mass ( extending to the group scale ) , a trend of increasing @xmath342 with temperature and mass is observed ( e.g. @xcite ) . = 5.5 cm = 5.5 cm    @xcite have also applied x - ray data to the @xmath342 test , using _ measurements at @xmath381 for 52 clusters in the redshift range @xmath382 , adopting similar priors on @xmath53 and @xmath373 and marginalizing over a prior on @xmath383 ( assumed constant with redshift ) . however , their data set was not restricted to dynamically relaxed systems , resulting in significantly weaker constraints ( table  [ tab : constraints : basic ] ) . @xcite employed _ chandra _ x - ray and ovro / bima sz observations of 38 clusters in the redshift range @xmath384 ( with no restriction dynamically relaxed systems ) , finding @xmath342 values consistent with previous x - ray work . ( their analysis did not take advantage of the relative normalization of the x - ray and sz signals to simultaneously provide a second distance constraint ; see section  [ sec : constraints : xsz ] . ) they adopted the simpler approach of assuming constant @xmath342 and marginalizing over its value , incorporating a _ wmap _ prior on the total density , @xmath385 . although this explicitly ignores the information available from the normalization of @xmath374 , their results clearly disfavor a dark matter dominated universe , preferring a low - density universe with dark energy ( table  [ tab : constraints : basic ] ) . the different dependence on distance of the gas density inferred from x - ray and sz observations of clusters can be exploited in a conceptually similar way to @xmath342 data ( section  [ sec : theory : xsz ] ) . the most recent contribution is that of @xcite , who measured distances to 38 clusters at redshifts @xmath384 . this cosmological test is intrinsically less sensitive to distance than the @xmath342 test , with the signal proportional only to @xmath386 ( equation  [ eq : daxraysz ] ) , and currently can constrain only one free parameter . assuming spatial flatness and fixing @xmath387 , @xcite obtained a constraint on the hubble parameter , @xmath388 , consistent with results from other data such as the hubble key project @xcite or the combination of @xmath342 and cmb data @xcite . we note that other works using the same method have typically found somewhat lower best fitting values ( @xmath3890.7 ; e.g. @xcite ) , but these discrepancies are not significant given the systematic uncertainties .      the cmb temperature power spectrum at multipoles @xmath390 encodes the thermal sz signature of unresolved clusters at all masses and redshifts ( section  [ sec : theory : tsz ] ) . although the primary cmb power decreases rapidly at these scales , extracting this cosmological information from the tsz spectrum has proved challenging due to uncertainties in , e.g. , the relevant observable  mass scaling relation at low masses and high redshifts ; the population of infrared and radio point sources ; the magnitude of the integrated kinetic sz effect ; and the form of the electron pressure profile at large cluster radii , where it is poorly constrained by current x - ray data ( @xcite ; section  [ sec : future_outskirts ] ) . the _ wmap _ , spt and act collaborations have all detected excess power at large multipoles , @xmath391 . subject to the systematic uncertainties mentioned , their results are broadly in agreement , and are consistent with estimates of @xmath58 obtained from studies of resolved clusters and the primary cmb ( @xcite ) . as discussed above , current cluster growth and @xmath342 data can constrain spatially flat models with constant @xmath313 , finding consistency with the cosmological constant model ( @xmath392 ) . constraints on constant @xmath313 models from the combination of these cluster data are shown in the left panel of figure  [ fig : physics : wevol ] . to go beyond this simple description of dark energy , it is necessary to include cosmological data from additional sources in the analysis . equation  [ eq : wwa ] provides a simple and commonly adopted model of evolving dark energy , in which the equation of state takes the value @xmath393 at @xmath110 and approaches @xmath394 at high redshift . constraints on this model were obtained from cluster growth data by @xcite and @xcite , assuming spatial flatness and in combination with external cmb , snia and bao ( and @xmath342 , in the case of @xcite ) data . in both studies , the results are consistent with the @xmath8cdm model ( @xmath395 and @xmath396 ; right panel of figure  [ fig : physics : wevol ] ) .    a slightly more general model due to @xcite , @xmath397 also makes a smooth transition from one value at the present day ( @xmath393 ) to another at early times ( @xmath398 ) , but has the advantage that the transition redshift , @xmath399 , can be marginalized over . ( when @xmath400 , equation  [ eq : physics : lcdm_wetmodel ] reduces to equation  [ eq : wwa ] with @xmath401 . ) @xcite used the data sets above to obtain @xmath402 and @xmath403 , again consistent with @xmath8cdm .    absent a concrete physical model for dark energy , phenomenological models of evolving dark energy can take any form ; for example , one possibility that has not yet been investigated using cluster data expands @xmath404 in principal components ( e.g. @xcite ) . although these ad - hoc descriptions of dark energy are perfectly straightforward to apply to measurements of the expansion history , their applicability to an analysis based on the cluster mass function or power spectrum is less clear . generically , descriptions of dark energy as a fluid with @xmath405 should include the effects of spatial variations in dark energy density ( e.g. @xcite ) . however , the results from cluster growth so far have made use of mass functions from @xmath8cdm simulations , accounting for the value of @xmath313 only in the expansion history ( @xcite ) or , at most , including the effect of density variations on the linear matter power spectrum ( @xcite ) . in principle , the cluster mass function must also be adjusted to account for the behavior of fluid dark energy on smaller scales and in higher density environments . encouragingly , preliminary work in this area suggests that the influence of dark energy on the mass function might be readily measurable , resulting in additional constraining power from clusters ( e.g. @xcite ) . such improvements , however , might come at the cost of requiring additional sophistication in the theoretical description of dark energy ( e.g. the dark energy sound speed and viscosity ; @xcite ) . the @xmath8cdm paradigm , while providing an excellent model for the large scale structure of the universe , incorporates little information on the physical nature of dark matter . it assumes only that dark matter is non - baryonic ; that it interacts weakly with baryonic matter and itself ; that it emits and absorbs no detectable electromagnetic radiation ; and that the dark matter particles move at sub - relativistic speeds . as clusters merge under the pull of gravity , their dark matter halos and x - ray emitting gas can become separated temporarily , as the gas experiences ram pressure and is slowed . the observed offsets between the dark matter and x - ray peaks in the bullet cluster ( 1e0657 - 558 ; @xcite ) and macsj0025.4 - 1222 ( @xcite ; figure  [ fig : bullet_macs0025 ] ) , both massive merging systems with relatively simple geometries , require conservatively that the scattering depth for the merging dark matter can not be greater than one . using gravitational lensing data to estimate the dark matter column densities through these clusters , @xcite and @xcite use the observed dark matter and x - ray peak separations to derive limits on the velocity - independent dark matter self - interaction cross - section per unit mass of @xmath406@xmath407g@xmath408 and @xmath409@xmath407g@xmath408 for 1e0657 - 558 and macsj0025.4 - 1222 , respectively . @xcite additionally use the non - detection of an offset between the lensing peaks and the galaxy centroids for the bullet cluster to refine this constraint to @xmath410@xmath407g@xmath408 . upcoming surveys ( section  [ sec : new_surveys ] ) , should provide hundreds of similar examples , removing the current systematic limitations set by small number statistics . in combination with multiwavelength follow - up observations and improved numerical simulations ( e.g. * ? ? ? * ) , this should allow the properties of dark matter in merging clusters to be studied in a robust , statistical manner . one of the most remarkable predictions of the cdm model is that the density profiles of relaxed dark matter halos , on all resolvable mass scales , can be approximated by a simple , universal profile ( equation  [ eq : nfw ] ) with an inner , density slope @xmath411 . in contrast , for dark matter models with significant self - interaction cross sections , halos are expected to exhibit flattened , quasi - isothermal cores @xcite . in the absence of significant rotational support , these cores are also expected to be approximately spherical .    using chandra x - ray data for a sample of massive , dynamically relaxed galaxy clusters , @xcite measure a mean central density slope of @xmath412 ( 95 per cent confidence limits ) , in good agreement with @xmath8cdm . detailed strong - plus - weak lensing analyses for a subset of these systems yields consistent constraints @xcite . on smaller scales , using stellar velocity dispersion profiles for the dominant cluster galaxies , @xcite ( see also @xcite ) infer a possible flattening of the central dark matter halos in abell 383 and ms2137.3 - 2353 . however , on these small scales , the impact of baryonic physics becomes important . @xcite use the lack of a dark matter core in x - ray and lensing data for the relaxed cluster ms1358 + 6245 , to place a limit on the velocity independent dark matter particle - scattering cross section @xmath413cm@xmath272g@xmath408 . @xcite use constraints on the ellipticity in the central regions of ms2137.3 - 2353 to infer @xmath415cm@xmath272g@xmath408 . to improve these constraints , combined multiwavelength observations for large samples of relaxed clusters , coupled with improved simulations modeling fully the interactions between dark matter and baryons , are required . certain dark matter candidates , including sterile neutrinos , possess a two - body radiative decay channel that produces a photon with energy @xmath416 , where @xmath417 is the dark matter particle mass ( e.g. * ? ? ? galaxy clusters have been the targets of searches for emission lines associated with such decays . the soft x - ray ( kev ) regime is particularly interesting , marking the lower limit of masses consistent with constraints from large scale structure formation . to date , all searches for monochromatic x - ray emission lines associated with non - baryonic matter in clusters ( as well as other dark matter rich objects ) have proved negative ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? gamma - ray observations of clusters with the fermi gamma - ray space telescope have been used to place interesting upper limits on dark matter annihilation , and the lifetimes of particles for a range of masses and decay final states @xcite . dark energy , though a key component of the standard cosmological model , provides by no means the only possible explanation for cosmic acceleration . various non - standard gravity models can also produce acceleration on cosmological scales ( for reviews , see @xcite ) . these include frameworks that consistently parametrize departures from general relativity ( gr ; @xcite ) ; full , alternative theories such as the dvali - gabadadze - porrati ( dgp ) braneworld gravity @xcite ; @xmath418 modifications of the einstein - hilbert action @xcite ; and modifications of gravity based on the mechanism of ghost condensation @xcite . a critical requirement for any modified gravity model is that it should mimic gr in the relatively small scale , high density regime where gr has been tested precisely . thus , in addition to investigating whether dark energy is well described by a cosmological constant , we are simultaneously interested in asking whether gr provides the correct description of gravity , and indeed whether dark energy is needed at all .    to discriminate among these possibilities , and between particular dark energy and modified gravity models , it is important to combine expansion history data with measurements of the growth and scale - dependence of cosmic structure . galaxy clusters provide some of our strongest constraints on structure growth . to utilize these constraints robustly , however , accurate predictions for the halo mass function are required . recently , a few mass functions for specific modified gravity models have been constructed and calibrated using n - body simulations . these include the self - accelerated branch @xcite and normal branch @xcite of dgp gravity , and an @xmath418 model @xcite . constraints on the latter model using the observed local cluster abundance and other data are presented by @xcite . an alternative to evaluating specific gravity theories is to adopt a convenient , parameterized description for the growth of structure . this can then be used to constrain departures from the predictions of @xmath8cdm+gr @xcite . at late - times , the linear growth rate can be simply parametrized as ( e.g. @xcite ) @xmath419 where @xmath420 is the density contrast and @xmath421 the growth index . conveniently , gr predicts a nearly constant and scale - independent value of @xmath422 for models consistent with current expansion data . as in the case of @xmath313 for dark energy models , constraining @xmath421 constitutes a phenomenological approach to studying gravity . @xcite report constraints on departures from gr on cosmic scales using this parameterization with cluster data . their results are simultaneously consistent with gr ( @xmath423 ) and @xmath8cdm ( @xmath392 ) at the 68 per cent confidence level ( left panel of figure  [ fig : physics : gamma_s8 ] ) . the mass of neutrinos directly influences the growth of cosmic structure , since any particle with non - zero mass at some point cools from a relativistic state , in which it effectively suppresses structure formation , to a non - relativistic state , in which it actively participates in the growth of structure ( details are reviewed in @xcite ) . in the standard scenerio where the neutrino species have approximately degenerate mass , the species - summed mass , @xmath424 , is sufficient to describe their cosmological effects . although current data lack the precision to directly detect the effect of neutrino mass on the time - dependent growth of clusters , cluster data do play a key role in cosmological constraints on neutrinos when combined with cmb observations . on its own , the cmb can place only a relatively weak upper bound on the mass , @xmath425 at 95% confidence for a @xmath8cdm model ( e.g. @xcite ; works discussed here used 5 years of _ wmap _ data , but their conclusions are not significantly changed by the 7-year update ) . incorporating cosmic distance measurements improves this to @xmath426 , with the results displaying a strong degeneracy between @xmath424 and the value of @xmath58 predicted from the amplitude of the primordial power spectrum , due to the integrated effect of neutrinos on the growth of structure . cluster data at low redshift provide a direct measurement of @xmath58 , breaking this degeneracy ( right panel of figure  [ fig : physics : gamma_s8 ] ) , improving the upper limit by a further factor of two , to @xmath427 ( @xcite ) . the degeneracy - breaking power of cluster observations also significantly improves the robustness of neutrino mass limits to the assumed cosmological model ( e.g. marginalizing over global curvature ; @xcite ) . in combination with _ planck _ data , expected near - term improvements in cluster mass measurements from high - quality lensing data could reduce the limit on @xmath428 to the point of distinguishing between the normal and inverted neutrino mass hierarchies . in the near future ( over the next 23 years ) , the completion of the spt , act and _ planck _ sz catalogs will extend our detailed , statistical knowledge of galaxy clusters out to @xmath357 . together , these projects expect to find @xmath429 new clusters , mostly at intermediate - to - high redshifts @xcite . used in combination with existing low - redshift x - ray and optical catalogs , they should provide significant improvements in our knowledge of cluster growth , and corresponding improvements in the constraints on dark energy and gravity models . as discussed in section  [ sec : obs_surveys_sz ] , challenges for these surveys will be defining and calibrating the optimal survey mass proxies , and understanding the impact of contamination from associated infrared sources and agn . looking further ahead , the development of experiments with improved spatial resolution , expanded frequency coverage and improved sensitivity , such as the cerro chajnantor atacama telescope ( ccat ) , will be of advantage .    at optical and near - infrared wavelengths , a suite of powerful , new ground - based surveys are about to come on - line . these include the panoramic survey telescope and rapid response system ( pan - starrs ) ; the dark energy survey ( des ) ; the kilo - degree survey ( kids ) and complementary vista kilo - degree infrared galaxy survey ( viking ) ; the subaru hyper suprime - cam survey ( hsc ) ; and , later , the large synoptic survey telescope ( lsst ) . these experiments offer significant potential for finding clusters , and will provide critical photometric redshift and lensing data . a primary challenge in constructing optical cluster catalogs will be the definition of robust mass proxies with minimal , well - understood scatter across the mass and redshift ranges of interest . planned , space - based survey missions such as the wide field infrared survey telescope ( wfirst ) and euclid also offer outstanding potential for cluster cosmology , complementing the ground - based surveys in providing lensing masses and photometric redshifts , and extending cluster search volumes out to higher redshifts .    in the near term , x - ray cluster samples constructed from serendipitous detections in _ chandra _ and _ xmm - newton _ observations offer the potential for important gains ( e.g. * ? ? ? * ; * ? ? ? the main advances at x - ray wavelengths , however , will be provided by the erosita telescope on the spektrum - roentgen - gamma mission . scheduled for launch in 2012 , erosita will perform a four year all - sky survey that should detect an estimated 50,000100,000 clusters with excellent purity and completeness . as discussed in section  [ sec : obs_surveys_xray ] , the unambiguous detection of clusters at x - ray wavelengths requires angular resolution to distinguish point sources from extended cluster emission . at modest redshifts and high fluxes , this should be straightforward for erosita ; at high redshifts ( @xmath357 ) and faint fluxes , however , separating cool - core clusters from agn will be challenging and will require follow - up observations . here , we note more than a dozen known examples of powerful agn surrounded by x - ray bright clusters at intermediate - to - high redshifts ( e.g. @xcite ) . looking further ahead , the development of improved x - ray mirrors with high spatial resolution across a wide field of view , such as those proposed for the wide field x - ray telescope , would be a major advance , rendering trivial the removal of contaminating agn emission and allowing surveys to take full advantage of the center - excised x - ray luminosity as a low - scatter mass proxy ( section  [ sec : future : coreevol ] ) .    for all of these surveys , accurate calibration will be important ; this is a particular challenge for space - based missions . in comparison to other cosmological probes , the demands on photo-@xmath137 calibration will be relatively modest . this is due both to the pronounced 4000 break in early - type galaxy spectra and the ability to combine measurements for many galaxies per cluster . fisher matrix studies by @xcite suggest targets of @xmath430 for bias error and @xmath431 for error in the scatter in surveys with @xmath0 clusters . for surveys with fewer counts , shot noise dominates , and photo-@xmath137 errors become less important .    extensive programs of follow - up observations using high resolution , high throughput telescopes will also be essential . x - ray observatories like _ chandra _ , _ xmm - newton _ , _ suzaku _ and _ astro - h _ are likely to remain cornerstones of this work , providing excellent , low - scatter mass proxy measurements for individual clusters . we note that such information need only be gathered for a fraction of the clusters in a survey to gain a significant boost in constraining power with respect to self - calibration alone ( e.g. * ? ? ? statistical calibration of the mean masses of clusters in flux and redshift intervals from weak lensing measurements will also be critical : in order for the intrinsic power of large surveys not to be impacted severely , calibration of the mean mass at the few per cent level is required @xcite . to achieve this accuracy , detailed redshift information for the lensed , background sources will be needed ; indeed , the use of full redshift probability density functions rather than simple color cuts may be required . detailed simulations will also be needed to probe the systematic limitations of these measurements , and to advise on the best analysis approaches ( section  [ sec : future_sims ] ) . inflation predicts a near scale invariant power spectrum and nearly gaussian distribution for the primordial curvature inhomogeneities that seed lss . for slow - roll , single - field inflation , departures from gaussianity are currently unobservable ( by at least four orders of magnitude ; @xcite ) . however , other multi - field and single - field inflation models predict observable non - gaussianity ( ng ) . examples include certain brane models such as the dirac - born - infeld ( dbi ) inflation ; single - field models with non - trivial kinetic terms ; ghost inflation ; models in which density perturbations are generated by another field such as the curvaton ; and models with varying inflation decay rate . ( see @xcite and @xcite for recent reviews . ) any detection of ng would provide critical information about the physical processes taking place during inflation . in particular , a convincing detection of ng of the ` local ' type ( referring to particular configurations in fourier space ) would rule out not only slow - roll but _ all _ classes of single field inflation models @xcite . current measurements of cmb anisotropies and lss are consistent with gaussianity @xcite . in principle , measurements of the clustering and abundance of galaxy clusters can also be used to place powerful , complementary constraints . galaxy clusters trace the rare , high - mass tail of density perturbations in the universe ( section  [ sec : theory ] ) and are uniquely sensitive to ng . n - body simulations have been used to study cluster formation under non - gaussian initial conditions ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? analytical mass functions have also been calculated using the press - schechter formalism @xcite or excursion set theory @xcite .    to date , cmb and lss studies have only been used to place scale - independent constraints on ng . however , certain models that produce relatively large ng ( see above ) also present strong scale - dependence of the signal , due to , e.g. , a changing ( effective or otherwise ) sound speed @xcite . the addition of measurements at the smaller , galaxy cluster scale should provide important constraints on such models @xcite . several recent works @xcite have investigated quantitatively the potential of future optical and/or x - ray cluster surveys for constraining ng .      to date , robust thermodynamic measurements have , in general , only been possible for the inner parts of clusters ( @xmath432 ) , where the x - ray emission is brightest and the sz signal strongest ; a large fraction of their volumes remain practically unexplored . precise , accurate measurements of the density , temperature , pressure and entropy out to the virial radii of clusters provide important insight into the physics of clusters and ( ongoing ) large - scale structure formation , and can improve the precision and robustness of cosmological constraints .    recently , the japanese - us _ suzaku _ satellite has opened a new window onto the outskirts of clusters . due to its lower instrumental background than flagship x - ray observatories like _ chandra _ and _ xmm - newton _ , which orbit beyond the earth s protective magnetic fields , _ suzaku _ can study the faint , outer regions of clusters more reliably . the primary challenge with _ suzaku _ is the relatively large point spread function of its mirrors , which limits detailed , spatially - resolved studies to systems at modest redshift ( @xmath433 ) .    over the past two years , a series of ground breaking measurements of the outskirts of clusters with _ suzaku _ have been reported @xcite . these have confirmed the presence of smoothly decreasing density and temperature profiles out to large radii , as was qualitatively expected from theoretical models and earlier data ( e.g. * ? ? ? * ; * ? ? ? interestingly , these observations have also suggested a possible flattening of the outer entropy profiles .    of particular interest are the advances made with the _ suzaku _ key project study of the perseus cluster ( @xmath434 ) , the nearest , massive galaxy cluster and brightest , extended extragalactic x - ray source . the temperature and metallicity profiles for the northwestern ( nw ) and eastern ( e ) arms of perseus measured with _ suzaku _ are shown in figure  [ fig : suzaku_kt ] @xcite . also plotted , for comparison , are the results from earlier , deep _ observations of the cluster core @xcite . the _ suzaku _ and _ chandra _ data show excellent agreement where they meet , and measure the thermodynamic structure of the icm over three decades in radius , out to @xmath233 . models of large scale structure formation show that gas is shock heated as it falls into clusters . entropy is an important tracer of this virialization process . numerical simulations predict that in the absence of non - gravitational processes such as radiative cooling and feedback , the entropy , @xmath435 , should follow a power - law with radius , @xmath436 , with @xmath437@xmath438 ( e.g. * ? ? ? . with the exception of a cold front seen along the eastern arm at @xmath439 , the entropy profile in perseus roughly follows the expected trend out to @xmath440 ( figure  [ fig : suzaku_thermo ] ) . beyond this radius , however , both arms flatten away from the predicted power - law shape . the pressure profile shows good agreement between the nw and e arms . within @xmath441 , the _ suzaku _ pressure results show good agreement with the predictions from numerical simulations @xcite , and previous measurements with _ xmm - newton _ @xcite . at larger radii , however , the observed pressure profile is shallower than a simple extrapolation of the _ xmm - newton _ results . = 6.5 cm = 6.5 cm    accurate estimates of the gas masses and total masses out to large radii are of particular importance for cosmological studies . the _ suzaku _ observations of the perseus cluster provide the first such measurements for a massive cluster . the best - fit nfw mass model determined from a hydrostatic analysis of the ( relatively relaxed ) nw arm has parameters in good agreement with the predictions from cosmological simulations ( @xmath442 , @xmath443 ) and provides a good description of the data @xcite . of particular interest is the cumulative @xmath342 profile , shown in the right panel of figure  [ fig : suzaku_thermo ] . within @xmath368 ( @xmath444 ) , the observed @xmath342 profile is consistent with previous _ chandra _ and sz measurements for other massive , relaxed clusters @xcite . from @xmath445@xmath446 ( i.e. excluding the central cooling core ) the @xmath342 profile is also consistent with the predictions from recent numerical simulations , incorporating a simple model of agn feedback @xcite . at @xmath447 , the enclosed @xmath342 value approximately matches the mean cosmic baryon fraction , as measured from the cmb @xcite . beyond @xmath448 , however , where the entropy also flattens away from the expected power - law behavior , the @xmath342 apparently _ exceeds _ the mean cosmic value . the most plausible explanation for the apparent excess of baryons at large radii is gas clumping : in x - rays , the directly measurable quantity from the intensity of the emission is the average of the square of the electron density , rather than the average electron density itself . if the density is not uniform , i.e. the gas is clumpy , then the average electron density estimated from the x - ray intensity will overestimate the truth ; this will lead to an overestimate of the gas mass and gas mass fraction , and will flatten the apparent entropy and pressure profiles .    the amount of gas clumping in perseus , as inferred from comparison of the observed and expected @xmath342 profiles , is shown in the lower panel of the right of figure  [ fig : suzaku_thermo ] . using this clumping profile to correct the pressure and entropy measurements gives the solid red curves in the left panel . the clumping - corrected entropy profile shows good agreement with the expected power - law form out to @xmath233 . likewise , the clumping - corrected pressure profile matches the form predicted by simulations @xcite . importantly , the _ suzaku _ results provide no evidence for the puzzling deficit of baryons at @xmath449 inferred from some previous studies of massive clusters using lower quality x - ray data ( at large radii ) and stronger modeling priors ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? the perseus data suggest that beyond the innermost core but within @xmath450 , x - ray measurements for massive , relaxed clusters can be used simply and robustly for cosmological work . at larger radii , the effects of gas clumping become increasingly important and must be accounted for .    in principle , the combination of x - ray and sz observations can also be used measure gas clumping , offering an important cross - check of the suzaku results . the origin of these density fluctuations is also an important , open question . although numerical simulations predict that the gas in the outskirts of clusters should be clumpy ( e.g. * ? ? ? * ) , the degree of inhomogeneity predicted depends in detail on a range of uncertain physical processes , including cooling , conduction , viscosity , and the impact of magnetic fields .    over the next few years , _ suzaku _ and sz studies of other bright , nearby clusters , complemented later by higher spectral resolution x - ray data from _ astro - h _ , can be expected to stimulate significant progress . in particular , @xmath342 measurements out to large radii for a statistical sample of nearby clusters , and measurements of the dispersion in this and other properties from region - to - region and system - to - system , should provide a robust low - redshift anchor for cosmological work and powerful constraints on astrophysical models . cluster cores  typically describing the central region of @xmath451100 radius  often host a variety of astrophysical processes including efficient radiative cooling , agn outbursts , modest star formation , and sloshing or other bulk motions ( see , e.g. , @xcite for a review ) . in the cores of even the most dynamically relaxed clusters , small - scale departures from hydrostatic equilibrium are often apparent ( @xcite ) . the astrophysics of cluster cores is especially important for cosmological studies at x - ray wavelengths . in particular , clusters for which radiative cooling in the core is efficient form an easily distinguishable sub - population with significantly enhanced central gas density and luminosity , also commonly accompanied by a drop in icm temperature @xcite . the sharp , central density peaks of these ` cool core ' clusters enhance their detectability at x - ray wavelengths . however , without spatial resolution @xmath286  arcsec or additional follow - up data , cool - core clusters at high redshift lying close to the survey flux limit can not be distinguished easily from non - cluster x - ray point sources . ( this selection bias has affected some studies in the literature . ) the prevalence and evolution of cool cores in the cluster population thus plays a significant role in determining the shape and evolution of the scatter in x - ray luminosity at fixed mass . confirmation of this fact can be found in the dramatic reduction in scatter obtained when emission from a central region of radius @xmath452 is excised from luminosity measurements in forming the @xmath208@xmath71 scaling relation , from @xmath453 to @xmath296 ( figure  [ fig : corescatter ] ; see also @xcite ) . ( while the excision radius of @xmath452 has become somewhat conventional , a similar reduction in scatter is evident when excising a fixed metric radius of similar size , e.g. 150 . ) @xcite found that roughly half of the intrinsic scatter can be attributed to radii @xmath454 , the typical scale of cool , dense cluster cores , with most of the remainder due to variations in the gas density profile at @xmath455 . the fraction of the total flux coming from @xmath456 can be as large as 50% , and correlates strongly with other observable signatures of cool cores such as central cooling time and cuspiness of the surface brightness profile ( @xcite ; see also @xcite ) . interestingly , the fractions of cool - core systems identified using this criterion are comparable in x - ray flux - selected samples , at least within @xmath457 ( @xmath269 ; see figure  [ fig : corescatter ] ; @xcite ) , although these fractions are biased relative to the full population due to selection effects . cool cores are also common in x - ray selected samples at high redshift ( @xcite ; @xmath458 ) . = 6.7 cm    although the number of massive clusters found in rass is too small to detect evolution or departures from log - normality in the luminosity  mass scatter , future , deeper x - ray surveys such as erosita will provide the larger samples necessary for these investigations . conversely , a better understanding of cluster cores will be required to fully exploit these surveys for cosmology . such studies will need to employ the full statistical apparatus described in section  [ sec : constraints : procedure ] to obtain unbiased results .    due to the weaker dependence of the sz signal on the central gas density ( section  [ sec : obs_observations ] ) , sz surveys are expected to be less sensitive to the presence of cool cores than x - ray surveys . this is consistent with the relatively low fraction of cool core systems observed among the newly discovered clusters in the _ planck _ early release cluster catalog @xcite . we note , however , that current sz surveys are not necessarily immune to biases associated with cool cores ; depending on the selection techniques employed , infrared emission associated with star formation in the central regions of cool core clusters may actually diminish the measured sz signal . unleashing the full statistical power of upcoming surveys will require careful control of theoretical uncertainties . for cluster abundance and evolution tests , uncertainties in scaling relation parameters are currently more important than uncertainties in the halo mass and bias functions @xcite . with aggressive mass calibration efforts to reduce the former , the mass function and bias errors will need to be limited to the percent level in order to avoid significant degradation in cosmological parameter constraints @xcite . simulation campaigns will be needed to address this challenge .    over the next decade , increased computing power will enable models with new capabilities ( e.g. multi - scale simulation in place of sub - grid models ) and will vastly expand the size of current simulation suites . by densely sampling a large control space of cosmological and astrophysical parameters , simulation ensembles can support survey analysis via functional interpolation , a method termed emulation @xcite . an initial application of this technique uses 38 @xmath459-particle n - body simulations to predict the non - linear matter power spectrum to @xmath460 accuracy ( for @xmath461 ) in a 5-dimensional space of cosmological parameters @xcite . precise n - body calibrations must be treated with some caution , as the baryons representing @xmath462 of the total mass undergo different small - scale dynamics than cdm . back - reaction effects of cooling and star formation could be important @xcite and should be systematically investigated . as discussed in sections  [ sec : physics : lcdm ] and [ sec : physics ] , studies of lss in non - standard cosmologies , including those for which dark energy may cluster or in which weak - field gravity is non - newtonian , should also be pursued . full solution of the galaxy formation problem from first principles remains challenging . direct simulation methods will continue to improve , though not likely to the point where multi - fluid simulations will offer definitive predictions for next - generation survey analysis . semi - analytic methods are increasingly informative ; for example , the frantic early merging that forms the central galaxies in clusters is becoming understood @xcite . but many elements required to predict the spectrophotometric properties of galaxies remain poorly understood , including stellar population synthesis and dust evolution @xcite . empirically tuned statistical approaches , essentially assigning galaxy properties to halos or sub - halos in a manner informed by observed clustering as a function of luminosity and/or color , will continue to provide a valuable complement to physical methods . upcoming surveys will require sophisticated , automated reduction and analysis pipelines for their large data streams , along with quality assurance to validate accuracy . simulations of sky expectations , often referred to as ` mock ' or ` synthetic ' surveys , will be required to provide realistic testbeds in which the answers are known . mock surveys provide key quality assurance support by : ( i ) incorporating line - of - sight projection effects on measured properties ; ( ii ) including distortions and noise from telescope / camera optics and other sources ; and ( iii ) validating image processing and data management pipelines . since the effort involved in producing such data is non - trivial , mechanisms to publish and enable their broader use within the community need to be pursued . with such significant improvements in data quality expected and such profound questions to be addressed , modeling considerations will become increasingly important . the key issues will be modeling and mitigating all important sources of bias and systematic error in the analyses , and using the information efficiently .    for tests based on the mass function and clustering of clusters and their associated mass  observable scaling relations , the impact of survey biases must be accounted for ( section  [ sec : theory : selection ] ) . typically this requires simultaneous modeling of the cluster population and scaling relations in a single likelihood function . such an approach also facilitates understanding the covariance between model parameters , and provides a structure within which to examine the impact of residual systematic uncertainties which often correlate with model parameters . to a large degree , the statistical frameworks required for such analyses have been developed ( section  [ sec : constraints : procedure ] ) and their application to existing data are discussed in this review . the application of judicious , blind ( see below ) cuts can improve the balance of statistical versus systematic uncertainties . for example , with the @xmath342 test ( section  [ sec : constraints : fgas ] ) , cosmological constraints are best derived from the most massive , dynamically relaxed clusters , which can be identified easily from short , snapshot _ chandra _ or _ xmm - newton _ x - ray observations . these clusters are also ideal targets for the xsz test ( section  [ sec : constraints : xsz ] ) . for the mass function and clustering tests , one must determine the optimal mass / flux / redshift range over which to compare the data and models , given the characteristics of the survey and follow - up data and one s understanding of mass  observable scaling relations . strong priors should be used with caution and their implications understood . in cluster cosmology , common assumptions include log - normal scatter in the mass  observable scaling relations , and negligible evolution of this scatter with redshift . when measuring masses from x - ray or optical dynamical data , the use of parameterized temperature , velocity and gas density models is also common . where such strong assumptions are not necessary , they should be avoided . where priors do become necessary , for example in parameterizing the impacts of known astrophysical effects , the validity of these assumptions should be checked empirically , where possible . the covariance between different measured quantities is an important consideration that is often wrongly neglected . the simplest example of this is simply the statistical correlation of quantities measured from the same observation ( e.g. due to poisson noise ) . a more subtle case follows from the definition of cluster radius for scaling relations in terms of the mean density enclosed : because the measured mass and radius covary , other quantities measured within that radius necessarily covary with mass , even if measured from independent observations . fortunately , within the framework of bayesian analysis , it is straightforward to account for such measurement correlations . for example , given a set of scaling relation parameters to be tested , the data likelihood can be integrated over all possible true values of , e.g. , the mass and temperature , with one of the terms in the integrand being the probability of the true mass and temperature values producing the observed values ( the sampling distribution ; see section  [ sec : constraints : procedure ] ) . the form of this probability density can be as general as is required  a multidimensional gaussian or log - gaussian in the simplest case  and in particular may have non - zero correlation . * see also @xcite ) discusses this general approach in the astrophysical context , and provides useful tools for bayesian linear regression .    in other situations , unnecessary covariance can be introduced by the model applied to the data . an example is the historically common practice of fitting parametric temperature and gas density profiles to x - ray data , and inferring overdensity radii and hydrostatic masses from these model profiles . here , the procedure introduces an artificial covariance between mass and temperature ( i.e. a prior on the scaling relation ) which could be avoided simply by modeling the temperature and gravitating mass profiles independently .    where binning will result in a loss of relevant information , it should be avoided . for example , the binning of x - ray surface brightness data to determine an integrated x - ray luminosity results in a significant loss of information ; the center - excised luminosity provides a tighter mass proxy ( section  [ sec : future : coreevol ] ) . hypothesis testing will remain a critical element of future analyses . wherever a model is fitted to data , one should check that it provides an adequate description ( i.e. that the goodness of fit is acceptable ) ; if it does nt , then the model can be ruled out . this is particularly important in the context of cosmological surveys and scaling relations that are subject to selection bias , since no straightforward , visual check of the goodness of fit is typically possible . where the simplest models fail to describe the data , one must evaluate carefully the additional degrees of freedom needed to alleviate the tension , in particular considering both astrophysical and cosmological possibilities ; simulations will play a critical role in motivating and validating alternatives to the simplest scaling models . a related situation arises in the combination of constraints from independent experiments : before combining , one should ensure that the model in question provides an adequate description of the data sets individually , and that the parameter values are mutually consistent , i.e. that their multi - dimensional confidence contours overlap . here it can be helpful to include the impact of all known systematics in the contours ; where contours do not overlap , the model is incomplete and/or unidentified systematic errors are present . historically , the combination of mutually _ inconsistent _ data sets has sometimes led to unphysically tight formal constraints . training sets can be useful to tune analyses . however , the biases introduced by this training must be understood and accounted for . the inclusion of sufficient redundancy into experiments ( i.e. having more than one independent way to make a measurement ) can also enhance significantly the robustness of conclusions . in cosmological studies , for example , we require more than a single way to measure both the growth and expansion histories , e.g. clusters and lensing , and snia and bao , respectively . a final modeling consideration is the impact of experimenter s bias : the subjective , subconscious bias towards a result expected by an experimenter . this is important where we are seeking to measure a quantity for which prior expectations exist , e.g. @xmath313 or the growth index , @xmath421 . evidence for experimenter s bias can be found readily in the literature , e.g. in historical measurements of the speed of light or the neutron lifetime ( see @xcite for a review ) . the best approach to overcoming experimenter s bias is blind analysis . a range of blind analysis techniques are commonly employed in the particle and nuclear physics communities . these are used to mask the ability of scientists to determine quantities of interest until all methods are finalized , the systematic errors that can be identified have been , and the final measurement is ready to be made . a technique with good applicability to cosmological studies is _ hiding the answer _ , wherein a fixed ( unknown to the experimenter ) offset is added to the parameter(s ) of interest , and only revealed once the final measurements are made . before un - blinding , it is advisable to think through how to proceed afterward , and what additional checks to employ . double - blind analyses , where two independent teams repeat the same process ( often used in biomedical research ) target additional sources of error . such techniques can be used powerfully , where resources allow . our article has summarized the status of cluster cosmology and the methods used to extract cosmological information from galaxy cluster observations . these methods have been applied successfully to samples of tens to thousands of massive clusters . much of this work has been pioneered at x - ray and optical wavelengths , but sz surveys and gravitational lensing studies are also now poised to play central roles . the availability of powerful , new cluster surveys will help to address profound mysteries such as the origin of cosmic acceleration , inflation and the nature of gravity . this work will require a multiwavelength approach , combining the strengths of the available techniques for finding clusters , calibrating their masses and obtaining low - scatter mass proxies . the analysis of these data will require dedicated efforts by large teams of researchers . the demands for follow - up observations will be high , and will require the support of time allocation committees . building on the progress made , we can expect clusters to remain at the forefront of cosmological work through the next decade . by combining cluster measurements with other , powerful probes such as snia , cmb , bao and cosmic shear , we can be optimistic of having the precision , complementarity and redundancy required to allow robust conclusions to be drawn . we are grateful to our colleagues for their many valuable insights . in particular , we thank roger blandford , hans bhringer , marusa brada , joanne cohn , carlos cunha , megan donahue , harald ebeling , andy fabian , pat henry , dragan huterer , andrey kravtsov , dan marrone , tim mckay , glenn morris , daisuke nagai , aaron roodman , eduardo rozo , robert schmidt , tim schrabback , neelima sehgal , aurora simionescu , alexey vikhlinin , mark voit , anja von der linden , risa wechsler and norbert werner . we particularly thank david rapetti for his detailed input to the material on gravity tests and non - gaussianities . swa was supported in part by the u.s . department of energy under contract number de - ac02 - 76sf00515 . aee acknowledges support from nasa astrophysics theory grant nnx10af61 g . abm was supported by an appointment to the nasa postdoctoral program at the goddard space flight center , administered by oak ridge associated universities through a contract with nasa . this research was supported in part by the national science foundation under grant no . nsf phy05 - 51164 .
Abstract:

Answer: 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 . keywords : cosmology , dark energy , dark matter , galaxy clusters , intracluster medium , large scale structure