The projected gravitational potential of the galaxy cluster MACS J1206 derived from galaxy kinematics
^{1} Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, Philosophenweg 12, 69120 Heidelberg, Germany
email: dennis.stock@stud.uniheidelberg.de
^{2} INAF–Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 34131 Trieste, Italy
^{3} Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen, Denmark
^{4} Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA
^{5} INAF–Osservatorio Astronomico di Capodimonte, via Moiariello 16, 80131 Napoli, Italy
^{6} Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, via Saragat 1, 44122 Ferrara, Italy
Received: 22 July 2015
Accepted: 7 October 2015
We reconstruct the radial profile of the projected gravitational potential of the galaxy cluster MACS J1206 from 592 spectroscopic measurements of velocities of cluster members. To accomplish this, we use a method we have developed recently based on the RichardsonLucy deprojection algorithm and an inversion of the sphericallysymmetric Jeans equation. We find that, within the uncertainties, our reconstruction agrees very well with a potential reconstruction from weak and strong gravitational lensing as well as with a potential obtained from Xray measurements. In addition, our reconstruction is in good agreement with several common analytic profiles of the lensing potential. Varying the anisotropy parameter in the Jeans equation, we find that isotropy parameters, which are either small, β ≲ 0.2, or decrease with radius, yield potential profiles that strongly disagree with that obtained from gravitational lensing. We achieve the best agreement between our potential profile and the profile from gravitational lensing if the anisotropy parameter rises steeply to β ≈ 0.6 within ≈ 0.5 Mpc and stays constant further out.
Key words: galaxies: clusters: general / galaxies: kinematics and dynamics / dark matter / gravitational lensing: strong / gravitational lensing: weak
© ESO, 2015
1. Introduction
Galaxy clusters offer several classes of observables reflecting their overall internal constitution: gravitational lensing effects in their weak and strong variants, Xray emission, the thermal SunyaevZel’dovich (tSZ) effect, and the kinematics of their member galaxies. We neglect radio emission, turbulence, and metal abundance in the intracluster gas and the population statistics of member galaxies because these effects are locally driven.
Gravitational lensing measures the gravitational tidal field, projected along the line of sight (l.o.s.) and thus directly probes the projected gravitational potential more precisely its curvature. The Xray emission and the thermal SunyaevZel’dovich effect depend on powers of the density and temperature of the intracluster medium. If equilibrium assumptions hold, hydrostatic and virial equilibrium foremost, the Xray emission, the thermal SunyaevZel’dovich effect and galaxy kinematics are also determined by the gravitational potential. Current discussions of the validity of these equilibrium assumptions can be found e.g. in Martino et al. (2014), von der Linden et al. (2014), Planck Collaboration Int. III (2013), Planck Collaboration XX (2014).
What is the gravitational potential that agrees best with all cluster observables? This question is relevant for different reasons. First, different observables trace the gravitational potential at different scales. The cluster core can be probed by the stellar kinematics of the brightest cluster galaxy, strong lensing and Xray emission probe the innermost regions, weak lensing and galaxy kinematics probe large scales, and the tSZ effect falls in between (see Umetsu et al. 2014; Merten et al. 2015; Newman et al. 2013). Aiming at a reliable reconstruction of cluster density profiles from their cores to their outskirts, combining all observables into a unique potential reconstruction offers the advantage of covering all relevant scales in a single step. Second, comparing cluster potential reconstructions based on lensing, on the one hand, and based on the rest of the observables, on the other hand, allows us to test the equilibrium assumptions or possible deviations therefrom. Third, lensing and the other observables do not necessarily see the same gravitational potential. While lensing is sensitive to the sum of the Bardeen potentials, the other observables only probe the spatial potential. In general relativity, the two Bardeen potentials agree in case of negligible anisotropic stress. Differences in potential reconstructions based on lensing compared to other observables may also hint at deviations from relativity (see also Sartoris et al. 2014; Barreira et al. 2015, in this context).
The gravitational potential has the major and important advantage that it is a locally measurable quantity, which is directly related to the observables listed above (see also Angrick & Bartelmann 2009, 2012).
We have recently developed methods for reconstructing the projected gravitational potential of galaxy clusters from their Xray emission, their tSZ effect and the kinematics of their member galaxies (Sarli et al. 2014; Konrad et al. 2013; Majer et al. 2013). These methods operate similarly, but with important differences in detail. That is, an observable is deprojected by means of the RichardsonLucy algorithm, requiring symmetry assumptions. The deprojected quantities are related to the threedimensional gravitational potential by relations derived from justifiable equilibrium assumptions. The gravitational potential can then be projected along the line of sight. These methods complement our techniques for joint cluster reconstruction from weak and strong gravitational lensing (Bartelmann et al. 1996; Cacciato et al. 2006; Merten et al. 2009; see also Coe et al. 2012; Merten et al. 2011; Merten 2014 for examples).
In this paper, we reconstruct the projected gravitational potential of the galaxy cluster MACS J1206.2−0847 based on galaxy kinematics, applying the technique developed in Sarli et al. (2014). This work is structured as follows: Sect. 2 briefly reviews the reconstruction method. In Sect. 3, we describe the selection and preparation of the data. Section 4 presents the results and compares the projected gravitational potential to that obtained from gravitational lensing. We summarise in Sect. 5 and discuss other methods in the Appendix.
2. Reconstruction method
The basic assumptions of our reconstruction method are the following: we treat the galaxy cluster as a spherically symmetric gas cloud of collisionless, pointlike test particles of mass m, i.e. galaxies, moving in the gravitational potential of their common dark matter halo. The system can then be described by the Jeans equation relating the radial velocity dispersion weighted by the galaxy number density ρ_{gal} to the gravitational potential φ: = Φ /m, (1)where the anisotropy parameter quantifies the ratio between the tangential and the radial velocity dispersions.
Expanding on the formal analogy with gas dynamics, we define an effective galaxy pressure . Using this definition and setting β to zero, (1) would turn into the equation of hydrostatic equilibrium for a gas. In addition, we introduce a polytropic relation between the effective galaxy pressure and the matter density, (2)This was justified for a simulated cluster in Sarli et al. (2014). In fact, the authors showed that the polytropic assumption is reasonable for several density profiles. Additionally, we shall see in our reconstruction results that the exact choice of the polytropic index only affects the final results very mildly.
The reconstruction algorithm described in detail in Sarli et al. (2014) proceeds along the following steps (see also Fig. 1):

