Issue 
A&A
Volume 547, November 2012



Article Number  A67  
Number of page(s)  14  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201219944  
Published online  30 October 2012 
The strongest gravitational lenses
II. Is the large Einstein radius of MACS J0717.5+3745 in conflict with ΛCDM?
^{1} Dipartimento di AstronomiaUniversità di Bologna, via Ranzani 1, 40127 Bologna, Italy
email: jcwaizmann@oabo.inaf.it
^{2} INAF – Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
^{3} INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{4} Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik, AlbertUeberleStr. 2, 69120 Heidelberg, Germany
Received: 3 July 2012
Accepted: 7 August 2012
Context. With the amount and quality of galaxy cluster data increasing, the question arises whether or not the standard cosmological model can be questioned on the basis of a single observed extreme galaxy cluster. Usually, the word extreme refers directly to cluster mass, which is not a direct observable and thus subject to substantial uncertainty. Hence, it is desirable to extend studies of extreme clusters to direct observables, such as the Einstein radius.
Aims. We aim to evaluate the occurrence probability of the large observed Einstein radius of MACS J0717.5+3745 within the standard ΛCDM cosmology. In particular, we want to model the distribution function of the single largest Einstein radius in a given cosmological volume and to study which underlying assumptions and effects have the strongest impact on the results.
Methods. We obtain this distribution by a Monte Carlo approach, based on the semianalytic modelling of the halo population on the past lightcone. After sampling the distribution, we fit the results with the general extreme value (GEV) distribution which we use for the subsequent analysis.
Results. We find that the distribution of the maximum Einstein radius is particularly sensitive to the precise choice of the halo mass function, lens triaxiality, the inner slope of the halo density profile and the massconcentration relation. Using the distributions so obtained, we study the occurrence probability of the large Einstein radius of MACS J0717.5+3745, finding that this system is not in tension with ΛCDM. We also find that the GEV distribution can be used to fit very accurately the sampled distributions and that all of them can be described by a (typeII) Fréchet distribution.
Conclusions. With a multitude of effects that strongly influence the distribution of the single largest Einstein radius, it is more than doubtful that the standard ΛCDM cosmology can be ruled out on the basis of a single observation. If, despite the large uncertainties in the underlying assumptions, one wanted to do so, a much larger Einstein radius (≳100″) than that of MACS J0717.5+3745 would have to be observed.
Key words: gravitational lensing: strong / methods: statistical / galaxies: clusters: general / cosmology: miscellaneous / galaxies: clusters: individual: MACS J0717.5+3745
© ESO, 2012
1. Introduction
Galaxy clusters are extreme objects from many points of view. They are the most massive gravitationally bound systems in the Universe and, hence, flag the rarest peaks of the initial density field. The gas contained in their gravitational potential wells is heated up to extremely high temperatures of the order of 10^{7} − 10^{8} K, resulting in the emission of Xray radiation. Furthermore, they can give rise to spectacular events of strong gravitational lensing. Individually and as a population, galaxy clusters contain rich information on the formation of structure in the Universe that will be recovered to greater extent in the near future by ongoing and upcoming surveys, like SPT (Carlstrom et al. 2011), eROSITA (Cappelluti et al. 2011) and EUCLID (Laureijs et al. 2011).
Recently, the interest in the extremest among the extreme, the most massive clusters, has substantially increased. This development was mainly triggered by the detection of very massive galaxy clusters at high redshifts, like XMMU J2235.3 − 2557 at z = 1.4 (Mullis et al. 2005; Rosati et al. 2009; Jee et al. 2009), ACTCL J0102 at z = 0.87 (Marriage et al. 2011; Menanteau et al. 2012) and SPTCL J2106 at z = 1.132 (Foley et al. 2011; Williamson et al. 2011). Several works studied the probability to find such objects in a standard ΛCDM cosmology (Holz & Perlmutter 2012; Baldi & Pettorino 2011; Cayón et al. 2011; Hotchkiss 2011; Mortonson et al. 2011; Chongchitnan & Silk 2012; Waizmann et al. 2012a). All these studies focused on the mass of galaxy clusters, which is unfortunately not a direct observable. The mass of a galaxy cluster, ill defined in the first place, is subject to substantial scatter and biases. Hence, it is desirable to study extremes in direct, better defined, observables, such as strong lensing signals.
A particular interesting case from this point of view is the extremely large critical curve of the Xray luminous galaxy cluster MACS J0717.5+3745 at redshift z = 0.546, which has been independently detected by the Massive Cluster Survey (MACS) (Ebeling et al. 2001, 2007) and as a host of a diffuse radio source (Edge et al. 2003). A stronglensing analysis revealed that the effective Einstein radius, with θ_{eff} = (55 ± 3)″ for an estimated source redshift of z ≃ 2.5, is the largest known at redshifts of z > 0.5 (Zitrin et al. 2009, 2011). It is unclear whether or not such a large Einstein radius is consistent with the ΛCDM cosmology (Zitrin et al. 2009).
In this work, we study the distributions of the single largest Einstein radius by means of semianalytic modelling of the halo distribution on the past lightcone, as introduced in Redlich et al. (2012). On the basis of the study of Oguri et al. (2003), the haloes are modelled using triaxial density profiles. By Monte Carlo (MC) sampling the distribution of the maximum for different underlying assumptions like the mass function, the allowed range of triaxiality, the inner slope and the massconcentration relation, we study the impact of different choices on the resulting distributions of the maximum. We use the results to assess the occurrence probability of the Einstein radius of MACS J0717.5+3745 in the redshift range of 0.5 ≤ z ≤ 1.0 and fit the generalized extreme value (GEV) distribution to the sampled distribution of the maximum Einstein radius.
This paper is structured as follows. In , we briefly review the basics of strong cluster lensing followed by an introduction of the semianalytic modelling of the distribution of Einstein radius in . In , we introduce extreme values statistics as far as it is relevant for the presented work before studying in detail the distribution of the largest Einstein radius for different underlying physical assumptions in . Afterwards, we perform a case study for MACS J0717.5+3745 in and close with a summary and the conclusions in .
Throughout this work we adopt the Wilkinson Microwave Anisotropy Probe 7year (WMAP7) parameters (Ω_{m0}, Ω_{Λ0}, Ω_{b0}, h, σ_{8}) = (0.727, 0.273, 0.0455, 0.704, 0.811) (Komatsu et al. 2011).
2. Strong lensing by galaxy clusters
The theory of gravitational lensing is well established and has become a very important tool for the study of the dark components of the Universe (for recent reviews, see e.g. Bartelmann 2010; Kneib & Natarajan 2011). In this work, we focus on strong lensing by galaxy clusters and in particular on the distribution of Einstein radii.
Einstein radii measure the size of the tangential critical curve, defined as the curve of vanishing tangential eigenvalue (1)with κ and γ denoting the convergence and the shear, respectively. In the axially symmetric case, the shear can be written as (2)where denotes the mean convergence within a circle of radius θ. One can then define the Einstein radius as the radius of a circle enclosing a mean convergence of unity (3)where θ_{E} denotes the Einstein radius. Recalling the relation between the convergence and the surface mass density κ = Σ / Σ_{crit}, one can formulate the definition of θ_{E} as (4)stating that the mean surface density within the Einstein radius equals the critical surface mass density Σ_{crit}.
The assumption of axially symmetric lenses is unsustainable for realistic lenses, and thus several definitions of the Einstein radius for the case of arbitrary lenses with irregular tangential critical curves exist.
Fig. 1 Distribution in mass and redshift of 4000 maxima of the effective (left panel) and the median Einstein radius (right panel) in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky based on the Tinker et al. (2008) mass function. The colour encodes the size of the individual largest Einstein radius from each simulation run. The green circles denote exemplary systems for which the two Einstein radius definitions selected different haloes as the maximum. 

