massive galaxy clusters provide several means for independently examining any viable cosmological model .
cluster samples used for this purpose are usually complete in terms of x - ray and redshift measurements such as the rdcs ( rosati et al .
1998 , borgani et al . 1999 ) and the macs survey ( ebeling , edge , & henry 2001 , ebeling et al . 2007 ) that we examine here with lensing .
ideally , controlled samples of clusters with reliable lensing masses would allow the most direct comparison with theory but these are presently still in their infancy with only a handful of clusters discovered in wide - field weak lensing searches ( taylor et al .
2004 , wittman et al .
2006 , massey et al .
2007 , miyazaki et al . 2007 ) , and may be often subject to projection biases ( e.g. , hennawi & spergel 2005 ) .
weak lensing surveys should have the advantage of being approximately volume - limited over the redshift range 0.1<@xmath3<1.0 , where the lensing distance ratio has a broad maximum .
ongoing measurements of the sunyaev - zeldovich ( sz ) effect ( which has already been detected in many tens of clusters , including some at relatively high redshifts @xmath11 , e.g. , laroque et al .
2006 ) in a large number of clusters over a wide redshift range should improve upon x - ray selection , due to the redshift - independence of the effect ( though clearly high redshift clusters are relatively less well resolved ) .
also , since the insensitivity of the effect to internal gas density variations will not be biased as with x - ray flux selection towards systems with luminous shocked gas , measurements should better correlate with cluster masses .
strong gravitational lensing ( sl ) is nearly always seen in sufficiently detailed observations of massive clusters .
einstein radii derived from such observations provide a reliable projected mass within the observed einstein radius , which depends only on fundamental constants , g and c , and knowledge of the lens and source distances .
constraining the inner mass profile from sl is much harder , requiring many sets of multiple images covering a wide range of source redshift , to overcome inherent lensing degeneracy between the gradient of the mass profile and the scaling of the bend angle with source distance , so that sl based mass profiles have been usefully constrained for only several clusters producing reliable mass maps ( e.g. , abell 370 , kneib et al .
1993 , richard et al . 2009a ; abell 901 , deb et al . 2010 ; abell 1689 , broadhurst et al . 2005 , coe et al . 2010 ; abell 1703 , limousin et al . 2008 ; cl0024 + 1654 , liesenborgs et al .
2008 , zitrin et al . 2009b
; ms 2137.3 - 2353 , gavazzi et al .
2003 , merten et al .
2009 ; rxj1347 , brada et al . 2008a , halkola et al .
2008 ; sdss j1004 + 4112 , sharon et al .
2005 , liesenborgs et al .
2009 ; `` the bullet cluster '' , brada et al .
2006 ) .    in quite a few of these clusters
there is a significant inconsistency between the einstein radius directly measured from sl analyses and the einstein radius predicted from model profiles fitted to weak - lensing measurements of the same clusters , often subject to background selection and dilution problems over a wide radial range ( e.g. , medezinski et al .
2010 ) , though more so towards the sl regime where the degree of contamination by cluster members is often poorly corrected .
the importance of a sl analysis of a significant sample is clear , especially in the statistical aspect , where the modelling method is similar for all examined clusters and thus supplies a coherent view for objective comparison ( e.g. , richard et al .
2009b ) .
the development of sl modelling - methods has increased in response to the improvement in data quality since the discoveries of the first giant arcs ( e.g. , lynds & petrosian 1986 , soucail et al .
1987 , kneib et al .
most methods can be classified as `` parametric '' if based on physical parameterisation , and as `` non - parametric '' if they are `` grid - based '' ( see also  4.4 in coe et al .
2008 , and references therein ) . still , many of these methods include too many parameters to be well - constrained by the number of initially known multiply - lensed systems .
here we use acs imaging to identify new multiply - lensed systems in order to constrain the mass distributions and define the critical curves of the sample , motivated by the successful minimalistic approach of broadhurst et al .
( 2005 ) to lens modelling . in following work
( zitrin et al .
2009b ) we presented an improved modelling method which involved only 6 free parameters enabling easier constraint , since the number of constraints has to be equal or larger to the number of parameters in order to get a reliable fit .
two of these parameters are primarily set to reasonable values so only 4 of these parameters have to be constrained initially , which sets a very reliable starting - point using obvious systems .
the mass distribution is therefore primarily well constrained , uncovering many multiple - images which can be then iteratively incorporated into the model , by using their redshift estimation and location in the image - plane .
the massive cluster survey ( macs ) has been aimed to create a complete sample of the very x - ray luminous clusters in the universe , successfully increasing the number of known such clusters to hundreds , particularly at @xmath12 ( ebeling , edge , & henry 2001 ) . from this a complete sample of 12 high-@xmath3 macs clusters ( @xmath0 ) was defined by ebeling et al .
( 2007 ) which have proved very interesting in several follow - up studies including deep x - ray , sz and hst imaging .
detection of a large - scale filament has been reported in the case of macs j0717.5 + 3745 by ebeling et al .
( 2004 ) , for which many multiply - lensed images have been recently identified by zitrin et al .
( 2009a ) , revealing this object to be the largest known lens , with an einstein radius equivalent to 55@xmath13 ( for a source at @xmath14 ) . in the case of macs
j1149.5 + 2223 ( zitrin & broadhurst 2009 ) , a background spiral galaxy at @xmath15 ( smith et al .
2009 ) has been shown to be multiply - lensed into several very large images , requiring a very shallow , unrelaxed central mass distribution .
another large multiply - lensed sub - mm source at a redshift of @xmath16 has been identified in macs j0454.1 - 0300 ( also referred to as ms 0451.6 - 0305 ; takata et al .
2003 , borys et al .
2004 , berciano alba et al .
2007 , 2009 ) , macs j0025.4 - 1222 was found to be a `` bullet cluster''-like ( brada et al .
2008b ) , and other macs clusters have been recently used for an extensive arc statistics study ( horesh et al .
2010 ) .    the x - ray data available for this sample ( see ebeling et al .
2007 ) along with the high - resolution hst / acs imaging and sz imaging ( e.g. , laroque et al . 2003 ) make these macs targets particularly useful for understanding the nature of the most massive clusters . here
we complete a full sl analysis of this 12 cluster sample with the available deep hst / acs imaging , principally to derive their einstein radii and also to help constrain the central mass distributions .
the effective einstein radii derived here are the corresponding radii of circles of equivalent areas to those enclosed within the critical curves , or similarly , the radii within which the averaged @xmath17 is equal to 1 . in addition , this work can also supply detailed magnification maps to help motivate deeper searches for high-@xmath3 galaxies behind these clusters .    another major motivation for pursuing accurate lensing - maps
is the increased precision of model predictions for cluster- size massive halos in the standard @xmath9cdm model and close variance ( see umetsu & broadhurst 2008 , e.g. , bullock et al .
2001 , spergel et al .
2003 , 2007 , tegmark et al .
2004 , hennawi et al .
2007 , neto et al . 2007 , duffy et al .
2008 , keselman , nusser & peebles 2009 , meneghetti et al .
2010 , fedeli et al .
2010 ) . in the standard hierarchical model , large massive bodies are due to collapse recently and thus
are not expected to be found in large numbers at high redshifts . on the other hand ,
the larger volume available with increasing distance means that in practice we can not expect to reside next to the most massive cluster ( see also zitrin et al .
the increasing number of accurately analysed clusters have revealed larger lenses than predicted by the standard @xmath18 model ( broadhurst & barkana 2008 , puchwein & hilbert 2009 ) .
this discrepancy is empirically supported by the surprisingly concentrated mass profiles measured for such clusters , when combining the inner strong lensing with the outer weak lensing signal ( gavazzi et al .
2003 , broadhurst et al .
2005 & 2008 , limousin et al .
2008 , donnarumma et al .
2009 , oguri et al .
2009 , umetsu et al .
2009 , zitrin et al .
2009b,2010 ) , or independently , when using the internal dynamics of cluster members ( lemze et al .
2009 ) and deep x - ray data ( lemze et al . 2008 ) .
geometrically , lensing is considered to be optimised at intermediate redshifts , where for a given mass the critical density for lensing is minimal , but this is partly offset by the late hierarchical growth of high - mass systems .
this trade - off results in estimates of the amplitude of strong lensing to favour the redshift range @xmath19 .
our recent analysis of two x - ray selected macs clusters reveals very large lenses ( e.g. , macs j1149.5 + 2223 , zitrin & broadhurst 2009 ; macs 0717.5 + 3745 , zitrin et al
. 2009a ) at high redshift , increasing the tension with @xmath9cdm predictions , since large einstein radius clusters at high-@xmath3 require earlier formation than implied by the standard @xmath9cdm model ( sadeh & rephaeli 2008 ) .
the largest lensing clusters have proven to be excellent targets for accessing the faint early universe due to their large magnification consistently providing the highest redshift galaxies ( ebbels et al .
1996 , franx et al .
1997 , frye & broadhurst 1998 , bouwens et al .
2004 , kneib et al .
2004 , bradley et al .
2008 , zheng et al .
2009 ) .
the statistics of large einstein radii provide an important opportunity to test the standard @xmath9cdm paradigm , as it probes both the high - mass end of the cluster mass function and central mass distributions of massive clusters ( oguri & blandford 2009 ) .
increasing numbers of theoretical and observational works have been done recently , deriving einstein radius distributions of various samples .
for example , cypriano et al .
( 2003 ) have derived the mass distributions of 24 x - ray selected abell clusters via weak - lensing . by fitting to sis profiles they obtain a mean einstein radius of @xmath20 .
more recently , hoekstra ( 2007 ) finds by weak - lensing analysis and sis fitting a mean einstein radius of @xmath21 for a sample of 20 x - ray luminous clusters .
okabe et al . ( 2009 ) find the mass distribution and einstein radii for a sample of @xmath22 locuss clusters by fitting the weak - lensing data to cis profiles , obtaining a mean value of @xmath23 ( see okabe et al .
2009 and references therein ) .
note that these values quoted here are calculated from the corresponding tables in these papers , and are not necessarily scaled or normalised to certain lens and source redshifts .
other recent work more relevant to our study summarises sl analysis for a sample of 20 , mostly relaxed _
undisturbed _ clusters ( richard et al .
they find that the einstein radii are distributed log - normally with a peak at @xmath24 and a corresponding @xmath25 enclosed mass ( within r@xmath26250 kpc ) , and show that the predicted distribution of einstein radii from @xmath9cdm cosmology falls short of the observed einstein radii by a factor of 2 .
it has been claimed that much larger einstein radii can be contemplated only with mass distributions which are highly prolate and aligned along the line of sight ( corless & king 2007 , oguri & blandford 2009 , see also hennawi et al .
2007 , meneghetti et al .
2007 , sereno , jetzer & lubini 2010 ) . analysing the 12 x - ray selected , high-@xmath3 macs clusters is important since they are predicted to be massive and should form efficient lenses , whose properties can be then compared to such studies , playing an important role in probing the @xmath9cdm scenario .
the paper is organised as follows : in  2 we describe the sample and observations ; the lensing analysis is described in  3 detailing each cluster separately .
our overall results are presented and discussed in  4 , followed by a conclusion (  5 ) . throughout the paper
we adopt the standard cosmology ( @xmath27 , @xmath28,@xmath29 ) .
the @xmath0 macs clusters have been imaged with the wide field channel ( wfc ) of the acs installed on hst , mainly in the framework of proposal i d 9722 ( pi : ebeling ) .
this mainly comprises two - filter observations in the f555w and f814w bands , taken during 2003 and 2004 with typical exposure times of @xmath30s in each filter .
we make use of some other hst data available in the hubble legacy archive ( hla ) , which we list in table [ sample ] .
in addition , available sextractor ( bertin & arnouts 1996 ) photometry catalogues were also downloaded from the hla , and were used to construct colour - magnitude diagrams for identifying the red cluster sequence galaxies belonging to each cluster , as a starting point for each lensing model , as detailed in  [ model ] .    [ cols="^,^,^,^,^,^,^,^,^",options="header " , ]      the 12 clusters have effective einstein radii in the range , 11@xmath13 to 55@xmath13 , with a median value of @xmath31 , ( and a mean value of @xmath32 ) assuming a fixed source redshift of @xmath33 , see figure [ hist ] .
this corresponds to a range of physical radii , 68 kpc to 353 kpc , with a median value of 180 kpc , when transforming to the redshift of each lens . in the hierarchical model large - scale perturbations
collapse recently and thus should be found relatively locally .
this sample includes very large einstein radii , exceeding even the most impressive lenses previously studied at lower redshifts as analysed by broadhurst & barkana ( 2008 ) .
more recently , richard et al . ( 2009b ) carefully studied the lensing properties of 20 mainly undisturbed clusters at lower-@xmath3 , selected to have sl features in hubble data and measured x - ray data , for which their mean einstein radius is @xmath24 , a factor of @xmath34two smaller than our sample .
this difference reflects the larger masses of the macs sub - sample that we studied , which is purely x - ray flux selected and restricted to @xmath0 .
this macs sample therefore does not suffer from a lensing - selection bias , but the effect of projection bias must be taken into account when evaluating 2d lensing observations with 3d mass profiles predicted by theory ( e.g. , hennawi et al .
2007 , oguri & blandford 2009 , sereno , jetzer & lubini 2010 ) , as discussed in  [ cdm ] below .
to @xmath35 .
the solid curve is the expected distribution of einstein radii as calculate in  [ cdm ] for the @xmath9cdm model , incorporating both the scatter on the @xmath36 relation and the projection bias from lensing .
we multiply the resulting theoretical @xmath37 curve by the width of the bins to normalise it .
the observed distribution is skewed to larger einstein radii than predicted .
the median values of both the observed and the theoretical distributions are plotted on the histogram in a blue dashed line , and a black dash - dotted line , respectively.,width=302 ]    cdm model ( black solid curve ) described in  [ cdm].the median values of both the observed and the theoretical cumulative distributions are plotted on the histogram in a blue dashed line , and a black dash - dotted line , respectively.,width=302 ]    , is tightly followed by the data , indicating that in general the central mass distributions of these clusters are well behaved.,width=302 ]     with the system scale size is apparent .
, width=302 ]    most of these clusters have resisted analysis by strong lensing despite the long availability of the hubble data .
our success in `` cracking '' these irregular lenses and defining their critical curves is encouraging for the application of our approach to complex unrelaxed clusters more generally .
for such clusters , where the deflection fields are not symmetric , the identification of multiple images must be aided by models that allow for flexibility in describing the mass distribution , as described in section [ model ] , free of the strictures imposed by the use of idealised elliptical potentials in conventional models .
the masses enclosed by the critical curves range from @xmath38 to @xmath39 with a median value of @xmath40 , and a mean value of @xmath41 as shown in table [ results_table ] .
these are calculated for each cluster by integrating the surface mass density distribution within the 2d model critical curve , scaled to a source redshift of @xmath42 .
we examine the relation between these einstein masses , @xmath43 enclosed within the critical curves , and the effective einstein radii derived above .
theoretically this should of course simply scale as @xmath44 for symmetric lenses . here
the lenses are not symmetric , but as can be seen in figure [ reme ] the quadratic relation is tightly followed , indicating that asymmetry is not very significant in terms of the mass distribution , and furthermore the measured or assumed redshifts of @xmath45 for the relevant systems is a reliable estimate .
in fact , we see clearly in that the 2d mass distributions are in general noticeably rounder than the critical curves which are very sensitive to substructure , as can be seen by comparing the 2d mass distributions with the critical curves , shown in figures [ curves0018 ] - [ contours2214 ] .
further encouragement for the accuracy of sl mass estimates is presented in recent work by meneghetti et al .
( 2010 ) , showing that sl methods based on parametric modelling are accurate at the level of few percent at predicting the projected inner mass .
richard et al .
( 2009b ) found a mean mass of @xmath25 enclosed mass ( within r@xmath26250 kpc ) for their sample of 20 clusters , while smith et al .
( 2005 ) examined the mass distribution of 10 x - ray selected clusters , and found a similar mean mass of @xmath46 ( see also richard et al .
these values are in good agreement with our sample median mass ( @xmath40 ) , and slightly lower than our sample mean mass which is boosted by few extremely massive clusters ( e.g. , macs j0717.5 + 3745 ) .
we examine the behaviour of the central mass - to - light ratio , @xmath47 , with the einstein radius and shown in figure [ ml ] .
the luminosity is summed over all cluster sequence galaxies identified as described in section [ model ] .
a clear correlation is obtained such that @xmath47 increases with the einstein radius , and with values quite typical of other massive clusters ( sarazin 1988 , rines et al .
2004 , medezinski et al .
this scaling towards higher @xmath48 provides confidence in our method . at the high mass end
our result is similar to the peak
m / l found in the detailed weak lensing profile studies of medezinski et al .
measuring @xmath47 in high-@xmath3 clusters is of interest also for the added insight on the amount of dm initially associated with individual galaxies as compared with the overall cluster dm component ( see also sarazin 1988 ) .
as can be seen in table [ results_table ] , most of the sample clusters provide large regions ( @xmath49 ) of highly magnified ( @xmath50 ) sky , useful for the search of the first stars and galaxies ( e.g. , zackrisson et al .
one has to be cautious when making such a statement since unlike the critical curves and the enclosed mass , the magnification is sensitive to the mass profile which requires source redshift estimation for several sources , which we have here for only a few clusters ( section [ model ] ) .
still , as explained in  2 , reasonable constraints are put on the mass profile slope by the image - plane minimisation , which generally has a broad minimum at the same point in the parameter space where the slope is also approximately correct . in addition , we are able to put a lower boundary on the high magnification area , by choosing the steeper profiles which naturally generate lower magnifications , or by assuming a boosted redshift for the main multiply - lensed system decreasing the lensing distance ratio for high-@xmath3 galaxies .
thus we constrain the area of high magnification for a high-@xmath3 source at the current limit of @xmath8 , showing that most of this sample clusters are excellent targets for high-@xmath3 galaxy searches due to these large high - magnification areas and lens power .
we estimate , based on previously analysed clusters ( e.g. , cl 0024 + 1654 , zitrin et al .
2009b ; macs j1149.5 + 2223 , zitrin & broadhurst 2009 ; macs j0717.5 + 3745 , zitrin et al .
2009a ) that in practice the areas of high magnification are about 10 - 20% higher than the lower limits presented in table [ results_table ] , with an upper limit of about 40 - 50% .
we generally follow broadhurst & barkana ( 2008 ; see also sadeh & rephaeli 2008 ) in constructing the theoretically predicted distribution of einstein radii . in particular , we adopt the nfw parameters measured by neto et al .
( 2007 ) for simulated halos , although we reduce @xmath51 by @xmath52 according to the results of duffy et al .
( 2008 ) who used the more recently measured values of the cosmological parameters .
studies based on large numerical simulations ( zhao et al .
2003,2009 , gao et al .
2008 ) have found for massive halos a rather weak decline of @xmath51 with increasing redshift ; we neglect this decline , which renders our results ( in terms of the disagreement with observations ) conservative .
we then correct the concentration parameter distribution using hennawi et al .
( 2007 ) to obtain the effective parameters for the population of clusters observed in projection .
note though that here we are considering a sample selected by x - ray flux ( not lensing cross - section ) , so there is no lensing bias , but there still is the projection bias ; this refers to the fact that when halos are seen in projection ( at randomly distributed angles ) , a different effective projected @xmath51 is measured than would be obtained from a spherical analysis of the halo density profile . in particular , projection increases the mean @xmath51 and adds substantially to the scatter compared to the distribution of the spherically - measured @xmath51 .
neto et al .
( 2007 ) split their halo sample into two groups ( `` relaxed '' and `` unrelaxed '' ) and analysed the statistics of each group separately , so we do the same , and in the end combine the two groups according to their relative numbers in the simulation by neto et al .
( 2007 ) , obtaining our predicted results for the total , combined population of halos .
this can then be compared with the observed sample , which does not select for relaxed clusters but includes both relaxed and unrelaxed ones .    in order to approximately simulate a flux - selected sample ,
we need to convert @xmath53 to halo mass .
we use the power - law relation from reiprich & bhringer ( 1999 ; with an updated hubble constant ) , along with the observed scatter of @xmath53 for a given mass , which we model as a lognormal distribution with a typical scatter that corresponds to a factor of 1.5 . note that this observed relation is based on @xmath54 , while the results from the above simulations refer to @xmath55 , so we must include the relation between @xmath54 and @xmath55 which is itself a function of @xmath51 .
an important caveat is that the observed @xmath56 relation was measured for relaxed clusters ( whose mass could be estimated from the x - ray emission assuming hydrostatic equilibrium ) .
we expect that unrelaxed clusters would tend to have an unusually high @xmath53 for a given mass ( e.g. , in a post - merger phase ) , which would insert more low - mass clusters within the sample than we assume , leading to smaller predicted einstein angles , increasing the discrepancy between the theory and the observations ( again making our results conservative ) .
on the other hand , it should be noted that much of the scatter in the @xmath56 relation probably arises from cool cores in the center ( since it is significantly reduced if the inner core is excised ) , implying a possible bias towards relaxed and more concentrated systems .
however , such a bias seems less probable in our case due to the high x - ray luminosity , the relatively high redshift , and the clear spread - out appearance of most of the clusters analysed here .
the observed properties of each halo ( flux , einstein radius distribution , etc . )
change slowly with redshift , so to simplify the calculation we approximate the clusters as all lying at a typical redshift @xmath57 . however , the halo abundance changes rapidly with redshift , so we calculate the integrated halo mass function for all halos observed at @xmath0 , in the fraction of the sky corresponding to the macs survey ( around a quarter of the sky ) .
we use the sheth & tormen ( 1999 ) formula , which accurately fits the halo mass function ( i.e. , the distribution of @xmath55 ) measured in simulations .
this halo mass function was convolved with with the @xmath51 distribution to obtain the distribution of @xmath54 , and then add the @xmath58 scatter in order to obtain the halo mass function corresponding to selection by a given minimum flux .
since the @xmath58 relation is not strictly valid for unrelaxed clusters , we only use the power - law and the scatter from this relation but allow some flexibility in the normalisation , choosing an effective minimum flux that corresponds to a predicted abundance of 12 clusters ( as in the observed sample ) .
this is reasonable also since we are interested here in testing the predicted @xmath59 distribution , not in testing the predicted cluster abundance .
note that the flux selection does not affect the distribution of the very highest @xmath59 values , since all clusters within the survey solid angle at @xmath0 that can produce the largest @xmath59 s will be well above the flux threshold .
given our predicted , x - ray selected mass function in terms of @xmath54 , we must convolve it with the @xmath59 distribution of each halo of mass @xmath54 . from the simulations we have the @xmath51 distribution ( and thus also the @xmath59 and @xmath54 distribution ) for a given @xmath55 , but we can invert this conditional probability using bayes theorem .
note that the @xmath55 to @xmath54 conversion should not really include the projection bias that is included in @xmath59 , but having two different @xmath51 distributions for each @xmath55 would makes things far more complicated .
our simplification of using a single @xmath51 distribution is reasonable since the projection bias typically affects @xmath54 ( for a given @xmath55 ) only by @xmath60 .
also , since @xmath51 affects both @xmath54 and @xmath59 in the same direction , the effect on our final result ( which derives @xmath59 for a given @xmath54 ) is much smaller .
our predicted @xmath59 distribution is compared to the observations in figure [ histr ] .
we compare the two cumulative distributions @xmath61 , i.e. , the total number of clusters expected above @xmath59 , where the theoretical distribution has been normalised ( as noted above ) to agree with the total of 12 clusters in the sample .
the observed distribution clearly lies at higher einstein radii compared with the predicted distribution .
one indication of this is the most discrepant cluster , j0717.5 + 3745 with @xmath62 .
the theoretical calculation yields a probability of only @xmath63 of finding such a large @xmath59 in this @xmath0 cluster sample , giving a significant ( though not extreme ) 2-@xmath64 discrepancy .
more generally , the observed distribution resembles the predicted one except offset towards higher angles by about a factor of 1.4 ( comparing , for example , the two medians ) .
the standard k - s statistic for comparing two distributions gives a probability of @xmath65 that the observed distribution is drawn from the predicted one .
this is also around a 2-@xmath64 discrepancy , largely independent of the above number since the k - s statistic focuses on the middle portion of the probability distributions and is insensitive to the edges .
the relatively low significance of the k - s discrepancy is an inevitable result of the poisson fluctuations expected with a sample of only 12 clusters .
the discrepancy that we find , of einstein radii that are too high by about a factor of 1.4 , appears smaller than the roughly factor of 2 difference found by broadhurst & barkana ( 2008 ) .
this may be explained at least partially by the difference between a lensing selected and an x - ray selected sample .
a more robust comparison would be possible with an x - ray selected cluster sample in which the clusters have reliable virial masses measured using weak lensing .
another caveat is that we have compared the observed einstein radii with the predictions based on pure dark matter simulations .
the possible effect of baryons on the halo density profile has been difficult to study with simulations , since hydrodynamic simulations that produce a large effect do so along with a central baryon concentration that disagrees with observations and is due to the `` overcooling '' problem in simulations .
analytical models suggest that even without producing a central baryon concentration , the complex coupled history of baryon and dark matter accretion may allow the baryons to significantly affect the final central density profile , potentially alleviating the discrepancy at least partially ( barkana & loeb 2010 ) .
on the other hand , recent simulations suggest that the discrepancy actually increases when agn feedback is taken into account to overcome the overcooling problem ( duffy et al .
2010 ) .      in this work
we determine the einstein radii and projected masses for the 12 clusters of the @xmath0 macs sample , and compare these to simulations .
we make use of our well - established modelling method to identify many sets of multiply - lensed images which are in turn used to constrain the models , and determine the einstein radii and projected masses for @xmath66 .
typical uncertainties in the source redshift estimates of @xmath67 impose only minor uncertainties on the comparison with @xmath9cdm .
firstly , such an uncertainty does not enable a unique determination of the profile slope , which we do not attempt to fully constrain here since both the einstein radius and the projected mass are indifferent to the mass profile slope and are determined by the data .
secondly , overestimating a source redshift would in practice increase the projected mass and the observed einstein radius for a source at @xmath68 , while maintaining the relation seen in figure [ reme ] and thus resulting only in growth of the discrepancy from @xmath9cdm simulations which are not dependent on the observed mass . on the other hand , underestimating a source redshift would in practice decrease the observed einstein radius and projected mass for a source at @xmath68 by less than 10% , thus insignificantly influencing the results .
the macs @xmath0 sample has proven to be of great value for many very different studies spanning the x - ray through radio spectrum . we have extracted additional useful information from this sample by solving the strong lensing for these clusters .
previous discoveries from this sample include the largest lensed images of a highly - magnified distant spiral galaxy ( macs j1149.5 + 2223 ; zitrin & broadhurst 2009 , smith et al .
2009 ) , the largest known lens ( macs j0717.5 + 3745 ; zitrin et al .
2009a ) , lensed sub - mm sources ( takata et al .
2003 , borys et al .
2004 , berciano alba et al .
2007 , 2009 ) , and another `` bullet cluster '' ( macs j0025.4 - 1222 ; brada et al .
2008b ) .
we presented mass models and the corresponding critical curves for the 12 high-@xmath3 x - ray luminous clusters of the sample via strong - lensing analysis in hst / acs images .
several of these clusters have only a few strongly - lensed images and corresponding low masses and einstein radii , while most are very massive and rich with multiply - lensed arcs due to large critical area and high lens power .
the equivalent einstein radii of this sample range from 11@xmath13 to 55@xmath13 with a median value of @xmath1 ( similar to the mean value ) , or 180 kpc .
the corresponding masses enclosed within these curves range from @xmath38 to @xmath39 with a median value of @xmath40 , and a mean value of @xmath41 .
we find that the enclosed mass follows tightly the quadratic relation with the equivalent einstein radius , indicating that the deviation from symmetry is not prominent in these clusters and that our measurements are highly accurate given the measured or assumed redshifts .
in addition , we have shown that @xmath47 increases as expected with the scale size of the system , or einstein radius .
we compare these results to the predictions of @xmath9cdm taking into account projection biases and find that these predictions fall short of our measured einstein radii by a factor of @xmath10 .
in addition , the standard k - s statistic gives a probability of only @xmath65 that the observed distribution is drawn from the predicted one , corresponding to a @xmath69 discrepancy .
it is apparent that the observed einstein radii and implied high masses are not likely at high redshifts in the context of the @xmath9cdm model when no lensing - selection biases are involved , showing again possible tension with the prediction of the standard model . finally , as a result of their unrelaxed mass distributions , most of these clusters cover a very large area ( @xmath70 ) of high magnification ( @xmath50 ) making them primary targets for high-@xmath3 galaxy searches , for which substantial hst time will shortly be forth - coming .    it should be acknowledged that theoretical predictions may clearly entail further uncertainties than those taken into account here , and a revised and more thorough analysis is needed in order to obtain higher significance results once source redshifts are available .
such comparisons should also be independently applied to other large samples , so that a significant number of clusters are compared to the @xmath9cdm model in order to statistically support any claimed discrepancy .
we thank the anonymous reviewer of this work for very useful comments .
az acknowledges eran ofek and salman rogers for their publicly available matlab scripts used in part of this work .
az thanks sharon sadeh for useful discussions , and elinor medezinski for sharing some subaru data .
this research is being supported by the israel science foundation ( isf ) grant 1400/10 .
rb is supported by the isf grant 823/09 .
acs was developed under nasa contract nas 5 - 32865 .
this research is based on observations provided in the hubble legacy archive which is a collaboration between the space telescope science institute ( stsci / nasa ) , the space telescope european coordinating facility ( st - ecf / esa ) and the canadian astronomy data centre ( cadc / nrc / csa ) .
part of this work is based on data collected at the subaru telescope , which is operated by the national astronomical society of japan .
barkana , r. & loeb , a. , 2010 , mnras , in press bentez , n. , et al . , 2009 ,
apj , 691 , 241 berciano alba , a. , garrett , m.a .
, koopmans , l.v.e . , wucknitz , o. , 2007 , a&a , 462 , 903 berciano alba , a. , koopmans , l.v.e . ,
garrett , m.a . ,
wucknitz , o. , limousin , m. , 2009 , arxiv:0910.1213 bertin , e. & arnouts , s. , 1996 , a&as , 117 , 393 borgani , s. , rosati , p. , tozzi , p. , norman , c. , 1999 , apj , 517 , 40 borys , c. , et al . , 2004 , mnras , 352 , 759 bouwens , r.j . , et al . , 2004 ,
apj , 616 , 79 brada , m. , et al . , 2008a , apj , 681 , 187 brada , m. , allen , s.w . , treu , t. , ebeling , h. , massey , r. , morris , r.g .
, von der linden , a. , applegate , d. , 2008b , apj , 687 , 959 brada , m. , et al . , 2006 , apj , 652 , 937 bradley , l.d . , et al . , 2008 , apj , 678 , 647 broadhurst , t. , et al .
, 2005 , apj , 621 , 53 broadhurst , t.j . &
barkana r. , 2008 , mnras , 390 , 1647 broadhurst , t , umetsu , k , medezinski , e. , oguri , m . , rephaeli , y. , 2008 , apj 685 , l9 bullock , j.s . ,
kolatt , t.s .
, sigad , y. , somerville , r.s . ,
kravtsov , a.v . , klypin , a.a . ,
primack , j.r . ,
dekel , a. 2001 , mnras , 321 , 559 clowe , d. , luppino , g.a . ,
kaiser , n. , gioia , i.m . ,
2000 , apj , 539 , 540 clowe , d. , luppino , g.a . ,
kaiser , n. , 2003 , a&a , 409 , 851 coe , d. , fuselier , e. , bentez , n. , broadhurst , t. , frye , b. , ford , h. , 2008 , apj , 681 , 814 coe , d. , bentez , n. , broadhurst , t. , moustakas , l. , ford , h. , 2010 , arxiv:1005.0398 corless , v.l . , king , l.j .
, 2007 , mnras , 380 , 149 deb , s. , goldberg , d.m . , heymans , c. , morandi , a. , 2010 , arxiv:0912.4260 donnarumma , a. , ettori , s. , meneghetti , m. , moscardini , l. , 2009 , mnras , 398 , 438 dressler , a. , gunn , j.e . , 1992 , apjs , 78 , 1 dressler , a. , smail , i. , poggianti , b.m . , butcher , h. , couch , w.j . , ellis , r.s . , oemler , a.jr . , 1999 , apjs , 122 , 51 duffy , a.r . , schaye , j. , kay , s.t . , & dalla vecchia , c. 2008 , mnras , 390 , 64 duffy , a.r . , schaye , j. , kay , s.t . , dalla vecchia , c. , battye , r.a . , booth , c.m . , 2010 , arxiv:1001.3447 ebbels , t.m.d . , le borgne , j .- f . , pello , r. , ellis , r.s . , kneib , j .-
p . , smail , i. , sanahuja , b. , 1996 , mnras , 281 , 75 ebeling , h. , edge , a.c . ,
henry , j.p .
, 2001 , apj , 553 , 668 ebeling , h. , barrett , e. , donovan , d. , 2004 , apj , 609l , 49 ebeling , h. , barrett , e. , donovan , d. , ma , c.j . ,
edge , a.c . , van speybroeck , l. , 2007 , apj , 661l , 33 ellingson , e. , yee , h.k.c . ,
abraham , r.g . ,
morris , s.l . ,
carlberg , r.g .
, 1998 , apjs , 116 , 247 ellis , s.c . , jones , l.r . , 2002 , mnras , 330 , 631 fedeli , c. , meneghetti , m. , gottloeber , s. , yepes , g. , 2010 , arxiv:1007.1551 franx , m. , illingworth , g.d . ,
kelson , d.d . , van dokkum , p.g . & tran , k. , 1997 , apj , 486 , 75 frye , b. & broadhurst , t. , 1998 , apj , 499 , 115 gao , l. , navarro , j.f . , cole , s. , frenk , c.s . ,
white , s.d.m .
, springel , v. , jenkins , a. , neto , a.f . , 2008 , mnras , 387 , 536 gavazzi , r. , fort , b. , mellier , y. , pello , r. , dantel - fort , m. , 2003 , a&a , 403 , 11 geach , j.e . , et al . , 2006 ,
apj , 649 , 661 halkola , a. , hildebrandt , h. , schrabback , t. , lombardi , m. , brada , m. , erben , t. , schneider , p. , wuttke , d. , 2008 , a&a , 481 , 65 hennawi , j.f . & spergel , d.n . , 2005 , apj , 624 , 59 hennawi j.f . , dalal n. , bode p. , ostriker j.p . , 2007 , apj , 654 , 714 hoekstra , h. , 2007 , mnras , 379 , 317 horesh , a. , maoz , d. , ebeling , h. , seidel , g. , bartelmann , m. , 2010 , mnras in press , arxiv:1004.2048 jones , t. , swinbank , m. , ellis , r. , richard , j. , stark , d. , 2009 , arxiv:0910.4488 keselman , j.a . ,
nusser , a. , peebles , p.j.e . , 2009 , arxiv:0912.4177 kneib , j.p . , mellier , y. , fort , b. , mathez , g. , 1993 , a&a , 273 , 367 kneib , j .-
, ellis , r.s . ,
smail , i. , couch , w.j . ,
sharples , r.m . , 1996 , apj , 471 , 643 kneib , j .-
p . , ellis , r.s . ,
santos , m.r .
, richard , j. , 2004 , apj , 607 , 697 laroque , s.j . , et al .
, 2003 , apj , 583 , 559 laroque , s.j . ,
bonamente , m. , carlstrom , j.e . ,
joy , m.k . ,
nagai , d. reese , e.d . ,
dawson , k.s .
, 2006 , apj , 652 , 917 lavery , r.j .
, 1996 , aj , 112 , 1812 lemze , d. , barkana , r. , broadhurst , t. , rephaeli , y. , 2008 , mnras , 386 , 1092 lemze , d. , broadhurst , t. , rephaeli , y. , barkana , r. , umetsu , k. , 2009 , apj , 701 , 1336 liesenborgs , j. , de rijcke , s. , dejonghe , h. , 2006 , mnras , 367 , 1209 liesenborgs , j. , de rijcke , s. , dejonghe , h. , bekaert , p. , 2008
, mnras , 389 , 415 liesenborgs , j. , de rijcke , s. , dejonghe , h. , bekaert , p. , 2009
, mnras , 397 , 341 limousin , m. , et al . , 2008 , a&a , 489 , 23 limousin , m. , et al . , 2009 , arxiv:0911.4125 luppino , g.a . , gioia , i.m . ,
hammer , f. , le fevre , o. , annis , j.a .
, 1999 , a&as , 136 , 117 lynds , r. & petrosian , v. , 1986 , baas , 18r , 1014 massey , r. , et al
. , 2007 , apjs , 172 , 239 medezinski , e. , broadhurst , t. , umetsu , k. , oguri , m. , rephaeli , y. , bentez , n. , 2010 , mnras , 405 , 257 meneghetti , m. , argazzi , r. , pace , f. , moscardini , l. , dolag , k. , bartelmann , m. , li , g. , oguri , m. , 2007 , a&a , 461 , 25 meneghetti , m. , rasia , e. , merten , j. , bellagamba , f. , ettori , s. , mazzotta , p. , dolag , k. , 2010 , arxiv:0912.1343 meneghetti , m. , fedeli , c. , pace , f. , gottloeber , s. , yepes , g. , 2010 , arxiv:1003.4544 merten , j. , cacciato , m. , meneghetti , m. , mignone , c. , bartelmann , m. , 2009 , a&a , 500 , 681 miyazaki , s. , hamana , t. , ellis , r.s . , kashikawa , n. , massey , r.j .
, taylor , j. , refregier , a. , 2007 , apj , 669 , 714 molnar , s.m . ,
hughes , j.p . , donahue , m. , joy , m. , 2002 , apj , 573l , 91 neto , a.f .
, et al . ,
2007 , mnras , 381 , 1450 oguri , m. , blandford , r.d . , 2009 , mnras , 392 , 930 oguri , m. , et al . , 2009 ,
apj , 699 , 1038 okabe , n. , takada , m. , umetsu , k. , futamase , t. , smith , g.p . , 2009 , arxiv:0903.1103o puchwein , e. & hilbert , s. , 2009 , 2009 , mnras , 398 , 1298 reiprich , t.h . & bhringer , h. , 1999 , an , 320 , 296 richard , j. , kneib j.p . , limousin , m. , edge . ,
a. , jullo , e. , 2009a , arxiv:0910.5553 richard , j. , et al . , 2009b , arxiv:0911.3302 ricker , p.m. & sarazin , c.l . , 2001 , apj , 561 , 621 rines , k. , geller , m.j
, diaferio , a. , kurtz , m.j . ,
jarrett , t.h .
, 2004 , aj , 128 , 1078 rosati , p. , della ceca , r. , norman , c. , giacconi , r. , 1998 , apj , 492l , 21 sadeh , s. & rephaeli , y. , 2008 , mnras , 388 , 1759 sarazin , c.l . , 1988 ,
x - ray emission from clusters of galaxies , ed .
sarazin sereno , m. , jetzer , ph . ,
lubini , m. , 2010 , arxiv:1001.1696 sharon , k. , et al . , 2005 ,
apj , 629 , 73 sheth , r.k . & tormen , g. , 1999 , mnras , 308 , 119 soucail , g. , fort , b. , mellier , y. , picat , j. p. , 1987
, a&a , 172 , 14 smail , i. , ellis , r.s . ,
fitchett , m.j . ,
edge , a.c .
, 1995 , mnras , 273 , 277 smail , i. , ellis , r.s . ,
dressler , a. , couch , w.j .
, oemler , a.jr . , sharples , r.m . , butcher , h. , 1997 , apj , 479 , 70 smith , g.p . , kneib , j .-
, smail , i. , mazzotta , p. , ebeling , h. , czoske , o. , 2005 , mnras , 359 , 417 smith , g.p . , et al . , 2009 , arxiv:0911.2003 spergel , d.n . , et al .
, 2003 , apjs , 148 , 175 spergel d.n . ,
et al . , 2007 , apjs , 170 , 377 tegmark , m. , et al . , 2004 , apj , 606 , 702 takata , t. , et al . , 2003 , pasj , 55 , 789 tanaka , m. , hoshi , t. , kodama , t. , kashikawa , n. , 2007 , mnras , 379 , 1546 taylor , a. n. , 2004 , mnras , 353 , 1176 umetsu , k. & broadhurst , t. , 2008 , apj , 684 , 177 umetsu , k. , medezinski , e. , broadhurst , t. , zitrin , a. , okabe , n. , hsieh , b.c . ,
molnar , s.m . , 2009 , arxiv:0908.0069 williams , l.l.r . , navarro , j.f . ,
bartelmann , m. , 1999 , apj , 527 , 535 wittman , d. , dellantonio , i.p .
, hughes , j.p . , margoniner , v.e . ,
tyson , j.a . ,
cohen , j.g . , norman , d. , 2006 , apj , 643 , 128 worrall , d.m . ,
birkinshaw , m. , 2003 , mnras , 340 , 1261 zackrisson , e. , et al . , 2010 , arxiv:1002.3368 zhao , d.h . , jing , y.p .
, mo , h.j . , brner , g. , 2003 , apj , 597 , l9 zhao , d.h . ,
jing , y.p .
, mo , h.j . , brner , g. , 2009 , apj , 707 , 354 zheng , w. , et al .
, 2009 , apj , 697 , 1907 zitrin , a. & broadhurst , t. , 2009 , apj , 703l , 132 zitrin , a. , broadhurst , t. , rephaeli , y. , sadeh , s. , 2009a , apj , 707l , 102 zitrin , a. , et al .
, 2009b , mnras , 396 , 1985 zitrin , a. , et al . , 2010 ,