1.
In order to obtain the effective pressure P, we first have to deproject the actual observable, i.e. the lineofsight projected velocity dispersion weighted by the galaxy number density, . This is accomplished via the RichardsonLucy deconvolution or deprojection, see Lucy (1974, 1994).

2.
Furthermore, we make use of the polytropic relation (2) to rewrite the Jeans Eq. (1) in terms of the effective pressure. This leads to a Volterra integral equation of the second kind for the gravitational potential Φ. After fixing the shape of the anisotropy profile, we can solve this equation in a quickly converging iteration process.

3.
Finally, we project Φ along the line of sight to find the projected gravitational potential Ψ.
This entire algorithm depends on three parameters and a function, which is carefully adapted during the reconstruction. Noise suppression within the RichardsonLucy deprojection requires regularisation, controlled by a smoothing scale L and an amplitude α. Furthermore, we have to choose the polytropic index γ in (2), usually being of order unity, as cluster analyses in Sarli et al. (2014) suggest. This also implies that the galaxy fluid can approximately be treated as an isothermal gas. The importance of the polytropic index is investigated in more detail in Sect. 4.1.
As we shall demonstrate later, the dominant parameter is the anisotropy profile β(r). However, if kinematic data are used exclusively for cluster reconstruction, there is in principle a degeneracy between the anisotropy parameter and the gravitational potential. We resolve this wellknown degeneracy by fixing the βprofile to obtain a nonparametric estimate for the gravitational potential. To avoid an arbitrary, unmotivated guess for β(r), we choose it such that the reconstructed potential agrees best with a reconstruction based on gravitationallensing data. Alternative methods for breaking the anisotropymass degeneracy and for recovering the gravitational potential with galaxy kinematics are discussed in the last section.
3. Data selection and preparation
The observational input data for the relaxed and massive galaxy cluster MACS J1206, Ebeling et al. (2001, 2009), (M_{lens} ≈ 5 × 10^{14}M_{⊙}h^{1}, Umetsu et al. 2012) investigated here consist of a catalogue of member galaxies (Biviano et al. 2013) observed in the context of the CLASH project (Postman et al. 2012) as part of a large spectroscopic campaign carried out with the Very Large Telescope (VLT; CLASHVLT Large Programme; Rosati et al. 2014). According to Biviano et al. (2013) and Annunziatella et al. (2014), the spatial incompleteness of the spectroscopic sample varies with position in the cluster by less than 20%.
Fig. 1 Scheme of the reconstruction algorithm. The observable is converted to a threedimensional quantity by RichardsonLucy deprojection. This is then turned into the threedimensional gravitational potential, solving the Jeans equation. The resulting potential is finally projected along the line of sight. 