Open with DEXTER 
In this work, we focus on two definitions: the first, introduced by Meneghetti et al. (2011) and known as the median Einstein radius, is defined as (5)where θ_{t} denotes the set of the tangential critical points and θ_{c} is the centre of the lens. The second definition is of geometrical nature and is usually referred to as effective Einstein radius. It is defined as (6)where A is the area enclosed by the critical curve. Of course, both definitions are identical to the original definition of θ_{E} in the axially symmetric case. In the next section, we will discuss in detail how the distributions of these two quantities can be modelled in a cosmological context.
3. Semianalytic modelling of the distribution of Einstein radii
For the semianalytic modelling of the distribution of Einstein radii, we closely follow the work of Redlich et al. (2012), which provides a more detailed discussion of the algorithm, including the computation of the lensing signal and the cosmological population of haloes. Therefore, we will in this section only briefly summarise those aspects of the semianalytic modelling that will be needed to follow the remainder of this work. However, in our paper, we conservatively neglect the impact of mergers on the extreme value distribution. The results of Redlich et al. (2012) indicate that the inclusion of mergers can be expected to significantly shift the distribution to larger Einstein radii. Thus, the results presented in our paper can be considered as conservative estimates.
3.1. Modelling triaxial haloes
The integral part for the modelling of the Einstein radius distribution is the inclusion of triaxiality as discussed, for instance, in Oguri & Blandford (2009, hereafter OB09), which is based on the work of Jing & Suto (2002, hereafter JS02). In their work, JS02 generalised the Navarro et al. (1996, hereafter NFW) profile to a triaxial model, where the axis ratios for a given mass, M, and redshift, z, can be sampled from the following empirically derived probability density functions, assuming the ordering \arraycolsep1.75ptwhere the latter relation holds for a / b ≥ r_{min} and is zero otherwise and (9)Here, M_{ ∗ } is the characteristic nonlinear mass scale and, according to JS02, the bestfitting parameter for the width of the axisratio distribution p(a / c) is σ_{s} = 0.113.
The concentration parameter c_{e} is defined as c_{e} ≡ R_{e} / R_{0}, where R_{e} is determined such that the mean density within the ellipsoid of the major axis radius R_{e} is Δ_{e}Ω(z)ρ_{crit}(z), with (10)where Δ_{vir}(z) is the overdensity of objects virialized at redshift z, which we approximate according to Nakamura & Suto (1997). In their work, JS02 found a lognormal distribution for the concentration, (11)with a dispersion of σ_{c} = 0.3. Following Oguri et al. (2003), we include a correlation between the axis ratio a / c and the mean concentration, where z_{c} is the collapse redshift. In order to avoid unrealistically small concentrations for particularly small axis ratios a_{sc}, we use the correction introduced by OB09, forcing f_{c} ≥ 0.3 in Eq. (13). Obviously, triaxial haloes with particularly small axis ratios a_{sc} (and hence also small concentrations c_{e}) are highly elongated objects whose lens potential is dominated by masses well outside the virial radius (see e.g. Oguri & Keeton 2004). There are two ways of dealing with the problem of these unrealistic scenarios. The first is based on truncating the density profile beyond the virial radius (see e.g. Baltz et al. 2009; OB09), the second approach suppresses particularly small axis ratios, a_{sc}, from the tail of the underlying axis ratio distribution (see ). Following JS02, we set the free parameter A_{e} = 1.1 for a standard ΛCDM model, unless stated otherwise. The expressions listed so far are valid for an inner slope of the density profile of α_{NFW} = 1.0. For the case of α_{NFW} = 1.5, we use the simple relation (Keeton & Madau 2001; JS02).
3.2. Preparatory considerations for the MC sampling
In order to obtain the extreme value distribution of the maxima of the Einstein radii, we have to sample the cluster populations for many mock realisations and to collect the largest Einstein radius from each realisation. In view of this, the choice of the redshift interval and the allowed halo mass range is decisive to keep computational costs under control.
Fig. 2 CDFs of the largest effective (left panels) and median (right panels) Einstein radius for a different number of maxima in the range between 125 and 4000, assuming the Tinker et al. (2008) mass function, M_{lim} = 2×10^{14} M_{⊙} / h and the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. The upper row shows the CDFs directly based on the MC simulations, and the lower row presents the corresponding fits of the GEV distribution. The small lower panels show the difference, ΔCDF, with respect to the highresolution run based on 4000 maxima. 