Open with DEXTER 
The following analysis is based on 592 member galaxies in total, selected via the method presented in Biviano et al. (2013). In order to arrive at the lineofsight projected velocity dispersion profile, we assign to every cluster member a projected radius from the cluster centre, which we take to be marked by the brightest cluster galaxy. This is achieved by multiplying the angular separation with the angulardiameter distance, calculated with a standard ΛCDMcosmology with recent cosmological parameters from Planck Collaboration XVI (2014).
We determine the lineofsight projected velocities from the measured redshifts as described in Harrison (1974). In our case, we can use the nonrelativistic relation between redshift z and lineofsight velocity v_{los}, since all involved velocities (see below) are clearly in the nonrelativistic regime. Thus, (3)As we are only interested in the galaxy velocities with respect to the cluster’s centreofmass, we have to take into account that the observed redshift includes four main contributions: the motion of the observer with respect to a local comoving observer, the motion of the cluster’s centreofmass with respect to a local comoving observer, the motion of each galaxy with respect to the cluster’s centreofmass, and finally the Hubble expansion. Harrison (1974) has shown that neglecting the motion of the cluster’s centre of mass leads to a multiplication of redshifts, (4)where z_{0} is the redshift due to the observer’s motion with respect to a local comoving observer, z_{cosm} represents the Hubble expansion, and z_{G} the motion of each galaxy with respect to a local comoving observer. To leading order, we can calculate the cluster average of Eq. (4) and establish that the average redshift of all member galaxies with respect to the local comoving frame vanishes because we assume an isotropic distribution of velocities along the l.o.s., (5)Furthermore, we can safely assume that the contribution by the Hubble expansion is the same for all cluster members, (6)Thus, we end up with (7)Inserting Eq. (7) into Eq. (4), we find the following expression for v_{los}, see also Fig. 2: (8)To arrive at the densityweighted, projected velocity dispersion profile , we adapt the bin width to that chosen for the potential from lensing for better comparability, i.e. a constant bin width of 0.1 Mpc. Next, we calculate the velocity dispersion within each bin and reconstruct the projected gravitational potential as described above by solving the radial Jeans equation.
Fig. 2 Velocities projected along the l.o.s. of cluster members. 