Open with DEXTER 
On the basis of the stronglensing analysis of 12 MACS clusters in the interval 0.5 ≤ z ≤ 0.6 by Zitrin et al. (2011), we focus our study of the distribution of the largest Einstein radii on clusters in the redshift range of 0.5 ≤ z ≤ 1, assuming a source redshift of z_{s} = 2.0. This choice already drastically reduces the number of haloes that have to be simulated. The remaining task is to identify a lower mass limit M_{lim}, such that the inferred sampled maxima distribution is not biased. To do so, we simulated 4000 maxima with M > 2×10^{14} M_{⊙} / h in 0.5 ≤ z ≤ 1 and present the results in for both the effective and the median Einstein radii. It can be directly inferred from the distribution of the maxima that M_{lim} = 2×10^{14} M_{⊙} / h is a sufficient lower mass limit, confirming the results of OB09. Thus, we will adopt this value throughout this work, unless stated otherwise.
The distribution of the maxima in mass and redshift, presented in , shows that the maxima stem from a wide range of masses. It is not unlikely that a rather low mass cluster gives rise to the largest Einstein radius. The fact that most of the maxima are found in the lower redshift range is a consequence of the selected lensing geometry determined by the choice of the source redshift z_{s} = 2.0. Since we are modelling triaxial haloes, this is a first indication that the orientation of the halo with respect to the observer, the lensing geometry and the concentration are more important than the mass. The green circles in denote systems for which the Einstein radius definitions select different haloes to exhibit the maximum Einstein radius. We note that the extremely large Einstein radii (black dots) are not affected by the choice of the Einstein radius definition.
Having fixed the mass and redshift range, the last step in optimizing the computational cost is to understand how many maxima actually have to be sampled in order to construct the cumulative distribution function (CDF) of the largest effective and median Einstein radii. To this end, we computed the respective CDFs for different sample sizes N_{samp} between 125 and 4000. Each CDF was computed at 50 linearly equidistant points between the largest and the smallest value. The results of these computations are presented in the upper row of , where the CDFs themselves are shown in the upper panels and the differences of the CDFs with respect to the highest resolution run (N_{samp} = 4000) are shown in the small lower panels. As expected, the noise of the CDFs decreases with increasing N_{samp}. For N_{samp} ≥ 1000, the difference with respect to the high resolution run is ≲ 0.02, corresponding to an overestimation of the occurrence probability of a given Einstein radius by less than two per cent. Hence, we will utilise N_{samp} = 1000 for all of our computations in the remainder of this work, unless stated otherwise.
4. Applying extreme value statistics to the distribution of the largest Einstein radii
Extreme value statistics (EVS) (for an introduction, see e.g. Gumbel 1958; Kotz & Nadarajah 2000; Coles 2001; Reiss & Thomas 2007) models the stochastic behaviour of the extremes and tries to give a quantitative answer to the question of how frequent unusual observations are. In the framework of EVS, there are two approaches to the modelling of rare events. The first one, also known as the Gnedenko approach (Fisher & Tippett 1928; Gnedenko 1943), models the distribution of the block maxima, while the second one, known as the Pareto approach (Pickands 1975), models the distribution of excesses over high thresholds. Since we are interested in the study of the distribution of the largest Einstein radii, we will discuss the first approach in more detail in the following.
Values of the location, scale and shape parameters, α_{eff}, β_{eff} and γ_{eff}, of the fitted GEV distributions of the maximum effective Einstein radius for different sample sizes N_{samp}, as shown in the lower left panel of .
4.1. The Gnedenko approach
This approach is concerned with the modelling of the block maxima M_{n} of independently identically distributed (i.i.d.) random variables X_{i}, which are defined as (14)It has been shown (Fisher & Tippett 1928; Gnedenko 1943) that, for n → ∞, the limiting CDF of the renormalised block maxima is given by one of the extreme value families: Gumbel (Type I), Fréchet (Type II) or Weibull (Type III). These three families can be unified (von Mises 1954; Jenkinson 1955) as a general extreme value (GEV) distribution (15)with the shape, scale and location parameters γ, β and α. In this generalisation, γ = 0 corresponds to the Type I, γ > 0 to Type II and γ < 0 to the Type III distributions. The corresponding probability density function (PDF) is given by g_{γ, β, α}(x) = dG_{γ, β, α}(x) / dx and reads for the case of γ ≠ 0, (16)From now on we will adopt the convention that capital initial letters denote the CDF (like G_{γ, β, α}(x)) and small initial letters denote the PDF (like g_{γ, β, α}(x)). The mode, the most likely value, of the GEV distribution reads (17)and the expected value is given by (18)where Γ denotes the Gamma function.
4.2. GEV and the distribution of Einstein radii
The MC approach to the distribution of the single largest Einstein radius, introduced in , provides simulated CDFs that are discrete by nature. Due to the complexity of the modelling of the Einstein radius distribution, it is not possible to find analytic relations for the GEV parameters, as can be done for halo masses (Davis et al. 2011; Waizmann et al. 2011). Hence, we use the limiting GEV distribution from Eq. (15) to fit the sampled distributions in order obtain analytic relations for the distribution of the largest Einstein radii.
Fig. 3 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and the median (right panels) Einstein radius for the mass functions of PressSchechter, Tinker, Sheth & Tormen and Crocce as labelled in the panels. All distributions are based on the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER 
The GEV distribution, given by Eq. (15), fits the MCsimulated distributions very well, as can be inferred from the lower row of . In the small panels below, we show the difference ΔCDF of the CDFs with respect to the fits based on the run with N_{samp} = 4000. It can be seen that the fitted functions deviate less from the highresolution run with respect to the MC data. The blue, dasheddotted line depicts the deviation of the fit from the corresponding data set for N_{samp} = 4000, showing that the GEVbased fits are capable of describing the MC simulated distributions of the Einstein radii very well. We also present in all fitted GEV parameters as well as the root mean square of the residuals for the case of the effective Einstein radius. In what follows, we will use the GEV fits for any subsequent analysis, like the calculation of PDFs, modes or quantiles.
The fitted shape parameters γ for all distributions discussed in this work are found to be in the range of 0.05 < γ < 0.2, which means that the distribution of the largest Einstein radii can in general be described by a Fréchet (Type II) distribution, indicating that the distribution is bounded from below. An exception to this will be discussed further in . The location parameter α is always very close to the mode, the most likely maximum, with the two values differing only by roughly one per cent. It is noteworthy that the location parameter α can be estimated very well with rather small sample sizes, whereas the shape parameter γ is subject to larger uncertainties. Even for only 125 samples, the difference of the mode with respect to the N_{samp} = 4000 case is less than 2 arcsec. This result is similar to the findings of Waizmann et al. (2011), who report the same behaviour for the case of halo masses.
5. The distribution of the largest Einstein radii
In this section, we study the impact of several underlying assumptions, like mass function, triaxiality, inner slope and massconcentration relation, on the distribution of the single largest Einstein radius.
5.1. The impact of the mass function
The choice of the mass function must have an effect on the distribution of the largest Einstein radii, since it alters the size of the halo population from which the maxima are drawn. This effect is particularly important for galaxy clusters since the exponentially suppressed tail of the mass function is naturally very sensitive to modification. As shown in , the maxima stem from a relatively broad range of masses. Hence, the larger the halo population in this mass range, the more likely it is to sample a particularly large Einstein radius.
In order to quantify the influence of different mass functions, we sampled the CDFs of the largest Einstein radii for four different mass functions, the Press & Schechter (1974, PS), Tinker et al. (2008), Sheth & Tormen (1999, ST) and the Crocce et al. (2010) mass functions. We decided to use the Tinker mass function as a reference because the halo masses are defined as spherical overdensities with respect to the mean background density, a definition that is closer to theory and actual observations than the friendoffriend masses that were used for the Crocce mass function. We added the Crocce et al. (2010) mass function to our analysis because it differs significantly at the highmass end from other simulations (Bhattacharya et al. 2011). The mass function is based on simulations with a box size much larger than the horizon scale, which gives more statistics at the highmass end at the price of leaving the realm of the Newtonian approximation. However, Green & Wald (2012) argue that the Newtonian approximation for Nbody simulations might also be valid on superhorizon scales.
The resulting extreme value distributions for the four different mass functions are presented in , on the basis of the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. The results show that the effect of different mass functions is substantial. The CDFs based on the ST and the Crocce mass functions exhibit the strongest difference with respect to the PS one. The Tinker results lie in between the two. The modes of the distributions, listed in , can differ by more than 10 arcsec for different choices of the mass function, which implies that the inferred occurrence probabilities for a given Einstein radius can be substantially different. From the shape parameters γ, it can be inferred that the different mass functions do not strongly affect the shape of the distributions. This is also confirmed by the fact that the CDFs shown in the upper row of seem just to be shifted in the θdirection.
Shape parameter γ and the mode of the CDFs of the largest effective and median Einstein radius for four different mass functions.
The distributions clearly reflect the different behaviour of the mass functions with mass, as shown in in the form of the ratio of the number of haloes more massive than M with respect to the Tinker case. It can be seen that the ST and Crocce mass functions lead to a substantial increase in the number of haloes, particularly at the highmass end, whereas the PS mass function results in far fewer haloes in the redshift interval of interest.
Considering the strong impact of the mass function on the extreme value distribution, it is certainly necessary to improve the accuracy of the mass function at the highmass end.
Fig. 4 Ratio of the number of haloes N( > M) for the different mass functions as labelled in the panel, with respect to the Tinker one, assuming a redshift interval of 0.5 ≤ z ≤ 1.0 and the full sky. 

Open with DEXTER 
5.2. The impact of triaxiality
The triaxiality of the lensing haloes has a substantial impact on the distribution of the maxima. For instance, a very elongated halo that is directed along the line of sight can lead to a highly concentrated, projected surface massdensity profile, which causes a large tangential critical curve (see e.g. Oguri et al. 2003; Dalal et al. 2004; Meneghetti et al. 2007, 2010). When sampling axis ratios from Eqs. (7) and (8), particularly small values of the sampled axis ratios will potentially propagate into extreme stronglensing events. Since the empirical fits from JS02 are only based on few data points in this regime, it is unclear how reliable the fitted axisratio distributions are.
In order to study the impact of this uncertainty, we introduced a cutoff in the distribution from Eq. (7) to remove extreme axis ratios. We cut off the distribution of the scaled axis ratios at different confidence levels n, according to (19)selecting values of 1σ, 2σ, 3σ and comparing them to the case without any cutoff. To do so, we MCsimulate the distributions of the largest Einstein radius for the different cutoffs based on 1000 maxima, assuming M_{lim} = 2×10^{14} M_{⊙} / h and the Tinker mass function. The results are presented in , which shows that the impact of the tail of the axis ratio distribution on the extreme value distribution is substantial. In comparison to the impact of the different mass functions that lead to a shift of the CDF, the cutoff value strongly affects the mode as well as the shape of the CDFs, as can be inferred from the shape parameters listed in . For the 1σ cutoff, the shape parameter becomes negative. Consequently, the CDF of the largest Einstein radius is then described by a Weibull (Type III) distribution, indicating that the underlying distribution is bounded from above. For decreasing values of the cutoff (less extreme axis ratios), the CDF steepens, which corresponds to the fact that a given observed Einstein radius appears to be less likely to exist.
In order to understand better what impact the triaxiality has on the sample of the largest Einstein radii, we study the distribution of the sampled haloes in scaled axis ratio a_{sc} and concentration c_{e} based on different selection criteria. We are interested in the question whether or not the largest Einstein radii always stem from very extreme axis ratios. In the lefthand panel of , we show the distribution of the full halo sample of a single realisation (black dots) together with the distribution of the minima of a_{sc} based on 1000 allsky realisations. As expected, the highly elongated haloes exhibit small values of the concentration parameter and typical values for the minima scatter around . Now we compare this distribution to the one based on the largest effective Einstein radii shown in the centre panel of . It is evident that the largest Einstein radii stem by no means exclusively from the haloes with a minimum of a_{sc}, but from a rather broad range of a_{sc}. Thus, the largest Einstein radii either stem from lowly concentrated very elongated haloes or from less elongated but higher concentrated ones. However, the largest maxima (the dark red to black dots in the centre panel of ) are confined to the region of small a_{sc} and small c_{e}. This fact explains the strong impact of the cutoff in the scaled axisratio distribution on the shape of the CDF of maxima because the cutoff effectively removes the largest maxima.
Due to the limited knowledge of the statistics of extremely small axis ratios, it is not possible to clearly define a proper choice of the cutoff (if present) until the triaxiality distributions of large halo samples (covering the largest cluster masses) are studied in numerical simulations. In the study of JS02 (see their Fig. 9), scaled axis ratios below ~ 0.2 were not found for any of the studied redshifts. The value of corresponds to the cutoff on the 3σ level. It should be noted that in general a_{sc} also depends on the underlying cosmology, as can be seen from Eq. (9). Adopting such a cutoff value would mean that the resulting CDFs of the maxima are very close to the one without any cutoff, as can be seen from . Of course, indicates that the impact of triaxiality should always be discussed together with that of the concentration.
Fig. 5 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and median (right panels) Einstein radius for cutoffs of the axisratio distribution as labelled in the figure. All distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER 
5.3. The impact of the inner slope and the c–M relation
In order to study the effect of the inner slope α_{NFW} of the density profile, we sampled the distribution of the maxima for two different values, α_{NFW} = (1.0,1.5), using the approach discussed at the end of . The resulting distributions are presented in for α_{NFW} = 1.0 (black solid line) and α_{NFW} = 1.5 (green dasheddotted line). The distribution for the steeper inner density profile is shifted to smaller Einstein radii, confirming the findings reported in Oguri & Keeton (2004) on p. 669. On average, steeper density profiles lead to slightly larger Einstein radii (see e.g. Oguri 2004). However, the distribution of the maxima is particularly sensitive to the triaxiality and the orientation of the halo along the line of sight. For the corresponding small axis ratios, the shallower density profile contributes stronger due to the projection effect^{1} if a very elongated, extended halo is aligned along the line of sight.
In addition to the previously discussed effects, the c–M relation naturally impacts significantly on the distribution of the largest Einstein radii. In order to study the impact, we mimic a variation in the c–M relation by computing the distributions for different values of the normalisation parameter, A_{e}, in Eq. (12) for the mean concentration. Here, assuming a fixed a_{sc}, smaller A_{e} correspond to smaller values of the mean concentration, and thus it is more likely that haloes have a smaller c_{e}. Following Oguri et al. (2003), we vary the value of A_{e} between 0.8 and 1.6, with A_{e} = 1.1 being the ΛCDM standard value, and present the resulting distributions in . It can be seen that the distributions with a larger value of A_{e} are shifted to smaller Einstein radii with respect to the standard choice of A_{e} = 1.1. Vice versa, the distribution for A_{e} = 0.8 is shifted to larger values, confirming the findings from the previous section and that the largest maxima stem from haloes with small a_{sc} and c_{e}. Thus, lowering the mean concentration with a smaller value of A_{e} naturally results in shifting the CDF of the maxima to larger values. In this sense, the results presented in reflect the high sensitivity of the maxima to particularly small axis ratios, which are connected to small concentrations (larger scale radii) and, hence, more extended haloes. When aligned along the line of sight, such a system will give rise to a particularly large Einstein radius. In fact, the righthand panel of shows that the vast majority of maxima of the Einstein radius are very well aligned along the line of sight.
Fig. 6 The lefthand and centre panels show the distribution in scaled axis ratio a_{sc} and concentration c_{e} of sampled haloes, according to different selection criteria. The lefthand panel shows the distribution of a full halo sample of a single realisation (black dots) and the sample of the minima in a_{sc} from 1000 realisations (red dots). The centre panel shows the distribution of the haloes that give rise to the largest effective Einstein ring (encoded in the colourbar) based on 1000 realisations. The small arrows denote the indicated cutoffs identical to . The righthand panel shows the distribution of the same maxima with respect to the scaled axis ratio and the alignment  cosθ  , where θ is the angle between the major axis of each halo with respect to the line of sight. 

Open with DEXTER 
Shape parameter γ and the mode of the CDFs for the largest effective and median Einstein radius for different cutoffs in the axisratio distribution.
Fig. 7 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and median (right panels) Einstein radius for different values of the inner slope of the density profile α_{NFW} = (1.0,1.5) as labelled in the figure. All distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER 
With the results of this section, we mainly want to emphasize that there are many uncertainties having a strong impact on the statistics of Einstein radii. In particular, we would like to underline the need to improve the statistics of triaxiality for extreme axis ratios and the closely related uncertainties stemming from projection effects.
Fig. 8 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and median (right panels) Einstein radius for different values of the normalisation parameter A_{e} of the massconcentration relation as labelled in the figure. All distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER 
5.4. Other important effects
The most important effect that will strongly influence the distribution of the maxima is the inclusion of dynamical mergers as discussed in Redlich et al. (2012). In particular, mergers perpendicular to the line of sight can cause a strong elongation of the critical curve. Closely related and to some extent equivalent, the inclusion of substructures is expected to have a similar effect. Also, the brightest cluster galaxy (BCG) can lead to an increase of the stronglensing effect (Puchwein et al. 2005), but to a lesser extent than the previously discussed effects. We decided to omit mergers in this study, since their inclusion is computationally expensive and, as the next section shows, they are not required to explain the presence of the largest known Einstein radius.
6. MACS J0717.5+3745 – a case study
After the statistical study of the extreme value distributions of the largest Einstein radius and the different effects that influence them, we use the inferred distribution to assess the occurrence probability of the Einstein radius of MACS J0717.5+3745, which has the largest currently known observed critical curve. In addition, we also discuss in the second part of this section, what mass MACS J0717.5+3745 would need to have in order to be considered to be in significant tension with the standard ΛCDM model.
6.1. The cluster MACS J0717.5+3745
The Xray luminous galaxy cluster MACS J0717.5+3745, independently observed at redshift z_{obs} = 0.546 by the MACS survey (Ebeling et al. 2001, 2007) and as a host of a diffuse radio source by Edge et al. (2003), is a quite remarkable system. The cluster is connected to a 4 Mpc long largescale filament (Ebeling et al. 2004) that leads into the cluster and exhibits ongoing merging activity (Ma et al. 2008). Furthermore, this cluster possesses the most powerful known radio halo (Bonafede et al. 2009) and has also been observed (LaRoque et al. 2003; Mroczkowski et al. 2012) via the Sunyaev Zeldovich effect (Sunyaev & Zeldovich 1972, 1980). A stronglensing analysis (Zitrin et al. 2009, 2011) of this highly interesting system revealed that, with θ_{eff} = 55 ± 3 arcsec for an estimated source redshift of z_{s} ~ 2.5, the effective Einstein radius is the largest known at redshifts z > 0.5. Also, the mass enclosed by this critical curve is with (7.4 ± 0.5)×10^{14} M_{⊙} very large, indicating that this cluster might qualify to be among the most massive known galaxy clusters at z > 0.5. The recent stronglensing analysis by Limousin et al. (2012) reports a higher redshift of z_{s} ~ 2.96 for the primary lensed system, which would lower the size of the Einstein radius as well as the overall mass estimate.
6.2. The probability of occurrence of the critical curve
In order to assess the occurrence probability of the observed effective Einstein radius of MACS J0717.5+3745, we compute the CDF of the largest Einstein radius for a redshift interval of 0.5 ≤ z ≤ 1.0, based on the Tinker mass function, a source redshift of z_{s} = 2.5 for the nominal MACS survey area (A_{s} = 22 735 deg^{2}), and the full sky (A_{s} = 41 153 deg^{2}). We decided to use the conservative cutoff , corresponding to the inclusion of 99 per cent of the possible axis ratios from the distribution in Eq. (7). In doing so, we cut off the distribution above the most likely minimum of . Thus the CDF will be steeper, resulting in a more conservative estimate of the occurrence probability of a given Einstein radius. In comparison to the previously discussed distributions that assumed a z_{s} = 2.0, the higher source redshift will shift the distribution to larger Einstein radii due to the modified lensing geometry, as discussed in Oguri & Blandford (2009).
Fig. 9 CDFs of the largest effective (left panel) and median (right panel) Einstein radius, assuming a redshift interval of 0.5 ≤ z ≤ 1.0, a source redshift of z_{s} = 2.5 and the nominal MACS survey area (upper panels) as well as the full sky (lower panels). Both distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h on the full sky. The dasheddotted lines, together with the yellow shaded area, illustrate the impact of the uncertainty in the WMAP7 value of σ_{8} on the CDFs and the grey shaded area denotes the uncertainty in the measurement of the Einstein radius. 