Open with DEXTER 
Because of the gauge freedom of the lensing potential, one is in particular allowed to add an arbitrary constant. By convention, we scale the potential such that Ψ(0) = 1 and use the normalisation Ψ(R_{cut}) = 0 for a given cutoff radius R_{cut} = 3 Mpc. This is even larger than the virial radius of MACS J1206, which is approximately 2 Mpc (Umetsu et al. 2012).
We perform a bootstrap analysis, i.e. for each bootstrap sample, we draw as many times with replacement from the original dataset as there are member galaxies to obtain error bars and to suppress the effect of outliers. Then, we reconstruct the gravitational potential for each individual bootstrap sample. This procedure is repeated 300 times. Finally, we calculate the mean of all potentials and their standard deviation. The bootstrap analysis also allows us to assess the effects of incompleteness and uncertain membership assignment by varying the sample of galaxies entering into the potential reconstruction.
4. Results
4.1. Effect of reconstruction parameters
As indicated above, we have to fix three reconstruction parameters to perform the reconstruction: the regularisation amplitude α, the smoothing scale L, and the polytropic index γ. Additionally, we have to model the anisotropy profile β(r).
Figure 3 shows how the variation of these reconstruction parameters affects the resulting potential Ψ. The amplitude of the anisotropy parameter β clearly has the dominant effect. In particular, an isotropic velocity distribution corresponding to β → 0 seems very implausible because of the strongly varying shape of the potential resulting from this assumption.
Although the assumption of a polytropic relation (2) may appear very bold, Fig. 3c shows that the particular choice of the polytropic index does not really matter. Also, the effects of the regularisation parameters α and L remain within the uncertainties, which is of similar order as shown in Fig. 6. When they are not varied, the parameters are set as follows: α = 0.3, L = 0.6, γ = 1.1, and β = 0.6. The value chosen for the polytropic index γ is motivated by the result of Sarli et al. (2014) that γ would usually be around unity. The regularisation parameters are chosen such that the agreement with the lensing reconstruction is best.
Fig. 3 Reconstructed radial profiles of the projected gravitational potential of MACS J1206, obtained with different parameter sets for the smoothing amplitude α (panel a)), smoothing scale L (panel b)), the polytropic index γ (panel c)), and the anisotropy parameter β (panel d)). 

Open with DEXTER 
4.2. Variable anisotropy profile
We also systematically investigate the effect of an anisotropy parameter varying with radius. To accomplish this, we compare radial profiles of gravitational potentials obtained with an anisotropy parameter increasing or decreasing linearly in discrete steps every 0.5 Mpc. In case of the decreasing anisotropyparameter profile, we choose β = 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, and reverse for the increasing case.
Figure 4 shows both cases together with the potential reconstruction from weak and strong gravitational lensing and a reconstruction with constant β = 0.6. First, one can observe that the qualitative form of the anisotropy profile affects the curvature of the potential: An increasing β profile leads to a potential with negative curvature. In contrast, a decreasing βprofile leads to a positive curvature and strong variations in the potential at radii beyond 1 Mpc. The reconstruction using a constant β leads to an almost straight line for the radial potential profile.
These results motivate the shape of the anisotropy profile used later in Eq. (9). Since the potential from lensing first shows a negative curvature at small radii and a linear evolution later, we choose an anisotropy profile increasing for small radii and turning constant afterwards to reproduce the linear behaviour.
Fig. 4 Radial profiles of the projected gravitational potential, reconstructed from galaxy kinematics with different profiles of the anisotropy parameter β. The other reconstruction parameters remain fixed at α = 0.3, L = 0.6 and γ = 1.1. 

Open with DEXTER 
4.3. Comparison with lensing and Xray data
As a first test, we compare the projected gravitational potential reconstructed from galaxy velocities with a reconstruction based on weak and strong gravitational lensing data of Merten et al. (2015) following constraints from Zitrin et al. (2012). After converting the 2Dlensing potential map (see Fig. 5) into a radial profile and using the same normalisation process as above, we can compare the result from gravitational lensing to the potential obtained from cluster kinematics (Fig. 6).
Fig. 5 Twodimensional map of the lensing potential of MACS1206 based on weak and strong lensing. 