Open with DEXTER 
Like mass function, the uncertainty in the normalisation of the matter power spectrum, σ_{8}, will also influence the distribution of the largest Einstein radius, since it influences the number of haloes from which the maxima will be drawn. In order to incorporate the uncertainty in the measured σ_{8}, we also computed the distributions for the upper and lower 1σ limits σ_{8} = 0.811 ± 0.023 of the WMAP7+BAO+SNSALT dataset (Komatsu et al. 2011). The results are shown in for the effective (lefthand panel) and the median (righthand panel) Einstein radius, where the dasheddotted lines show the CDFs for the upper and lower allowed value of σ_{8}. The upper and the lower panels of show the distributions for the nominal MACS survey area and for the full sky, respectively. The grey shaded area illustrates the uncertainty due to the accuracy of the measurement of the Einstein radius itself and the yellow shaded region depicts the uncertainty due to the allowed range of σ_{8}. For the MACS survey area and the effective Einstein radius θ_{eff}, we find an occurrence probability of ~(16−32) per cent based on the uncertainty in θ_{eff} alone. When the additional uncertainty of the precision of σ_{8} is included, this range widens to ~(11−42) per cent. The CDF for the median Einstein radius leads to smaller occurrence probabilities of ~(4−9) per cent and ~(3−12) per cent; however, the median Einstein radius is more sensitive to the individual structure of the system. Thus, we decide to base our study on the more conservative choice of the effective Einstein radius. When the survey area is extended to the full sky, the CDFs are shifted to larger values of the Einstein radius. As a result, the occurrence probabilities for a given observation will increase. In the case of MACS J0717.5+3745, we find ~(18−61) per cent for the effective Einstein radius and ~(4−20) per cent for median Einstein radius when taking the uncertainty of σ_{8} into account.
From the ranges of occurrence probabilities, it can be directly inferred that the large critical curve of MACS J0717.5+ 3745 cannot be considered in tension with the ΛCDM model. This finding is supported by the results of the previous sections, which showed that the uncertainty of the mass function, particularly at the highmass end, and the uncertainty in the shape of galaxy clusters allow a wide range of distributions. The recently reported higher redshift for the primary lensed system (Limousin et al. 2012) and hence a smaller inferred Einstein radius would further strengthen our conclusions. Because MACS J0717.5+3745 is one of the most dynamically active known clusters and because Redlich et al. (2012) show that the distributions of the largest Einstein radius will be significantly shifted to larger values if dynamical mergers are accounted for, it can with certainty be deduced that the large Einstein radius of MACS J0717.5+3745 is consistent with the standard ΛCDM model.
With this in mind, these inferred occurrence probabilities should just be considered a rough estimate. Because of the uncertainties in modelling the distribution of Einstein radii, an observed critical curve should exhibit a much larger extent in order to be taken as clearly in tension with the ΛCDM model. To quantify this statement, we take the CDF of θ_{eff} for MACS J0717.5+3745 and calculate the values of θ_{eff} for which the CDF takes values that correspond to confidence levels nσ with n ∈ [ 1,5] , as shown in . In order to lie beyond the 3σ level, corresponding to an occurrence probability of ~0.3 per cent, and to account for the uncertainty stemming from σ_{8}, θ_{eff} should be larger than ~115 arcsec. This number will also be strongly affected by the redshift of the source population in the sense that sources at lower (higher) redshifts will result in a smaller (larger) value of . Of course, the inclusion of dynamical mergers will increase this limit to even larger values.
Fig. 10 CDF of the largest effective Einstein radius, assuming a redshift interval of 0.5 ≤ z ≤ 1.0, the full sky and a source redshift of z_{s} = 2.5. The dasheddotted lines illustrate the uncertainty in the CDF due to the imprecision of the WMAP7 value of σ_{8}. The vertical shaded regions denote the values of θ_{eff} for which the CDF takes a value equal to nσ with n ∈ [ 1,5] , as labelled in the panel. The red filled circle denotes the measured value of θ_{eff} for MACS J0717.5 and the short red arrow points at the corresponding value of the CDF. 

Open with DEXTER 
6.3. The probability of occurrence of the mass
Apart from the larger critical curve, MACS J0717.5+3745 is also considered to be one of the most massive clusters in the redshift range of 0.5 ≤ z ≤ 1.0. The stronglensing analysis by Zitrin et al. (2009) revealed that the mass enclosed by the critical curve is (7.4 ± 0.5)×10^{14} M_{⊙} and that the mass enclosed within the larger critical curve for the multiplylensed dropout galaxy at z ~ 4 is found to be ~1×10^{15} M_{⊙}. These values are the masses in the innermost regions of the cluster, and thus the total mass can be considered to be a multiple of these values. Therefore, it is also interesting to study the expectation for the most massive galaxy cluster in the redshift range of interest of 0.5 ≤ z ≤ 1.0.
Fig. 11 CDF of the most massive halo that is expected to be found in the redshift interval of 0.5 ≤ z ≤ 1.0 in an survey area of A_{s} = 22 735 deg^{2}. The dasheddotted lines illustrate the uncertainty in the CDF due to the imprecision of the WMAP7 value of σ_{8}. The vertical shaded regions denote the values of the mass for which the CDF takes a value equal to nσ with n ∈ [ 1,5] , as labelled in the panel. We also added for reference the mass of ACTCL J0102 (red errorbar), the most massive known cluster in the given redshift range, where the two small red arrows point at the values of the CDF stemming from the upper allowed mass for the distributions based on σ_{8} (lower arrow) and σ_{8} − Δσ_{8} (upper arrow). 

Open with DEXTER 
EVS can also be applied to study the distribution of the most massive halo in a given volume (Chongchitnan & Silk 2012; Davis et al. 2011; Harrison & Coles 2011, 2012; Waizmann et al. 2011, 2012a,b); we will follow the procedure shown in Waizmann et al. (2012a) to compute the distribution function. The results are presented in , where the CDF of the most massive halo in 0.5 ≤ z ≤ 1.0 is shown for the nominal MACS survey area of A_{s} = 22 735 deg^{2}. For reference, we added ACTCL J01024915 (Menanteau et al. 2012), which is currently the most massive known cluster in the range of 0.5 ≤ z ≤ 1.0. The red filled circle with errorbars illustrates the mass and the upper/lower mass limits after the correction for the Eddington bias (Eddington 1913) according to Mortonson et al. (2011). The two small red arrows point at the values of the CDF stemming from the upper allowed mass for the distributions based on σ_{8} (lower arrow) and σ_{8} − Δσ_{8} (upper arrow). As already discussed in Waizmann et al. (2012a), this cluster can be considered to be in agreement with ΛCDM. The vertical, blue shaded areas indicate the values of the maximum halo mass for which the CDF corresponds to the values equal to the labelled confidence level. In order to let MACS J0717.5+3745 be in tension with ΛCDM, its mass should lie at least above the value indicated by the 3σ area, which would correspond to an occurrence probability of ≲ 0.3 per cent.
Therefore, the mass that should be at least exceeded is M_{3σ} ≃ 4 × 10^{15} M_{⊙}, which is per se already a very high mass for a cluster and, since cluster masses are subject to significant uncertainties, the lower mass limit should lie above this value. Furthermore, for a statistical analysis similar to the one of ACTCL J01024915, one also has to shift the observed mass to a smaller value to account for the Eddington bias (Mortonson et al. 2011).
Considering the complex dynamical state of MACS J0717.5+3745, the embedding in a large scale filament and the very high mass that would be required, it seems to be more than doubtful that, despite the large mass enclosed in the critical curve (Zitrin et al. 2009), the total mass of this system can be used to exclude ΛCDM. When summing up our findings for both the Einstein radius and the mass, we find that the characteristics of MACS J0717.5+3745 are not unlikely to be found in a ΛCDM framework. These results are substantially different from the findings of Zitrin et al. (2009), who report that the probability to find such a system is of the order of ~10^{7}. The main reason for this difference is that in order to sample the distribution of the maxima, a larger number of universes ( ~ 1000) have to be simulated, which is feasible with a semianalytic approach but where a Nbody based approach (Broadhurst & Barkana 2008) falls short. The latter can be used to infer the statistical characteristic of large Einstein radii in general but not of the single largest observation. Furthermore, one biases the results by a posteriori choosing the redshift interval and the assigned Einstein radius, since it can not be known before at which redshift the most massive cluster will be realised and what Einstein radius it will have (see e.g. Hotchkiss 2011; Waizmann et al. 2012a).
7. Summary and conclusions
In this work, we have presented a study of the distribution of the single largest Einstein radius at redshifts 0.5 ≤ z ≤ 1.0, based on the MC simulation of triaxial halo populations in mock universes extending the work of OB09. The details of the implementation of our semianalytic method can be found in Redlich et al. (2012). Our work can be divided into three distinct parts: first, a preparatory study; second, a study of the effects that impact on the distribution of the maximum Einstein radius; and third, a case study for MACS J0717.5+3745.
In the first preparatory part, we showed that ~1000 mock universes are sufficient to sample the distribution of the maxima and that the resulting distribution can be very well fitted by the functional shape of the generalized extreme value distribution. In general, we find that the distribution of the largest Einstein radii can be well described by a Fréchet (Type II) distribution, indicating that the distribution is bound from below. Furthermore, we confirm the findings of OB09 that the sample of maxima is distributed in a wide range in the massredshift plane, indicating that the single largest Einstein radius has its origin by no means necessarily in the most massive haloes. This indicates that the triaxiality, together with the halo orientation with respect to the observer, has a stronger impact than the mass of the cluster itself. We also report that different definitions of the Einstein radius, like the effective or the median one, can lead to the selection of different haloes for the largest Einstein radius, particularly for smaller values of the maximum Einstein radius. However, both definitions lead to identical results for the largest realizations.
In the second part, we studied the influence of different underlying assumptions and effects on the resulting extreme value distributions. The results of this part can be summarised as follows.

Mass function. We sampled the extreme value distribution for four different mass functions, comprising the Press & Schechter (1974), Sheth & Tormen (1999), Tinker et al. (2008) and Crocce et al. (2010) mass functions. Modifying the size of the halo sample in each mock universe from which the maximum Einstein radius will be sampled, the different mass functions lead mainly to a shift to smaller or larger Einstein radii, while the impact on the shape of the distribution is less pronounced. Of course, the effect of a different normalisation of the matter power spectrum, σ_{8}, is similar in nature.

Triaxiality. We studied the impact of triaxiality by introducing different cutoffs in the underlying axisratio distribution, with the result that the extreme value distribution of the largest Einstein radius is very sensitive to the presence of small axis ratios and hence, very elongated objects. The different cutoffs lead not only to shifts of the resulting extreme value distributions, but also substantially influence the shape of the distribution. The more elongated objects are allowed to exist, the higher will be the tail of the extreme value distributions towards very large values of the Einstein radius.