Open with DEXTER 
The reconstruction parameters are chosen as listed above. For the radial profile of the anisotropy parameter β, we assume the following radially dependent profile, which is motivated in Sect. 4.2. This profile turns out to be in qualitative agreement with that in Biviano et al. (2013), i.e. (9)Additionally, we compare it with a reconstruction based on Xray emission of the intracluster gas based on Donahue et al. (2014) using Chandra data up to 0.72 Mpc, which is a typical range of validity using this method for galaxy clusters.
Bestfit parameters for different models.
The growing difference between the lensing and kinematic profile outside ~ 2 Mpc can be explained by the limited range of validity of each reconstruction method. The small difference at the innermost radii is due to increasing baryonic effects causing the equilibrium assumption ultimately to break down. The Xray data agree well within the error boundaries.
Fig. 6 Radial profiles of the gravitational potential of MACS J1206, reconstructed from combined strong and weak gravitational lensing (blue), Xray emission (yellow), and from galaxy kinematics (red) with the parameters α = 0.3, L = 0.6, γ = 1.1 and β as in (9). 

Open with DEXTER 
4.4. Comparison with common lensing potentials
Furthermore, we can compare the radial profile of the projected gravitational potential obtained from galaxy kinematics with radial profiles common in gravitationallensing studies, such as the singular (SIS) and the nonsingular (NIS) isothermal spheres and the NFW profile (Navarro et al. 1996; Narayan & Bartelmann 1996); see Fig. 7, i.e.: (10)with the function h(x), depending on the radius in units of the scale radius x: = r/r_{s} (see Golse & Kneib (2002), i.e. (11)where s_{c}, r_{s} and A...F are free parameters determined by fitting the models (10) to the potential obtained from galaxy kinematics. The parameter values are summarised in Table 1.
We can compare the results for the scales s_{c} in case of the softened isothermal sphere and r_{s} in case of NFW with Umetsu et al. (2012). They obtain The agreement is quantified by a goodnessoffit parameter Ω, which we define in analogy to the χ^{2}function appropriate for uncorrelated measurements. Since our data points are not independent but correlated by the bootstrap method, i.e. some data points appear multiple times, we cannot interpret Ω as a χ^{2}function. Because we just want to single out the best model, we do not go into a more elaborate analysis for correlated measurements here. For the ith data point x_{i} with standard deviation σ_{i} and the corresponding model prediction f(x_{i}), we write (14)As one can already see by eye, the softened isothermal sphere and the NFWmodel describe the data equally well. More quantitatively, the goodnessoffit parameters for all lensing model are (15)So far, all analytic lensing potentials are compatible with the reconstruction from galaxy kinematics, but the error bars are substantial.
Fig. 7 Comparison of three analytic profiles of gravitationallensing potentials with the gravitationalpotential profile of MACS J1206 reconstructed from galaxy kinematics using α = 0.3, L = 0.6, γ = 1.1 and β as in Eq. (9). 

Open with DEXTER 
5. Summary
We have recently developed a method for reconstructing the projected gravitational potential of galaxy clusters from the kinematics of their member galaxies. In our study, we applied this method for the first time to a real galaxy cluster (MACS J1206). The reconstruction assumes that the threedimensional gravitational potential of this cluster is approximately spherically symmetric. The potential reconstruction is based on 592 measured galaxy velocities. Our results can be summarised as follows: 1) The reconstruction algorithm requires four parameters to be set: the anisotropy parameter β; the polytropic index γ of the effective galaxy pressure; the amplitude α of the regularisation term in the RichardsonLucy deconvolution and the length scale L of the regularisation. We found that the reconstructed potential is almost insensitive to the socalled nuisance parameters γ, α, and L, while it does depend very sensitively on the anisotropy parameter β. 2) Within the bootstrap error bars, the reconstructed potential profile is indistinguishable from the potential profile obtained from the combination of weak and strong gravitational lensing and from Xray analysis. Typical analytic lensing profiles, such as the singular and nonsingular isothermal spheres and the NFW profile, also agree well with our reconstruction from galaxy kinematics. 3) Small anisotropy parameters β ≲ 0.2 or anisotropyparameter profiles decreasing with radius yield potential profiles differing strongly from the profile obtained from gravitational lensing. Best agreement with gravitational lensing is achieved with an anisotropy parameter increasing to β ≈ 0.6 within 0.5 Mpc and staying constant further out.
Our results show that the method returns convincing results on one wellstudied galaxy cluster, which we use here as a test case. Even though the uncertainties are still large, the anisotropy parameter can be rather well constrained by comparing our results with those obtained from gravitational lensing. We see three main future applications of this method: (1) in joint reconstructions of galaxycluster potentials compatible with all cluster observables; (2) in constraints of the anisotropy parameter; and (3) in tests of fundamental assumptions, such as hydrostatic or virial equilibrium, and possibly also of general relativity, if applied to large and wellmeasured cluster samples.
Acknowledgments
This work was supported in part by the Collaborative Research Center TR 33 “The Dark Universe” of the Deutsche Forschungsgemeinschaft. This work is based on data collected at the ESO VLT (prog. ID 186.A0798).
References
 Angrick, C., & Bartelmann, M. 2009, A&A, 494, 461 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Angrick, C., & Bartelmann, M. 2012, A&A, 538, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Annunziatella, M., Biviano, A., Mercurio, A., et al. 2014, A&A, 571, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barreira, A., Li, B., Jennings, E., et al. 2015, MNRAS, 454, 4085 [NASA ADS] [CrossRef] [Google Scholar]
 Bartelmann, M., Narayan, R., Seitz, S., & Schneider, P. 1996, ApJ, 464, L115 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J., & Mamon, G. A. 1982, MNRAS, 200, 361 [NASA ADS] [CrossRef] [Google Scholar]
 Biviano, A., Rosati, P., Balestra, I., et al. 2013, A&A, 558, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cacciato, M., Bartelmann, M., Meneghetti, M., & Moscardini, L. 2006, A&A, 458, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Coe, D., Umetsu, K., Zitrin, A., et al. 2012, ApJ, 757, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Dejonghe, H., & Merritt, D. 1992, ApJ, 391, 531 [NASA ADS] [CrossRef] [Google Scholar]
 Donahue, M., Voit, G. M., Mahdavi, A., et al. 2014, ApJ, 794, 136 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Edge, A. C., & Henry, J. P. 2001, ApJ, 553, 668 [NASA ADS] [CrossRef] [Google Scholar]
 Ebeling, H., Ma, C. J., Kneib, J.P., et al. 2009, MNRAS, 395, 1213 [NASA ADS] [CrossRef] [Google Scholar]
 Golse, G., & Kneib, J.P. 2002, A&A, 390, 821 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Harrison, E. R. 1974, ApJ, 191, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Konrad, S., Majer, C. L., Meyer, S., Sarli, E., & Bartelmann, M. 2013, A&A, 553, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Łokas, E. L. 2002, MNRAS, 333, 697 [NASA ADS] [CrossRef] [Google Scholar]
 Lucy, L. B. 1974, AJ, 79, 745 [NASA ADS] [CrossRef] [Google Scholar]
 Lucy, L. B. 1994, A&A, 289, 983 [NASA ADS] [Google Scholar]
 Majer, C. L., Meyer, S., Konrad, S., Sarli, E., & Bartelmann, M. 2013, unpublished [arXiv:1304.6522] [Google Scholar]
 Mamon, G. A., & Boué, G. 2010, MNRAS, 401, 2433 [NASA ADS] [CrossRef] [Google Scholar]
 Mamon, G. A., Biviano, A., & Boué, G. 2013, MNRAS, 429, 3079 [NASA ADS] [CrossRef] [Google Scholar]
 Martino, R., Mazzotta, P., Bourdin, H., et al. 2014, MNRAS, 443, 2342 [NASA ADS] [CrossRef] [Google Scholar]
 Merten, J. 2014, ArXiv eprints [arXiv:1412.5186] [Google Scholar]
 Merten, J., Cacciato, M., Meneghetti, M., Mignone, C., & Bartelmann, M. 2009, A&A, 500, 681 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Merten, J., Coe, D., Dupke, R., et al. 2011, MNRAS, 417, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Merten, J., Meneghetti, M., Postman, M., et al. 2015, ApJ, 806, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Bartelmann, M. 1996, ArXiv eprints [arXiv:astroph/9606001] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 [NASA ADS] [CrossRef] [Google Scholar]
 Newman, A. B., Treu, T., Ellis, R. S., et al. 2013, ApJ, 765, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration Int. III. 2013, A&A, 550, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XX. 2014, A&A, 571, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Postman, M., Coe, D., Benítez, N., et al. 2012, ApJS, 199, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Richardson, T., & Fairbairn, M. 2013, MNRAS, 432, 3361 [NASA ADS] [CrossRef] [Google Scholar]
 Rosati, P., Balestra, I., Grillo, C., et al. 2014, The Messenger, 158, 48 [NASA ADS] [Google Scholar]
 Sarli, E., Meyer, S., Meneghetti, M., et al. 2014, A&A, 570, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sartoris, B., Biviano, A., Rosati, P., et al. 2014, ApJ, 783, L11 [NASA ADS] [CrossRef] [Google Scholar]
 Umetsu, K., Medezinski, E., Nonino, M., et al. 2012, ApJ, 755, 56 [NASA ADS] [CrossRef] [Google Scholar]
 Umetsu, K., Medezinski, E., Nonino, M., et al. 2014, ApJ, 795, 163 [NASA ADS] [CrossRef] [Google Scholar]
 van der Marel, R. P., Magorrian, J., Carlberg, R. G., Yee, H. K. C., & Ellingson, E. 2000, AJ, 119, 2038 [NASA ADS] [CrossRef] [Google Scholar]
 von der Linden, A., Mantz, A., Allen, S. W., et al. 2014, MNRAS, 443, 1973 [NASA ADS] [CrossRef] [Google Scholar]
 Wolf, J., Martinez, G. D., Bullock, J. S., et al. 2010, MNRAS, 406, 1220 [NASA ADS] [Google Scholar]
 Zitrin, A., Rosati, P., Nonino, M., et al. 2012, ApJ, 749, 97 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Discussion of other methods