Inner slope and c–M relation. Both underlying assumptions impact on the resulting distributions and exhibit a particular behaviour, which is caused by projection effects. Since the extreme value distribution will be naturally based on elongated haloes that are oriented along the line of sight, the stronglensing signal can be significantly enhanced from the outer regions of the halo. It remains unclear if and where a radius cutoff of the density profile, e.g. at the virial radius (see e.g. Oguri & Keeton 2004; Baltz et al. 2009), should be imposed. This adds an additional uncertainty to the results that have been listed so far.
With our study, we could show that a multitude of underlying assumptions strongly influence the extreme value distribution of the largest Einstein radius. Many of those require more detailed studies, e.g. the triaxiality of dark matter haloes. Another effect having a strong impact on the extreme value distribution is the presence of dynamical mergers as shown in the work of Redlich et al. (2012), and it will be studied by the authors in further detail in future work. In view of this complexity, it is unlikely that the extreme value distribution of the Einstein radius can be used for consistency test of ΛCDM. However, due to its enhanced sensitivity to the underlying assumptions, it could very well be used to learn more about these assumptions.
In the last part of this work, we used the previously studied extreme value distributions to assess the probability of occurrence for the largest known Einstein radius of MACS J0717.5+3745 (Zitrin et al. 2009). Accounting only for the uncertainty in σ_{8}, we find for the observed effective Einstein radius of θ_{eff} = 55 ± 3 arcsec an occurrence probability of ~(11−42) per cent for the MACS survey area and of ~(18−61) per cent on the full sky, indicating that this observation can not be considered in conflict with ΛCDM. This conclusion is supported by the fact that the probability range would widen further if we would account for the uncertainties in the underlying assumptions, for instance, mass function and triaxiality, rendering any claim of tension with ΛCDMuntenable. Furthermore, we neglected the impact of dynamical merging for which MACS J0717.5+3745 is a prime example, which again would make extremely large critical curves more likely to be found.
However, apart from our results for the large Einstein radius, MACS J0717.5+3745 is a candidate for the most massive known galaxy cluster in the redshift range of 0.5 ≤ z ≤ 1.0, as indicated by the mass enclosed by the critical curve of ~1×10^{15} M_{⊙} (Zitrin et al. 2009). Since this is the mass that is contained only in the innermost region, the overall cluster mass is expected to be significantly larger^{2}. A more thorough mass estimate for MACS J0717.5+3745 is expected to be provided by the Cluster Lensing And Supernova survey with Hubble (CLASH) (Postman et al. 2012). Inspired by this result, we calculated the total mass a galaxy cluster would need to have in the redshift range of 0.5 ≤ z ≤ 1.0 in order to exhibit significant tension with ΛCDM. To do so, we utilised the extreme value statistical approach for halo masses used in Waizmann et al. (2011) and found that for a 3σ deviation from the ΛCDM model, the cluster would need to have a mass of at least M_{3σ} ≃ 4×10^{15} M_{⊙}. This value needs to be even larger in order to account for the correction for the Eddington bias (see e.g. Mortonson et al. 2011) and the unavoidable uncertainties in the mass determination. Whether the mass for MACS J0717.5+3745 will reach such high values remains to be seen.
As a closing remark, we conclude that it seems to be more than doubtful that the single largest observed Einstein radius can be used as a basis for ΛCDM falsification experiments. However, we expect nevertheless useful insights into the underlying assumptions that enter the modelling of the Einstein radius distribution. In the future, we intend to perform further studies along these lines.
It should be mentioned that the recent strong lensing analysis of MACS J0717.5+3745 by Limousin et al. (2012) finds a redshift for the primary lensed galaxy of z_{s} ~ 2.96 instead of z_{s} ~ 2.5 used by Zitrin et al. (2009). This change would lead to a shift in the surface mass normalization, κ, and hence to a smaller mass estimate.
Acknowledgments
We would like to thank Adi Zitrin and Marceau Limousin for the very helpful discussions. J.C.W. acknowledges financial contributions from the contracts ASIINAF I/023/05/0, ASIINAF I/088/06/0, ASI I/016/07/0 COFIS, ASI EuclidDUNE I/064/08/0, ASIUni BolognaAstronomy Dept. EuclidNIS I/039/10/0, and PRIN MIUR 2008 Dark energy and cosmology with large galaxy surveys. M.R. thanks the Sydney Institute for Astronomy for the hospitality and the German Academic Exchange Service (DAAD) for the financial support. Furthermore, M.R.’s work was supported in part by contract research Internationale Spitzenforschung II1 of the BadenWürttemberg Stiftung. M.B. is supported in part by the TransregioSonderforschungsbereich The Dark Universe of the German Science Foundation.
References
 Baldi, M., & Pettorino, V. 2011, MNRAS, 412, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Baltz, E. A., Marshall, P., & Oguri, M. 2009, Cosmol. Astropart. Phys., 1, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Bartelmann, M. 2010, Class. Quant. Grav., 27, 233001 [NASA ADS] [CrossRef] [Google Scholar]
 Bhattacharya, S., Heitmann, K., White, M., et al. 2011, ApJ, 732, 122 [NASA ADS] [CrossRef] [Google Scholar]
 Bonafede, A., Feretti, L., Giovannini, G., et al. 2009, A&A, 503, 707 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Broadhurst, T. J., & Barkana, R. 2008, MNRAS, 390, 1647 [NASA ADS] [Google Scholar]
 Cappelluti, N., Predehl, P., Böhringer, H., et al. 2011, Mem. Soc. Astron. It. Supp., 17, 159 [Google Scholar]
 Carlstrom, J. E., Ade, P. A. R., Aird, K. A., et al. 2011, PASP, 123, 568 [NASA ADS] [CrossRef] [Google Scholar]
 Cayón, L., Gordon, C., & Silk, J. 2011, MNRAS, 415, 849 [NASA ADS] [CrossRef] [Google Scholar]
 Chongchitnan, S., & Silk, J. 2012, Phys. Rev. D, 85, 063508 [NASA ADS] [CrossRef] [Google Scholar]
 Coles, S. 2001, An Introduction to Statistical Modeling of Extreme Values (Springer) [Google Scholar]
 Crocce, M., Fosalba, P., Castander, F. J., & Gaztañaga, E. 2010, MNRAS, 403, 1353 [NASA ADS] [CrossRef] [Google Scholar]
 Dalal, N., Holder, G., & Hennawi, J. F. 2004, ApJ, 609, 50 [NASA ADS] [CrossRef] [Google Scholar]
 Davis, O., Devriendt, J., Colombi, S., Silk, J., & Pichon, C. 2011, MNRAS, 413, 2087 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Edge, A. C., & Henry, J. P. 2001, ApJ, 553, 668 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Barrett, E., & Donovan, D. 2004, ApJ, 609, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Barrett, E., Donovan, D., et al. 2007, ApJ, 661, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Eddington, A. S. 1913, MNRAS, 73, 359 [NASA ADS] [CrossRef] [Google Scholar]
 Edge, A. C., Ebeling, H., Bremer, M., et al. 2003, MNRAS, 339, 913 [NASA ADS] [CrossRef] [Google Scholar]
 Fisher, R., & Tippett, L. 1928, Proc. Cambridge Phil. Soc., 24, 180 [NASA ADS] [CrossRef] [Google Scholar]
 Foley, R. J., Andersson, K., Bazin, G., et al. 2011, Phys. Rev. D, 731, 86 [Google Scholar]
 Gnedenko, B. 1943, Ann. Math., 44, 423 [CrossRef] [Google Scholar]
 Green, S. R., & Wald, R. M. 2012, Phys. Rev. D, 85, 063512 [NASA ADS] [CrossRef] [Google Scholar]
 Gumbel, E. 1958, Statistics of Extremes (New York: Columbia University Press, reprinted in 2004 New York: Dover) [Google Scholar]
 Harrison, I., & Coles, P. 2011, MNRAS, 418, L20 [NASA ADS] [Google Scholar]
 Harrison, I., & Coles, P. 2012, MNRAS, 421, L19 [NASA ADS] [Google Scholar]
 Holz, D. E., & Perlmutter, S. 2012, ApJ, 755, L36 [NASA ADS] [CrossRef] [Google Scholar]
 Hotchkiss, S. 2011, J. Cosmol. Astropart. Phys., 7, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Jee, M. J., Rosati, P., Ford, H. C., et al. 2009, ApJ, 704, 672 [NASA ADS] [CrossRef] [Google Scholar]
 Jenkinson, A. F. 1955, Quarterly J. Roy. Metereol. Soc., 81, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Jing, Y. P., & Suto, Y. 2002, ApJ, 574, 538 [NASA ADS] [CrossRef] [Google Scholar]
 Keeton, C. R., & Madau, P. 2001, ApJ, 549, L25 [NASA ADS] [CrossRef] [Google Scholar]
 Kneib, J.P., & Natarajan, P. 2011, A&ARv, 19, 47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kotz, S., & Nadarajah, S. 2000, Extreme Value Distributions – Theory and Applications (London: Imperial College Press) [Google Scholar]
 LaRoque, S. J., Joy, M., Carlstrom, J. E., et al. 2003, ApJ, 583, 559 [NASA ADS] [CrossRef] [Google Scholar]
 Laureijs, R., Amiaux, J., Arduini, S., et al. 2011 [arXiv:1110.3193] [Google Scholar]
 Limousin, M., Ebeling, H., Richard, J., et al. 2012, A&A, 544, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ma, C.J., Ebeling, H., Donovan, D., & Barrett, E. 2008, ApJ, 684, 160 [NASA ADS] [CrossRef] [Google Scholar]
 Marriage, T. A., Acquaviva, V., Ade, P. A. R., et al. 2011, ApJ, 737, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Menanteau, F., Hughes, J. P., Sifón, C., et al. 2012, ApJ, 748, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Meneghetti, M., Argazzi, R., Pace, F., et al. 2007, A&A, 461, 25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meneghetti, M., Fedeli, C., Pace, F., Gottlöber, S., & Yepes, G. 2010, A&A, 519, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meneghetti, M., Fedeli, C., Zitrin, A., et al. 2011, A&A, 530, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mortonson, M. J., Hu, W., & Huterer, D. 2011, Phys. Rev. D, 83, 023015 [NASA ADS] [CrossRef] [Google Scholar]
 Mroczkowski, T., Dicker, S., Sayers, J., et al. 2012, ApJ, submitted [arXiv:1205.0052] [Google Scholar]
 Mullis, C. R., Rosati, P., Lamer, G., et al. 2005, ApJ, 623, L85 [NASA ADS] [CrossRef] [Google Scholar]
 Nakamura, T. T., & Suto, Y. 1997, Progr. Theor. Phys., 97, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M. 2004, Ph.D. Thesis, The University of Tokyo [Google Scholar]
 Oguri, M., & Blandford, R. D. 2009, MNRAS, 392, 930 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M., & Keeton, C. R. 2004, ApJ, 610, 663 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M., Lee, J., & Suto, Y. 2003, ApJ, 599, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Pickands, J. 1975, Annals of Statistics, 3, 119 [CrossRef] [MathSciNet] [Google Scholar]
 Postman, M., Coe, D., Benítez, N., et al. 2012, ApJS, 199, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 [NASA ADS] [CrossRef] [Google Scholar]
 Puchwein, E., Bartelmann, M., Dolag, K., & Meneghetti, M. 2005, A&A, 442, 405 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Redlich, M., Bartelmann, M., Waizmann, J.C., & Fedeli, C. 2012, A&A, 547, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reiss, R.D., & Thomas, M. 2007, Statistical Analysis of Extreme Values, 3rd edn. (Basel: Birkhauser Verlag) [Google Scholar]
 Rosati, P., Tozzi, P., Gobat, R., et al. 2009, A&A, 508, 583 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comm. Astrophys. Space Phys., 4, 173 [NASA ADS] [EDP Sciences] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, I. B. 1980, ARA&A, 18, 537 [NASA ADS] [CrossRef] [Google Scholar]
 Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709 [NASA ADS] [CrossRef] [Google Scholar]
 von Mises, R. 1954, Americ. Math. Soc., Vol. II, 271 [Google Scholar]
 Waizmann, J.C., Ettori, S., & Moscardini, L. 2011, MNRAS, 418, 456 [NASA ADS] [CrossRef] [Google Scholar]
 Waizmann, J.C., Ettori, S., & Moscardini, L. 2012a, MNRAS, 420, 1754 [NASA ADS] [CrossRef] [Google Scholar]
 Waizmann, J.C., Ettori, S., & Moscardini, L. 2012b, MNRAS, 2810 [Google Scholar]
 Williamson, R., Benson, B. A., High, F. W., et al. 2011, ApJ, 738, 139 [NASA ADS] [CrossRef] [Google Scholar]
 Zitrin, A., Broadhurst, T., Rephaeli, Y., & Sadeh, S. 2009, ApJ, 707, L102 [NASA ADS] [CrossRef] [Google Scholar]
 Zitrin, A., Broadhurst, T., Barkana, R., Rephaeli, Y., & Benítez, N. 2011, MNRAS, 410, 1939 [NASA ADS] [Google Scholar]