On top of the conditions leading to the Jeans equation, the method presented above implicitly rests on the assumption that all test particles (i.e. galaxies) have the same mass, that their number density is proportional to their mass density and that one can establish a polytropic relation between their mass density ρ and the effective galaxy pressure . Fixing the velocity anisotropy profile β(r) enables us to constrain the gravitational potential of the cluster via the radial Jeans equation using the observed densityweighted galaxy velocity dispersion. Our method does not assume the mass to follow light a priori, i.e. ρ_{DM} does not necessarily need to follow ρ_{gal}. Since the resulting potential is nonparametric, one can test the validity of particular parametric models as done above.
Another common deprojection method is the Abel inversion using derivatives of observables. However, because of the fluctuating nature of most observables, which cause strongly fluctuating derivatives, we decided to use the RichardsonLucy deprojection instead. This deprojection only involves integrals.
There are many other methods to determine the gravitational potential of a cluster, or its mass, or to constrain the velocity anisotropy. For example, Binney & Mamon (1982) use the lineofsight velocity dispersions as observables together with the surface brightness in spherical galaxies. Assuming a constant masstolight ratio, they can replace the unknown density in the Jeans equation by the luminosity, and can thus determine the velocity anisotropy profile β(r). However, in the context of galaxy clusters it is not clear that the same relation between mass and luminosity should hold as well.
Dejonghe & Merritt (1992) determine the potential Φ by taking into account higher than secondorder velocity moments of the collisionless Boltzmann equation (CBE). Given the velocity dispersions, and making use of a linear relation between the observables and Φ, they expand the gravitational potential in terms of basis functions with unspecified coefficients. Using the positivity constraint of the distribution function and higher order moments of the CBE, they are able to determine the coefficients and thus the potential. However, from a practical point of view, it becomes increasingly difficult to constrain higherorder moments of the velocity distribution from their observable lineofsight projections, as the projection integrals turn out to be substantially more complex with increasing order (see Richardson & Fairbairn 2013, for instance).
A method very similar to ours is proposed by Mamon & Boué (2010). Starting again from the Jeans equation and fixing β(r), they determine both the number density and the velocity dispersion observationally. Besides their using the Abel inversion, the only difference from our method is that they determine the number density of galaxies appearing in the Jeans equation by counting. Hence, they do not need our assumed polytropic relation between density and effective pressure. In future work, one could use this to test directly to what degree the polytropic assumption is appropriate.
Wolf et al. (2010) and van der Marel et al. (2000) are considering a general parametrisation of β(r) and of the mass M(r) and constrain all free parameters via a maximumlikelihood analysis of the velocity dispersion predicted by the Jeans equation compared to the observed dispersion. van der Marel et al. (2000) use a constant velocity anisotropy whereas Wolf et al. (2010) take a varying parametric profile.
Another way to break the anisotropy mass degeneracy is taking the fourth order velocity moments of the CBE into account (Łokas (2002); see also Richardson & Fairbairn (2013) in this context). By expressing the velocity distribution function in terms of two integrals of motion, the energy E and the angular momentum L, and further assuming that the distribution is separable in E and L with β = const., the two fourthorder moments are reduced to one equation involving the radial velocity to the fourth power, β, the radial velocity dispersion and Φ. Thus, by measuring the velocity dispersion and , one can infer the gravitational potential by solving this fourthorder equation together with the Jeans equation.
The MAMPOSSt method introduced by Mamon et al. (2013) breaks the degeneracy assuming parametric forms for the gravitational potential, the velocity anisotropy, and the distribution of 3Dvelocities. By performing a maximumlikelihood fit of the galaxy distribution in the projected phase space, all parameters can be determined. However, it has recently been shown (Richardson & Fairbairn (2013) that their assumption of a Gaussian 3Dvelocity distribution is incompatible with the equilibrium assumption underlying the CBE. One would thus have to choose a different 3Dvelocity distribution.
All Tables
All Figures
Fig. 1 Scheme of the reconstruction algorithm. The observable is converted to a threedimensional quantity by RichardsonLucy deprojection. This is then turned into the threedimensional gravitational potential, solving the Jeans equation. The resulting potential is finally projected along the line of sight. 

Open with DEXTER  
In the text 
Fig. 2 Velocities projected along the l.o.s. of cluster members. 

Open with DEXTER  
In the text 
Fig. 3 Reconstructed radial profiles of the projected gravitational potential of MACS J1206, obtained with different parameter sets for the smoothing amplitude α (panel a)), smoothing scale L (panel b)), the polytropic index γ (panel c)), and the anisotropy parameter β (panel d)). 

Open with DEXTER  
In the text 
Fig. 4 Radial profiles of the projected gravitational potential, reconstructed from galaxy kinematics with different profiles of the anisotropy parameter β. The other reconstruction parameters remain fixed at α = 0.3, L = 0.6 and γ = 1.1. 

Open with DEXTER  
In the text 
Fig. 5 Twodimensional map of the lensing potential of MACS1206 based on weak and strong lensing. 

Open with DEXTER  
In the text 
Fig. 6 Radial profiles of the gravitational potential of MACS J1206, reconstructed from combined strong and weak gravitational lensing (blue), Xray emission (yellow), and from galaxy kinematics (red) with the parameters α = 0.3, L = 0.6, γ = 1.1 and β as in (9). 

Open with DEXTER  
In the text 
Fig. 7 Comparison of three analytic profiles of gravitationallensing potentials with the gravitationalpotential profile of MACS J1206 reconstructed from galaxy kinematics using α = 0.3, L = 0.6, γ = 1.1 and β as in Eq. (9). 

Open with DEXTER  
In the text 