All Tables
Values of the location, scale and shape parameters, α_{eff}, β_{eff} and γ_{eff}, of the fitted GEV distributions of the maximum effective Einstein radius for different sample sizes N_{samp}, as shown in the lower left panel of .
Shape parameter γ and the mode of the CDFs of the largest effective and median Einstein radius for four different mass functions.
Shape parameter γ and the mode of the CDFs for the largest effective and median Einstein radius for different cutoffs in the axisratio distribution.
All Figures
Fig. 1 Distribution in mass and redshift of 4000 maxima of the effective (left panel) and the median Einstein radius (right panel) in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky based on the Tinker et al. (2008) mass function. The colour encodes the size of the individual largest Einstein radius from each simulation run. The green circles denote exemplary systems for which the two Einstein radius definitions selected different haloes as the maximum. 

Open with DEXTER  
In the text 
Fig. 2 CDFs of the largest effective (left panels) and median (right panels) Einstein radius for a different number of maxima in the range between 125 and 4000, assuming the Tinker et al. (2008) mass function, M_{lim} = 2×10^{14} M_{⊙} / h and the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. The upper row shows the CDFs directly based on the MC simulations, and the lower row presents the corresponding fits of the GEV distribution. The small lower panels show the difference, ΔCDF, with respect to the highresolution run based on 4000 maxima. 

Open with DEXTER  
In the text 
Fig. 3 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and the median (right panels) Einstein radius for the mass functions of PressSchechter, Tinker, Sheth & Tormen and Crocce as labelled in the panels. All distributions are based on the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER  
In the text 
Fig. 4 Ratio of the number of haloes N( > M) for the different mass functions as labelled in the panel, with respect to the Tinker one, assuming a redshift interval of 0.5 ≤ z ≤ 1.0 and the full sky. 

Open with DEXTER  
In the text 
Fig. 5 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and median (right panels) Einstein radius for cutoffs of the axisratio distribution as labelled in the figure. All distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER  
In the text 
Fig. 6 The lefthand and centre panels show the distribution in scaled axis ratio a_{sc} and concentration c_{e} of sampled haloes, according to different selection criteria. The lefthand panel shows the distribution of a full halo sample of a single realisation (black dots) and the sample of the minima in a_{sc} from 1000 realisations (red dots). The centre panel shows the distribution of the haloes that give rise to the largest effective Einstein ring (encoded in the colourbar) based on 1000 realisations. The small arrows denote the indicated cutoffs identical to . The righthand panel shows the distribution of the same maxima with respect to the scaled axis ratio and the alignment  cosθ  , where θ is the angle between the major axis of each halo with respect to the line of sight. 

Open with DEXTER  
In the text 
Fig. 7 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and median (right panels) Einstein radius for different values of the inner slope of the density profile α_{NFW} = (1.0,1.5) as labelled in the figure. All distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER  
In the text 
Fig. 8 CDFs (upper panels) and PDFs (lower panels) of the largest effective (left panels) and median (right panels) Einstein radius for different values of the normalisation parameter A_{e} of the massconcentration relation as labelled in the figure. All distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h in the redshift interval of 0.5 ≤ z ≤ 1.0 on the full sky. 

Open with DEXTER  
In the text 
Fig. 9 CDFs of the largest effective (left panel) and median (right panel) Einstein radius, assuming a redshift interval of 0.5 ≤ z ≤ 1.0, a source redshift of z_{s} = 2.5 and the nominal MACS survey area (upper panels) as well as the full sky (lower panels). Both distributions are based on the Tinker et al. (2008) mass function and the simulation of 1000 maxima with M_{lim} = 2×10^{14} M_{⊙} / h on the full sky. The dasheddotted lines, together with the yellow shaded area, illustrate the impact of the uncertainty in the WMAP7 value of σ_{8} on the CDFs and the grey shaded area denotes the uncertainty in the measurement of the Einstein radius. 

Open with DEXTER  
In the text 
Fig. 10 CDF of the largest effective Einstein radius, assuming a redshift interval of 0.5 ≤ z ≤ 1.0, the full sky and a source redshift of z_{s} = 2.5. The dasheddotted lines illustrate the uncertainty in the CDF due to the imprecision of the WMAP7 value of σ_{8}. The vertical shaded regions denote the values of θ_{eff} for which the CDF takes a value equal to nσ with n ∈ [ 1,5] , as labelled in the panel. The red filled circle denotes the measured value of θ_{eff} for MACS J0717.5 and the short red arrow points at the corresponding value of the CDF. 

Open with DEXTER  
In the text 
Fig. 11 CDF of the most massive halo that is expected to be found in the redshift interval of 0.5 ≤ z ≤ 1.0 in an survey area of A_{s} = 22 735 deg^{2}. The dasheddotted lines illustrate the uncertainty in the CDF due to the imprecision of the WMAP7 value of σ_{8}. The vertical shaded regions denote the values of the mass for which the CDF takes a value equal to nσ with n ∈ [ 1,5] , as labelled in the panel. We also added for reference the mass of ACTCL J0102 (red errorbar), the most massive known cluster in the given redshift range, where the two small red arrows point at the values of the CDF stemming from the upper allowed mass for the distributions based on σ_{8} (lower arrow) and σ_{8} − Δσ_{8} (upper arrow). 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.